武汉工程大学数字信号处理实验二时域离散系统及系统响应0001.docx

实验二时域离散系统及系统响应 一、实验目的 1、 掌握求解离散时间系统冲激响应和阶跃响应的方法； 2、 进一步理解卷积定理，掌握应用线性卷积求解离散时间系统响应的基本方法; 3、 掌握离散系统的响应特点。 二、实验内容 1、请分别用impz和dstep函数求解下面离散时间系统的冲激响应和阶跃响应 (1)系统的差分方程为: y(n) =0.8y(n 一 1) _0.64y(n 2) 0.866x(n) a=[1,-0.8,0.64]; b=[0.866,0,0]; n=20; hn=impz(b,a, n); %冲激响应 gn=dstep(b,a ,n); %阶跃响应 subplot(2,1,1),stem(h n, filled ); %显示冲激响应曲线 title(系统的单位冲激响应’)； ylabel( h( n) );xlabel( n); axis([0, n,1.1*mi n(h n),1.1*max(h n)]); subplot(2,1,2),stem(g n, filled ); %显示阶跃响应曲线 title(系统的单位阶跃响应); ylabel( g(n) );xlabel( n); axis([0, n,1.1*mi n(g n),1.1*max(g n)]); 系统的单位冲激响应 系统的单位阶跃响应 1 -0.5z， 1 -z」Z (2)系统的系统函数为: H(z)二 a=[1,-1,1]; b=[1,-0.5,0]; n=20; hn=impz(b,a, n); %冲激响应 );%显示冲激响应曲线gn=dstep(b,a ,n); );%显示冲激响应曲线 subplot(2,1,1),stem(h n, filled title(系统的单位冲激响应’)； ylabel( h(n) );xlabel( n); );% );%显示阶跃响应曲线 title(系统的单位阶跃响应); ylabel( g(n) );xlabel( n); axis([0, n,1.1*mi n(g n),1.1*max(g n)]); 系统的单位冲激响应 系统的单位阶跃响应 1.5 1 0.5 0 -0.5 10 12 14 16 18 20 n 2、运行例题2.3,理解卷积过程和程序中每一句的意义。 elf; nf1=0:20; %f1y ?e± ???da? f1=0.8.A nf1; lf1=le ngth(f1); nf2=0:10; % f2 卩?e± ???da? lf2=length(nf2); %e ?f2 卩?e± ???oa ?卩?3O?e f2=o nes(1,lf2); Imax=max(lf2,lf1); if lf2lf1 n f2=0; nf1=lf2-lf1; elseif lf2lf1 n f1=0; nf2=lf1-lf2; else n f2=0;lf1=0; end lt=lmax; u=[zeros(1,lt),f2,zeros(1, nf2),zeros(1,lt)]; t1= (-lt+1:2*lt); f1=[zeros(1,2*lt),f1,zeros(1, nf1)]; hf1= fliplr(f1); N=le ngth(hf1); y=zeros(1,3*lt); for k=0:2*lt p=[zeros(1,k),hf1(1:N-k)]; y1=u.*p yk=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 on pause(1); %i￡ ?u 1???o end -201-1010203040500-101020304050IT诂一-1010203040500「 -20 -20 1 -10 10 20 30 40 50 0 -10 10 20 30 40 50 IT诂一 -10 10 20 30 40 50 0「 -20 -1 -20 5「 3、利用第1题求得的系统冲激响应求解系统在激励 x(n)二u(n -3)下的响应 a=[1,-0.8,0.64]; b=[0.866,0,0]; yi=[0 0]; xi=0; xic=filtic(b,a,yi,xi); n=20; hn=impz(b,a,n); %3?? ? io| subplot(2,1,1),stem(h n, filled ); %?望？3??O ? io |? u ?? title( ? 3 卩?^￥ ??3??O ?i 0; | yl