实验2:离散LSI系统的时域分析

实验二:离散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');

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4