matlab 符号运算的化简问题

shouce

首先来个简单

>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
+ d) w9 h3 i6 i$ s# b) A1 k>> X=-x*cos(q)-y*sin(q);
! }8 k& u/ t7 Q1 O4 Q6 Y>> simplify(X)
6 B5 d* U6 a" X6 ?# _9 ~3 x3 B
ans =
5 S' m. q+ e( ^  W7 W6 `$ c
+ e$ b# q  H5 C$ u* c3 Y7 G- a*cos(q) - b*sin(q) - r*cos(q - t)   

+ M+ B2 n) t; C8 R5 k' M& u复杂点   就不化简了
5 F# _1 a7 Z* y >> syms t r a b k q A
+ F. r. p0 J1 G2 ?. Y>>  x=r*cos(t)+a;
! ~- ~; f$ S0 V7 A: C7 R5 }>> y=r*sin(t)+b;
>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
  [% B, A" F* H! ]+ B/ G% R>> M=[x;y;A];
>> NM=N*M
3 e9 Z/ u( m$ }' W8 ~
NM =
: V" |" G% b: r' I$ 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)
9 t3 D) c  X" ?( i
>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
0 G  m2 k6 ~& c8 x6 V! Z>> simplify(x1)

ans =
3 }) x/ ~; e+ @$ |9 B0 Q1 }4 ^
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
" C' p1 H. D0 E1 l$ d. l) t
A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
7 t( K9 h( U" c  e* }
. `2 W/ I3 @2 J6 P2 \' h>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

  V" _- g/ j0 [9 K% C* Fans =
% B5 ~3 [; V# V
' o4 z) m2 |2 V1 s; Tcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

, t! {; p( d8 |, K/ o* b

明月山河

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

shouce

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

% I" {% D" P0 D/ z8 X' ?大侠   能举个例子吗   我这个不懂

明月山河

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

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

3 I+ v) {" z3 A% G) _% r3.带有附加条件的化简
4 F2 ~6 ?- d* A>simplify(expr,{x*y=1});
& D* ?+ p4 u6 n  O>siderel:=x^2+y^2=1;
" i. G7 F: o* o  g>simplify(expr,{siderel},[y,x]);
0 o2 H/ P" J4 t# m#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
( P" T3 {+ s; e" T( N4.展开与组合
expand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起


$ X! g$ U. D- E, Q% t6 d1 R
' l: t+ \1 ~7 _3 Amatlab是用的maple的符号运算内核，使用规则一样。

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
3 @1 t) H9 m3 h+ i0 d& u  }$ E2 B
- i$ z- O8 ?6 J- F1 f3.带有附加条件的化简

. y' {+ Y1 D9 h" H% L多谢了   等下去试试
0 v6 {8 V- ^$ N$ I1 j

shouce

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

3.带有附加条件的化简

8 e" x* a# H% l试了下还是不行   可能matlab和maple 使用不一样
9 e: Y% \: I1 N>> syms x
>> combine(sin(x)^2+cos(x)^2)   

& q$ j5 b2 D( b, x$ zans =
) ^( G# X0 Z& Y$ k8 w2 x
) {& m9 s' r- [. J; X* j; y1 L, U* Rcos(x)^2 + sin(x)^2          这个不化简单
  i3 l7 ]9 q: ^" v( h6 M
: t9 x+ \* t& Z" N# q! F>> simplify(sin(x)^2+cos(x)^2)

ans =
! c% [7 O/ P" A1 v  f. F9 o
& w3 g+ D( K- P- t- P2 L& v1


  y6 _. h, D2 t8 D6 F>> syms x
: E+ i2 a6 c' k: B8 E$ j7 h( E# r>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
& ^' w( E' e1 ?& v. B
( i7 H' X8 `  w7 l$ r4 s7 TA =

( b8 A# P# J. h4 Z5 ?9 Y$ jlog(3*x) + cos(x)^2 + sin(x)^2
' H" e3 E& f4 E8 u! h2 w
/ u1 A2 R1 s8 Q>> simplify(A)
( ~+ m: O: @$ A* I3 o1 ]0 i
ans =
; O* p0 H5 G' o! h4 v; {! e7 |
log(3*x) + 1

1 R" Z7 l1 P- }( _

>> simplify(A,x)
2 Q* f2 y: |6 \1 h# A. m" ^% X错误使用 sym/simplify (line 99)
Too many input arguments.
! D% a- q& z  [
$ l& s' k# Q0 i# t1 p) Z3 g

393810298

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

shouce

>>  expand(x1)
ans =
  H: w5 @% b( ]0 s: S  h2 p/ \, k, X* nA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> 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)')
; e1 {8 V  T; i! @9 HKMM =
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 =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
4 \% ?" ~# _+ L5 \7 C
要先展开然后合并