实验二：离散LSI系统的时域分析
一、实验目的：
1加深对离散系统的差分方程、单位脉冲响应、单位阶跃响应的理解。 2.初步了解用MATLAB语言进行离散时间系统时域分析的基本方法。. 二、实验内容：
1、已知描述某离散LSI系统的差分方程为2y(n)?3y(n?1)?y(n?2)?x(n?1)，分别用impz和dstep函数、filtic和filter函数两种方法求解系统的单位序列响应和单位阶跃响应。 用impz和dstep函数 程序如下：
a=[1,-3/2,1/2]; b=[0,1/2,0]; N=32; n=0:N-1;
hn=impz(b,a,n); gn=dstep(b,a,n);
subplot(1,2,1);stem(n,hn,'k'); title('系统的单位序列响应'); ylabel('h(n)');xlabel('n');
axis([0,N,1.1*min(hn),1.1*max(hn)]); subplot(1,2,2);stem(n,gn,'k'); title('系统的单位阶跃响应'); ylabel('g(n)');xlabel('n');
axis([0,N,1.1*min(gn),1.1*max(gn)
系统的单位序列响应系统的单位阶跃响应10.90.80.7302520h(n)0.50.40.30.20.10010n2030g(n)0.6151050010n2030
x01=0;y01=0; a=[1,-3/2,1/2]; b=[0,1/2,0]; N=32;n=0:N-1;
xi=filtic(b,a,0); x1=[n==0];
hn=filter(b,a,x1,xi); x2=[n>=0];
gn=filter(b,a,x2,xi);
subplot(1,2,1);stem(n,hn,'k'); title('系统的单位序列响应'); ylabel('h(n)');xlabel('n');
axis([0,N,1.1*min(hn),1.1*max(hn)]); subplot(1,2,2);stem(n,gn,'k'); title('系统的单位阶跃响应'); ylabel('g(n)');xlabel('n');
axis([0,N,1.1*min(gn),1.1*max(gn)]);
系统的单位序列响应系统的单位阶跃响应10.90.80.7302520h(n)0.50.40.30.20.10010n2030g(n)0.6151050010n2030
2、编写程序描绘下列序列的卷积波形： n1=0:10;
N1=length(n1); f1=[n1>=2];
subplot(2,2,1);stem(n1,f1,'filled'); title('f1(n)'); n2=0:10;
N2=length(n2); f2=ones(1,N2);
subplot(2,2,2);stem(n2,f2,'filled'); title('f2(n)'); y=conv(f1,f2);
subplot(2,1,2);stem(y,'filled');