数字信号处理实验指导书(3)

三、.实验内容：

1、阅读并上机验证实验原理部分的例题程序，理解每一条语句的含义。

改变例题中的有关参数（如信号的频率、周期、幅度、显示时间的取值范围、采样点数等），观察对信号波形的影响。

2、编写程序，产生以下离散序列： （1）f(n)=δ(n) (-3

(0.1+j1.6∏)n

(0

（4）f(n)=3sin(nП/4) (0

3、一个连续的周期性方波信号频率为200Hz，信号幅度在-1～+1V之间，要求在图形窗口上显示其两个周期的波形。以4kHz的频率对连续信号进行采样，编写程序生成连续信号和其采样获得的离散信号波形。

四、.实验预习：

1、预先阅读附录部分的MATLAB基础介绍，认真阅读实验原理，明确本次实验任务，读懂例题程序，了解实验方法。

2、根据实验内容预先编写实验程序。

3、预习思考题：产生单位脉冲序列和单位阶跃序列各有几种方法？如何使用？ 五、实验报告

1、列写调试通过的实验程序，打印实验程序产生的曲线图形。

2、思考题：通过例题程序，你发现采样频率Fs、采样点数N、采样时间间隔dt在程序编写中有怎样的联系？使用时需注意什么问题？

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

一、.实验目的：

1、加深对离散系统的差分方程、单位脉冲响应、单位阶跃响应和卷积分析方法的理解。

2、初步了解用MATLAB语言进行离散时间系统时域分析的基本方法。

3、掌握求解离散时间系统的单位脉冲响应、单位阶跃响应、线性卷积以及差分方程的程序的编写方法，了解常用子函数的调用格式

二、实验原理：

1、离散LSI系统的响应与激励

由离散时间系统的时域分析方法可知，一个离散LSI系统的响应与激励可以用如下框图表示：

x[n]y[n]Discrete-timesystmeNM

其输入、输出关系可用以下差分方程描述：

?ak?0ky[n?k]??bkx[n?m]

k?02、用函数impz和dstep求解离散系统的单位脉冲响应和单位阶跃响应。

例2-1 已知描述某因果系统的差分方程为6y(n)+2y(n-2)=x(n)+3x(n-1)+3x(n-2)+x(n-3) 满足初始条件y(-1)=0，x(-1)=0，求系统的单位脉冲响应和单位阶跃响应。

解： 将y(n)项的系数a0进行归一化，得到

y(n)+1/3y(n-2)=1/6x(n)+1/2x(n-1)+1/2x(n-2)+1/6x(n-3)

分析上式可知，这是一个3阶系统，列出其bk和ak系数： a0=1， a,1=0， a,2=1/3， a,3=0 b0=1/6，b,1=1/2， b,2=1/2， b,3=1/6

程序清单如下： a=[1,0,1/3,0]; b=[1/6,1/2,1/2,1/6]; 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)]); 程序运行结果如图2-1所示：

系统的单位序列响应1.20.51.10.410.90.30.8h(n)系统的单位阶跃响应0.2g(n)0.70.60.10.500.40.30.2010n2030010n2030-0.1

图2-1

3、用函数filtic和filter求解离散系统的单位序列响应和单位阶跃响应。

例2-2 已知描述某因果系统的差分方程为6y(n)-2y(n-4)=x(n)-3x(n-2)+3x(n-4)-x(n-6)，满足初始条件y(-1)=0，x(-1)=0，求系统的单位脉冲响应和单位阶跃响应。时间轴上N取32点。

注意：原式非标准形式，必须化为标准形式后再列出系数b,a。 程序清单如下： x01=0;y01=0; a=[1,0,0,0,-1/3,0,0]; b=[1/6,0,-1/2,0,1/2,0,-1/6]; 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)]); 程序运行结果如图2-2所示：

系统的单位序列响应0.60.40.2)n(h0-0.2-0.40102030n系统的单位阶跃响应0.20.10-0.1-0.2-0.30102030ng(n)