matlab 符号运算的化简问题

shouce

首先来个简单

>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
- s' M1 Q" t( ~) U) I>> simplify(X)
: x" L, g9 c+ l" W7 Q5 Q
+ V/ t( U0 P1 c: rans =
* L8 v7 ^5 y  y/ Z) g9 L
1 O" o. D" W" d7 l& x. x- u; F- f- a*cos(q) - b*sin(q) - r*cos(q - t)   
+ @  z& @" s) c5 A* b
; R" M  F2 h: s  K1 l( ]1 \复杂点   就不化简了
- Z( B! E5 P5 p+ |$ O$ _- g& J% F >> syms t r a b k q A
! j& E- X) Y# b: ^# S: b2 I" x8 w>>  x=r*cos(t)+a;
2 ?9 r' V3 e7 {  t! f# ~- F$ Z>> y=r*sin(t)+b;
3 e8 E; W, j+ ?3 Q* R) }! F>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
>> M=[x;y;A];
>> NM=N*M
& j% K+ A6 H& [- E, E1 u* r
NM =
" a: r( H* H( M6 E
9 ^2 J0 h' I: ?/ m) 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)
6 N4 `- w3 p, F4 v( w, |
8 @  @6 P* v( E' ]1 M$ s& H>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)
- M  I3 Y/ ]4 L0 k
$ [4 @( F8 j2 y$ D7 Q9 P- tans =
4 U5 B& \/ ?- H
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)

6 _. b, y7 ?# h5 _>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)
$ U$ U+ W* n9 d
, F/ I) q' w+ q+ m8 wans =
$ D( V2 q9 r+ y& a1 R/ ]' f
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
/ O; Y4 e4 K# a

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
3 F0 A, c, z/ q2 ]simplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

shouce 发表于 2016-2-20 18:33
/ o: n! I0 |; C4 ~# \$ ?+ e: |大侠   能举个例子吗   我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

! F7 h5 a1 w8 E/ I' @+ G2 b7 ^; r3.带有附加条件的化简
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
5 P) w' E% ?/ {& }& Q4 h; y* Y* y#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
' ^& P' z& |6 J! A* g4.展开与组合
# A' O3 |( s0 B8 g$ H( S' g, |expand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起



matlab是用的maple的符号运算内核，使用规则一样。
) f9 k8 B, ?0 K8 L5 W1 l$ c
5 i& N: K; W, W# l- F

shouce

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

3.带有附加条件的化简

多谢了   等下去试试
4 `: i' t* T3 z2 ]- J, y

shouce

明月山河 发表于 2016-2-20 18:38
2 D; V- O# G+ z  D- G& Fhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

0 a+ P6 q5 T8 F3.带有附加条件的化简

0 s! L; E& n( U/ C5 C6 K试了下还是不行   可能matlab和maple 使用不一样
>> syms x
  K. q  w2 `0 `* F! G& b>> combine(sin(x)^2+cos(x)^2)   
/ V# G# ]( e0 u1 s
ans =
$ k/ C1 m+ ~/ k; P" x4 B# R
cos(x)^2 + sin(x)^2          这个不化简单
: o9 W! Y' B- J5 J2 z# i
# a* i4 g) J: {8 p% a>> simplify(sin(x)^2+cos(x)^2)

: g6 ]$ Q; k4 _7 `ans =

; D# n" }/ A% J3 F9 b- ]. I' I1

' j! e8 n  e/ P
>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

1 {! Z5 i0 k7 ?- fA =
2 i+ U, o2 T" b: [
# m* R' ?8 ^& u3 Xlog(3*x) + cos(x)^2 + sin(x)^2

1 E! {3 g1 Q6 @% a>> simplify(A)
8 I9 b3 M6 {6 t
, \4 n2 D% z3 L0 ?9 V2 b* ^( Tans =

4 O8 U- u+ @7 ilog(3*x) + 1
" ^8 p- R% j/ Y% k7 D8 ~) ]% B- M" e
. T9 U, n! a6 y, @+ Y! v# \8 }

7 d* V0 a: x2 I  C>> simplify(A,x)
错误使用 sym/simplify (line 99)
Too many input arguments.
. V6 a% @& \/ p, u0 A

393810298

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

shouce

>>  expand(x1)
/ J$ B! a1 _* D+ Oans =
- q8 ~, \0 B5 I, d+ M9 GA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
0 f2 J4 y% j6 ^0 r>> 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)')
! d* W3 X2 I2 IKMM =
( L' {+ ^9 l! T/ p. E- G6 A9 p( VA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
" x9 Y& c& J4 ~; r) p/ Y>> simplify(KMM)
7 v) W  U, L, `. ians =
- L% y. e$ R, o, n3 e( bA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并