Matlab工程應(yīng)用基礎(chǔ)-3-2

1、第3章 矩陣、數(shù)組和符號(hào)運(yùn)算二、符號(hào)及運(yùn)算掌握內(nèi)容：（1）了解 MATLAB 6.0 的符號(hào)變量，掌握 MATLAB 符號(hào)表達(dá)式、符號(hào)矩陣的兩種創(chuàng)建方法。 （2）掌握 MATLAB 符號(hào)數(shù)學(xué)函數(shù)的創(chuàng)建。 （3）掌握符號(hào)矩陣的基本運(yùn)算及MATLAB 關(guān)于不同精度的控制方法。 （4）掌握符號(hào)微積分內(nèi)容，包括求函數(shù)的極限、對(duì)符號(hào)表達(dá)式求導(dǎo)數(shù)和微分、符號(hào)積分、符號(hào)求和、傅立葉變換及其逆變換等。（5）掌握各種符號(hào)方程的求解方法和函數(shù)命令。（6）了解 MATLAB 可視化的符號(hào)函數(shù)分析界面及使用。 （7）初步了解 MAPLE 的符號(hào)資源。第3章 矩陣、數(shù)組和符號(hào)運(yùn)算抽象運(yùn)算：公式推導(dǎo)、因式分解、求解代數(shù)方

2、程或微分方程的精確解 符號(hào)數(shù)學(xué)工具箱 1）通過(guò)基本符號(hào)數(shù)學(xué)工具箱的專用函數(shù)；符號(hào)表達(dá)式和符號(hào)矩陣的操作； 多項(xiàng)式的化簡(jiǎn)、展開(kāi)和代入； 線性代數(shù)； 微積分； 符號(hào)方程的求解； 特殊的數(shù)學(xué)函數(shù)。2）通過(guò) maple.m、mpa.m 兩個(gè)專門(mén)設(shè)計(jì)的 M 文件進(jìn)行符號(hào)運(yùn)算；3) 通過(guò) MATLAB 中的函數(shù)計(jì)算器（Function Caculator）。 第3章 矩陣、數(shù)組和符號(hào)運(yùn)算1、符號(hào)變量的創(chuàng)建a. sym 函數(shù) S=sym(arg) ，從表達(dá)式 arg 創(chuàng)建一個(gè) sym 對(duì)象 Sx=sym(x)x = sym(x,real)x = sym(x,unreal) 附加屬性x = sym(x,pos

3、itive)pi = sym(pi)delta = sym(1/10) S = sym(A, flag) ，將數(shù)值或矩陣轉(zhuǎn)化為符號(hào)形式其中 flag 選項(xiàng)有四項(xiàng)參數(shù)f, r, e 和 d，r為缺省項(xiàng)。f：代表十六進(jìn)制浮點(diǎn)形式；r：代表有理數(shù)形式；e：估計(jì)誤差；d：表示十進(jìn)制小數(shù)。第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 A=2/5,4/0.78,sqrt(23)/3;0.33,0.3333,log(4) 輸入數(shù)值矩陣AA = 0.4000 5.1282 1.5986 0.3300 0.3333 1.3863 FA=sym(A) 將數(shù)值矩陣A轉(zhuǎn)化為符號(hào)矩陣FA FA = 2/5, 200/39, sqrt(

4、23/9) 33/100, 3333/10000, 6243314768165359*2(-52)不管數(shù)值矩陣的元素是以分?jǐn)?shù)或是浮點(diǎn)數(shù)表示，轉(zhuǎn)換后的符號(hào)矩陣都將以最接近有理式的形式給出。b. syms 函數(shù) syms arg1 arg2 . syms a b c x y2、符號(hào)表達(dá)式和矩陣的創(chuàng)建a.字符串直接輸入創(chuàng)建 符號(hào)表達(dá)式和符號(hào)方程對(duì)空格很敏感。因此，在創(chuàng)建符號(hào)表達(dá)式或符號(hào)方程時(shí)，不要在字符間任意加空格符；符號(hào)計(jì)算中出現(xiàn)的數(shù)字也是當(dāng)作符號(hào)處理的； f=a*x2+b*x+cf =a*x2+b*x+c f=a*x2+b*x+c=0f =a*x2+b*x+c=0第3章 矩陣、數(shù)組和符號(hào)運(yùn)算第3

5、章 矩陣、數(shù)組和符號(hào)運(yùn)算這種方法輸入符號(hào)矩陣與字符串矩陣的輸入相似。但要保證在同一列中各元素字符串有同樣的長(zhǎng)度，在較短的字符串前后用空格符填充；這種方法要求符號(hào)矩陣每一行的兩端都有方括號(hào)，而字符串矩陣僅在首尾有方括號(hào)。 B=4+x x2 x ;x3 5*x-3 x*aB = 4+x x2 x x3 5*x-3 x*a 第3章 矩陣、數(shù)組和符號(hào)運(yùn)算b.由 sym 命令創(chuàng)建 f=sym(a*x2+b*x+c)f =a*x2+b*x+c f1=sym(a*x2+b*x+c=0)f1 =a*x2+b*x+c=0 A=sym(4+x, x2, x;x3, 5*x-3, x*a)A = 4+x, x2,

6、x x3, 5*x-3, x*a第3章 矩陣、數(shù)組和符號(hào)運(yùn)算c.由 syms 命令創(chuàng)建 syms x a b c f=a*x2+b*x+cf =a*x2+b*x+c syms x a B=4+x x2 x;x3 5*x-3 x*aB = 4+x, x2, x x3, 5*x-3, x*a不能創(chuàng)建符號(hào)方程 第3章 矩陣、數(shù)組和符號(hào)運(yùn)算3、數(shù)字矩陣和符號(hào)矩陣的轉(zhuǎn)換 MATLAB 中的數(shù)值型、字符型和符號(hào)型三種數(shù)據(jù)類型中數(shù)值變量級(jí)別最低，字符變量級(jí)別居中，符號(hào)變量級(jí)別最高；三種變量參與的混合運(yùn)算，系統(tǒng)將會(huì)把所有參與運(yùn)算的變量自動(dòng)統(tǒng)一轉(zhuǎn)換為變量等級(jí)最高的類型，然后進(jìn)行計(jì)算；可以通過(guò)命令來(lái)完成對(duì)不同數(shù)據(jù)

7、類型之間的轉(zhuǎn)換，大致可以分為三種情況：轉(zhuǎn)換為數(shù)值變量: double, str2num,numeric轉(zhuǎn)化為符號(hào)變量: sym轉(zhuǎn)化為字符變量: int2str, num2str第3章 矩陣、數(shù)組和符號(hào)運(yùn)算4、MATLAB 關(guān)于不同精度的控制 針對(duì)浮點(diǎn)運(yùn)算的數(shù)值算法計(jì)算速度最快，占用計(jì)算機(jī)內(nèi)存最少的算法，與 C、FORTRAN 語(yǔ)言中的浮點(diǎn)運(yùn)算算法完全相同。在機(jī)器內(nèi)的表達(dá)和計(jì)算都是一個(gè)被“ 截?cái)唷钡? 位浮點(diǎn)近似值。針對(duì)精確運(yùn)算的符號(hào)算法計(jì)算時(shí)間最長(zhǎng)，內(nèi)存占用最多，精度也最高。任意精度的算法運(yùn)算時(shí)間、內(nèi)存占用和計(jì)算精度均介于以上兩種運(yùn)算之間。采用函數(shù) digits 來(lái)控制十進(jìn)制結(jié)果的有效位數(shù)。d

8、igits 的缺省值為 32，大約對(duì)應(yīng)于浮點(diǎn)精度。符號(hào)數(shù)學(xué)工具箱中，用 vpa 函數(shù)執(zhí)行任意精度運(yùn)算。第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 syms x f=x-cos(x)f = x-cos(x) f1=subs(f,x,pi) %將符號(hào)表達(dá)式中的變量替換為另一變量f1 = pi+1 digits(25) vpa(f1)ans = 4.141592653589793238462643 vpa(f1,6)ans = 4.14159 numeric(f1)ans = 4.1416 double(f1)ans = 4.14165、符號(hào)矩陣的運(yùn)算基本運(yùn)算 四則運(yùn)算兩個(gè)符號(hào)矩陣的大小相等方可進(jìn)行加減運(yùn)算，符號(hào)

9、矩陣和符號(hào)標(biāo)量的加減運(yùn)算按照數(shù)組運(yùn)算規(guī)則進(jìn)行； 兩個(gè)符號(hào)矩陣只有內(nèi)積相等時(shí)才可以進(jìn)行乘法運(yùn)算； 符號(hào)的乘方運(yùn)算 Sp，若 S 為符號(hào)表達(dá)式，p 可以為符號(hào)表達(dá)式或數(shù)值表達(dá)式；若 S 為符號(hào)矩陣，則 p 必須是整數(shù)。第3章 矩陣、數(shù)組和符號(hào)運(yùn)算第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 a=sym(1/x, 1/(x+1); 1/(x+2), 1/(x+3)a = 1/x, 1/(x+1) 1/(x+2), 1/(x+3) b=sym(x, 1; x+2, 0)b = x, 1 x+2, 0 b-aans = x-1/x, 1-1/(x+1) x+2-1/(x+2), -1/(x+3) abans = -6*

10、x-2*x3-7*x2, 3/2*x2+x+1/2*x3 6+2*x3+10*x2+14*x, -1/2*x3-2*x2-3/2*x a.bans = x2, x+1 (x+2)2, 0 a2ans = 1/x2+1/(x+1)/(x+2), 1/x/(x+1)+1/(x+1)/(x+3) 1/(x+2)/x+1/(x+3)/(x+2), 1/(x+1)/(x+2)+1/(x+3)2 exp(b)ans = exp(x), exp(1) exp(x+2), 1第3章 矩陣、數(shù)組和符號(hào)運(yùn)算矩陣運(yùn)算 aans = 1/conj(x), 1/(2+conj(x) 1/(1+conj(x), 1/(3

11、+conj(x) inv(a)ans = 1/2*x*(x+1)*(x+2), -1/2*x*(x+3)*(x+2) -1/2*x*(x+3)*(x+1), 1/2*(x+3)*(x+1)*(x+2) det(a)ans =2/x/(x+3)/(x+1)/(x+2) rank(a)ans =2 eig(b)ans = 1/2*x+1/2*(x2+4*x+8)(1/2) 1/2*x-1/2*(x2+4*x+8)(1/2) triu(a)ans = 1/x, 1/(x+1) 0, 1/(x+3) diag(a)ans = 1/x 1/(x+3) tril(a)ans = 1/x, 0 1/(x+2

12、), 1/(x+3)第3章 矩陣、數(shù)組和符號(hào)運(yùn)算b.符號(hào)矩陣的簡(jiǎn)化因式分解 syms x factor(x9-1)ans =(x-1)*(x2+x+1)*(x6+x3+1)符號(hào)矩陣展開(kāi) syms x y expand(x+1)3)ans =x3+3*x2+3*x+1 expand(sin(x+y)ans =sin(x)*cos(y)+cos(x)*sin(y)同類項(xiàng)合并 syms x y collect(x2*y+y*x-x2-2*x)ans =(y-1)*x2+(y-2)*x分式通分 syms x y n,d=numden(x/y+y/x)n =x2+y2d =y*x 符號(hào)簡(jiǎn)化 syms x

13、 simple(cos(x)2-sin(x)2)simplify:2*cos(x)2-1radsimp:cos(x)2-sin(x)2combine(trig):cos(2*x)factor:(cos(x)-sin(x)*(cos(x)+sin(x)expand:cos(x)2-sin(x)2combine:cos(2*x)第3章 矩陣、數(shù)組和符號(hào)運(yùn)算convert(exp):(1/2*exp(i*x)+1/2/exp(i*x)2+1/4*(exp(i*x)-1/exp(i*x)2convert(sincos):cos(x)2-sin(x)2convert(tan):(1-tan(1/2*x)

14、2)2/(1+tan(1/2*x)2)2-4*tan(1/2*x)2/(1+tan(1/2*x)2)2collect(x):cos(x)2-sin(x)2ans =cos(2*x) simplify(cos(x)2-sin(x)2)ans =2*cos(x)2-1 simplify(sin(x)2+cos(x)2)ans =1第3章 矩陣、數(shù)組和符號(hào)運(yùn)算第3章 矩陣、數(shù)組和符號(hào)運(yùn)算6、符號(hào)微積分Matlab自變量確定原則：除i、j外，字母位置最接近x的小寫(xiě)字母為自變量；如果表達(dá)式中沒(méi)有變量，x會(huì)被視為默認(rèn)的變量。由函數(shù)findsym可一找到默認(rèn)變量a. 符號(hào)極限limit(F,x,a)limi

15、t(F,a)，變量為由findsym定義的默認(rèn)變量limit(F,x,a,right) /limit(F,x,a,left)Limit(F), a=0 syms x a t h limit(sin(x)/x)ans =1 limit(1+2*t/x)(3*x),x,inf)ans =exp(6*t) limit(1/x,x,0,right)ans =infb.符號(hào)積分int(S) 不定積分int(S,v)int(S,a,b) 定積分int(S,v,a,b) syms x x1 alpha u t; A=cos(x*t),sin(x*t);-sin(x*t),cos(x*t)A = cos(x*

16、t), sin(x*t) -sin(x*t), cos(x*t) int(A,t)ans = 1/x*sin(x*t), -cos(x*t)/x cos(x*t)/x, 1/x*sin(x*t) int(x1*log(1+x1),0,1)ans =1/4 第3章 矩陣、數(shù)組和符號(hào)運(yùn)算符號(hào)合計(jì)函數(shù)Symsum(S)Symsum(S,v)Symsum(S,a,b) syms k n; simplify(symsum(k)ans =1/2*k2-1/2*k simplify(symsum(k2,0,n)ans =1/3*n3+1/2*n2+1/6*n simplify(symsum(k2,0,10)

17、ans =385 c. 符號(hào)微分和差分微分和差分Diff(S)Diff(S,v)Diff(S,n)Diff(S,v,n) syms x diff(sin(x2)ans = 2*cos(x2)*x diff(t6,6)ans = 720第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 syms x t; A=cos(x*t),sin(x*t);-sin(x*t),cos(x*t)A = cos(x*t), sin(x*t) -sin(x*t), cos(x*t) diff(A,x) % diff(A,x)ans = -sin(x*t)*t, cos(x*t)*t -cos(x*t)*t, -sin(x*t)*t d

18、iff(A,t,2)ans = -cos(x*t)*x2, -sin(x*t)*x2 sin(x*t)*x2, -cos(x*t)*x2第3章 矩陣、數(shù)組和符號(hào)運(yùn)算梯度函數(shù)DF=gradient(F) 一維梯度FX,FY=gradient(F) 二維梯度FX,FY=gradient(F,H)FX,FY=gradient(F,HX,HY)FX,FY,FZ=gradient(F) 三維梯度FX,FY,FZ=gradient(F,HX,HY,HZ)H指定間距，默認(rèn)為1 x,y=meshgrid(-2:0.2:2,-2:0.2:2); z=x.*exp(-x.2-y.2); px,py=gradien

19、t(z,0.2,0.2); contour(z) hold on quiver(px,py) %畫(huà)矢量圖 hold offd.傅立葉變換和傅立葉逆變換傅立葉快速離散變換 MATLAB 提供了 fft(內(nèi)置函數(shù))、ifft、fft2、ifft2、fftn、ifftn、fftshift、ifftshift 等函數(shù)，用來(lái)計(jì)算矩陣的離散快速傅立葉變換。 函數(shù) fft 和 ifft 函數(shù) fft 最完整的調(diào)用格式為： Y=fft(X,dim) 或 Y= fft(X,n,dim) 數(shù)據(jù)長(zhǎng)度n是 2 次冪時(shí)，可以采用基-2 算法進(jìn)行快速計(jì)算。輸入?yún)?shù) X 可以是向量、矩陣。 dim 指定變換的實(shí)施方向。當(dāng)

20、X 是矩陣時(shí)， 1 指明變換按列進(jìn)行（默認(rèn)），2 指明變換按行進(jìn)行。 第3章 矩陣、數(shù)組和符號(hào)運(yùn)算第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 X=1,2,3;4,5,6;7,8,9X = 1 2 3 4 5 6 7 8 9 Y=fft(X)Y = 12.0000 15.0000 18.0000 -4.5000 + 2.5981i -4.5000 + 2.5981i -4.5000 + 2.5981i -4.5000 - 2.5981i -4.5000 - 2.5981i -4.5000 - 2.5981i函數(shù) fft2 和 ifft2 函數(shù) fft2 和 ifft2 是對(duì)數(shù)據(jù)做二維快速傅立葉變換和逆傅立葉變

21、換。數(shù)據(jù)的二維傅立葉變換 fft2(X)相當(dāng)于 fft(fft(X)，即先對(duì) X 的列做一維傅立葉變換，然后對(duì)變換結(jié)果的行做一維傅立葉變換。其調(diào)用格式為： Y=fft2(X,mrows,ncols)和函數(shù) fft2 和 ifft2 類似，函數(shù) fftn 和 ifftn 對(duì)數(shù)據(jù)做多維快速傅立葉變換。函數(shù) fftshift 和 ifftshift 函數(shù) fftshift 用于把傅立葉變換結(jié)果（ 頻域數(shù)據(jù)）中的直流分量（ 頻率為 0 處的值）移到中間位置。其調(diào)用格式為： Y=fftshift(X) X=rand(3,3)X = 0.9501 0.4860 0.4565 0.2311 0.8913 0

22、.0185 0.6068 0.7621 0.8214 Y=fft(X)Y = 1.7881 2.1394 1.2964 0.5311 + 0.3254i -0.3407 - 0.1119i 0.0365 + 0.6953i 0.5311 - 0.3254i -0.3407 + 0.1119i 0.0365 - 0.6953i第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 Y1=fftshift(X)Y1 = 0.8214 0.6068 0.7621 0.4565 0.9501 0.4860 0.0185 0.2311 0.8913傅立葉積分變換及其反變換 離散傅立葉變換（ DFT）作用于有限數(shù)據(jù)采樣，傅立葉變換

23、作用于連續(xù)函數(shù) 。傅立葉變換調(diào)用格式為： F = fourier(f)：求表達(dá)式 f 的傅立葉變換。缺省的自變量為 x，缺省的返回值是關(guān)于 w 的函數(shù)。F = fourier(f,v)：返回函數(shù) F 是關(guān)于符號(hào)表達(dá)式對(duì)象 v 的函數(shù)，而不是缺省的 w。F=fourier(f,u,v)：對(duì)關(guān)于 u 的函數(shù) f 進(jìn)行變換，返回函數(shù) F 是 關(guān)于 v 的函數(shù)。傅立葉逆變換調(diào)用格式為： f = ifourier(F)：符號(hào)表達(dá)式對(duì)象F的傅立葉逆變換。缺省的自變量為 w，缺省返回是關(guān)于 x 的函數(shù)。f = ifourier(F,u)：返回函數(shù) f 是關(guān)于符號(hào)表達(dá)式對(duì)象 u 的函數(shù)，而不是缺省的 x的函數(shù)。f = ifourier(F,v,u)：對(duì)關(guān)于 v 的函數(shù) F 進(jìn)行變換，返回關(guān)于 u 的函數(shù) f。第3章 矩陣、數(shù)組和符號(hào)運(yùn)算 syms x t w u; fourier(exp(-x2)ans = pi(1/2)*exp(-1/4*w2) fourier(exp(-x2),u)ans =pi(1/2)*exp(-1/4*u2) fourier(exp