1、信號、系統(tǒng)與信號處理實驗實驗報告實驗名稱：離散系統(tǒng)頻率響應和零極點分布一、 實驗目的通過MATLAB仿真簡單的離散時間系統(tǒng)，研究其時域特性，加深對離散系統(tǒng)的沖激響應，頻率響應分析和零、極點分布的概念的理解。二、 實驗內容與要求考一個LTI離散時間系統(tǒng)的輸入輸出差分方程為y(n)-1.6y(n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)（1）編程求出此系統(tǒng)的單位沖激響應序列，并畫出其波形。（2）若輸入序列x(n)=(n)+2(n-1)+3(n-2)+4(n-3)+5(n-4),編程求此系統(tǒng)輸出序列y(n),并畫出其波形。（3）編程得到系統(tǒng)頻響的幅度響應和相位響應，并畫圖。（4

2、）編程得到系統(tǒng)的零極點分布圖，分析系統(tǒng)的因果性和穩(wěn)定性。系統(tǒng)2：y(n)=0.45x(n)+0.5x(n-1)+0.45x(n-2)+0.53y(n-1)-0.46y(n-2)輸入x(n)=cos(20n/256)+cos(200n/256)0n299（5）編程得到系統(tǒng)頻響的幅度響應和相位響應，并畫圖。（6）編程得到系統(tǒng)的零極點分布圖，分析系統(tǒng)的因果性和穩(wěn)定性。三、 實驗程序與結果1、N=40;num=0.5 0.1 0;den=1 -1.6 1.28;y=impz(num,den,N);stem(y);2、n=0:49;x=1 2 3 4 5 zeros(1,45);num=0.5 0.1;

3、den=1 -1.6 1.28;y=filter(num,den,x);stem(y);3、fs=1000;num=0.5 0.1;den=1 -1.6 1.28;h,f=freqz(num,den,256,fs);mag=abs(h);ph=angle(h);ph=ph*180/pi;subplot(2,1,1),plot(f,mag);xlabel(頻率);ylabel(幅度);subplot(2,1,2),plot(f,ph);xlabel(頻率);ylabel(相位);4、num=0.5 0.1;den=1 -1.6 1.28;z,p,k=tf2zp(num,den);zplane(z

4、,p);5、fs=1000;num=0.45 0.5 0.45;den=1 -0.53 0.46;h,f=freqz(num,den,256,fs);mag=abs(h);ph=angle(h);ph=ph*180/pi;subplot(2,1,1),plot(f,mag);xlabel(頻率);ylabel(幅度);subplot(2,1,2),plot(f,ph);xlabel(頻率);ylabel(相位);6、num=0.45 0.5 0.45;den=1 -0.53 0.46;z,p,k=tf2zp(num,den);zplane(z,p);四、仿真結果分析對于系統(tǒng)y(n)-1.6y(

5、n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)，由圖4可知，零點在單位圓內，所以是因果的；極點在單位圓外，所以是不穩(wěn)定的。對于系統(tǒng)y(n)=0.45x(n)+0.5x(n-1)+0.45x(n-2)+0.53y(n-1)-0.46y(n-2)，由圖6可知，零點在單位圓內，所以是因果的；極點在單位圓內，所以是穩(wěn)定的。5、 實驗問題解答與體會1、分析系統(tǒng)零極點分布與系統(tǒng)頻率響應的關系極點會使調節(jié)時間變短，是系統(tǒng)反應更快，但是也會使系統(tǒng)的穩(wěn)定性變差，零點一般是使得穩(wěn)定性增加，但是會使調節(jié)時間變長；極點主要影響頻率響應的峰值，極點愈靠近單位圓，峰值愈尖銳；零點主要影響頻率特性的谷值，零點愈靠近單位圓，谷值愈深（當零點在單位圓上時，頻率特性為零）。(1)沖激響應波形是指指數(shù)衰減還是指數(shù)增長或等幅振蕩，主要取決于極點位于s左半平面還是右半平面或在虛軸上。(2)沖激響應波形衰減或增長快慢，主要取決于極點離虛軸的遠近。(3)沖激響應波形振蕩的快慢，主要取決于極點離實軸的遠近。零點分布只影響沖激響應函數(shù)的幅度和相位，不影響響應模式。對于結構不穩(wěn)定系統(tǒng)，改變系統(tǒng)結構后，只要適當選配參數(shù)就可使系統(tǒng)穩(wěn)定。2、實驗心得通過本次實驗，學會了通過MATLAB仿真