matlab 符号运算的化简问题

shouce

首先来个简单
# J$ z+ k: x; N* K$ r/ i
. A4 ~+ p! m% O: Z: ~2 }9 {/ j>> syms t r a b q
4 Y+ u% `$ b) }, `) q6 I>> x=r*cos(t)+a;
* M+ B- m& y7 w4 C7 m>> y=r*sin(t)+b;
$ c' B- k/ Q8 j>> X=-x*cos(q)-y*sin(q);
% |% [5 E  ?" y+ J& Q3 w: P, \>> simplify(X)

ans =
6 d" A+ a5 Q8 B, h3 E
- a*cos(q) - b*sin(q) - r*cos(q - t)   
' S+ ^9 d7 n) T. }
复杂点   就不化简了
# c5 U6 I% @$ W2 V >> syms t r a b k q A
+ B% V  u( b* n  [7 @' K- _" a6 P>>  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)];
; h! V( q' |6 Z>> M=[x;y;A];
- W2 K( M' u6 P4 V2 q>> NM=N*M

NM =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
* c  j4 G  e3 f# x) U8 p6 ^* D cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

; K1 w+ a7 b9 `; k+ Z2 v! U>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)

7 h6 E' J, T  v; ^6 N0 f* ]6 c! p; yans =
, m% i* L6 ^' K8 E0 ^5 w6 e8 A
. p2 ?% _3 f* ~; GA*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)
# D; E( E+ A% n- w
; T9 i9 k, h3 Q. K$ }# m- W+ w" i0 o>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
: w3 a& E) p5 m8 j' O>> simplify(y1)

! c/ f' y! K1 \1 [  W3 Pans =
; z/ A) q" ^* a8 z% k+ p* u
. v. _5 V  ^  Lcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
* z3 v2 p0 W# a) s6 \" Y- U3 tsimplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

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

, O' d) ?/ ]5 e1 I6 Khttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简
4 Y. t1 k& ?: i0 F% u>simplify(expr,{x*y=1});
1 |; Q  S9 C% _. d6 ~>siderel:=x^2+y^2=1;
7 t9 `2 W# w, |( z5 g6 s+ X>simplify(expr,{siderel},[y,x]);
5 b6 X6 E# j3 I% ]/ z5 B2 t#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
6 n! J% ?# U4 r& H# n( b' K; L5 a4.展开与组合
expand &combine
6 b8 V" f4 p, xcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起


. i) J$ p7 `5 f5 P
matlab是用的maple的符号运算内核，使用规则一样。
, t4 z7 l7 m4 a8 t+ Z
0 d, Q9 _% F8 w2 }" _- A& B" J+ Y

shouce

明月山河 发表于 2016-2-20 18:38
6 O' n7 v8 k- J4 Vhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

* }# x4 _' {. _2 d3.带有附加条件的化简

多谢了   等下去试试
8 g) V- y2 t4 K. o

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
# c: f/ C) G/ N2 v! u5 O7 }
2 p0 K/ I* J! @! e3.带有附加条件的化简

5 n# d6 j, n3 k! v$ F试了下还是不行   可能matlab和maple 使用不一样
8 X. ?+ t+ U/ B# R>> syms x
>> combine(sin(x)^2+cos(x)^2)   
, P9 H& w; b3 Q: R: T3 d5 _" M
# q( J$ X& R4 `! a) oans =
  Y$ [# d0 a: A1 v- ^
cos(x)^2 + sin(x)^2          这个不化简单

>> simplify(sin(x)^2+cos(x)^2)

ans =

1
* Q/ U- v& Q/ L3 d( J' X
. i! d+ S% @* }5 ?2 {  ^8 C  V
>> syms x
" G8 e! t/ Y3 t/ x( d- F* ]3 y>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =
" a  `5 u' D4 D1 T# A/ q# n1 G
) _) ]& p# _; b6 Llog(3*x) + cos(x)^2 + sin(x)^2
( r5 k& m/ r  C
! R# \( U$ [- G6 p6 w>> simplify(A)

ans =
" Z) `# |- `0 c" T$ [
( d: x: A% _8 j0 ~log(3*x) + 1



>> simplify(A,x)
  R3 u9 V+ v& D错误使用 sym/simplify (line 99)
Too many input arguments.
8 O) I. f* W( u% w) ]: ?

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)')
" U2 h! o! W/ l4 ]9 s) aKMM =
$ C  U/ L, L! d5 a" wA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
% R9 G) {# B( u  s  _' C' f0 P>> simplify(KMM)
ans =
  y* F8 R# Q  u$ r( B' I. M! kA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

5 X  ~9 y* z- }. K9 |" S要先展开然后合并
- o( ~( w9 f0 ^1 D, T