matlab 符號運(yùn)算的化簡問題

shouce

首先來個簡單

! J2 M% l0 d8 U  L1 }( x>> syms t r a b q
, ]9 r; ~- a) K" Q' w: c% q>> x=r*cos(t)+a;
; t5 M! \& Q# _* z/ _' v>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
& F2 H" W& \7 a" s>> simplify(X)

ans =
- ?4 x" q/ z5 }: T& {# e, S
- a*cos(q) - b*sin(q) - r*cos(q - t)   

復(fù)雜點(diǎn)   就不化簡了
8 s. @- ^9 G. u* W, N) P >> syms t r a b k q A
0 }' i& c7 o# `* k# m- u% \( c% s>>  x=r*cos(t)+a;
7 N) n7 L0 [0 m# `6 w>> y=r*sin(t)+b;
2 V! C" `7 }  x# @4 M>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
>> NM=N*M
. m7 k; ]8 R2 z' Z! @
) d( j- }3 p" Z$ [& D5 n3 S& CNM =
9 j$ l0 s3 M9 d, ?# [5 i) X6 w
7 b$ K& R% P8 b) S) q 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)
, K: m' z  k9 [- K/ V9 C
>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
% m9 @9 O6 X2 K9 Q/ k4 }( G+ P8 a; z>> simplify(x1)

ans =
5 w2 Z" r, z! m0 c& x5 @
" m( `8 `! b0 Q8 X& RA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化簡無效果    按正常應(yīng)該是  
4 g6 X& E$ y- _8 @: h* T* a8 p
, s- r& u% c$ K  L0 T& g# o9 ~- kA*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);
/ `' S  U# k8 h$ w, q9 h& x>> simplify(y1)
4 e* r3 w- B5 W5 m! D; F# V
% D+ n& k8 ~8 U" N4 sans =
3 R$ v0 e9 D7 p& K5 H: N9 ]
! L$ F( h% X" I' q' i( H- \cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
5 F8 j' L- p1 D; X
  V* |) ?  A; \8 {

明月山河

simplifty函數(shù)要有控制參數(shù)才會按照你想像的運(yùn)算。

shouce

明月山河 發(fā)表于 2016-2-20 18:27
3 @, o# j/ n* z" {4 `7 _simplifty函數(shù)要有控制參數(shù)才會按照你想像的運(yùn)算。

+ Z) ]& b% _8 m" b. F+ I; P& L大俠   能舉個例子嗎   我這個不懂

明月山河

shouce 發(fā)表于 2016-2-20 18:33
) j, x* v' ~0 p. K大俠   能舉個例子嗎   我這個不懂

5 u7 B. C! m- u2 T7 Ghttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

0 U6 H& N% r. L9 n3.帶有附加條件的化簡
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
7 ^% Z# g2 W' q( b! p5 i>simplify(expr,{siderel},[y,x]);
' K, N# d& f0 I6 i5 Q. d4 b#maple在表達(dá)式中作代換y^2=1-x^2，然后再試圖對x^2做代換，由于找不到，就停止了
4.展開與組合
expand &combine
combine是與expand相反的命令，它可以按照數(shù)學(xué)規(guī)則將表達(dá)式中的某些項組合在一起
# U3 o2 P2 {  _4 j- Q
7 g! C/ K7 w. n0 }1 E$ W
* m( k0 Q1 o+ l
, R2 H/ c/ E5 _# l$ J6 h1 M  l( Wmatlab是用的maple的符號運(yùn)算內(nèi)核，使用規(guī)則一樣。
6 k# d6 p! W" k
0 Z" d1 D0 c. B  q7 Z  v

shouce

明月山河 發(fā)表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.帶有附加條件的化簡

多謝了   等下去試試
% u* T1 T8 I1 j4 e/ @

shouce

明月山河 發(fā)表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
5 T. ?  ]6 z* I* ?+ ~* @8 t# V
2 i  O3 K& `9 x, E) e' Z  g, u3.帶有附加條件的化簡

試了下還是不行   可能matlab和maple 使用不一樣
>> syms x
( @, a- [& }& S% U9 Y7 R3 j' d  \>> combine(sin(x)^2+cos(x)^2)   

ans =
" f& H& N5 G3 _8 W7 l6 i" V
cos(x)^2 + sin(x)^2          這個不化簡單

>> simplify(sin(x)^2+cos(x)^2)
. q0 Z7 F, h9 l, U4 X, k
. P$ [. w2 b4 p. d3 B1 L  zans =

6 J1 A1 |2 x$ }) x1


>> syms x
: k+ p+ o/ s2 L  M$ E! d' Z+ I- O( E7 a>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
0 Z; O/ F; P' F+ g& f7 X" }
A =

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

>> simplify(A)

/ Y  C- H0 r. [ans =

7 W4 X$ \0 \  J/ Dlog(3*x) + 1


5 o1 d, ?. {; E3 P. \
) ^/ u0 I1 k: }7 L>> simplify(A,x)
錯誤使用 sym/simplify (line 99)
Too many input arguments.

+ G# v( s4 S' n5 z

393810298

matlab大學(xué)里學(xué)完了，我就在也沒有用過了，已經(jīng)看不懂了

shouce

>>  expand(x1)
+ u8 X, R+ }/ R. pans =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
  F1 T" T2 w1 e4 h* a, W: e>> 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 c  L7 X. s% }1 V& v5 R; EKMM =
8 @4 m3 m( |7 P. G1 z/ d5 ]$ C$ u& gA*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)

要先展開然后合并