matlab 符號運(yùn)算的化簡問題

shouce

首先來個(gè)簡單

>> syms t r a b q
>> x=r*cos(t)+a;
: T/ E8 m* w/ J7 y# [>> y=r*sin(t)+b;
>> X=-x*cos(q)-y*sin(q);
3 W" o$ a4 a. e$ P9 x>> simplify(X)
* J* c% V: j% F+ M7 P& I8 ]. G
7 A% }: d* ~8 ^/ r" aans =

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

復(fù)雜點(diǎn)   就不化簡了
" q) L& F1 O! c' G >> syms t r a b k q A
>>  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];
  G# W/ s$ q. s7 p0 T>> NM=N*M

# g4 l" ]9 @+ [6 P* i4 n- zNM =
1 \4 p9 Q, p9 m# Q, p& Q% O
A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
, @% O! c0 n3 b- M: U# z( T2 K 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));
  n4 X! e: B* f>> simplify(x1)
- A8 }& P: o3 u, Z" h
6 i0 g+ q& g3 K: H' |. uans =
# q$ ?3 K8 E- t7 ^
! G( U% ~( l8 U2 W9 J+ ~A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化簡無效果    按正常應(yīng)該是  
$ _) z/ y0 B8 m  s% d6 c$ g
" J& t) l' ^6 r3 K( N. `A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)

+ T9 E' E: g9 d7 y* V6 g>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
>> simplify(y1)

ans =
0 e/ [+ c: M; o: h, T/ |
cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

明月山河

simplifty函數(shù)要有控制參數(shù)才會(huì)按照你想像的運(yùn)算。

shouce

明月山河 發(fā)表于 2016-2-20 18:27
simplifty函數(shù)要有控制參數(shù)才會(huì)按照你想像的運(yùn)算。

大俠   能舉個(gè)例子嗎   我這個(gè)不懂

明月山河

shouce 發(fā)表于 2016-2-20 18:33
大俠   能舉個(gè)例子嗎   我這個(gè)不懂

http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

  |2 y8 Y( n1 j4 R; p3.帶有附加條件的化簡
: C; Q) _% k7 G6 L>simplify(expr,{x*y=1});
/ M2 r, _$ B. u4 }% `; I+ u+ P9 _" t" l>siderel:=x^2+y^2=1;
# N5 q$ a9 G* K( ?) }' |* h>simplify(expr,{siderel},[y,x]);
+ _5 P# ^1 n+ v3 a  U5 a0 T# q3 d#maple在表達(dá)式中作代換y^2=1-x^2，然后再試圖對x^2做代換，由于找不到，就停止了
2 z+ m. }  j4 u" b% g- V4.展開與組合
expand &combine
5 U# T/ {/ Q8 s9 g( \& w. ccombine是與expand相反的命令，它可以按照數(shù)學(xué)規(guī)則將表達(dá)式中的某些項(xiàng)組合在一起
% {5 c6 E7 c4 _5 U! x) c

0 v; d% V: t& H2 S4 g% T; x
- K* R% |/ w* l& d3 s+ Tmatlab是用的maple的符號運(yùn)算內(nèi)核，使用規(guī)則一樣。
8 c' @$ e6 x7 D- ^8 j" I/ ]+ t

shouce

明月山河 發(fā)表于 2016-2-20 18:38
+ r/ S. Q8 O9 h6 _http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
- s8 b# T6 p. c9 B3 N  S! [/ d& [' |) k
3.帶有附加條件的化簡

多謝了   等下去試試

shouce

明月山河 發(fā)表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.帶有附加條件的化簡

0 f( x; g; l8 o( }" b  R( x試了下還是不行   可能matlab和maple 使用不一樣
>> syms x
>> combine(sin(x)^2+cos(x)^2)   

ans =
& D/ ]" g' k7 e5 ^  C+ m
cos(x)^2 + sin(x)^2          這個(gè)不化簡單

  ^  Z7 g" Y4 Q. Z0 Y9 Z>> simplify(sin(x)^2+cos(x)^2)

, f& a1 K5 K2 x8 Ians =
$ i. U2 [8 _# M! q7 a
5 v9 H, c+ |" F1
( X* c8 L$ \0 k& Y, x6 ^8 ?

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

A =
! G7 N, m2 B( K" w/ M% s
log(3*x) + cos(x)^2 + sin(x)^2

>> simplify(A)

7 H/ a' b- ]" G) N' P$ U; ians =

6 R. k$ e# w7 b6 t  K4 g- Dlog(3*x) + 1
9 ]: Q1 E4 L! {


/ y1 l+ S& H! n6 k>> simplify(A,x)
錯(cuò)誤使用 sym/simplify (line 99)
7 a& u6 g, z9 y) q) SToo many input arguments.

) G, U6 P6 t1 p9 `) S3 N0 `8 ^

393810298

matlab大學(xué)里學(xué)完了，我就在也沒有用過了，已經(jīng)看不懂了

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)
% H3 y. m1 W1 C# O>> 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)')
, }, x* `& a* T% I  V/ \+ H$ uKMM =
( Y: F1 _4 ]- Q% U; r/ VA*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
>> simplify(KMM)
ans =
A*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
5 D' J7 L$ S6 ~9 J
要先展開然后合并