云南师范大学dsp(数字信号处理)试验报告资料.doc

本科学生实验报告 学号 114090395 姓名 李华福 学院 物电学院 专业、班级 11电子 实验课程名称 数字信号处理（实验） 教师及职称 李宏宁老师 开课学期 2013 至 2014 学年 下 学期 填报时间 2014 年 5 月 7 日 云南师范大学教务处编印 实验序号 8 实验名称 离散系统分析 实验时间 2014年5月7日 实验室 数字信号处理实验室（313） 一．实验预习 1．实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频率特性及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频率响应和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 2．实验原理、实验流程或装置示意图 实验原理： MATLAB提供了许多可以用于分析线性时不变离散系统的函数，主要包含有系统函数、系统时域响应、系统频率响应等分析函数。 （1）离散系统的时域响应 描述离散LTI系统的输入x[k]与输出y[k]关系的差分方程为： (2.2.1) 对上式两边同取z变换可得离散系统的系统函数H(z)为 (2.2.2) 使用向量b=[]和向量分别保存分子多项式和分母多项式的系数，这些系数都是从按z的降幂排列。 离散系统的单位冲激响应h[k]的计算 [h,k]=impz(b,a)计算系统的单位脉冲响应的相应的时间轴向量，其中b和a分别是系统函数H(z)的分子多项式和分母多项式的系数矩阵；也可以简写为h=impz(b,a)。 [h,k]=impz(b,a,n)计算n点单位脉冲响应；也可以简写为h=impz(b,a,n)。 Impz(b,a)自动绘制单位脉冲响应的图形。 2）离散系统响应y[k]的计算 y=filter(b,a,x)计算系统在输入x作用下的零状态响应为y。 y=filter(b,a,x，zi)计算系统在输入x和初始状态作用下的完全响应y。zi是由系统的初始状态经过filtic函数转换而得到的初始条件：zi=filtic(b,a,Y0),其中Y0为系统的初始状态，Y0=[y[-1],y[-2],y[-3],…]。 也可利用lsim函数求解离散系统的系统响应，该函数在连续系统分析试验中介绍过，在此不再重复。 [例2.2.1]已知系统函数为，求 离散系统的单位脉冲响应h[k]; 输入x[k]=u[k]，求系统的零状态响应y[k]; 输入就x[k]=u[k]，初始条件y[-1]=1,y[-2]=2，求系统的完全响应y[k]。 【解】 计算前40个点的单位脉冲响应h[k] N=40; a=[1,0.4,-0.12];b=[1,2]; y=impz(b,a,N); stem(y) xlabel(k);title(h[k]) 运行结果如下图2.1.1所示。 图2.2.1 系统的单位脉冲响应 计算前100个点的零状态响应y[k] N=40; a=[1,0.4,-0.12];b=[1,2]; x=ones(1,N); y=filter(b,a,x) 计算前100个时刻的完全响应y[k] Filter(b,a,x,zi)中的初始值zi不是y[-1]=1,y[-2]=2，它可以由filtic函数计算。 N=100; a=[1,0.4,-0.12];b=[1,2]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); （2）离散系统的系统函数零极点分析 LTI离散系统的系统函数H(z)可以表示为零极点形式，即 (2.2.3) 可以使用MATLAB提供的roots函数计算离散系统的零极点，使用zplane函数绘制离散系统的零极点分布图。在利用这些函数时，要求H(z)的分子多项式和分母多项式的系数的个数相等，若不相等需要补零。 [例2.2.2]已知系统函数为计算该系统函数的的零极点，并画出系统函数零极点分布图。 【解】 b=[1,2,0];a=[1,0.4,-0.12]; z=roots(b); p=roots(a); zplane(b,a) 程序运行输出零点为 z = 0 -2 输出的极点为 p = -0.6000 0.2000 绘出的零极点分布图如图2.2.2所示。可看见系统函数含有两个零点，两个极点。极点都位于单位圆内，因此该离散系统稳定。 图2.2.2 系统函数的零极点分布图 （3）离散系统的频率响应 当离散系统因果LTI系统的系统函数H(z)的极点全部位于z平面单位圆