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

实 验 报 告

课程名称 数字信号处理 实验成绩 指导教师 王 丽 霞

院系 信息工程学院

班级 09通信一班 学号 090110009 姓名 张莹 日期 2011.12.22

实验2 离散LSI系统的时域分析 一、实验内容

1、列写调试通过的实验程序以及实验程序产生的卷积波形。 （1）f1(n=u(n),f2(n)=u(n-2), (0≤n<10) 程序清单：

n1=0:10; N1= length(n1); f1=ones(1,N1);

subplot(2,2,1);stem(n1,f1,'filled'); title('f1(n)');

n2=0:10;nt=length(n2); f2=ones(1,nt);

[y,ny]=convu(f1,n1,f2,n2); nys=0;nyf=20;

y=conv(f2,f1);ny=nys:nyf subplot(2,2,1);stem(n1,f1); subplot(2,2,2);stem(n2,f2); subplot(2,1,2);stem(ny,y);

卷积波形：

10.80.60.40.20051010.80.60.40.20051015105002468101214161820

（2）x(n)=sin(n/2),h(n)=(0.5)n (-3≤n≤4П) 程序清单： n1=-3:4*pi; a=1/2; b=0.5; f1=sin(a*n1);

subplot(2,2,1);stem(n1,f1,'filled'); title('f1(n)'); f2=b.^n1;

subplot(2,2,2);stem(n1,f2,'filled'); title('f2(n)'); y=conv(f1,f2);

subplot(2,1,2);stem(y,'filled');

卷积波形：

f1(n)10.50-0.5-1-520100-10-2086420-5f2(n)05101505101505101520253035

2、已知某离散LSI系统的单位序列响应为h(n)=3δ(n-3)+0.5δ(n-4)+0.2δ(n-5)+0.7δ(n-6)-0.8δ(n-7) 求输入为x(n)=e-0.5nu(n)时的系统响应。 程序清单： N=16; n=0:N-1; x=exp(-0.5*n)

subplot(2,2,1);stem(n,x); a=[1,0,0,0,0,0,0,0]; b=[0,0 ,0,3,0.5,0.2,0.7,-0.8]; hn=impz(b,a,n);

subplot(2,2,2);stem(n,hn); y=conv(x,hn); subplot(2,1,2);stem(y);

10.80.63210.40.200510150-10510153210-105101520253035

3、已知描述某离散LSI系统的差分方程为y(n)=0.7y(n-1)+2x(n)-x(n-2)，求输入为x(n)=u(n-3)时的系统响应。

程序清单： N=16; n=0:N-1;

x=[zeros(1,3),ones(1,(N-3))]; subplot(2,2,1);stem(n,x); a=[1,-7/10,0]; b=[2,0,-1];

hn=impz(b,a,n);

subplot(2,2,2);stem(n,hn) y=conv(x,hn);

subplot(2,1,2);stem(y);

10.80.60.40.2005101521.510.50-0.505101543210-105101520253035

二、列出本实验提出的有关MATLAB函数在调用时应注意的问题。

答：1）.MATLAB提供的求卷积函数conv默认两个序列的序号均从n=0开始，卷积结果y对应的序列的序号也从n=0开始。当两个序列不是从0开始时，必须对conv函数稍加扩展。 2）.用函数filtic和filter求解离散系统的单位序列响应和单位阶跃响应时，原式为非标准形式，必须化为标准形式后再列出系数矩阵a，b 。

本文来源：https://www.bwwdw.com/article/5a8t.html