數(shù)字信號處理上機題15,16,17,18.doc

生物醫(yī)學(xué)工程111班 耿慧超 61034110161、課本31頁第15題：已知系統(tǒng)的差分方程和輸入信號y（n）-0.5y(n-1)=x(n)+2x(n-2)，X(n)=1,2,3,4,2,1;用遞推法求零狀態(tài)響應(yīng)。輸入代碼如下：a1=-0.5;b1=0;b2=2;B=1,b1,b2;A=1,a1;xn=1,2,3,4,2,1; yn=filter(B,A,xn); n=0:length(yn)-1;subplot(1,1,1);stem(n,yn,.)title(a);xlabel(n);ylabel(y(n)最終得到的波形圖為：結(jié)果分析：因為是求零狀態(tài)響應(yīng)，所以y（-1）=0.當(dāng)n=0時，y（0）=x（0）+2x（-2）=1當(dāng)n=1時，y（1）=0.5y（0）+x（1）+2x（-1）=2.5當(dāng)n=2時，y（2）=0.5y（1）+x（2）+2x（0）=6.25當(dāng)n=3時，y（3）=0.5y（2）+x（3）+2x（1）=11.125當(dāng)n=4時，y（4）=0.5y（3）+x（4）+2x（2）=13.5625當(dāng)n=5時，y（5）=0.5y（4）+x（5）+2x（3）=15.78125因為后面一個數(shù)都是在前面的基礎(chǔ)上加正值，所以結(jié)果會越來越大，但是增加的幅度變小。2、課本31頁第16題：已知兩個系統(tǒng)的差分方程分別為（1）Y(n)=0.6y(n-1)-0.08y(n-2)+x(n)（2）y(n)=0.7y(n-1)-0.1y(n-2)+2x(n)-x(n-2) 分別求出所描述的系統(tǒng)的單位脈沖響應(yīng)和單位階躍響應(yīng)。（1）輸入代碼如下：ys=0;xn=1,zeros(1,30);B=1;A=1,-0.6,0.08;xi=filtic(B,A,ys);yn=filter(B,A,xn,xi);n=0:length(yn)-1;subplot(2,1,1);stem(n,yn,.)title(a);xlabel(n);ylabel(y(n)ys=0;xn=1,ones(1,30);B=1;A=1,-0.6,0.08;xi=filtic(B,A,ys);yn=filter(B,A,xn,xi);n=0:length(yn)-1;subplot(2,1,2);stem(n,yn,.)title(a);xlabel(n);ylabel( h(n)最終得到的波形圖為：結(jié)果分析：單位沖激響應(yīng)：y(0)=0.6y(-1)-0.08y(-2)+x(0)=1 y(1)=0.6y(0)-0.08y(-1)+x(1)=0.6 y (2)=0.6y(1)-0.08y(0)+x(2)=0.28 y(3)=0.6y(2)-0.08y(1)+x(3)=0.12 y(4)=0.6y(3)-0.08y(2)+x(4)=0.0496 y(5)=0.6y(4)-0.08y(3)+x(5)=0.02016 單位階躍響應(yīng)：y(0)=0.6y(-1)-0.08y(-2)+x(0)=1 y(1)=0.6y(0)-0.08y(-1)+x(1)=1.6 y (2)=0.6y(1)-0.08y(0)+x(2)=1.88 y(3)=0.6y(2)-0.08y(1)+x(3)=2 y(4)=0.6y(3)-0.08y(2)+x(4)=2.0496 y(5)=0.6y(4)-0.08y(3)+x(5)=2.06976由結(jié)果可以看出與圖形相符，說明程序是正確的，在n=5之前y（n）所取的值變化大，之后就幾乎不變這是由單位采樣序列（n）或單位階躍序列u（n）本身的性質(zhì)決定的。（2）輸入代碼如下：ys=0;xn=1,zeros(1,30);B=2,0,-1;A=1,-0.7,0.1;xi=filtic(B,A,ys);yn=filter(B,A,xn,xi);n=0:length(yn)-1;subplot(2,1,1);stem(n,yn,.)title(a);xlabel(n);ylabel(y(n)ys=0;xn=1,ones(1,30);B=2,0,-1;A=1,-0.7,0.1;xi=filtic(B,A,ys);yn=filter(B,A,xn,xi);n=0:length(yn)-1;subplot(2,1,2);stem(n,yn,.)title(a);xlabel(n);ylabel(h(n)最終得到的波形圖為：結(jié)果分析：單位沖激響應(yīng)：y(0)=0.7y(-1)-0.1y(-2)+2x(0)-x(-2)=2 y(1)=0.7y(0)-0.1y(-1)+2x(1)-x(-1)=1.4 y (2)=0.7y(1)-0.1y(0)+2x(2)-x(0)=-0.22 y(3)=0.7y(2)-0.1y(1)+2x(3)-x(1)=-0.294 y(4)=0.7y(3)-0.1y(2)+2x(4)-x(2)=-0.146 y(5)=0.7y(4)-0.1y(3)+2x(5)-x(3)=-0.0728 單位階躍響應(yīng)：y(0)=0.7y(-1)-0.1y(-2)+2x(0)-x(-2)=2 y(1)=0.7y(0)-0.1y(-1)+2x(1)-x(-1)=3.4 y (2)=0.7y(1)-0.1y(0)+2x(2)-x(0)=3.18 y(3)=0.7y(2)-0.1y(1)+2x(3)-x(1)=2.886 y(4)=0.7y(3)-0.1y(2)+2x(4)-x(2)=2.7022 y(5)=0.7y(4)-0.1y(3)+2x(5)-x(3)=2.60294（2）和（1）題型相同，結(jié)果分析類似。17、已知系統(tǒng)的差分方程為y(n)=a1y(n1)a2y(n2)+bx(n)其中，a1=0.8，a2=0.64，b=0.866。（1）編寫求解系統(tǒng)單位脈沖響應(yīng)h(n)(0n49)的程序，并畫出h(n)(0n49)； (2）編寫求解系統(tǒng)零狀態(tài)單位階躍響應(yīng)s(n)(0n100）的程序，并畫出s(n)(0n100）輸入代碼如下:a1=-0.8;a2=0.64;b=0.886;B=b;A=1,a1,a2;xn=1,zeros(1,49);yn=filter(B,A,xn);n=0:length(yn)-1;subplot(2,1,1);stem(n,yn,.)title(a);xlabel(n);ylabel(y(n)a1=-0.8;a2=0.64;b=0.886;B=b;A=1,a1,a2;xn=1,ones(1,100);yn=filter(B,A,xn);n=0:length(yn)-1;subplot(2,1,2);stem(n,yn,.)title(a);xlabel(n);ylabel(y(n)最終得到的波形圖為：結(jié)果分析：計算所得結(jié)果與圖形相符，無論是單位脈沖響應(yīng)還是單位階躍響應(yīng)都是在n=20之前y（n）的波動比較大，之后就幾乎不變。18、在題18*圖中，有四個分系統(tǒng)T1、T2、T3和T4，四個分系統(tǒng)分別用下面的單位脈沖響應(yīng)或者差分方程描述： T1：h1(n)=1/2n n=0，1，2，3，4，5 0 其它 T2：h2(n)=1 n=0，1，2，3，4，5 0 其它 T3：y3(n)=1/4x(n)+1/2x(n-1)+1/4x(n-2) T4：y(n)=0.9y(n-1)-0.81y(n-2)+v(n)+v(n-1)編寫程序計算整個系統(tǒng)的單位脈沖響應(yīng)h(n)，0n99。設(shè) x(n)=(n)， 則 v(n)=h1(n)*h2(n)+h3(n)該式調(diào)用conv 函數(shù)計算。h(n)=T4v(n)該式調(diào)用filter 函數(shù)計算。調(diào)用MATLAB 函數(shù)conv 和filter 計算該系統(tǒng)的系統(tǒng)響應(yīng)的程序ex118.m如下：h1n=1,1/2,1/4,1/8,1/16,1/32;h2n=ones(1,6);h3n=1/4,1/2,1/4,zeros(1,97);h12n=conv(h1n,h2n);h12n=h12n,zeros(1,89);vn=h12n+h3n;B4=