matlab 符号运算的化简问题

shouce

首先来个简单

>> syms t r a b q
>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
( L. a. L* N, g+ X>> simplify(X)

9 F5 T! V) x- t5 U0 g" l6 {ans =

- a*cos(q) - b*sin(q) - r*cos(q - t)   
' p, D+ |$ i) A% H" L+ |: C) }
1 [; L0 c: W. o* b# }8 h& ?" w" _6 }复杂点   就不化简了
( F% U& m' b6 T2 Q" X >> syms t r a b k q A
2 F. n# M% P# e>>  x=r*cos(t)+a;
>> y=r*sin(t)+b;
" |% r7 j  J. q. u$ E9 p- W>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
3 f: L% ]$ c% ~- c! ]>> M=[x;y;A];
' r+ ^: P8 |1 y& {>> NM=N*M

NM =
% q. ^; t& g1 i5 f7 [! a$ ?
! C0 ^0 M4 F& u8 u4 X% ~ 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));
/ h; N0 m1 @. x8 f9 P* J>> simplify(x1)

: d/ [8 p' ?: R' B4 cans =
  Y6 O+ t* w6 f3 r* I
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  

/ F) A  g  B1 {  z3 P% x, r/ pA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

" M7 h! [% u' j. e0 o# z- H+ w. h5 Y>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
; o3 I2 T/ m9 v$ F) V) ~/ y>> simplify(y1)
# t4 R/ x# l) T8 b+ _
ans =

cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
( x$ _* S, [* U$ T( n: e5 w

shouce

>>  expand(x1)
% m: b0 P' \4 _$ D( q" D7 z6 rans =
0 c$ U0 P! D; I& ]1 |: g. C. B4 o5 bA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
% ~1 H5 g& ^  x' e" q; F0 `>> 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 =
2 e: \$ o2 J; ]A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
% i: J9 ?1 t% E- z>> simplify(KMM)
! G5 R5 u8 k  v0 q  J* g) Xans =
: N& X4 \/ ^6 u! AA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)

要先展开然后合并

393810298

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

shouce

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

3.带有附加条件的化简

; N4 S0 ]( H* y试了下还是不行   可能matlab和maple 使用不一样
& [/ C9 t9 j3 w# k( f+ A0 o8 p0 c>> syms x
>> combine(sin(x)^2+cos(x)^2)   

- c" [2 x' `* Dans =

cos(x)^2 + sin(x)^2          这个不化简单
0 X( u; n" |8 K2 I' _  G6 u6 b
3 r2 l' k: e5 q6 j+ u4 m>> simplify(sin(x)^2+cos(x)^2)

& J0 f  G( F1 z& Tans =
3 Q% m8 J. _5 I3 U8 R8 _7 P# n* _
  w6 g; t# l' d5 R! Y1
3 r9 A5 v+ ]" ?& _! z! W

>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

/ I8 f9 T& h) N2 CA =
" [: ]( r$ D9 B+ T. @  {
- p" y, }3 q' blog(3*x) + cos(x)^2 + sin(x)^2
: p, c6 `; [2 H9 W2 |
>> simplify(A)

ans =

4 v. h* ?1 H) nlog(3*x) + 1



1 l7 Z( X( d- e. ]9 M, \>> simplify(A,x)
错误使用 sym/simplify (line 99)
/ |2 F8 q& ?9 lToo many input arguments.

shouce

明月山河 发表于 2016-2-20 18:38
& H% `, f5 i' J4 Z6 x3 ?4 _) m1 Hhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
* b2 E' O$ O5 Q+ B9 R) `* D# v. l
& y- J1 F: F  t3.带有附加条件的化简

/ e! u0 _8 ^* @: B+ I# v/ Q# @. Y多谢了   等下去试试
9 t8 _$ ?6 w: V+ f

明月山河

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

' ~; n. y/ {2 \! ~http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
1 {6 |9 }: p; M) {
3.带有附加条件的化简
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
' J1 C; y! C. _  U4 x$ z3 |>simplify(expr,{siderel},[y,x]);
#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
5 n. x# _; m! b3 J4.展开与组合
" o3 J& Y/ K  i0 }* F/ a$ `expand &combine
6 c( P6 o# d8 v$ E5 b2 `' Mcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起


$ ~$ l/ j' C" B
matlab是用的maple的符号运算内核，使用规则一样。
; ^4 k) J: K; H+ l
0 ~" Q; S3 ]0 M) i1 ^4 ~, G

shouce

明月山河 发表于 2016-2-20 18:27
' p7 `1 K8 g+ c& B) p# _8 ]simplifty函数要有控制参数才会按照你想像的运算。

. t' W9 R$ x3 O3 K) V大侠   能举个例子吗   我这个不懂

明月山河

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