matlab 符號運算的化簡問題

shouce

首先來個簡單

3 m3 ?& x9 R. w3 b& I; R>> syms t r a b q
>> x=r*cos(t)+a;
& X9 D. @- b  ?; u+ [>> y=r*sin(t)+b;
) ?; B. ?) R; t8 Y) n>> X=-x*cos(q)-y*sin(q);
>> simplify(X)
5 n# C) v  l5 U6 @. [, J
4 |$ [) Z+ J# Cans =
: T* d$ q+ Z: e$ B. j/ a
- a*cos(q) - b*sin(q) - r*cos(q - t)   

復(fù)雜點   就不化簡了
>> syms t r a b k q A
/ l) W0 q; |9 S- D0 z( L# E>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
0 W( P' D# E, O7 t>>  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 =
. B+ Y7 F( T7 ~- T- J
, b2 s2 `; q+ l A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
& M) [* e6 a9 q. q% E$ i cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
" t- X, z9 h+ [' f  i' T
; r6 S& N; ~  z* ^$ K% q>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
5 n6 [: _& S; ~0 I6 d0 {>> simplify(x1)

ans =
/ B" K) d9 U- S7 s' e7 J0 M7 s0 d
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化簡無效果    按正常應(yīng)該是  

  V5 }" D) {- ~6 O+ c3 b8 QA*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);
7 M" W8 D1 V, C. p>> simplify(y1)
9 `* I2 }& _6 ^- G9 R
ans =
# A6 T. C* i- S( Z: G! H* t
) ]! [6 B) x% Y8 {& t( i4 ]. F1 {cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
0 i6 A4 r) h: P

明月山河

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

shouce

明月山河 發(fā)表于 2016-2-20 18:27
simplifty函數(shù)要有控制參數(shù)才會按照你想像的運算。

5 z2 U1 E: ^, l1 c大俠   能舉個例子嗎   我這個不懂

明月山河

shouce 發(fā)表于 2016-2-20 18:33
大俠   能舉個例子嗎   我這個不懂

. H  J$ z! w- A5 d) f4 Lhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
8 E7 [6 |9 L- z. u9 t( i
: I$ r7 \# x8 c( S+ m/ H4 p3.帶有附加條件的化簡
5 s, U' ^$ H  w! ~/ U* K3 T>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
; N+ I; o! T# r, b. k#maple在表達式中作代換y^2=1-x^2，然后再試圖對x^2做代換，由于找不到，就停止了
. j: M2 w% i8 h, s$ a& s; o# F4.展開與組合
expand &combine
combine是與expand相反的命令，它可以按照數(shù)學(xué)規(guī)則將表達式中的某些項組合在一起

( K+ W- H/ c/ T
9 V$ G& `" Q6 e  e. W6 m1 F
matlab是用的maple的符號運算內(nèi)核，使用規(guī)則一樣。

shouce

明月山河 發(fā)表于 2016-2-20 18:38
0 n$ Z4 l3 A+ Y7 x  Ahttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
& ?3 d5 t: F( G3 Y- y* I
$ u3 r; W7 A7 p5 E/ [' \" f5 u9 W3.帶有附加條件的化簡

多謝了   等下去試試

shouce

明月山河 發(fā)表于 2016-2-20 18:38
" d; G" H3 ~! Q5 Zhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

* Q6 E# U  P5 W' e; M3 X& T3.帶有附加條件的化簡

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

' g& Y% U( V9 F6 d9 Gans =
2 b9 z0 W6 _$ L( v1 V. W
) ~2 l9 `2 z0 O3 }, t7 n. _5 z2 Wcos(x)^2 + sin(x)^2          這個不化簡單
7 N+ R, W0 y( w9 A; l; D
* p  C& e, D' B2 k4 s, v& f>> simplify(sin(x)^2+cos(x)^2)
  c  U& l# \9 Y7 P/ D
* K" L0 H$ |: v/ Oans =

, @8 r' q( G5 `6 G* ^1
6 o: E. P! i3 [

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

* S$ r$ a4 n2 G+ o( Z! gA =
% b- l  I8 T0 t! b$ P; K
log(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)
+ ^/ K/ `5 q$ B3 r4 Z0 M( v
ans =

6 F+ E/ [1 R. ~  Y2 w- Alog(3*x) + 1
8 A4 K+ \* h, _' _" N& c* V


% q: p0 C: V: c/ M, b8 s9 q6 ^>> simplify(A,x)
7 i( c" g# C: D& N錯誤使用 sym/simplify (line 99)
Too many input arguments.

: B, B1 J; V2 B7 ^# t. ~

393810298

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

shouce

>>  expand(x1)
% R6 o2 L- d  ?) X2 tans =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
+ r8 m7 X$ N6 K2 c9 X0 I* |>> 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)')
& c+ C4 T% Z2 V6 ]! q8 L. mKMM =
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 =
1 U. ~& [' P; c$ y. i9 W0 m5 RA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展開然后合并