matlab 符号运算的化简问题

shouce

首先来个简单

& w, u3 B6 A. l0 d* r' r0 J: U>> syms t r a b q
* M" y% u5 a/ V, O>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

ans =
. L- {. w* L. ?( b
! B( N8 E& k; P7 l F! t R8 g! r- a*cos(q) - b*sin(q) - r*cos(q - t)

2 \4 d5 n0 s' @2 u8 W6 t复杂点 就不化简了
>> syms t r a b k q A
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
+ }* d" b. p7 B- y# ?( L1 d>> N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
) F- H! ]+ }7 t3 d7 Z>> M=[x;y;A];
& P5 |; [) B( x! i5 z! N( i& L k>> NM=N*M
6 N% m1 p( D+ P! B; W- j& W
4 Q# L0 R- T% Q" B/ j N- B: vNM =
' a) ^& v& q! l: n. B
4 a; f4 e7 i5 E& A8 jA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
3 }5 Z( w2 M3 F6 q0 s8 V8 ccos(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));
( F$ p r. j3 F; z" J, y>> simplify(x1)
7 `2 D3 N* G% |9 I) X
+ p% `1 ~* l& L; \ans =
0 w* s+ j5 p1 n7 U# k6 I) @- [4 n
! u* {) q7 F3 D3 z9 W% BA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))化简无效果 按正常应该是
$ ~ m/ ~; P, r" V8 M4 U( K
3 a* ^. f2 v; Y6 l* A7 VA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

: X T- s' Z# j E>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)
; x2 F7 l7 b! _" G0 j }
9 B$ @* k4 i% L- Tans =
% D! S+ i: n/ c5 w* \( G% Y0 V9 x
5 q- e" a) @3 W: ~5 v. Dcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
8 {/ P; A4 U, z1 F/ a8 qsimplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

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

: J) _: n- d2 vhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

) p4 F2 F5 n9 G8 n3.带有附加条件的化简
4 _( A: _1 z" d8 w- q>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
$ _$ n& x+ j6 U! y7 k' J/ A9 T2 L>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
expand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
5 h: ?; @- f: W
( O5 \. ~( \2 [9 l( z/ T* |

matlab是用的maple的符号运算内核，使用规则一样。
+ K$ t; H' K& \; A5 A& C
3 p [; ~( L& z2 y4 d) x

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
l( J! z. r1 [: Q+ h
3.带有附加条件的化简

2 I# N! p' T0 G* }+ m多谢了 等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
" d5 `; z7 T3 J# hhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简

! s7 O B, N: B试了下还是不行 可能matlab和maple 使用不一样
>> syms x
* V$ A+ s2 [! j4 R>>combine(sin(x)^2+cos(x)^2)
% p `5 O1 r7 ^' e1 l0 T7 v
! m2 F( f2 `) n5 v6 vans =
, c: I& C. ]* K6 y) I
cos(x)^2 + sin(x)^2这个不化简单
. d k/ @* ~+ f. S
: a. D& w0 R: }# K. B>> simplify(sin(x)^2+cos(x)^2)
( X- f! y7 E9 I* o/ R( @1 n
w H/ O' s2 A# |ans =

2 ?2 ^5 R4 H+ U' ?3 b1
- e# K4 D- U$ W
& w' ?* K% {) _4 z( Y% G
8 ]" m# c* j, U0 [ i>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

" Z7 [% q4 R8 I2 [" x0 sA =

log(3*x) + cos(x)^2 + sin(x)^2

" y- H8 I# R; D- P0 X! g>> simplify(A)

ans =
- C* h* E1 ^$ Y3 t, b
log(3*x) + 1
, V1 J7 m% A1 |. e( L
$ Q# n4 Y6 z4 ?$ z( A( _) }

* p7 c% }* Q; b+ a>> simplify(A,x)
: O) P- q: P# M5 D错误使用 sym/simplify (line 99)
2 K+ i# T5 D' d8 D0 \Too many input arguments.
! t- ~) @0 Y: a" }/ F% ~) z' x2 ?

393810298

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

shouce

>> expand(x1)
% ]% }0 t2 e4 T3 n: t* k1 Ians =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
1 B, i* P7 }5 G>> 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)')
) g. K8 z F' V r; T5 Q5 PKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> simplify(KMM)
, _" g; s% m' Cans =
/ z' K0 z! `) ^1 F6 A/ GA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

' ?( `6 }' s5 ^8 a. `要先展开然后合并
9 \3 R& J3 @0 O* B- O