matlab 符号运算的化简问题

shouce

首先来个简单

>> syms t r a b q
* U0 `0 k/ k, S1 l% u; y>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
9 z0 g* F% i( L$ _>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

ans =

1 }& l: W5 @' b0 m1 ]8 n- a*cos(q) - b*sin(q) - r*cos(q - t)   

# {" Z; K* t1 y复杂点   就不化简了
>> 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)];
+ U0 a* C4 U+ @>> M=[x;y;A];
( T2 E% j2 Z: J# {% y' v8 _>> NM=N*M
5 \' \  M& u; o1 q5 z
) g$ b  ^  Y% n5 l7 mNM =
4 q+ M2 Q6 t- v+ V2 n. j9 [: T
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
; o+ b4 h" b# I$ K 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 =
$ N! M. c3 @, f# e* N
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
  D3 M* D+ c! l. R- |8 M: d
  f5 W! }: ^5 R- b. u( WA*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)
2 O; G% D" s/ g# W( N0 N: R9 N* W
ans =
- H+ G1 G$ U6 R5 F' K
5 N) T; Q- C( y/ ]cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
0 C' u/ A% \1 x; _7 E
& I: C0 ?; L* ]; J8 Z3 h

明月山河

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

shouce

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

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

明月山河

shouce 发表于 2016-2-20 18:33
! O$ x' R1 v3 Y) K: x- A6 }大侠   能举个例子吗   我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
+ L. {- w9 Q, S( \$ F6 ]
3.带有附加条件的化简
0 C3 P: H1 J" T! Q. z3 y; Y>simplify(expr,{x*y=1});
* c1 ^: x' {+ E# C>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
' D  h6 m) I  p  m: @4 `#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
# F( K# ?! F% _& l; Y7 z2 Qexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起

8 z# U" j% x9 R  y

6 q3 ^3 P, x- w+ _0 o: {matlab是用的maple的符号运算内核，使用规则一样。

& u& W" g' I  A4 A' M4 O

shouce

明月山河 发表于 2016-2-20 18:38
' ~: ?: {7 w. I: n% P  G: I4 X: R6 khttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
6 M' i" p4 Q" X$ W9 e" Q
' ]* G  Z) X: z; ?) I% B3.带有附加条件的化简

1 W0 R4 ]$ F5 F2 ?: o多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
6 e+ m. g  \8 O) E  thttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
2 d) L# x5 V* u- [$ Z, O
: u4 [0 I  L: f3.带有附加条件的化简

/ G2 @& \( J; E1 F9 F试了下还是不行   可能matlab和maple 使用不一样
>> syms x
# X5 b1 _- c5 r1 e>> combine(sin(x)^2+cos(x)^2)   

# [* @$ b: X1 i! @ans =
& h4 `4 p# p& G9 i3 m* u4 o
cos(x)^2 + sin(x)^2          这个不化简单
' ~7 l% L7 |8 R
) N$ H" u6 I. w, n0 W9 e>> simplify(sin(x)^2+cos(x)^2)

/ d5 [) w9 q* P* Ians =

/ H4 p) j0 o: ?" A4 g8 P/ N* i' G1


>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
5 I, ?8 l2 g- i7 p8 G' Z
4 k0 u! o! a2 w+ Q' Q! XA =

log(3*x) + cos(x)^2 + sin(x)^2
& \2 }# O: r9 L2 y) p# _' C
>> simplify(A)

+ k1 c0 n0 T& l1 c- Qans =

7 u9 Q! b6 A" E2 k: H( ?' Q' M; Llog(3*x) + 1


. K5 {- a( J/ P1 T3 Q5 _3 z+ |% w
, w6 P8 E$ b& o>> simplify(A,x)
错误使用 sym/simplify (line 99)
Too many input arguments.

393810298

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

shouce

>>  expand(x1)
) u/ x' o. l+ A7 b% H8 Kans =
" Q. u. Q5 I* T3 F7 l3 r& kA*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)')
( _! I& _' L. {: I% OKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
* \2 W5 j7 {2 a2 c, K; A5 s  `- q>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
; {% L) ^, q. g
要先展开然后合并