序列逆Z變換Matlab實(shí)現(xiàn)課件

、序列逆乙變換的Matlab實(shí)現(xiàn)函數(shù)residues:適合計(jì)算離散系統(tǒng)有理函數(shù)的留數(shù)和極點(diǎn),可以用于求解序列的逆乙變換。X(z)=b+b1z-+…+bnzB(z2.98)a,+a,z+“+azA(z)A-11-Pz-0函數(shù)residues基本調(diào)用方式>>[r,p,c=residue(b,a)輸入?yún)?shù):b=b0,b1,…,bM]為分子多項(xiàng)式的系數(shù),a=a0,排列為分母多項(xiàng)式的系數(shù),這些多項(xiàng)式都按z的降冪輸出參數(shù):r是極點(diǎn)的留數(shù),p是極點(diǎn),c是無窮項(xiàng)多項(xiàng)式的系數(shù)項(xiàng),僅當(dāng)MN時(shí)存在。例:計(jì)算逆乙變換例:計(jì)算X(2)=2-32+1的逆Z變換。解:有理分式X2)分子和分母0+zX(z)多項(xiàng)式都按z的降冪排列2z2-3z+12-3x1+>b=[0,1];a=[2,3,1;%多項(xiàng)式的系數(shù)Tr;p,c}=residue(b,a);%求留數(shù)、極點(diǎn)和系數(shù)項(xiàng)啦disp(留數(shù):");disp(r);%顯示輸出參數(shù)disp(極點(diǎn):);disp(p);X(2)的部分分式形式為dsp(系數(shù)項(xiàng):)disp(c);程序運(yùn)行結(jié)果為留數(shù):1-1逆z變換為極點(diǎn):1.00000.5000系數(shù)項(xiàng):x(n)=t(n)-(0.5)"u(n)DFT物理意義的Matlab實(shí)現(xiàn)·序列的N點(diǎn)DFT的物理意義:對(duì)X(ejo)在[O,2m]上進(jìn)行N點(diǎn)的等間隔取樣函數(shù)ft用于快速計(jì)算離散傅里葉變換,調(diào)用方式為y=fft(x,N)y=ft(x利用FFT算法計(jì)算序列x的離散傅里葉變換當(dāng)x為矩陣時(shí),y為矩陣x每一列的FFT當(dāng)x長度為2的整數(shù)次冪時(shí),函數(shù)ft采用基-2的FFT算法,否則采用混合基算法。y=ft(,N)采用N點(diǎn)FFT當(dāng)序列x長度小于N時(shí),函數(shù)ft動(dòng)對(duì)序列尾部補(bǔ)零,構(gòu)成N點(diǎn)數(shù)當(dāng)κ長度大于N時(shí),函數(shù)ft自動(dòng)截取序列前面N點(diǎn)數(shù)據(jù)進(jìn)行FFT旳數(shù)ift用于快速計(jì)算冋量或矩陣的離散傅里葉逆變換,與函數(shù)ft的調(diào)用規(guī)則基本相同調(diào)用方式為>>y=ifft(x)>>y=ifft(x,N)例:利用FFT實(shí)現(xiàn)線性卷積例:利用FFT實(shí)現(xiàn)線性卷積。已知序列x(n)=R4(m),求(1)用conv函數(shù)求x(m)與x(m)的線性卷積ym),并繪出圖形;(2)用FFT求x(m)與x(m)的4點(diǎn)循環(huán)卷積y1(m),并繪出圖形(3)用FFT求x(m)