实验二 时域离散系统及系统响应

一、实验目的

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