matlab 符号运算的化简问题

shouce

首先来个简单
& D: p; z- a9 ]* X5 |' u2 z: ]
>> syms t r a b q
% r. _% T  ]( v9 b6 j# m>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
9 ~0 N) w, ]& U! B$ ^>> X=-x*cos(q)-y*sin(q);
>> simplify(X)
1 G7 }; z6 e/ }1 k2 x7 U
ans =

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

复杂点   就不化简了
; r1 d' A+ X0 c" i( S* j) N >> syms t r a b k q A
>>  x=r*cos(t)+a;
1 G7 c* T4 _  P' |: h( l. \>> y=r*sin(t)+b;
* J1 w' f' v7 n/ K>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
' ~1 d8 A# F' M9 G- t* F>> M=[x;y;A];
>> NM=N*M

9 j$ B# @# V# fNM =

" j/ E, W; q. X% a$ _+ \0 U 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)

>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
4 d' Z7 e) ~9 k/ A9 [# ]>> simplify(x1)

ans =
! s2 {7 q1 z9 B1 t" u6 k8 m
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
2 _0 r; _# z' Z9 w9 L
; Y) C& H- o1 |% w& PA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
) _- q$ M7 ^' Z* V3 l
, F- L2 B, b: j! Y- x4 ^, E: J6 Q7 Q>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)
! ], f% G! l2 ?) ?( e
ans =

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

; W1 j* I: [+ o+ O

明月山河

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

shouce

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

% f; `* C& V% D0 A' u& Q  U大侠   能举个例子吗   我这个不懂

明月山河

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

: n! X3 n7 j  ~1 l6 B& m0 n. ?http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
: z( g7 l  M- c( a
3.带有附加条件的化简
>simplify(expr,{x*y=1});
9 }2 N9 r+ W; I# W>siderel:=x^2+y^2=1;
' C# w/ Z, L8 @7 X0 B>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
expand &combine
' C8 r9 h9 o. U3 g3 Zcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起



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

9 F7 d* d, X0 E

shouce

明月山河 发表于 2016-2-20 18:38
9 v# N' j7 f& F" U4 W7 chttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简

多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
) i5 W( S$ G' vhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简

5 f3 s- c+ \9 t4 f' F! C试了下还是不行   可能matlab和maple 使用不一样
# R+ J: I) F# o>> syms x
>> combine(sin(x)^2+cos(x)^2)   

9 }! j8 J3 I9 W$ j' }% a' vans =
$ B7 c6 G: b* l, d  o* @
cos(x)^2 + sin(x)^2          这个不化简单
; h& }/ K( e1 D! j. h; k
>> simplify(sin(x)^2+cos(x)^2)

4 |6 H# l2 J  z. J* A  cans =
# n! I; G4 Q  {4 z3 p/ ?
1


>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
' ^2 X) O. x/ X, z
A =
, z# P# |- F; a/ b9 [
log(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)
" V- \5 A+ E, \* W# p
ans =
- {; [, J  n" |: w
log(3*x) + 1
: y6 @7 }! q* }- M4 ]

* }" p, c2 K4 ]4 E
1 h& C! a( {( T% I1 v" \>> simplify(A,x)
4 U" A1 N/ a" ?! Z) L2 }" g错误使用 sym/simplify (line 99)
Too many input arguments.
! _+ q0 ^. K' @1 u3 D

393810298

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

shouce

>>  expand(x1)
ans =
% h3 p. I% N. @, L, G% GA*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)
1 e0 Q* u1 y6 o( j- ]>> simplify(KMM)
9 N4 q7 V# X& O2 j$ Aans =
3 r, l* J9 k( w6 QA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并
: p6 ]  o, j7 Z; U2 m* l( K$ Y