实验六离散系统分析.docVIP

本科学生验证性实验报告 学号 114090413 姓名 罗朝斌 学院 物理与电子信息学院专业、班级11光电子 实验课程名称 离散信号分析 教师及职称 杨卫平 （教授） 开课学期 2013 至 2014 学年 第二 学期 填报时间 2014 年 4 月 21 日 云南师范大学教务处编印 实验序号 实验六 实验名称 离散系统分析 实验时间 2014-4-21 实验室 同析3栋（313）数字信号处理实验室 一．实验预习 一实验目的 深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频率特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB分析离散系统的时域响应、频率响应和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 二实验原理、实验流程或装置示意图 MATLAB提供了许多可用于分析线性时不变离散系统的函数，主要包含有系统函数、系统时域响应、系统频率响应等分析函数。 离散系统的时域响应? 描述离散LTI系统输入x[k]与输出y[k]关系的差分方程为 （2.2.1） 对上式两边同取z变换可得离散系统的系统函数为 （2.2.2） 使用向量和向量分别保存分子多项式和分母多项式的系数，这些系数都是从z0按z的降幂排列。 （1）离散系统的单位脉冲响应的计算 计算系统的单位脉冲响应和相应的时间轴向量，其中b和a分别是系统函数的分子多项式和分母多项式的系数矩阵；也可简写为。 计算n点单位脉冲响应；也可简写为。 自动绘制单位脉冲响应的图形。 （2）离散系统响应的计算 计算系统在输入x作用下的零状态响应y。 计算系统在输入x和初始状态作用下的完全响应y。zi是由系统的初始状态经过filtic函数转换而得到的初始条件：，其中Y0为系统的初始状态，。 也可以利用lsim函数求解离散系统的系统响应，该函数在连续系统分析实验中介绍过，在此不再重复。 【例2.2.1】 已知系统函数为，求 （1）离散系统的单位脉冲响应； （2）输入，求系统的零状态响应； （3）输入，初始条件，，求系统的完全响应。 [解]（1）计算前40个点的单位脉冲响应 N=40; a=[1,0.4,-0.12];b=[1,2]; y=impz(b,a,N); stem(y) xlabel(‘k’);title(‘h[k]’) 运行结果如图2.2.1所示。 图2.2.1 系统的单位脉冲响应 （2）计算前100个点的零状态响应 N=100; b=[1,2];a=[1,0.4,-0.12]; x=ones(1,N); y=filter(b,a,x) 前十个数据： 1 2.6 2.08 2.48 2.2576 2.39456 2.31308 2.362112 22(3)计算前100个时刻的完全响应 中的初始值zi不是，，它可以由filtic函数计算。 N=100; b=[1,2];a=[1,0.4,-0.12]; x=ones(1,N); zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); 前十个数据： 0.84 2.784 1.9872 2.5392 2.222784 2.4155904 22.3696716 222．离散系统的系统函数零极点分析 LTI离散系统的系统函数可以表示为零极点形式，即 （2.2.3） 可以使用MATLAB提供的roots函数计算离散系统的零极点，使用zplane函数绘制离散系统的零极点分布图。在利用这些函数时，要求的分子多项式和分母多项式的系数的个数相等，若不等则需要补零。 【例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系统的系统函数的极点全部位于z平面单位圆内时，系统的频率特性可由求出，即 （2