matlab 符號運算的化簡問題

shouce

首先來個簡單
; {9 \# y' |: ^# W
>> syms t r a b q
+ I' j( d$ y/ c>> x=r*cos(t)+a;
>> y=r*sin(t)+b;
* D1 l  k% \7 e  z" D>> X=-x*cos(q)-y*sin(q);
0 F- J; W- J4 m% w: u9 O>> simplify(X)
0 E! A) @% b1 N; A! O
ans =
, O2 \7 P- p( Y* {1 \
: u% W' f5 h4 T! I) I! u- a*cos(q) - b*sin(q) - r*cos(q - t)   
5 [) I, I. F: m' l+ n/ `0 j
復(fù)雜點   就不化簡了
>> syms t r a b k q A
>>  x=r*cos(t)+a;
2 ]" O' }& G( V* p4 H: z: w>> y=r*sin(t)+b;
3 _2 ~  `3 Z! G+ J. w>>  N=[-cos(k*q) -sin(k*q) cos(q);-sin(k*q) cos(k*q) sin(q)];
, X  T5 J: U  i# a/ B6 x>> M=[x;y;A];
>> NM=N*M

" n2 |+ L' f3 a7 i2 q( |7 [NM =

$ ~5 _0 i! X; c) o' Z A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))
- D- K, o6 I  O- `, e cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)

( Q3 v. z5 r! {1 g>> x1=A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t));
>> simplify(x1)

1 L4 {' a2 I5 d) Y/ W% M6 j! |. _& Fans =

A*cos(q) - cos(k*q)*(a + r*cos(t)) - sin(k*q)*(b + r*sin(t))       化簡無效果    按正常應(yīng)該是  
$ z: n! T! q: \
A*cos(q) - a*cos(k*q) - b*sin(k*q)-r*cos(k*q-t)
' y3 Z# s! L) S3 l& }* _
>> y1=cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q);
- @4 Z" w1 I0 {>> simplify(y1)

" D3 u. f/ M' Y8 E$ zans =

cos(k*q)*(b + r*sin(t)) - sin(k*q)*(a + r*cos(t)) + A*sin(q)
/ f9 D' n$ T' @4 ^, g& x5 F4 ]* Y: z

明月山河

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

shouce

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

大俠   能舉個例子嗎   我這個不懂

明月山河

shouce 發(fā)表于 2016-2-20 18:33
. ^6 ^" O; E4 ^$ C# Y9 z3 }大俠   能舉個例子嗎   我這個不懂

! |' A  u! U7 V9 ~. _; r: fhttp://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

3.帶有附加條件的化簡
>simplify(expr,{x*y=1});
>siderel:=x^2+y^2=1;
7 |, |  r5 Z6 E6 o. e4 I! a>simplify(expr,{siderel},[y,x]);
#maple在表達(dá)式中作代換y^2=1-x^2，然后再試圖對x^2做代換,，由于找不到,，就停止了
& y3 [2 I& N' \4.展開與組合
( e4 o8 @" \9 j0 Z* Gexpand &combine
, {6 {# ]6 v; ^7 U, A- U/ {combine是與expand相反的命令，它可以按照數(shù)學(xué)規(guī)則將表達(dá)式中的某些項組合在一起
, t' G3 ~2 P; k8 W

5 y/ Y4 y* \: w5 o$ K9 A( n
) A/ \) w8 I7 k, c4 |matlab是用的maple的符號運算內(nèi)核,，使用規(guī)則一樣,。
6 u' N' @, T; D2 J) S" `4 i
& `& n* m, s6 W/ D) S1 K/ c

shouce

明月山河 發(fā)表于 2016-2-20 18:38
0 w  G+ B1 A) U7 t* |http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html

# b: s! u! Y) M4 g; Y3.帶有附加條件的化簡

多謝了   等下去試試

shouce

明月山河 發(fā)表于 2016-2-20 18:38
http://blog.sina.com.cn/s/blog_636a8b120100jnyt.html
  S8 F7 x- O' n' @
3.帶有附加條件的化簡

試了下還是不行   可能matlab和maple 使用不一樣
>> syms x
>> combine(sin(x)^2+cos(x)^2)   
! l6 ]- p. S$ l4 w# Y3 ~
ans =
" x, F7 Z4 @! J
1 [9 s; _7 Y* i2 r" s# kcos(x)^2 + sin(x)^2          這個不化簡單

5 F. i% a  X9 {>> simplify(sin(x)^2+cos(x)^2)

8 @6 s1 _+ g' x4 ]  F$ w& [ans =

1

3 j0 \; ?' R; i1 h9 Z  |+ b
+ S2 f" P9 ~* j' q$ M' V8 {8 @' j& n>> syms x
>> A=sym('ln(3*x)+sin(x)^2+cos(x)^2')

A =

log(3*x) + cos(x)^2 + sin(x)^2

0 ~2 Z; E8 T' F/ @( t1 P2 e6 ~+ X>> simplify(A)

ans =

3 \% x- Y  B! Q; s, B* \! Clog(3*x) + 1



>> simplify(A,x)
錯誤使用 sym/simplify (line 99)
Too many input arguments.
( n4 ]3 Q2 D1 J5 K$ G8 a8 l

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)
9 Q  _: S) C9 J! t! v# g% q8 P0 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)')
7 ~3 P4 K4 \; TKMM =
A*cos(q) - b*sin(k*q) - a*cos(k*q) - r*cos(k*q)*cos(t) - r*sin(k*q)*sin(t)
, x0 U% Z/ J. N' ^9 D% ?>> simplify(KMM)
ans =
) \' I! I2 e8 I* @, P9 T  MA*cos(q) - a*cos(k*q) - b*sin(k*q) - r*cos(t - k*q)
+ e4 F0 q  Z8 ]+ P% v
; D" ~# X+ ~, \' Z要先展開然后合并
; q% O2 }4 D7 N; R4 v9 g