matlab 符号运算的化简问题

shouce

首先来个简单
1 }) d" C4 B5 q) S2 P8 P
' t+ B% `. J) b>> syms t r a b q
>> x=r*cos(t)+a;
6 F4 ~6 T, k, B" n( i; d  p>> y=r*sin(t)+b;
: e. T2 @  ]) g4 U  W7 f>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

! u, q. `4 b4 dans =

- a*cos(q) - b*sin(q) - r*cos(q - t)   

+ R0 v* b; ^+ y1 I复杂点   就不化简了
8 P/ N" s4 T; W  b! x >> syms t r a b k q A
>>  x=r*cos(t)+a;
, n/ }! U# Y5 e! Y8 k>> y=r*sin(t)+b;
  j  _$ N' G/ J! Z>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
4 a! b% ?9 V. z: C+ }/ K>> M=[x;y;A];
>> NM=N*M
% r& Y1 y5 I8 [! l/ P9 V+ t
! f2 H2 ]1 N" WNM =

$ m; ~# a' [4 g6 {) W* v A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
+ h% d# ?4 H& o) @% u cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
; C" J1 o# y' H9 u6 v# E5 C
& U3 S" [- B+ ^" e# a>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
8 N, E7 z/ x" P$ E9 p>> simplify(x1)

ans =
- \6 P6 W) x2 W
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  

A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

3 ^8 b9 x" D. D" p2 ?>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

: s: r" D' o6 w, _8 k) z1 r9 xans =

1 n6 D7 m1 ]. U* K) f! l" U* \7 ucos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

/ ]' d/ R& c7 ?. r" `1 f

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
' @9 U: n$ G& V; n9 }1 R8 @simplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

shouce 发表于 2016-2-20 18:33
" d. j' H* M9 T大侠   能举个例子吗   我这个不懂

! o5 E0 n4 s/ s. K3 q8 I1 Ohttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简
% [& L5 {( T/ I>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
3 [6 m+ L0 F& I3 g+ C4.展开与组合
expand &combine
1 u2 ]8 v+ C% ]6 p3 I0 V6 qcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起



matlab是用的maple的符号运算内核，使用规则一样。

6 i/ {/ K9 k$ L5 l: y6 }/ _

shouce

明月山河 发表于 2016-2-20 18:38
4 E- W" ?% D+ \6 N4 I; [http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
. B# t: S6 A* s: ?! |+ {! c
8 ^4 d8 h# j, l, G3 E0 P& B3.带有附加条件的化简

% Z  j/ X  j( b多谢了   等下去试试
) i, c+ j& K" _. }9 H) ?

shouce

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

3.带有附加条件的化简

! q5 S$ G9 V. b1 m& D+ z( k试了下还是不行   可能matlab和maple 使用不一样
) N( U- V( p. z' K- o3 g>> syms x
>> combine(sin(x)^2+cos(x)^2)   

ans =

$ W. K0 @3 O2 A: t/ L7 mcos(x)^2 + sin(x)^2          这个不化简单

3 c3 E9 ^6 d) `6 H>> simplify(sin(x)^2+cos(x)^2)

ans =

1
6 y+ c4 Z9 p* `0 v1 K+ _
6 S! R; i( R) k- o2 Z% o: ?
1 C! G8 ]5 B, M% |! [5 H>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
. b" m+ h/ t9 M/ |* P6 k6 V+ d( g0 @
) C% Q$ u, b& m* AA =

9 H2 p6 D9 h5 K+ a9 H8 c* K3 }log(3*x) + cos(x)^2 + sin(x)^2

4 n2 A' A+ x7 d! B- X7 }  W7 x6 |>> simplify(A)

' n0 S/ y4 O, ^9 P; `  fans =

log(3*x) + 1
8 [2 D+ [+ c9 }7 L9 O


/ J7 A* G2 S. A& l0 W1 A0 w/ v>> simplify(A,x)
错误使用 sym/simplify (line 99)
Too many input arguments.

; Q; U( t& {: x

393810298

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

shouce

>>  expand(x1)
ans =
A*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)')
/ |/ L3 R1 q' v2 \7 PKMM =
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)

% e2 C/ U5 O* S4 u9 e要先展开然后合并
" v  F& Q2 T( I1 B( \