matlab 符号运算的化简问题

shouce

首先来个简单

>> syms t r a b q
  S+ X  o9 C2 Z) _! }# u>> x=r*cos(t)+a;
8 F) x+ {+ k+ T7 T>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

+ ]. Y) s. e/ v% o- {& \. N1 a; vans =

- a*cos(q) - b*sin(q) - r*cos(q - t)   
$ [1 H  ]( e' T' S7 N6 `: L# N
1 A  t0 W  X' q. E# g复杂点   就不化简了
2 l& Q5 L  N1 a* P/ R" H% X >> syms t r a b k q A
>>  x=r*cos(t)+a;
2 X/ j9 d* [( y! S+ h. v% |: ~' ~>> y=r*sin(t)+b;
>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
4 e6 v7 n" e* ^: L$ @+ F) ?$ a9 I; F>> M=[x;y;A];
>> NM=N*M
( |  B/ X+ P6 r0 \; v
NM =

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)

* a% Y/ B- l. w% q0 D9 F$ ~- r>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
9 ]1 [7 X: }" m$ A- [! q>> simplify(x1)
6 e! w& \3 L6 @; z4 e( S% ~
5 h9 Z; u3 p6 _% G' P/ Aans =

! D/ V$ d$ p4 m0 x3 o+ _A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
: i' _; j6 L8 @: K7 P
2 v3 J) ]! Q. T" t2 \A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

6 }" V/ G: j: ?. Z* K# R>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

$ M+ F  R, M) k% ^$ J1 ^, qans =

cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
6 F# F) B9 _2 }, E" G. A4 M

明月山河

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

shouce

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

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

明月山河

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

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

3.带有附加条件的化简
>simplify(expr,{x*y=1});
# \6 W7 a; y. I1 d' m' X* l; b>siderel:=x^2+y^2=1;
7 D- H# s7 `5 i) i9 @' |>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
expand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
: _$ R4 E: R' f. E
! _. M! t* u9 X& k

( ^$ Q1 \2 ?: \6 Lmatlab是用的maple的符号运算内核，使用规则一样。

9 o& \" I% a" f: e* J& a

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
: [" I: B2 W7 w* m- k
. J1 H4 ]; G% J! ^3 n* ^3.带有附加条件的化简

多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
) Z7 |, `$ N8 r1 X+ v/ q: c$ o2 Z! Qhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
3 `% T$ N* o7 x
/ {6 [5 f* m7 p* \6 a! A1 |# ]3.带有附加条件的化简

" ~1 ]: J. O5 x试了下还是不行   可能matlab和maple 使用不一样
& ^( {5 @" R! m1 Z" o' z) z* {% ]>> syms x
4 t7 s" r  n) I' _6 r/ ^# j. M>> combine(sin(x)^2+cos(x)^2)   
/ ?9 i4 N4 K1 Y! d/ ^$ Q
! c# }* h( |$ }; B( ]" n( V! gans =
! `1 T1 s7 y# j# s" Q. \
$ I, H1 D) r. |' B" wcos(x)^2 + sin(x)^2          这个不化简单
8 [8 E  ]4 i" v# D* f" j
>> simplify(sin(x)^2+cos(x)^2)
3 {+ K% u, `9 o! u) v0 B
/ }4 D. E+ G% E% ?/ q8 Lans =

; S+ k7 J  R1 }: j2 B0 l$ C1
( T9 d: C$ O% Y
& m* }) R# J' u+ D% r  Y. d
3 }% E" Z) m" b: |( W>> syms x
# T5 h* a/ T& p' l1 R" s1 [; p>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
: ~9 }% K- x9 _
) k' L; O5 M4 {) ~/ d; G; u* UA =

4 W# o: l5 B* T1 R& Tlog(3*x) + cos(x)^2 + sin(x)^2
5 F+ `7 S8 P" h" l$ D  M' x
>> simplify(A)
8 M& m. l: t( }# E
ans =
8 q9 n% l7 }9 g, I' V
  V& r2 a: U4 x/ \  [; Ylog(3*x) + 1

; b1 i% I/ k$ i$ d5 w
" x/ I3 `, h1 r) Z$ F5 t* G
( l8 E7 _& p* E3 j1 Z! R: m. t>> simplify(A,x)
$ z$ K, N4 g2 ~. W( C6 ?# l错误使用 sym/simplify (line 99)
Too many input arguments.
8 i, \1 [2 P& `( ~2 D% x
3 |' K* l- Q0 Z+ c

393810298

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

shouce

>>  expand(x1)
; Y: T# }* U. I" s6 Ians =
+ D- F8 J! ]/ t$ [; V" q, Q; YA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
$ [2 S& s. ?  A7 k, C- D6 s>> 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)')
9 l7 R8 N9 n& F8 j6 NKMM =
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)
( O+ S7 v4 A" kans =
7 w6 ^1 W! \* Q; E* _$ EA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并
- N  f$ o; J& k# {1 j# d