matlab 符号运算的化简问题

shouce

首先来个简单
9 k. ]8 A$ _5 x% q% j
>> syms t r a b q
>> x=r*cos(t)+a;
* I3 o# F: m2 ]7 O2 K>> y=r*sin(t)+b;
* x( N' t) T! Q: n% b* n1 @' n& K>> X=-x*cos(q)-y*sin(q);
( z, g; r0 O' b5 a: D) d  D) s>> simplify(X)
- l- |1 A$ A. _% _' R9 v' f
5 q  y$ \8 _  s' Wans =
4 ^: e5 d/ n  b( p: J
, E4 e) Y  }% ?, W( d: ?# ^4 I- a*cos(q) - b*sin(q) - r*cos(q - t)   

6 _& ?: a7 I0 g, ^% s复杂点   就不化简了
>> syms t r a b k q A
>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
+ f+ h5 J* m" Q3 D, m) A6 V+ L>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
) X& q* \, [5 \8 D/ u+ w* r>> M=[x;y;A];
>> NM=N*M

6 g/ j# y1 Q( M' \3 ]3 u5 vNM =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
* ]8 N  r- X+ M5 W" N 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));
) Z8 @4 J( `7 M9 k, A9 G>> simplify(x1)
# \& [/ I* I' ~/ a) v' E$ |
% R. X  C- J$ v4 p5 qans =

  C( S# u4 G9 Q8 r0 j5 Q4 R! n7 s7 dA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
$ i" M/ L) d% N8 h8 q" w/ f
A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
' ~# p6 f# b4 _, t3 Z
' u  o2 Y" l' N6 {! s>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
& N! c4 F# i  T8 q' ]>> simplify(y1)

- [9 o6 X' h; Ians =
, d( o1 q" ]& Y" H8 N; U/ D
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
6 k5 k; L; ~5 Q: D" L: n

明月山河

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

shouce

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

$ |* ^5 ?) D& O+ w. X, C) o大侠   能举个例子吗   我这个不懂

明月山河

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

6 t* i; k* w% A9 G+ uhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
+ }' ^" o2 O8 f2 j" V& Z( p2 e
3.带有附加条件的化简
" F6 ?9 l9 q& e>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
% O) _, z" h8 G>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
- Q: N/ j& F4 u2 w/ i& f- W4.展开与组合
1 z! q$ Z( C! u. y  B; p$ \! d. f9 Vexpand &combine
" {4 K8 G5 t) \  vcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
$ X" \( j( Z9 i( U

( r- f% w  n- o7 }" W" Q3 x
" K2 z' O1 A5 z- O/ W$ e7 i- G; Xmatlab是用的maple的符号运算内核，使用规则一样。
8 B. q$ X' f2 v
! h2 k% r2 J% \( N1 @  M4 x

shouce

明月山河 发表于 2016-2-20 18:38
) H2 I! H; @% z9 a# h8 g6 ?http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
2 m( U1 I+ K" t5 Z$ I6 w
3.带有附加条件的化简

8 p' X9 C# a" c- c  r- `多谢了   等下去试试
& c0 {- |' k7 H3 U( g

shouce

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

- G8 q  O  r& R2 w2 u' R9 \, ^1 ?3.带有附加条件的化简

$ c0 T# [: Z3 p% \试了下还是不行   可能matlab和maple 使用不一样
>> syms x
6 B7 a8 g/ H# a, p>> combine(sin(x)^2+cos(x)^2)   

ans =
: T3 K: r" A! Q3 o
cos(x)^2 + sin(x)^2          这个不化简单
- q  \' L( `( U8 j: Q/ [0 m
; F3 w& |# A# \# Y/ o8 m>> simplify(sin(x)^2+cos(x)^2)
3 a4 _: c  L$ K. {3 T; I7 A
ans =
* h  t% G( D7 }( A6 W5 {
$ z9 ?+ `- _: g% [2 V! r1


>> syms x
0 L* e+ L0 m* }>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

4 H+ n3 K  I0 H2 s& v$ m( rA =
; T9 B2 a/ Z2 T7 L- _! Q
log(3*x) + cos(x)^2 + sin(x)^2

& Z- q, ]& ?4 G3 D: ^>> simplify(A)

8 t- L) W3 w, pans =

* {8 B2 j6 l6 K& v3 a$ x" [log(3*x) + 1

& y" Z) {" h* d& b" T6 p
1 h3 G, ^* L2 C
4 T' N) @# I' k- n2 @8 {) G>> simplify(A,x)
错误使用 sym/simplify (line 99)
Too many input arguments.
1 H& ^/ n- G/ z% I$ u
; C' e0 H, w8 }" i

393810298

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

shouce

>>  expand(x1)
; o6 Q/ C( X: Q9 H9 _ans =
  b. K; ]6 l# F7 `A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
: Z$ i2 a& }4 M- P>> 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)')
0 F+ N9 x4 M% c9 [KMM =
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)
0 _  ^$ ~) ]! S! j" s  J7 y: j- I" [ans =
# ?# j4 s' J. rA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

" c7 ]% F1 y( s, i要先展开然后合并
! o; X& M, R. R9 X* h' P