matlab 符号运算的化简问题

shouce

首先来个简单
; f: [$ l8 q; P# P
>> syms t r a b q
, \2 l+ Q5 `; S/ y8 x>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)
; H& b  p, x. ~
ans =
/ l# t6 y1 `- b8 `" p" e0 A0 b
- a*cos(q) - b*sin(q) - r*cos(q - t)   

6 P$ J% f9 j  t  p- d2 _复杂点   就不化简了
>> syms t r a b k q A
1 U" u& W, X# s>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
3 l9 G' H* C- V4 J1 G>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
/ T  D( }( L7 d8 i5 N>> M=[x;y;A];
; P3 q$ t  Y+ }>> NM=N*M

1 C  y6 [7 Y$ oNM =

) }; y9 h2 Z# b 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)

, x7 z( [; U5 H/ z" f>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)

ans =

7 o0 l& ~2 q$ q3 F" SA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
1 j8 Z( o2 v- M4 h  d5 N+ e! U6 B
A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
/ Z7 e. p% c* r3 B# g
>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)
) x' x; U+ q" L2 j2 w1 W: T0 U
- h# u; E- _& {' Z: D/ {2 b% ]ans =

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

1 A1 C% i- R  ]' z* n

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
simplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

shouce 发表于 2016-2-20 18:33
大侠   能举个例子吗   我这个不懂

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

3.带有附加条件的化简
+ r. h" u- Y2 o8 q>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
: ^3 r/ B& L3 h4 K3 _) s* ^; P2 M>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4 D9 S% n( B1 }! ^+ L; n4.展开与组合
expand &combine
- f1 b' J3 x$ [5 Zcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起

* E; r  }$ ?2 o* ^$ s8 o. h

matlab是用的maple的符号运算内核，使用规则一样。
4 u8 V3 C9 X/ l- q/ l: C

shouce

明月山河 发表于 2016-2-20 18:38
1 O; J# T( C1 @$ W5 }  {: F) R& c0 `http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

2 V+ O5 n0 U+ a  w, x$ w3.带有附加条件的化简

8 M) e. Y7 j+ u" Q; g% M多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
3 _1 X, x: `0 K1 R! f5 x
3.带有附加条件的化简

; {* m! d# {; L; g/ u试了下还是不行   可能matlab和maple 使用不一样
) L3 C, q5 _) D# V1 [>> syms x
: o$ R, c, c- ?! I3 }) R+ S: r* B>> combine(sin(x)^2+cos(x)^2)   

ans =

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

- q0 [' M+ p, `* g3 g' Y>> simplify(sin(x)^2+cos(x)^2)
; H- h  ^: d) W5 B
ans =

' P- c  A# |8 c5 I1


8 @7 Z& r7 D( e>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

( x1 P9 @, f! c6 VA =
( F) [- x0 ~! S) m
2 ^8 |! h+ Z+ d. _8 @; dlog(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)

; {! C) X) b$ y$ m5 E7 N  Oans =
' m& d" O6 l/ w% n$ I
& S1 ^2 U" G* l0 t' [log(3*x) + 1
: I) }9 n3 y( r9 r


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

393810298

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

shouce

>>  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)')
0 `% Q9 B0 h6 e3 S& W/ A  pKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
! E6 [1 m- T/ h/ q2 _>> simplify(KMM)
9 ?: `1 Q1 J* C- [& qans =
7 x' o; ?/ u; M$ \; H: P% PA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
: e, m) y* n3 N) b( u* _
( ^/ {/ t5 K- r; ?. R, }要先展开然后合并
* f! m6 Y% W6 W! v  m! P