數(shù)字信號(hào)處理實(shí)驗(yàn)實(shí)驗(yàn)二.doc

數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)名稱：離散系統(tǒng)頻率響應(yīng)和零極點(diǎn)分布學(xué)生姓名： z 學(xué)生學(xué)號(hào)： 學(xué)生班級(jí)： 上課時(shí)間： 周二上午 指導(dǎo)老師： 一、 實(shí)驗(yàn)?zāi)康耐ㄟ^(guò)MATLAB仿真簡(jiǎn)單的離散時(shí)間系統(tǒng)，研究其頻域特性，加深對(duì)離散系統(tǒng)的沖激響應(yīng)，頻率響應(yīng)分析和零極點(diǎn)分布概念的理解。二、 實(shí)驗(yàn)原理若用x(n)和y(n)分別表示一個(gè)線性時(shí)不變（LTI）離散時(shí)間系統(tǒng)，則其輸入，輸出關(guān)系可用以下常系數(shù)線性差分方程描述。記LTI離散時(shí)間系統(tǒng)單位沖激響應(yīng)為h(n),h(n)反映了系統(tǒng)固有的特征，它是離散系統(tǒng)的一個(gè)重要參數(shù)。任意LTI系統(tǒng)都可由系統(tǒng)單位沖激響應(yīng)h(n)表示，相應(yīng)地在頻域響應(yīng)H（w）表示，它是h(n)的傅里葉變換。使系統(tǒng)函數(shù)的分母多項(xiàng)式等于零的z值，稱為系統(tǒng)的極點(diǎn)；同理，使系統(tǒng)函數(shù)的分子多項(xiàng)式等于零的z值，稱為系統(tǒng)的零點(diǎn)。通過(guò)系統(tǒng)的零極點(diǎn)增益表達(dá)式，可以判斷一個(gè)LTI離散時(shí)間系統(tǒng)的穩(wěn)定性。對(duì)一個(gè)因果的離散時(shí)間系統(tǒng)，若所有的極點(diǎn)都位于單位圓內(nèi)，則系統(tǒng)是穩(wěn)定的。同理，由零極點(diǎn)分布圖可大致估計(jì)出系統(tǒng)的頻率響應(yīng)：（1） 單位圓附近的零點(diǎn)對(duì)幅度響應(yīng)的谷點(diǎn)的位置與深度有明顯影響，當(dāng)零點(diǎn)位于單位圓上時(shí)，谷點(diǎn)為零。零點(diǎn)可在單位圓外。（2） 單位圓附近的極點(diǎn)對(duì)幅度響應(yīng)的峰點(diǎn)位置和高度有明顯影響。三、實(shí)驗(yàn)內(nèi)容一個(gè)LTI離散時(shí)間系統(tǒng)的輸入輸出差分方程為y(n)-1.6y(n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)（1） 編程求此系統(tǒng)的單位沖激響應(yīng)序列，并畫(huà)出其波形。MATLAB程序如下：N=100;num=0.5 0.1 0;den=1 -1.6 1.28;y=impz(num,den,N);stem(y);xlabel(時(shí)間序號(hào));ylabel(信號(hào)幅度);title(沖激響應(yīng));grid;執(zhí)行結(jié)果如下圖形所示：(2) 如輸入序列x（n）=,編程求此系統(tǒng)輸出序列y（n）,并畫(huà)出其圖形。MATLAB程序如下：n=100;x=1 2 3 4 5 zeros(1,n-5) ;subplot(2,1,1);stem(x);xlabel(時(shí)序信號(hào));ylabel(信號(hào)幅度);title(x(n);grid;y=filter(num,den,x);subplot(2,1,2);stem(y);xlabel(時(shí)序信號(hào));ylabel(信號(hào)幅度);title(y(n);grid; 執(zhí)行結(jié)果如下圖形所示：(3) 編程得到系統(tǒng)頻率響應(yīng)的幅度響應(yīng)和相位響應(yīng)，并畫(huà)圖。MATLAB程序如下：fs=1000;b=0.5 0.1 ;a=1 -1.6 1.28;h,f=freqz(b,a,256,fs);mag=abs(h);ph=angle(h);ph=ph*180/pi;subplot(2,1,1);plot(f,mag);xlabel(頻率(Hz);ylabel(幅度);subplot(2,1,2);plot(f,ph);xlabel(頻率(Hz);ylabel(幅度); 執(zhí)行結(jié)果如下圖形所示：(4) 編程得到系統(tǒng)的零極點(diǎn)分布圖，分析系統(tǒng)的因果性和穩(wěn)定性。MATLAB程序如下：.b=0.5 0.1 ;a=1 -1.6 1.28;z,p,k=tf2zp(b,a);zplane(z,p); xlabel(實(shí)部);ylabel(虛部);title(系統(tǒng)的零極點(diǎn)分布圖); 執(zhí)行結(jié)果如下圖形所示：. MATLAB程序如下：b=0.5 0.1 ;a=1 -1.6 1.28;zplane(b,a);xlabel(實(shí)部);ylabel(虛部);title(系統(tǒng)的零極點(diǎn)分布圖); 執(zhí)行結(jié)果如下圖形所示：答：由上面圖形可知系統(tǒng)的極點(diǎn)落在單位圓外，所以本系統(tǒng)是不穩(wěn)定的。三、 記錄調(diào)試運(yùn)行情況及所遇問(wèn)題的解決方法（1） 采用filter函數(shù)時(shí)，輸出y與輸入x的長(zhǎng)度是相同的，所以在（2）中n=100,x必須補(bǔ)95個(gè)0才能與y長(zhǎng)度相同，所以程序如下x=1 2 3 4 5 zeros(1,n-5) ; （2）在MATLAB中編寫(xiě)程序要注意不能用中文標(biāo)點(diǎn)，否則程序編譯時(shí)會(huì)出現(xiàn)錯(cuò)誤。（3）在問(wèn)題（2）中，x(n)的波形是n=1時(shí)，幅度=1，n=2時(shí)，幅度=2.實(shí)際情況應(yīng)該是n=0時(shí)，幅度=1，n=1時(shí)，幅度=2