[精选]《数字信号处理》实验指导书(完整)

《[精选]《数字信号处理》实验指导书(完整)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《数字信号处理》实验指导书通信教研室安阳工学院二零零九年三月第1章系统响应及系统稳定性1.1实验目的•学会运用MATLAB求解离散时间系统的零状态响应；•学会运川MATLAB求解离散时间系统的单位取样响应；•学会运用MATLAB求解离散时间系统的卷积和。1.2实验原理及实例分析1.2.1离散时间系统的响应离散时间LTI系统可用线性常系数差分方程来描述，即NM工吋(〃-0=工bjX(n-j)(1-1)z=Oj=0其中，a.(z=0,1,N)和巧(j=0,1,...»M)为实常数。MATLAB中函数filter可对式(13・1)的差分方程

2、在指定时间范围内的输入序列所产生的响应进行求解。函数filter的语句格式为y=filter(b,a,x)其中，x为输入的离散序列；y为输出的离散序列；y的长度Lx的长度一样；b与a分别为差分方程右端与左端的系数向量。【实例1・1】已知某LTI系统的差分方程为3y(n)-4y(z?-1)+2y(n一2)=x(n)+2x(n一1)试用MATLAB命令绘出当激励信号为兀(/!)=(1/2)5®)时,该系统的零状态响应。解：MATLAB源程序为»a=[3-42];»b=[l2];»n=0:30;»x=(l/2).An;»y=filter(b

3、,a,x);»stem(n,y/fiir),gridon>>xlabclCnJ’titleC系统响应y(n)J程序运行结果如图1・1所示。系统响应y(n)图1・1实例1・1系统的零状态响应1.2.2离散时间系统的单位取样响应系统的单位取样响应定义为系统在55)激励下系统的零状态响应，用/?(7?)表示。MATLAB求解单位取样响应可利用函数filter,并将激励设为前面所定义的impDT函数。例如，求解实例1・1中系统的单位取样响应吋，MATLAB源程序为»a=[3-42];»b=[l2];»n=0:30;»x=impDT(n);»h

4、=filter(b,a,x);>>stem(n,h,'fill'),gridon»xlabcl(,n,),titlc(,系统单位取样响应h(n)J程序运行结果如图1・2所示。图1-2实例的系统单位取样响应MATLAB另一种求单位取样响应的方法是利用控制系统工具箱捉供的函数impz来实现。impz函数的常用语句格式为impz(b,a,N)其中，参数N通常为正整数，代表计算单位取样响应的样值个数。【实例1・2】已知某LTI系统的差分方程为3y(n)一4y(/i-1)+2y(n-2)=x(n)+2x(n一1)利用MATLAB的impz函数绘

5、出该系统的单位取样响应。解：MATLAB源程序为»a=[3-42];»b=[l2];»n=0:30;»impz(b,a,30),gridon»title(,系统单位取样响应h(n)*)程序运行结果如图1・3所示，比较图1・2和图1・3,不难发现结果相同。o6■o系统卑位取样响应h(n)510152025n(samples)图1-3系统单位取样响应28o.642-0bbb€prn=dE<24001.2.3离散时间信号的卷积和运算由丁系统的零状态响应是激励为系统的单位取样响应的卷积，因此卷积运算在离散时间信号处理领域被广泛应用。离散时间信

6、号的卷积定义为8y(n)=x(n)*h(n)=为x(m)h(n-m)(1-2)W=-00可见，离散时间信号的卷积运算是求和运算，因而常称为“卷积和MATLAB求离散时间信号卷积和的命令为conv,-K语句格式为y=conv(x,h)其中，x与h表示离散时间信号值的向量；y为卷积结果。用MATLAB进行卷积和运算时，无法实现无限的累加，只能计算时限信号的卷积。例如，利用MALAB的conv命令求两个长为4的矩形序列的卷积和，即g(n)=[u⑷一u(n-4)]*[u(n)-u(n-4)],其结果应是长为7(4+4-1=7)的三角序列。用向

7、量[1111]表示矩形序列，MATLAB源程序为»xl=[l111];»x2=[l111];»g=conv(xl,x2)1234321如果要绘出图形来，则利用stem命令，即»n=l:7;»stem(n,g/fill,),gridon,xlabelC『)程序运行结果如图1・4所示。图1・4卷积结果图对于给定函数的卷积和，我们应计算卷积结果的起始点及其长度。两个时限序列的卷积和长度一般等于两个序列长度的和减lo【实例1・3】已知某系统的单位収样响应为/?(H)=0.8,,[w(n)-w(n-8)]，试用MATLAB求当激励信号为x(n)

8、=u(n)-u(n-4)时，系统的零状态响应。解：MATLAB中可通过卷积求解零状态响应,即x(n)*h(n).由题意可知,描述/2(砒向量的长度至少为8,描述兀⑺)向量的长度至少为4,因此为了图形完整美观，我们将/?(