matlab 符号运算的化简问题

shouce · 发表于 2016-2-20 17:49:24

首先来个简单

>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

ans =

- a*cos(q) - b*sin(q) - r*cos(q - t)

复杂点就不化简了
>> syms t r a b k q A
>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
>> NM=N*M

NM =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)

ans =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))    化简无效果按正常应该是

A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

ans =

cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

明月山河 · 发表于 2016-2-20 18:27:03

simplifty函数要有控制参数才会按照你想像的运算。

shouce · 发表于 2016-2-20 18:33:16

明月山河发表于 2016-2-20 18:27
% t* l$ I# W# G: n; G1 m3 ~4 o+ ysimplifty函数要有控制参数才会按照你想像的运算。

大侠能举个例子吗我这个不懂

明月山河 · 发表于 2016-2-20 18:38:11

shouce 发表于 2016-2-20 18:33
( w, M/ ?% p+ t2 I大侠能举个例子吗我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
expand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起

matlab是用的maple的符号运算内核，使用规则一样。

shouce · 发表于 2016-2-20 18:42:54

明月山河发表于 2016-2-20 18:38
- x4 e8 L2 S& t# k+ e2 }' U6 }http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
. j3 r1 l* p3 U1 \# D' ?- H
+ v' ~- C" A6 O/ t. N+ u: c1 d3.带有附加条件的化简

多谢了等下去试试

shouce · 发表于 2016-2-20 19:09:49

明月山河发表于 2016-2-20 18:38
8 A4 v" ]1 O2 P( o0 W2 Z& Ihttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
' P9 s2 V2 M) d- R- X* v- Q& D3 v7 u2 _- n1 z* [ A7 F; P
3.带有附加条件的化简

试了下还是不行可能matlab和maple 使用不一样
>> syms x
>> combine(sin(x)^2+cos(x)^2)

ans =

cos(x)^2 + sin(x)^2 这个不化简单

>> simplify(sin(x)^2+cos(x)^2)

ans =

1

>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =

log(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)

ans =

log(3*x) + 1

>> simplify(A,x)
错误使用 sym/simplify (line 99)
Too many input arguments.

393810298 · 发表于 2016-2-20 19:28:09

matlab大学里学完了，我就在也没有用过了，已经看不懂了

shouce · 发表于 2016-3-23 23:44:22

>> expand(x1)
ans =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> KMM=sym('A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)')
KMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并

		自动登录	找回密码
密码			注册会员

matlab 符号运算的化简问题

相关帖子

点评