matlab 符号运算的化简问题

shouce

首先来个简单

; ^- V- U* i0 e# n' y>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)
6 J/ V4 @. M& q* Y
9 Q+ z' U0 U7 Vans =
; W4 P/ Q( I' x3 H
. w% }4 O, w& [) w  y- a*cos(q) - b*sin(q) - r*cos(q - t)   

/ ^! g2 U- R5 s& {# N7 b# r复杂点   就不化简了
>> syms t r a b k q A
>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
* R" s$ Z; }& G, K8 [/ f# L* U>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
>> NM=N*M

; x" Y: U5 H! ^3 n5 ^" Z/ F3 ANM =

) b& }, q  X: q7 C, F' ?% y 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));
$ X- C5 R6 |& a5 a>> simplify(x1)

ans =
% x7 H$ v, h8 f. R5 p
: Z9 w! P7 q4 e. i2 HA*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)
+ h$ J( a6 K2 h/ Y& r- j" n. C6 h
) V4 ]* D& a0 U/ T>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
$ X. R/ o0 u) X$ M5 v>> simplify(y1)
2 [+ N# n/ y$ k3 g. Q3 V+ w$ U, L
8 N/ ~" I/ U1 D  \- Qans =

cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
# v% R: W& j, e5 h0 Y& a

明月山河

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

shouce

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

4 |% z' @' V9 j( o' P2 E大侠   能举个例子吗   我这个不懂

明月山河

shouce 发表于 2016-2-20 18:33
: p. h' G+ M! x( O7 t, d# d大侠   能举个例子吗   我这个不懂

% q8 q8 U( w3 I; Yhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
6 ?% B: V% h% [% x7 O* s
3.带有附加条件的化简
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
4 X4 i, j! V% N3 b#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
5 b+ ]3 W) n5 q; R4 @7 _$ o+ rexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起


& F! P# U9 U* x$ d, \$ t0 N
matlab是用的maple的符号运算内核，使用规则一样。
6 e8 o, I/ w& q0 `

shouce

明月山河 发表于 2016-2-20 18:38
7 Q( b! K& \1 g' d! \& {9 |http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
5 i3 i9 {& @5 T: q& n
3.带有附加条件的化简

多谢了   等下去试试
: B! T, P  C+ x0 s& u

shouce

明月山河 发表于 2016-2-20 18:38
& Y# i$ o* T. S2 T) nhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
) X9 ^1 P7 g, e0 l6 V. i
; N& [; l# C$ ~3.带有附加条件的化简

试了下还是不行   可能matlab和maple 使用不一样
>> syms x
>> combine(sin(x)^2+cos(x)^2)   

$ M* R5 k9 L/ T# d- Xans =
1 b" g1 l, S8 I. u( P
7 k  g4 \0 e3 ]+ h2 hcos(x)^2 + sin(x)^2          这个不化简单

) C( m  t; K; c7 [) t>> simplify(sin(x)^2+cos(x)^2)
( m. ?' ^. ?8 M' K" Y
* g# ?' w; ^# f' S* v) Q( v1 Z3 m* oans =
4 {/ p5 i8 j5 w) _, I- X; C
) @$ S# w5 ~( l1 l& A" i1
! t" J- }) F8 e( ]$ X2 O
& z3 V% |5 _& e
>> syms x
+ T$ Z& Z& }) P' a>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =
9 I6 m9 V4 i0 ]0 A
, l: a" S% p# r/ g, `& @+ X: Glog(3*x) + cos(x)^2 + sin(x)^2
+ a) o1 Z: N% `% m5 N! T
# o* c$ F+ M" r- s; Y5 x>> simplify(A)
1 q+ S2 l! u" R- n( d( d% P
ans =

log(3*x) + 1

/ K# ^* W7 W, C! L4 d7 g$ O/ I

5 W! T: e/ k, a/ X! G# G>> simplify(A,x)
3 g, \% H" X9 c5 {  e2 u错误使用 sym/simplify (line 99)
Too many input arguments.

# ]$ x+ Y! I3 J, ~

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)
/ P+ t( e2 g6 u, v+ T7 u/ c>> 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)')
) a6 T/ h6 |- u2 a8 [- k: JKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
& s" \. e7 U" ?# I) q; B$ R/ _>> simplify(KMM)
) j' [1 O$ c6 `+ M/ u4 Lans =
4 P, j. c% G+ OA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
% d7 V/ F, F: [# v
' J3 }  |0 {# I' [* }. P. i要先展开然后合并