利用MATLABFDATOOL观察零极点.docVIP

一、实验目的 1.利用MATLAB的fdatool观察传递函数H(z)的零极点分布、幅度响应、脉冲响应及阶跃响应，并观察零极点分布与系统稳定性之间的关系。 2.利用MATLAB的fdatool验证有关稳定性的三个重要结论。 3.利用MATLAB绘制传递函数的三维图像，更直观的查看传递函数在实轴与虚轴中的变化情况。 二、实验条件 PC机，MATLAB7.0 实验内容 实验前准备： 打开MATLAB中的fsatool，将相关参数设置好。 在Numerator中输入H(z)中分子的系数（注意：系数按降幂顺序输入），在Denominator中输入H(z)中分母的系数。点击Import Filter观察图像。 （一）观察不同传递函数H（z)之间相关响应图像之间的差别 a. 由幅度响应的图像可发现，该系统为高通滤波器。 b. 由幅度响应的图像可发现，该系统为高通滤波器。 c. 由幅度响应的图像可发现，该系统为低通滤波器。 d. 由幅度响应的图像可发现，该系统为低通滤波器。 e. 由幅度响应的图像可发现，该系统为低通滤波器。 f. 由幅度响应的图像可发现，该系统为高通滤波器。 g. 由幅度响应的图像可发现，该系统为高通滤波器。 h. 由幅度响应的图像可发现，该系统为带通滤波器。 实验结论：通过观察极点在坐标系中的位置，可以得出以下两个结论 1.极点越靠近原点，系统收敛的越快。2.极点在z轴的左半平面，则相应在正负轴间交替。 为了更直观的观察零极点的位置与响应之间的关系，可以直接在零极点图中拖拽极点改变极点的位置，直接观察冲击响应或阶跃响应的变化。以为例来进行验证。 将极点向原点拉近，可以发现脉冲响应的收敛变快。 将极点拉向z轴的左半平面，可以发现脉冲响应开始正负交替。 利用验证零点与极点之间距离与脉冲响应之间的关系。 实验结论：由图中的横坐标可以看出，当零点与极点越接近，脉冲响应衰减的越快。 （二）绘制三维的传递函数图像 程序代码： clear all t=-4:0.05:4; x=cos(pi/2*t); y=sin(pi/2*t); z=x+1i.*y; H=(0.5-0.1.*z.^-1-0.04.*z.^-2)./(1+1.4.*z.^-1+0.85.*z.^-2); stem3(x,y,H) 运行程序后出现图像如下： 图为XYZ轴的三维图像 图为通过旋转，XZ轴的图像 由图可以看出，x,y分别表示实轴与 虚轴，并且在复平面形成了单位圆。 而Z轴则表示了传递函数H（z)在复 平面上的变化情况。并且根据x,y的 范围，表示了传递函数H（z)在复平 面单位圆上的变化情况。 图为XY轴的三维图像 图为该函数通过滤波器工具所显示的频响特性和零极点分布图 实验中遇到的问题及解决： 在输入传递函数数值的时候（例如）应按照降幂次序输入 0 0 1，而不是输入1。 验证零点与极点之间的距离与脉冲响应之间的关系时，使用系统移动极点时，出现的脉冲响应图片中会有星号。 后面询问后知道这是因为零点不共轭所产生的星号，而这种系统在物理上是不可实现的系统。 （3）利用fdatool上方的窗口直接观看脉冲和阶跃响应时，调整零极点的位置图像会跟着相应变化。但是频繁调整后会发现显示出来的响应图片会有延迟导致前一刻的图片与现在的图片重叠在一起，不方便我们观察。而点击会弹出单独的窗口来显示我们所需要的图片，并且图片中没有延迟，较为清晰。 四、实验结论和讨论 1.掌握了MATLAB中的fdatool工具，实现了模拟滤波器的功能，不用写程序就能观察到不同传递函数的零极点图、幅度响应和脉冲响应等相关图像，十分便捷快速。 2.通过改变零极点图中的零点或极点，观察对