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

1、实验二：离散LSI系统的时域分析1、 实验目的：1加深对离散系统的差分方程、单位脉冲响应、单位阶跃响应的理解。2.初步了解用MATLAB语言进行离散时间系统时域分析的基本方法。.二、实验内容：1、 已知描述某离散LSI系统的差分方程为，分别用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'

2、;); 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) x01=0;y01=0; a=1,-3/2,1/2; b=0,1/2,0; N=32;n=0

3、: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('系统的单位阶跃响应');

4、ylabel('g(n)');xlabel('n'); axis(0,N,1.1*min(gn),1.1*max(gn);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

5、(n)'); y=conv(f1,f2);subplot(2,1,2);stem(y,'filled');N=32;nt=1;n=-3:4*pi; f1=sin(n/2) subplot(2,2,1);stem(n,f1,'filled'); title('f1(n)'); n2=-3:4*pi; f2=0.5.n2 subplot(2,2,2);stem(n2,f2,'filled'); title('f2(n)'); y=conv(f1,f2); subplot(2,1,2);stem(y,'f

6、illed');3、已知某离散LSI系统的单位序列响应为求输入为时的系统响应。程序如下： n1=3:7; h(n1)=3,0.5,0.2,0.1,-0.8; subplot(2,2,1);stem(n1,h(n1),'filled'); title('h(n1)'); n=10; n2=1:n; g(n2)=exp(-0.5*n2); subplot(2,2,2);stem(n2,g(n2),'filled'); title('g(n2)'); y=conv(h(n1),g(n2); subplot(2,1,2);stem

7、(y,'filled');title('y')图形为：4、 已知描述某离散LSI系统的差分方程为，求输入为时的系统响应。程序如下：a=1,-0.7,0;b=2,0,-1;N=25;n=0:N-1;xi=filtic(b,a,0);x=n>=3;gn=filter(b,a,x,xi);stem(n,gn,'filled');title('ϵͳÏìÓ¦');ylabel('g(n)');xlabel('n');axis(0,N,1.1*min(gn),1.1*max(gn);图形为：三、实验总结：（1）通过本次试验加深了对离散系统的差分方程、单位脉冲响应、单位阶跃响应和卷积分析方法的理解。用MATLAB编程能很简单的实现系统的卷积和响应。MATLAB中丰富的函数库为我们编程提供很大的便利。 （