第5章Matlab符號(hào)運(yùn)算

1、1,第5章 MATLAB符號(hào)運(yùn)算, matlab 不僅具有數(shù)值運(yùn)算功能，還開(kāi)發(fā)了在matlab環(huán)境下實(shí)現(xiàn)符號(hào)計(jì)算的工具包Symbolic Math Toolbox,主講：劉忠偉,2/60,符號(hào)計(jì)算的特點(diǎn)：,一、運(yùn)算以推理解析的方式進(jìn)行，因此不受計(jì)算誤差積累問(wèn)題困擾； 二、符號(hào)計(jì)算，或給出完全正確的封閉解，或給出任意精度的數(shù)值解（當(dāng)封閉解不存在時(shí)）； 三、符號(hào)計(jì)算指令的調(diào)用比較簡(jiǎn)單，經(jīng)典教科書(shū)公式相近； 四、計(jì)算所需時(shí)間較長(zhǎng)，有時(shí)難以忍受。,3/60,符號(hào)運(yùn)算的功能,符號(hào)表達(dá)式、符號(hào)矩陣的創(chuàng)建 符號(hào)線(xiàn)性代數(shù) 因式分解、展開(kāi)和簡(jiǎn)化 符號(hào)代數(shù)方程求解 符號(hào)微積分 符號(hào)微分方程,4/60,一、符號(hào)運(yùn)算

2、的基本操作,什么是符號(hào)運(yùn)算 與數(shù)值運(yùn)算的區(qū)別 數(shù)值運(yùn)算中必須先對(duì)變量賦值，然后才能參與運(yùn)算。 符號(hào)運(yùn)算無(wú)須事先對(duì)獨(dú)立變量賦值，運(yùn)算結(jié)果以標(biāo)準(zhǔn)的符號(hào)形式表達(dá)。,5/60,2. 符號(hào)常量,sym(A)或 sym(A，flag) flag為可選參數(shù)，有r、d、e、f r用有理數(shù)格式表達(dá)符號(hào)量（如分式，指數(shù)式，開(kāi)方式） 默認(rèn)形式 d十進(jìn)制數(shù)格式（32位） 將數(shù)值量轉(zhuǎn)換成符號(hào)量,6/60,f1=sym(log(2) f1 = log(2) %符 號(hào)常量 f2=sym(3*4-2)/5+1) %表達(dá)式形式的符 %號(hào)常量 f2 = (3*4-2)/5+1,7/60,3、符號(hào)變量,（1）使用函數(shù)sym() 格

3、式 sym(x)等 舉例： (2)使用命令syms 格式 syms arg1 arg2 舉例：syms a b c (可以定義多個(gè)符號(hào)變量，中間用空隔分隔變量）,8/60,4、符號(hào)表達(dá)式,符號(hào)表達(dá)式：由符號(hào)對(duì)象參與運(yùn)算的表達(dá)式 syms x y z r s x2+2*x+1,9/60,5、符號(hào)矩陣,符號(hào)矩陣的元素可以是符號(hào)常量、符號(hào)變量和符號(hào)表達(dá)式。 可用函數(shù)sym直接創(chuàng)建,10/60,符號(hào)矩陣的創(chuàng)建 數(shù)值矩陣A=1,2;3,4 A=a,b;c,d 不識(shí)別 用matlab函數(shù)sym創(chuàng)建矩陣（symbolic 的縮寫(xiě)) 命令格式：A=sym( ) 符號(hào)矩陣內(nèi)容同數(shù)值矩陣 需用sym指令定義 需用

4、 標(biāo)識(shí),11/60,例如：A = sym(a , 2*b ; 3*a , 0) A = a, 2*b 3*a, 0 這就完成了一個(gè)符號(hào)矩陣的創(chuàng)建。 注意：符號(hào)矩陣的每一行的兩端都有方括號(hào)，這是與 matlab數(shù)值矩陣的一個(gè)重要區(qū)別。,12/60,將數(shù)值矩陣轉(zhuǎn)化為符號(hào)矩陣 函數(shù)調(diào)用格式：sym(A) A=1/3,2.5;1/0.7,2/5 A = 0.3333 2.5000 1.4286 0.4000 sym(A) ans = 1/3, 5/2 10/7, 2/5,符號(hào)矩陣與數(shù)值矩陣的轉(zhuǎn)換,13/60,將符號(hào)矩陣轉(zhuǎn)化為數(shù)值矩陣 函數(shù)調(diào)用格式： numeric(A) A = 1/3, 5/2 10

5、/7, 2/5 numeric(A) ans = 0.3333 2.5000 1.4286 0.4000,14/60,用函數(shù)size( )可得到符號(hào)矩陣的大小,注：函數(shù)返回?cái)?shù)值或向量，而不是符號(hào)表達(dá)式。,A=sin(1) cos(2);s=size(A)s= 1 2,15/60,二、符號(hào)運(yùn)算,加 減 乘 除,16/60,例1：f= 2*x2+3*x-5; g= x2+x-7; syms x f=2*x2+3*x-5; g= x2+x-7; h=f+g h = 3*x2+4*x-12 例2：f=cos(x);g= sin(2*x); syms x f=cos(x);g=sin(2*x); f/g

6、+f*g ans = cos(x)/sin(x)+cos(x)*sin(x),17/60, syms a b c d e f g h; A=a b;c d;B=e f;g h; A+B ans = a+e, b+f c+g, d+h,18/60,A*B A/B A.*B,syms a b c d e f g h; A=a b;c d;B=e f;g h; C1=A.*B C1 = a*e, b*f c*g, d*h,19/60,syms a b c d e f g h; A=a b;c d;B=e f;g h; C2=A*B C2 = a*e+b*g, a*f+b*h c*e+d*g, c*f

7、+d*h,20/60, syms x fx gx; fx=3*x5-x4+2*x3+x2+3; gx=1/3*x3+x2-3*x-1; fx*gx ans = (3*x5-x4+2*x3+x2+3)*(1/3*x3+x2-3*x-1),21/60,fx/gx ans = (3*x5-x4+2*x3+x2+3)/(1/3*x3+x2-3*x-1),22/60,三、獨(dú)立變量與表達(dá)式化簡(jiǎn),findsym 查找符號(hào)表達(dá)式中或矩陣中哪一個(gè)變量是獨(dú)立變量。,23/60,三、獨(dú)立變量與表達(dá)式化簡(jiǎn),【例】對(duì)獨(dú)立自由符號(hào)變量的自動(dòng)辨認(rèn)。 syms a b x X Y;k=sym(3);z=sym(c*sqrt(

8、delta)+y*sin(theta); EXPR=a*z*X+(b*x2+k)*Y; findsym(EXPR) ans = X, Y, a, b, c, delta, theta, x, y findsym(EXPR,1) ans = x findsym(EXPR,2) ans = x,y,24/60,【例】findsym確定自由變量是對(duì)整個(gè)矩陣進(jìn)行的。 syms a b t u v x y;A=a+b*x,sin(t)+u;x*exp(-t),log(y)+v findsym(A,1) A = a+b*x, sin(t)+u x*exp(-t), log(y)+v ans = x,25/

9、60,合并同類(lèi)項(xiàng)collect 多項(xiàng)式展開(kāi)expand 因式分解factor 不定化簡(jiǎn)simple 通分numden 書(shū)寫(xiě)格式美化pretty,表 達(dá) 式 化 簡(jiǎn),26/60,【例】按不同的方式合并同冪項(xiàng)。 EXPR=sym(x2+x*exp(-t)+1)*(x+exp(-t); expr1=collect(EXPR) expr2=collect(EXPR,exp(-t) expr1 = x3+2*exp(-t)*x2+(1+exp(-t)2)*x+exp(-t) expr2 = x*exp(-t)2+(2*x2+1)*exp(-t)+(x2+1)*x,合并同類(lèi)項(xiàng)collect,27/60,多

10、項(xiàng)式展開(kāi)（expand）,【例】 syms x y a b c t E1=expand(x-2)*(x-4)*(y-t) E1 = x2*y-x2*t-6*x*y+6*x*t+8*y-8*t,28/60,【例】factor指令的使用 syms a x; f1=x4-5*x3+5*x2+5*x-6; factor(f1) ans = (x-1)*(x-2)*(x-3)*(x+1) f2=x2-a2;factor(f2) ans = (x-a)*(x+a),29/60,【例】符號(hào)表達(dá)式的化簡(jiǎn),MATLAB提供的對(duì)符號(hào)表達(dá)式化簡(jiǎn)的函數(shù)有： simplify(s)：應(yīng)用函數(shù)規(guī)則對(duì)s進(jìn)行化簡(jiǎn)。 simp

11、le(s)：調(diào)用MATLAB的其他函數(shù)對(duì)表達(dá)式進(jìn)行綜合化簡(jiǎn)，并顯示化簡(jiǎn)過(guò)程。,30/60,diff(f) 對(duì)缺省變量求微分 diff(f,v) 對(duì)指定變量v求微分 diff(f,v,n) 對(duì)指定變量v求n階微分 int(f) 對(duì)f表達(dá)式的缺省變量求積分 int(f,v) 對(duì)f表達(dá)式的v變量求積分 int(f,v,a,b) 對(duì)f表達(dá)式的v變量在（a,b) 區(qū)間求定積分,四. 符號(hào)微積分與積分變換,31/60,四. 符號(hào)微積分與積分變換,1 符號(hào)極限 limit函數(shù)的調(diào)用格式為： (1) limit(f,x,a)：求符號(hào)函數(shù)f(x)的極限值。即計(jì)算當(dāng)變量x趨近于常數(shù)a時(shí)，f(x)函數(shù)的極限值。 (

12、2) limit(f,a)：求符號(hào)函數(shù)f(x)的極限值。由于沒(méi)有指定符號(hào)函數(shù)f(x)的自變量，則使用該格式時(shí)，符號(hào)函數(shù)f(x)的變量為函數(shù)findsym(f)確定的默認(rèn)自變量，即變量x趨近于a。,32/60,四. 符號(hào)微積分與積分變換,(3) limit(f)：求符號(hào)函數(shù)f(x)的極限值。符號(hào)函數(shù)f(x)的變量為函數(shù)findsym(f)確定的默認(rèn)變量；沒(méi)有指定變量的目標(biāo)值時(shí)，系統(tǒng)默認(rèn)變量趨近于0，即a=0的情況。 (4) limit(f,x,a,right)：求符號(hào)函數(shù)f的極限值。right表示變量x從右邊趨近于a。 (5) limit(f,x,a,left)：求符號(hào)函數(shù)f的極限值。left表

13、示變量x從左邊趨近于a。,33/60,四. 符號(hào)微積分與積分變換,2 符號(hào)微分（導(dǎo)數(shù)） diff函數(shù)用于對(duì)符號(hào)表達(dá)式求導(dǎo)數(shù)。該函數(shù)的一般調(diào)用格式為： diff(s)：沒(méi)有指定變量和導(dǎo)數(shù)階數(shù)，則系統(tǒng)按findsym函數(shù)指示的默認(rèn)變量對(duì)符號(hào)表達(dá)式s求一階導(dǎo)數(shù)。 diff(s,v)：以v為自變量，對(duì)符號(hào)表達(dá)式s求一階導(dǎo)數(shù)。 diff(s,n)：按findsym函數(shù)指示的默認(rèn)變量對(duì)符號(hào)表達(dá)式s求n階導(dǎo)數(shù)，n為正整數(shù)。 diff(s,v,n)：以v為自變量，對(duì)符號(hào)表達(dá)式s求n階導(dǎo)數(shù)。,34/60, s1=sym(6*x3-4*x2+b*x-5); s2=sym(sin(a); diff(s1) ans

14、= 18*x2-8*x+b diff(s2) ans = cos(a) diff(s1,b) ans = X diff(s1,2) ans = 36*x-8,35/60,四. 符號(hào)微積分與積分變換,3 符號(hào)積分 符號(hào)積分由函數(shù)int來(lái)實(shí)現(xiàn)。該函數(shù)的一般調(diào)用格式為： int(s)：沒(méi)有指定積分變量和積分階數(shù)時(shí)，系統(tǒng)按findsym函數(shù)指示的默認(rèn)變量對(duì)被積函數(shù)或符號(hào)表達(dá)式s求不定積分。 int(s,v)：以v為自變量，對(duì)被積函數(shù)或符號(hào)表達(dá)式s求不定積分。 int(s,v,a,b)：求定積分運(yùn)算。a,b分別表示定積分的下限和上限。該函數(shù)求被積函數(shù)在區(qū)間a,b上的定積分。a和b可以是兩個(gè)具體的數(shù)，也可

15、以是一個(gè)符號(hào)表達(dá)式，還可以是無(wú)窮(inf)。當(dāng)函數(shù)f關(guān)于變量x在閉區(qū)間a,b上可積時(shí)，函數(shù)返回一個(gè)定積分結(jié)果。當(dāng)a,b中有一個(gè)是inf時(shí)，函數(shù)返回一個(gè)廣義積分。當(dāng)a,b中有一個(gè)符號(hào)表達(dá)式時(shí)，函數(shù)返回一個(gè)符號(hào)函數(shù)。,36/60,例.計(jì)算二重不定積分, syms x y F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y),37/60,4.符號(hào)代數(shù)方程求解,matlab符號(hào)運(yùn)算能夠解一般的線(xiàn)性方程、非線(xiàn)性方程及一般的代數(shù)方程、代數(shù)方程組。當(dāng)方程組不存在符號(hào)解時(shí)，又無(wú)其他自由參數(shù)，則給出數(shù)值解。 命令格式： solve(f) 求一個(gè)方程的解 Solve(f1,f2, fn) 求n個(gè)方程的解,38/60,例1. f = ax2+bx+c 求解 f=sym(a*x2+b*x+c); solve(f) 對(duì)缺省變量x求解 ans = 1/2/a*(-b+(b2-4*a*c)(1/2) 1/2/a*(-b-(b2-4*a*c)(1/2),計(jì)算機(jī) 格式,一般格式,39/60,例. 符號(hào)方程cos(x)=sin(x) 求解 f1=solve(cos(x)=sin(x), f1 = 1/4*pi,