matlab 符号运算的化简问题

shouce

首先来个简单
4 ~8 A( a% b9 |- Q+ E
4 U4 {; G7 d# s* Z) Z) M6 T>> syms t r a b q
0 w2 @4 i, `& k/ i, Y3 H- Y/ Y. @9 `2 }>> x=r*cos(t)+a;
% w# z; e' B6 K8 ~$ F% l>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)
* {: u3 x. t' Y' f# g: ~' J1 R
8 d" V" B& B+ f; Q; q1 e) p; Tans =

- a*cos(q) - b*sin(q) - r*cos(q - t)   
- ~1 Q% `) `% T" r0 V0 c5 b, _
复杂点   就不化简了
. C/ C4 N5 u* Z  c4 `: Z >> syms t r a b k q A
>>  x=r*cos(t)+a;
: P+ w, y6 W# Q# G( D" X* y5 i>> 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];
>> NM=N*M
6 Q) g5 n( P5 N( S4 H$ V8 ]
NM =
3 }5 R2 f; {% X' W8 t! J
1 i8 l3 j, Y% n4 P1 g! n 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));
>> simplify(x1)
  a7 i, C  E3 Y: I
% J0 b& T8 j1 _! Tans =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
0 T5 ~6 y. C( a* ?( [' |5 o
A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
5 @. T2 V) p+ y, g5 \$ q* w
>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
) [/ ^( ^/ h  R9 C>> simplify(y1)

, N# L( E  [; M% \& ?ans =
( \* C2 g* q$ q, e: x7 p
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
6 ]4 k8 |3 o7 E$ W. {$ h; M
$ P& g) T2 {. H& @8 B+ B4 u

明月山河

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

shouce

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

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

明月山河

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

# X2 w; ~  w6 E, |3 n! fhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
  r  h6 ~6 p8 D# s% X  H
3.带有附加条件的化简
; v$ }) F0 E/ b& t3 s. B0 i>simplify(expr,{x*y=1});
3 u2 Q5 A/ X9 r# f- ]  o- V3 y2 g9 a>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
9 |; o+ ]( v8 W#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
/ g* y2 L* ~3 U5 Q  {: Vexpand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
' L6 t3 A# H  `1 O5 l3 z
; S3 n0 Y# g$ r/ J/ H8 ]; x# V
3 M3 S5 n6 f- E7 g* K& y
matlab是用的maple的符号运算内核，使用规则一样。

shouce

明月山河 发表于 2016-2-20 18:38
1 q0 Z" d/ v# Z  l: ]http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

5 s! C4 t2 Z! i3 z6 y1 R6 c3.带有附加条件的化简

多谢了   等下去试试
" }- b2 R" D# ^

shouce

明月山河 发表于 2016-2-20 18:38
8 F/ v/ ~$ F, M& H4 s: f. ehttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
. M+ i3 S8 \5 I7 M- t$ p+ c' j
. b! B) s* ~) [9 e2 `( J3.带有附加条件的化简

试了下还是不行   可能matlab和maple 使用不一样
; U! p2 O9 Y9 Z6 m6 m6 T+ @; I>> syms x
' d" p; @( |6 Z0 |" c>> combine(sin(x)^2+cos(x)^2)   
0 C: _4 o+ I  H9 L- Z6 ?% @
- Y: p1 ]6 S' Gans =

' G$ H3 B1 L  G; y( }6 }  {cos(x)^2 + sin(x)^2          这个不化简单
  f3 l2 s1 j, J. r, R8 }3 }
>> simplify(sin(x)^2+cos(x)^2)
9 u$ ]! o& b6 e. ?" F
5 h+ y& l! z. b" dans =
6 |  R3 b0 Q: W, C
% T0 y3 X% H1 `2 k1

$ \1 t8 O+ }# {! _/ B
+ E7 `" W& ?# U* @; d>> syms x
) A6 Z  q. R. ?0 R! X7 f1 P>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =
- v) w' D/ t0 I: ?
log(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)
1 T# x2 Q# W7 f
! j6 f* a! r! C$ l0 |) j% x( p; ^: Cans =

log(3*x) + 1



>> simplify(A,x)
5 {2 l- j4 r1 m. a' u2 j# `, y: ~错误使用 sym/simplify (line 99)
" y5 E- m+ @3 IToo many input arguments.

$ _9 i* N6 y: r+ n$ q! S

393810298

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

shouce

>>  expand(x1)
ans =
8 w. w- E9 {) V1 T9 D3 i1 NA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
$ |8 E0 a3 M% u! i# s- m>> 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)')
6 k3 U8 Q7 J# y- }, b7 O. ]. B$ |KMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
0 I+ U* p/ \5 X: B. k! J>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

5 ?3 Q. u* O. G; y/ P- z要先展开然后合并
  q* P8 s- ]3 ~5 K7 r" s