數(shù)字信號處理實驗二時域離散系統(tǒng)及系統(tǒng)響應(yīng)

1、實驗二 實驗二 時域離散系統(tǒng)及系統(tǒng)響應(yīng) 時域離散系統(tǒng)及系統(tǒng)響應(yīng)一、 一、實驗?zāi)康?實驗?zāi)康?、掌握求解離散時間系統(tǒng)沖激響應(yīng)和階躍響應(yīng)的方法； 2、進一步理解卷積定理，掌握應(yīng)用線性卷積求解離散時間系統(tǒng)響應(yīng)的基本方法；3、掌握離散系統(tǒng)的響應(yīng)特點。二、實驗內(nèi)容 二、實驗內(nèi)容 1、請分別用impz 和dstep函數(shù)求解下面離散時間系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)。（1）系統(tǒng)的差分方程為： ) ( 866 . 0 ) 2 ( 64 . 0 ) 1 ( 8

2、 . 0 ) ( n x n y n y n y ? ? ? ? ?（2）系統(tǒng)的系統(tǒng)函數(shù)為： 2 1115 . 0 1 ) ( ? ??? ?? ? z zz z H2、運行例題2.3，理解卷積過程和程序中每一句的意義。3、利用第 1 題求得的系統(tǒng)沖激響應(yīng)求解系統(tǒng)在激勵 下的響應(yīng)。 ) 3 ( ) ( ? ? n u n x三、 三、實現(xiàn)步驟 實現(xiàn)步驟4、請分別用impz 和dstep函數(shù)求解下面離散時間系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)。（1）

3、系統(tǒng)的差分方程為： ) ( 866 . 0 ) 2 ( 64 . 0 ) 1 ( 8 . 0 ) ( n x n y n y n y ? ? ? ? ?% y(n) = 0.8y(n -1) - 0.64y(n - 2) + 0.866x(n)% y(n) - 0.8y(n -1) + 0.64y(n - 2) = 0.866x(n)% a0=1 ,a1= -0.8 ,a2 =0.64 % b0=0.866 ,b1= 0 ,b2 =0a

4、=[1,-0.8,0.64];b=[0.866,0,0];n=20;hn=impz(b,a,n); %沖激響應(yīng)gn=dstep(b,a,n); %階躍響應(yīng)subplot(2,1,1),stem(hn,'filled'); %顯示沖激響應(yīng)曲線title('系統(tǒng)的單位沖激響應(yīng)');ylabel('h(n)');xlabel('n'); axis([0,n,1.1*min(hn)

5、,1.1*max(hn)]);subplot(2,1,2),stem(gn,'filled'); %顯示階躍響應(yīng)曲線title('系統(tǒng)的單位階躍響應(yīng)');ylabel('g(n)');xlabel('n');axis([0,n,1.1*min(gn),1.1*max(gn)]);5、運行例題2.3，理解卷積過程和程序中每一句的意義。clf;nf1=0:20; %f1的時間向

6、量f1=0.8.^nf1;lf1=length(f1);nf2=0:10; % f2的時間向量lf2=length(nf2); %取f2的時間向量的長度f2=ones(1,lf2);lmax=max(lf2,lf1);if lf2>lf1 nf2=0;nf1=lf2-lf1;elseif lf2<lf1 nf1=0;nf2=lf1-lf2;else nf2=0;lf1=0;endlt=lmax;

7、 u=[zeros(1,lt),f2,zeros(1,nf2),zeros(1,lt)]; 圖9t1=(-lt+1:2*lt);f1=[zeros(1,2*lt),f1,zeros(1,nf1)];hf1=fliplr(f1);N=length(hf1);y=zeros(1,3*lt);for k=0:2*ltp=[zeros(1,k),hf1(1:N-

8、k)];y1=u.*pyk=sum(y1);y(k+lt+1)=yk;subplot(4,1,1);stem(t1,u);subplot(4,1,2);stem(t1,p);subplot(4,1,3);stem(t1,y1);subplot(4,1,4);stem(k,yk);axis([-20,50,0,5]);hold onpause(1); %停頓1秒鐘end6、利用第 1 題求得的系統(tǒng)沖激響應(yīng)求解系統(tǒng)在激勵 下的響應(yīng)。 ) 3