數(shù)字信號處理上機題15,16,17,18.doc_第1頁
數(shù)字信號處理上機題15,16,17,18.doc_第2頁
數(shù)字信號處理上機題15,16,17,18.doc_第3頁
數(shù)字信號處理上機題15,16,17,18.doc_第4頁
數(shù)字信號處理上機題15,16,17,18.doc_第5頁
已閱讀5頁,還剩4頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

生物醫(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=

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論