实验五零极点实验及其频率响应.docVIP

实验五 零极点实验及其频率响应 【实验目的】 通过matlab仿真简单的离散时间系统，研究其频率域特性，加深对离散系统的冲激响应、频率响应分析和零、极点分布概念理解。 【实验原理】 【实验内容及要求】 设离散时间系统对应的输入输出差分方程为 y n -0.4y n-1 -0.5y n-2 0.2x n +0.1x n-1 求系统的频率响应。 画出系统的零极点图。 解答：（1） fs 1000; b [0.2 0.1]; a [1 -0.4 -0.5]; [h,f] freqz b,a,256,fs ; mag abs h ; ph angle h ; ph ph*180/pi; subplot 2,1,1 ; plot f,mag ; grid on xlabel 频率（HZ） ; xlabel 幅度 ; subplot 2,1,2 ; plot f,ph ; grid on xlabel 频率 HZ ; xlabel 相位 ; 2 b [0.2 0.1]; a [1 -0.4 -0.5]; zplane b,a ; grid on 3 b [0.2 0.1]; a [1 -0.4 -0.5]; [z,p,k] tf2zp b,a zplane b,a ; grid on 结果：z -0.5000 p 0.9348 -0.5348 k 0.2000 一个LTI离散时间系统的输入输出差分方程为 Y n -1.6y n-1 +1.28y n-2 0.5x n +0.1x n-1 （1）编程得到系统频率响应的幅度响应和相位响应，并画图。 （2）编程得到系统的零极点分布图。 解：（1） fs 1000; b [0.5 0.1]; a [1 -1.6 1.28]; [h,f] freqz b,a,256,fs ; mag abs h ; ph angle h ; ph ph*180/pi; subplot 2,1,1 ; plot f,mag ; grid on xlabel 频率（HZ） ; xlabel 幅度 ; subplot 2,1,2 ; plot f,ph ; grid on xlabel 频率 HZ ; xlabel 相位 ; （2） b [0.5 0.1 0]; a [1 -1.6 1.28]; subplot 2,1,1 ; zplane b,a ; grid on [z,p,k] tf2zp b,a subplot 2,1,2 ; zplane b,a ; grid on 结果：z 0 -0.2000 p 0.8000 + 0.8000i 0.8000 - 0.8000i k 0.5000 系统零极点分布与系统频率响应的关系：单位圆附近的零点对幅度响应的谷点的位置与深度有明显影响，当零点位于单位圆上时，谷点为；单位圆附近的极点对幅度响应的峰点位置和高度有明显影响。