matlab 符号运算的化简问题

shouce

首先来个简单
5 K1 ^& `6 d& v- }/ B( o  P
9 F3 T7 _/ `( E" \3 v>> syms t r a b q
+ ?  y7 G# C8 \) j! T+ {: C>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
>> simplify(X)

4 ]7 x  O' L- C9 wans =

4 F/ l8 d! D' \' c/ C- a*cos(q) - b*sin(q) - r*cos(q - t)   
! ?' B; ^# n4 b+ [+ G1 Z
复杂点   就不化简了
* |6 f0 o4 f. {$ T, h >> syms t r a b k q A
! T  m* v! T+ a4 w) a. v>>  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];
$ I' @+ Q" A% {' J>> NM=N*M

NM =
+ f& L4 X( i8 _
: g; h! b: T2 Z& l! q0 {- L 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));
. ]: S/ a3 Q# W! _1 p( w5 P4 ]( U>> simplify(x1)

5 d3 H3 z; l0 b2 {ans =

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)
7 m7 l' Q; h# X( b* b) V- a
; w+ G& G- l9 v7 B2 O>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
+ V- P- Q- W7 [; C>> simplify(y1)

ans =
6 W  W, a" N. Q) l( T! o
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
5 x7 [: R9 J- E4 {# N

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
  e* T' ]. y- k0 F$ [simplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

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

" s2 F, t! O2 F- yhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
7 l; `. s3 w- l5 U>simplify(expr,{siderel},[y,x]);
/ y, Z  ^* s, Q# o8 Q$ A#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
3 Q' G( r+ {: q8 u: Q  B  n4.展开与组合
expand &combine
9 B3 G% E& x: P  A+ [1 B  G$ \! ycombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
! D2 H/ i( d& m0 k


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

1 k5 B; H* M' ?& B5 K1 j" P

shouce

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

, |4 Q' G; A- x$ K% I- W5 b7 |3.带有附加条件的化简

多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
9 r* v' Y* W3 \http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

2 r" V) {; ^5 V3.带有附加条件的化简

1 N% X( L8 W/ J试了下还是不行   可能matlab和maple 使用不一样
>> syms x
/ t% @7 f+ B- z( V0 b5 e: B>> combine(sin(x)^2+cos(x)^2)   

2 i& ^* _6 d# {5 @+ h2 k7 Q" Oans =
: q0 }& \2 c7 b
cos(x)^2 + sin(x)^2          这个不化简单
: h+ L3 i2 o# g! D! [5 p
>> simplify(sin(x)^2+cos(x)^2)
. ^$ ?* v9 E6 z1 K& s# M
5 T/ U4 W+ v& h" H* k) c8 wans =
+ ]7 s9 R0 B! v' q$ c& M+ S2 ~' g. }
1
' E! K+ I6 u# U1 W6 V0 s! J! w  w% T
# v* @! ?: n7 G) `
>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
- G( o) U/ j2 S+ y1 V- W
4 {2 @1 t4 d$ A5 R2 E' AA =
4 }9 g2 W- H4 W' p" k
log(3*x) + cos(x)^2 + sin(x)^2
0 C" m/ w0 [/ a8 M3 a8 N
>> simplify(A)

0 b! e2 r! d  U* D2 hans =
& a3 J: a' _4 ?) i' I' h
log(3*x) + 1
! W5 {8 W  B4 v+ x  g  v6 M


>> simplify(A,x)
3 D  l! ^3 u+ O* `# q( X错误使用 sym/simplify (line 99)
Too many input arguments.
0 ~) A2 x; ~  Q7 A' v
/ L, a. r( f" h/ h

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)
( S% {  L0 w  o, _1 ]>> 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)')
5 F- ^4 x8 V. t6 _  _KMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
+ Y) P9 J2 {# M1 q>> simplify(KMM)
; K7 G* Q& A& q* ~2 d' m% Wans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
6 Y- `1 U. f, Z! X" i# i. X
: `- X: m6 T  L! T( `- L: F要先展开然后合并