matlab 符号运算的化简问题

shouce

首先来个简单
0 m1 f6 Y5 k6 L: G" }& I( q
% W) ]* ~; e$ i>> syms t r a b q
: ]# [. a7 n* r( a: W>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
9 I6 M: P7 w% M7 f+ {6 Z" V>> X=-x*cos(q)-y*sin(q);
  B6 `6 g# ?: D( q>> simplify(X)
% b' X. a! L& s2 V! O) p
ans =
4 O. a, D2 C7 Q  x9 n. J. ^2 a
- a*cos(q) - b*sin(q) - r*cos(q - t)   

* F3 j# ^! o& D/ k: s$ ~' ^: M复杂点   就不化简了
0 q- L7 i0 W0 s" F >> syms t r a b k q A
2 g# E" B- |) t: x: t>>  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];
>> NM=N*M
- a& O! l- z6 I) d" T, T5 o
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)
' {' t" `! H; x0 K$ I& U* \
>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)
  p* E* }  {  h- ?8 y
ans =
- |0 R5 c( B# q
! U, f( q: j3 q" jA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
) `5 i* r, c+ K. R
5 c0 d8 w6 c$ XA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
) h) P$ A; l+ m. ^0 o; ^
) S7 j  G0 B# |7 Q1 p# B>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)
2 j: J3 {2 V  a7 {
. l% E% u, }" k7 lans =
) s5 h: [. u% r6 V  u
* A) T: f0 u) W4 g0 kcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
1 n; g5 Z+ Z1 N6 R. U

明月山河

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

shouce

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

5 a7 q( L: |+ ^大侠   能举个例子吗   我这个不懂

明月山河

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

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
6 Q# S& w% \: {" ~
3.带有附加条件的化简
>simplify(expr,{x*y=1});
3 ?) x) S# J6 I" M6 ?>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
! a$ S: n. ~7 H6 _) s( A#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
' h0 F- ]* w1 m; }* H* v* ~4.展开与组合
expand &combine
combine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起



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

; d( S2 g& E  Q% S

shouce

明月山河 发表于 2016-2-20 18:38
3 \8 ^" `7 s7 T: j/ e, Zhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
' x5 q0 p0 E* |" E: C  o3 w# D7 T
3.带有附加条件的化简

多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
; `* ?% V1 \; g* N% V6 x+ k
3.带有附加条件的化简

% R. ?' ^. X" f9 y; ]. \, ?) g试了下还是不行   可能matlab和maple 使用不一样
>> syms x
>> combine(sin(x)^2+cos(x)^2)   
4 y9 R- Z0 k! F, O
) X4 Y% E1 s/ U1 A0 M8 [" B$ [0 qans =
  ?# O  J, y5 `& L9 Z" {
5 p! G8 j1 y, o, H$ J; Mcos(x)^2 + sin(x)^2          这个不化简单
, o- F; G/ R" ~  R
1 k- O. }7 e: z% e9 u0 o>> simplify(sin(x)^2+cos(x)^2)

# b1 |0 K2 g* w5 X$ K5 s# W. bans =

% W7 K* T2 @" x4 G1

# W4 D( E: |6 X/ p
>> syms x
% t% o5 u3 M. M  j7 ?0 R>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =
0 E$ B4 `& O: s: K
: R: I  K- k6 b: k& z3 b% Hlog(3*x) + cos(x)^2 + sin(x)^2
( L0 v% I) ?$ S; `
4 m9 ?. z2 L+ O- `4 p4 z>> simplify(A)
! x9 V0 k& u  a" H5 Y% F. F
+ H& J9 j% g7 |4 @7 c7 l, Aans =

* G/ r6 y# X6 A2 G: Plog(3*x) + 1

* y& O4 Z) G) P4 ]
4 e3 J( d8 i) g# ~
>> simplify(A,x)
错误使用 sym/simplify (line 99)
9 V5 l4 l1 j: R* f. B4 TToo many input arguments.
$ f& |& M9 u2 S& [& J) F

393810298

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

shouce

>>  expand(x1)
; a1 o# _1 j- p: T* {ans =
4 A: q0 O& i1 t6 X+ f; AA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
$ r! X. W* K& z5 L+ B% p>> 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)')
2 m5 x7 I# I( p( j& u  kKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
7 D% ?+ }( g# O' a# `) a) _>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
* }4 a; \! V, E' ]- ]
* j9 y" O1 X& X* |% |要先展开然后合并