数字信号处理--实验四极零点分布对系统频率响应的影响.docVIP

学生实验报告 开课学院及实验室： 电子楼317 2013 年 4 月 15 日 学院 机械与电气工程学院 年级、专业、班 姓名 学号 实验课程名称 数字信号处理实验 成绩 实验项目名称 实验四　极零点分布对系统频率响应的影响 指导老师 一、实验目的 学习用分析零极点分布的几何方法分析研究系统频率响应，理解零点、极点对系统频率响应的影响。二、实验原理 如果知道信号的Z变换以及系统的系统函数，可以得到它们的零极点分布，由零极点分布可以很方便地对它们的频率响应进行定性分析。 按照教材分析结果，信号的幅度特性由零点矢量长度之积除以极点矢量的长度之积，当频率从0变化到时，观察零点矢量长度和极点矢量长度的变化，重点观察那些矢量长度较短的情况。另外，由分析知道，极点主要影响频率响应的峰值，极点愈靠近单位圆，峰值愈尖锐；零点主要影相频率特性的谷值，零点愈靠近单位圆，谷值愈深，如果零点在单位圆上，那么频率特性为零，根据这些规律可以定性画出频率响应的幅度特性。 峰值频率和谷值频率可以近似用响应的极点和零点的相角表示，例如极点，峰值频率近似为，极点愈靠近单位圆，估计法结果愈准确。 本实验借助计算机分析信号和系统的频率响应，目的是掌握用极、零点分布的几何分析法分析频率响应，实验时需要将代入信号的Z变换和系统函数中，再在之间，等间隔选择若干点，并计算它的频率响应。 [提示：本实验可以采取两种变成方法：①先求出系统函数，再调用MATLAB函数freqz计算比绘制幅频特性和相频特性曲线；②先求出系统的函数的封闭表达式，再编程序计算其在给定离散频率点上的值，最后调用函数abs，求出模值并打印曲线。] 三、使用仪器、材料 1、硬件：计算机 2、软件：Matlab 四、实验步骤 1.　假设系统用下面差分方程描述： 假设a=0.7，0.8，0.9，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。 2.　假设系统用下面差分方程描述： 假设a=0.7，0.8，0.9，分别在三种情况下分析系统的频率特性，并打印幅度特性曲线。 3.　假设系统函数用下式描述： 试分析它的频率特性，要求打印其幅度特性曲线，并求出峰值频率和谷值频率。 五、实验过程原始记录（数据、图表、计算等） 1.%a=0.7 B=1; a=0.7;A=[1,-a]; %设置系统函数系数向量A和B subplot(2,2,1);zplane(B,A); %绘制零极点分布图 [H,w]=freqz(B,A,whole); %计算频率响应 subplot(2,2,2);plot(w/pi,abs(H)); %绘制幅频响应曲线 grid on; %网格效果 xlabel(\omega/\pi);ylabel(|H(e^j^\omega)|); subplot(2,2,4);plot(w/pi,angle(H)); %绘制相频响应曲线 xlabel(\omega/\pi);ylabel(\phi(\omega)); grid on; %网格效果 实验图像： %a=0.8 B=1; a=0.8;A=[1,-a]; %设置系统函数系数向量A和B subplot(2,2,1); zplane(B,A); %绘制零极点分布图 [H,w]=freqz(B,A, whole); %计算频率响应 subplot(2,2,2); plot(w/pi,abs(H)); grid on; %绘制幅频响应曲线 xlabel(\omega/\pi); ylabel(|H(e^j^\omega)|); subplot(2,2,4); plot(w/pi,angle(H)); %绘制相频响应曲线 xlabel(\omega/\pi); ylabel(\phi(\omega)); grid on; 实验图像： %a=0.9 B=1;a=0.9;A=[1,-a]; %设置系统函数系数向量A和B subplot(2,2,1); zplane(B,A); %绘制零极点分布图 [H,w]=freqz(B,A, whole); %计算频率响应 subplot(2,2,2); plot(w/pi,abs(H)); grid on; %绘制幅频响应曲线 xlabel(\omega/\pi); ylabel(|H(e^j^\omega)|); subplot(2,2,4); plot(w/pi,angle(H)); %绘制相频响应曲线 xlabel(\omega/\pi); ylabel(\phi(\omega)); g