matlab 符号运算的化简问题

shouce

首先来个简单

' N  H' k* N- Z+ S; w0 ?>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
+ v4 A# g8 S* H7 k* n  D1 z% L4 C>> simplify(X)

4 U  \. H5 ]# M5 Q6 qans =
4 c1 u2 A& e0 r$ q; O
9 N/ c! H+ ?3 `6 t, Y- a*cos(q) - b*sin(q) - r*cos(q - t)   
' t( v* M! a/ r8 K) F5 y7 L
6 T2 t* @  }7 ]复杂点   就不化简了
>> syms t r a b k q A
9 ~. |7 j3 y: `' f, H! r! S>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
+ L6 |/ Y) p% t; M>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
0 b: l$ |! ^) s0 c3 t5 W$ M>> M=[x;y;A];
>> NM=N*M
" i1 M5 j, P0 u' C, j
$ N' h: q6 V& kNM =

' z0 `1 t; G7 m2 ~% d/ g 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)
; H) s+ H9 v3 F: Z3 a' H5 n( V" s
ans =

8 m9 m! t7 L6 gA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  
7 q1 |' F/ ]' V6 |) C
; G( X6 z5 a- {7 b4 Q% k# NA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
! }" V$ f. p  j
>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
$ E) y0 q; K' O$ ^) R  ^+ s>> simplify(y1)
* W, |+ t& y6 ~0 @" D; [
* W; C! G0 K- S) S! h5 h. Aans =
9 |7 q( e' u% `) O
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
0 v& p) q" H) w. [& ^! {

明月山河

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

shouce

明月山河 发表于 2016-2-20 18:27
: O9 p  [0 p5 ^! Rsimplifty函数要有控制参数才会按照你想像的运算。

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

明月山河

shouce 发表于 2016-2-20 18:33
# w% a# d# |7 i, O大侠   能举个例子吗   我这个不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
8 T. T% ]# U9 p6 t
! l  `  u" ~! M% v3.带有附加条件的化简
>simplify(expr,{x*y=1});
, W" c- c0 A/ F+ v3 `6 R- U0 E- `: o  q>siderel:=x^2+y^2=1;
>simplify(expr,{siderel},[y,x]);
8 M& P9 k0 R% i: |! |0 W; O#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
9 N6 G# B3 }' {% i9 N4.展开与组合
! V6 C/ t0 `) W+ S  u" Pexpand &combine
5 r" V4 R4 C9 F$ dcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起


+ a3 q1 D0 m  K0 F8 Q5 P  _+ j
6 `2 }& R& x$ ]$ r+ I* `0 j* }6 I* \matlab是用的maple的符号运算内核，使用规则一样。

, I% [% ~, @& e8 d- I; e

shouce

明月山河 发表于 2016-2-20 18:38
3 M4 b1 p' C# ?) r5 L% \http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
1 Z0 Q; y2 i- Z  ]% W: y
3.带有附加条件的化简

多谢了   等下去试试

shouce

明月山河 发表于 2016-2-20 18:38
8 s1 C& |; Q$ n& [  G* V+ Ghttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
8 L2 w. S2 m9 o" v, i
4 C; ?' M5 W/ ?: d2 h0 {3.带有附加条件的化简

试了下还是不行   可能matlab和maple 使用不一样
>> syms x
3 l. Y, D- z  n  `; \1 L" x>> combine(sin(x)^2+cos(x)^2)   
% g; O9 i, I1 U+ l% h
/ J4 e( n& q$ w. Uans =
, J- m& r/ [7 G0 Z- a. T
3 K# n# _9 q: \# _' s  L& C9 ~cos(x)^2 + sin(x)^2          这个不化简单
9 V6 y% h; G# R6 P' @8 X
>> simplify(sin(x)^2+cos(x)^2)
5 v9 h4 F( l+ C& Q
2 Y8 ~. ^5 W+ x/ ians =

1
! ]/ Z* O; b% a( y" A* ]' {( n
8 D: A: N* d/ h7 E/ G* p( Y0 T
7 T" g& c# C0 g5 N1 ^+ s>> syms x
  {6 e! o) j) Y6 N>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')
" N2 ~( x' A, |7 M8 M9 A' K( B
A =

log(3*x) + cos(x)^2 + sin(x)^2
! f$ L7 w% ^, Z- D+ g% F/ w4 J6 r
, g! V- R5 J+ e) n; Q4 X8 [>> simplify(A)

ans =

log(3*x) + 1

0 Z7 p" c# a( |1 f: ?% J

! V3 [; Q$ S- L3 o  n>> simplify(A,x)
错误使用 sym/simplify (line 99)
1 i$ _% X& C0 ?: KToo many input arguments.

; O5 @6 p# Y  O; I3 A3 }

393810298

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

shouce

>>  expand(x1)
6 {% W  g5 Y+ x) r: `' jans =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
  G8 y5 \& k* y8 i4 U>> 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)')
( N$ |4 L6 c# @. v" aKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
# i7 j4 m. N- I4 C>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
- n+ W. a3 n  Y. @3 f0 h, J
2 {" t" R3 h( S9 x( u& H要先展开然后合并