matlab 符号运算的化简问题

shouce

首先来个简单
) A: b6 l; j( U$ d
>> syms t r a b q
& E1 A( i+ m3 F; ^3 E>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
% [  N" q8 U+ |. a>> X=-x*cos(q)-y*sin(q);
0 w: e! x1 Z' g% Q>> simplify(X)

; Y/ s) f. u4 i/ S- h; {ans =

- a*cos(q) - b*sin(q) - r*cos(q - t)   

复杂点   就不化简了
>> syms t r a b k q A
7 T( I- U$ _6 T. {>>  x=r*cos(t)+a;
% G0 f# ~! B$ q& u>> y=r*sin(t)+b;
>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
/ e1 @9 c0 d0 ]) A1 B>> M=[x;y;A];
/ F0 M6 f7 O. j. y) M8 H. A>> NM=N*M

& p6 r+ E6 v! @# Z& T% i7 \NM =
6 H) `( Z/ ^3 H' f0 \
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)
8 |# [" e  R1 c! S$ q
>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)

ans =

: P" U3 X; q. a* e( |) L; X3 p% c3 kA*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化简无效果    按正常应该是  

* [2 l. r5 r3 r2 IA*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
5 q3 _% p; N2 A3 t/ G
; G, t. t  a% I/ Z! p3 x. D0 b9 ^3 Z>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

ans =
; e3 ]* I: P: l& W6 P
0 d# v- K( z8 K: K* Zcos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

% j# G2 x1 d/ e& a: O7 V

明月山河

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

shouce

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

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

明月山河

shouce 发表于 2016-2-20 18:33
- M  ~3 E7 w9 {: L) [! M9 E大侠   能举个例子吗   我这个不懂

8 s' o- y- U! C+ w7 fhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
# u% |- Y- }, F* y
3.带有附加条件的化简
  U( G+ P& k$ Q# V>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
- z: o  E3 O0 g# a2 f) W- N>simplify(expr,{siderel},[y,x]);
$ H+ [6 F+ x+ v1 i$ i; p#maple在表达式中作代换y^2=1-x^2，然后再试图对x^2做代换，由于找不到，就停止了
4.展开与组合
0 ^3 C9 M9 i5 T2 B' r, q3 e0 aexpand &combine
5 P* V. L/ W  v! T" U: G$ w) `/ Lcombine是与expand相反的命令，它可以按照数学规则将表达式中的某些项组合在一起
! r( e) M' B& U* ~$ j% ^( [

: ]! X0 _0 k( }3 E5 n9 t% _
matlab是用的maple的符号运算内核，使用规则一样。

shouce

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

3.带有附加条件的化简

3 B3 U+ S. o  S# _" U8 I多谢了   等下去试试
# X% `: a1 B( X& J

shouce

明月山河 发表于 2016-2-20 18:38
  J  H- ]/ ^7 T  `: o0 P/ ]: Whttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.带有附加条件的化简

$ B  n6 }: _, i( z# b, b! S; y# Y' |试了下还是不行   可能matlab和maple 使用不一样
- w6 ^7 U$ D% H6 ]+ V, d5 e; j- y>> syms x
>> combine(sin(x)^2+cos(x)^2)   
: y/ J( u1 `% Y  z7 B3 @
1 R; _6 o+ v" ^6 d6 v1 hans =

) g" {9 _+ ~8 n/ t5 |8 Q6 Y8 m. @cos(x)^2 + sin(x)^2          这个不化简单

>> simplify(sin(x)^2+cos(x)^2)

ans =

* w9 w/ ]  A$ m2 k% q# Z* w* L8 I3 |1
3 h" e3 G5 z3 [7 }7 I; }3 J. a. y

( z0 J! }" \% |7 |>> syms x
6 \2 D* {$ v' k>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =
. F" {3 O5 b! F* L& l
log(3*x) + cos(x)^2 + sin(x)^2
- z$ _1 r8 Z% r% O
>> simplify(A)
: }0 \: \3 l) O( f; H  U& T
6 h' {2 e2 `+ s4 Vans =

log(3*x) + 1
/ e( G0 n# P( y! _2 Y# t

- v( N- I  `+ h5 W
! J/ L# V$ Z) A* ~>> simplify(A,x)
3 y# f: S8 V6 X9 ~# ^" Z! M错误使用 sym/simplify (line 99)
Too many input arguments.

393810298

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

shouce

>>  expand(x1)
" K& K, Y  g. n4 [ans =
" L9 \5 ~: l+ N$ PA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> 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 =
/ G2 B9 F" H, x7 k6 jA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
/ S( `" ]/ s( E3 O8 T>> simplify(KMM)
ans =
5 {# _% m& C4 L& b0 G7 M: U' z9 FA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
9 o+ Q3 J2 @/ K+ J8 h
$ _6 C6 M( Y1 M7 @5 t; B要先展开然后合并