实验二 时域离散系统及系统响应
一、实验目的
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(hn,'filled'); %显示冲激响应曲线 title('系统的单位冲激响应');
ylabel('h(n)');xlabel('n');
axis([0,n,1.1*min(hn),1.1*max(hn)]);
subplot(2,1,2),stem(gn,'filled'); %显示阶跃响应曲线 title('系统的单位阶跃响应');
ylabel('g(n)');xlabel('n');
axis([0,n,1.1*min(gn),1.1*max(gn)]);
系统的单位冲激响应0.80.60.40.20-0.2-0.40246101214n系统的单位阶跃响应8161820h(n)g(n)1.61.41.210246810n1214161820
1?0.5z?1（2）系统的系统函数为：H(z)?
1?z?1?z?2a=[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(hn,'filled'); %显示冲激响应曲线 title('系统的单位冲激响应');
ylabel('h(n)');xlabel('n');
axis([0,n,1.1*min(hn),1.1*max(hn)]);
subplot(2,1,2),stem(gn,'filled'); %显示阶跃响应曲线 title('系统的单位阶跃响应');
ylabel('g(n)');xlabel('n');
axis([0,n,1.1*min(gn),1.1*max(gn)]);
系统的单位冲激响应10.5h(n)0-0.5-10246101214n系统的单位阶跃响应81618201.51g(n)0.50-0.50246810n1214161820
2、运行例题2.3，理解卷积过程和程序中每一句的意义。
clf;
nf1=0:20; ?μ?ê±???òá? f1=0.8.^nf1; lf1=length(f1);
nf2=0:10; % f2μ?ê±???òá?
lf2=length(nf2); %è?f2μ?ê±???òá?μ?3¤?è f2=ones(1,lf2); lmax=max(lf2,lf1);
if lf2>lf1 nf2=0;nf1=lf2-lf1;
elseif lf2 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=length(hf1); y=zeros(1,3*lt); for k=0:2*lt