matlab 符号运算的化简问题

shouce

首先来个简单
0 k6 h! ?5 s- E3 d$ E
2 T# h* Y* ^* t>> syms t r a b q
>> x=r*cos(t)+a;
6 \9 l" M" e1 g1 O' A, d>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

ans =
) W* A' {5 v, ]8 _, l
+ z( \& X) n; O. r$ S2 w6 I- a*cos(q) - b*sin(q) - r*cos(q - t)

复杂点 就不化简了
>> syms t r a b k q A
4 H8 x$ L1 _! X, x# w3 B, R>> 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)];
>> M=[x;y;A];
8 F5 }: x6 _2 V3 e>> NM=N*M
F' V+ x6 j6 h8 G& p9 k
8 H' Y- o3 V, d, _. \NM =

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));
* P3 d& i' h4 w4 n! p; w>> simplify(x1)

6 E. _1 x+ q0 m/ `8 Cans =

( H4 Q6 i' G7 y' E" \1 Y i4 B R" v$ RA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))化简无效果 按正常应该是
3 C1 [. `5 H9 B" e/ A
. v& M1 t& Z- ?( J' U9 W2 Q9 pA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

. P, @# k: }! @>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
" @6 A/ F6 }8 y0 J( u0 s8 B>> simplify(y1)
' ^1 }" J$ \* \7 R
0 ~5 S3 ^) z9 u& x: gans =
$ D* f2 r+ W( o% E
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

2 P$ T* x" U0 h: Z W3 v

明月山河

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

shouce

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

# B" X l }0 z3 U大侠 能举个例子吗 我这个不懂

明月山河

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

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
$ t; H6 \. C, n+ d# y2 l0 m
3.带有附加条件的化简
9 b' G$ y2 z* K' y- v5 |# Q>simplify(expr,{x*y=1});
2 m. m+ T. ?9 \9 O>siderel:=x^2+y^2=1;
3 ^2 R( ~9 R6 |>simplify(expr,{siderel},[y,x]);
# C v/ ?3 q/ Q#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
+ e2 R0 h! u! B$ }0 x7 B0 f ^! j4.展开与组合
, r6 z9 o4 I9 `5 |- L# Hexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起



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

( h$ c) A$ C5 J! a3 n" K

shouce

明月山河 发表于 2016-2-20 18:38
3 ^4 x" `9 p2 k5 qhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简

) x* \0 [3 ^% r8 x5 V多谢了 等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
" M$ a& B5 A0 N! L- Uhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
9 @1 [) X3 w; p$ B$ v
3.带有附加条件的化简

- S0 m+ ~! I* V m. }3 j9 [. Z. |试了下还是不行 可能matlab和maple 使用不一样
>> syms x
>>combine(sin(x)^2+cos(x)^2)
3 m& l$ a2 E" U2 T' l9 m) Z
' U7 r9 M& {+ @/ P+ Uans =

cos(x)^2 + sin(x)^2这个不化简单
! X/ x# B& {/ U% A7 p* I7 S# J
>> simplify(sin(x)^2+cos(x)^2)
4 B6 ]& E3 m9 t2 a4 y# b+ u
ans =

& Y/ P$ d. b1 s" X1


: O \; @' K9 B: @6 I; {, L>> syms x
; s3 U9 i0 c8 j9 T( X>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

7 s- }' W% R/ z n1 V3 p1 uA =

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

>> simplify(A)
6 J/ I5 @, _4 Y$ n
p1 J5 ?4 X1 k) {& b4 u5 Qans =

log(3*x) + 1
1 z( a3 W' }: D3 I- j
) V/ I& i& a7 _5 x8 k! r* [ i
9 ?9 g& x) {3 Z9 o0 K1 |9 a
>> simplify(A,x)
: p7 r. L' n, t6 T& ~0 Y错误使用 sym/simplify (line 99)
Too many input arguments.
( j+ \( D3 G$ F/ Z4 G
3 }' i2 ?9 @! N6 N: a. p

393810298

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

shouce

>> expand(x1)
ans =
4 F, _7 W8 F4 \5 y: e/ ~A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
3 d: N: s& G1 C# T* `) X4 }>> 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)
>> simplify(KMM)
* s* m8 A& ]) J, W% z! W9 N2 K% Pans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

+ h% I! T4 T% I5 p要先展开然后合并
- N( g% @; o% u o" e- g2 x7 r) Z