matlab 符号运算的化简问题

shouce

首先来个简单

3 q6 ?- D! \& }- I, @, y>> syms t r a b q
) Z6 J' E$ w& ]0 B/ v' e>> x=r*cos(t)+a;
1 z) o/ q- W% E! |. X" [) L2 F; b+ S>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
; A0 f: G: k' W2 D4 Q>> simplify(X)

6 Y* j3 g' `! [; Z! J% p" `- T4 Q5 d5 Kans =

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

6 B# R. q* s0 c8 a& R5 D复杂点   就不化简了
9 ?1 p, K1 |/ h4 e) D2 z2 v >> syms t r a b k q A
+ i% N# E5 U4 c+ }6 X>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
8 b% \" C" L. f! S>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
: i5 T* u! `$ ~% q; J>> M=[x;y;A];
>> NM=N*M

NM =
& G7 M2 m$ c: C* Y7 |6 r+ d
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)
" P2 \9 H% l, W  s2 H
>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
. C: P, G2 N8 X1 q2 H+ T>> simplify(x1)

: k  P$ A4 ^" n1 dans =

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

! U6 K: p0 a* \7 T1 j8 C7 @A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

( o: U( D" S$ p! l6 g: G& `( v>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
. x7 H! T/ H' T7 i>> simplify(y1)

" @: X! A! s3 E! J7 W& }; tans =

. W4 ~" X" ]$ ?( r5 k! S3 F" Zcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

) N1 m5 p/ n: p/ T+ m3 U

明月山河

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

shouce

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

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

明月山河

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

/ s. Y8 }( q9 q! Mhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

. ]. M1 ]' V$ Q3.带有附加条件的化简
>simplify(expr,{x*y=1});
% `- B6 ?2 N; a>siderel:=x^2+y^2=1;
" V. }# k& I( p9 F+ ^>simplify(expr,{siderel},[y,x]);
1 ?; S2 K. M' H#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
0 X- o6 w- p) eexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
3 j+ a/ \' J! P


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

shouce

明月山河 发表于 2016-2-20 18:38
) j( D) b: b5 I7 Q# _http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

9 U0 w4 f7 ~8 S7 ]5 o; }9 l3.带有附加条件的化简

) j# [' L/ l: P$ E" X多谢了   等下去试试
5 j- W0 @& b, [' V2 Q& j

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

( V4 u; y8 U5 L) I3.带有附加条件的化简

试了下还是不行   可能matlab和maple 使用不一样
>> syms x
>> combine(sin(x)^2+cos(x)^2)   
/ k; \# Z' r; w. }% r
ans =
$ p# Y/ a: L; n/ g; I' d+ A
cos(x)^2 + sin(x)^2          这个不化简单
1 W  I! B- I1 b9 }
>> simplify(sin(x)^2+cos(x)^2)
; v# h; s. P/ G) j0 a3 o# V
ans =

1
5 S! P+ \* j; H& D$ Q, r
  d* `" s% p$ ?9 L! o
: ^1 O- ]1 y9 W% g>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
  ^- V2 ?% s" f5 r+ R0 S& c
A =

/ k% _) D3 m# ]6 U6 ]log(3*x) + cos(x)^2 + sin(x)^2

1 u+ g3 @/ S$ D1 V  Z7 L. }* N1 `>> simplify(A)

  I" T! e; [5 [ans =

log(3*x) + 1
! F, L0 |, J0 N7 ~& o! r
& e3 D! T6 |$ X1 F# q0 _

>> simplify(A,x)
错误使用 sym/simplify (line 99)
1 d% g+ ~+ M' v6 Q* }+ R7 i' _( qToo many input arguments.

0 a1 b  p6 t" Q) A: a

393810298

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

shouce

>>  expand(x1)
ans =
( d* D2 ~- G* K! f# DA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
7 f* l9 j  V6 t8 j! x7 h>> 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)')
4 }2 `& e, Q, x, B/ a0 oKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
" z& w0 J* w" W7 C* X>> simplify(KMM)
7 C8 W& O: i. O# n( J4 e, {ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
4 H; @2 R! G9 x6 ~5 i, E8 x' w  D4 G
要先展开然后合并