matlab 符号运算的化简问题

shouce

首先来个简单

0 S1 o% z5 J- l8 {- \8 w) v>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
+ G$ R. d) q# B>> simplify(X)
. r7 c9 Z$ V6 K! f2 b. Z/ x. U
( m% J$ a2 K4 d2 [" Eans =

/ C4 h# L9 P; ]9 w0 L& p- a*cos(q) - b*sin(q) - r*cos(q - t)   

复杂点   就不化简了
>> syms t r a b k q A
) ^: V/ D* K# k; d( n>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
6 F/ h# `' H" ?2 ~>> NM=N*M

0 i& O0 ]/ Z0 G. l8 k3 ]- UNM =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
# W7 b* |8 q6 o9 z8 M cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
+ C7 [  @) @. ]9 ]" ^; C7 B. t
' P* e' R: l8 w>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
- e$ G8 o3 Y! {: u) c7 B$ k>> simplify(x1)

7 }; F) ^! T3 {0 Uans =

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)

>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

  o4 @4 @3 M, a2 r& {! L5 Aans =
' K6 B( F  ^, ]8 w/ C
4 A' ~$ R! [6 q9 Bcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
+ s) r: F+ c- r8 |% a/ G& K! N% p

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
2 e# o- f# Z, h2 nsimplifty函数要有控制参数才会按照你想像的运算。

( F9 D& k# j1 }$ g大侠   能举个例子吗   我这个不懂

明月山河

shouce 发表于 2016-2-20 18:33
8 i+ _# Y5 E9 R$ s% Z7 ^+ E大侠   能举个例子吗   我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
2 q: r- t, m: Y0 k6 W( [% L) R
3.带有附加条件的化简
1 h" q! o$ u4 o/ B' e7 x8 u>simplify(expr,{x*y=1});
2 X+ D+ H6 A7 g1 u: D6 I) V>siderel:=x^2+y^2=1;
5 P1 K) {  r, u! N9 n$ ~; e>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
expand &combine
$ `& u" T/ I8 U) l) a1 Ycombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起

$ ]' S. Q) D$ x5 h% L0 r
4 l1 i; r* T! C( |
- g* T% h% x4 w9 m0 G6 E: O% S* Fmatlab是用的maple的符号运算内核，使用规则一样。
7 q6 c- F, {2 k! U
( D! j# R4 V* `! c2 Z. N% j7 ~3 G# u

shouce

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

' e, z1 ]6 n. P4 y( B& _6 K3.带有附加条件的化简

多谢了   等下去试试
/ C9 K' H6 ~, B& H; A( o

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
6 d& f  }- z4 r. Z! m( d# X1 v
& c& v! ^3 c& F; F2 F9 b9 z; `& q3.带有附加条件的化简

试了下还是不行   可能matlab和maple 使用不一样
8 O* p  `4 k+ N& I6 Y: k' U+ V' e>> syms x
>> combine(sin(x)^2+cos(x)^2)   

ans =

cos(x)^2 + sin(x)^2          这个不化简单
' C( ?, h: H/ O3 K) ?
/ p9 M! w4 N6 y>> simplify(sin(x)^2+cos(x)^2)
) D* t, [" A( F' p, A' y1 |
ans =
: j+ F( f" p' G4 `1 u, m
4 s1 u" E& r' b0 d' Z) @1

' i2 I. R7 X! h, O) M
>> syms x
. Z& O9 H, d( ~>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
  f# l( P, \% h: L5 U' t' W
% ?5 S2 ^4 ~8 ?) `A =

& V: s+ \- l8 T' l  b, Slog(3*x) + cos(x)^2 + sin(x)^2
. O+ h( r/ {- l3 h2 K
>> simplify(A)

- i# j1 ~: b/ t/ {9 U2 }* ~ans =

log(3*x) + 1
1 h2 B1 l: _0 Q$ ?


>> simplify(A,x)
2 @+ k1 @! j0 @错误使用 sym/simplify (line 99)
Too many input arguments.

" V6 ~. J9 z) g

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)')
KMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
$ ^/ k2 Q/ @1 m" d/ r0 P  H: V>> simplify(KMM)
3 W2 a4 C% g. |! Gans =
/ e7 J* d4 [5 ?8 pA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并
! u# B+ u% W2 Y7 ]1 M  i