第10章离散系统的频率响应.ppt

一、实验目的　　(1)加深对离散系统的频率响应特性基本概念的理解。　　(2)了解离散系统的零极点与频响特性之间关系。　　(3)熟悉MATLAB中进行离散系统分析频响特性的常用子函数，掌握离散系统幅频响应和相频响应的求解方法。 二、实验涉及的MATLAB子函数　　1.freqz　　功能：用于求解离散时间系统的频率响应函数H(ejw)。　　调用格式：　　［h，w］＝freqz(b，a，n)；可得到数字滤波器的n点复频响应值，这n个点均匀地分布在［0，p］上，并将这n个频点的频率记录在w中，相应的频响值记录在h中。缺省时n＝512。 　　［h，f］＝freqz(b，a，n，Fs)；用于对H(ejw)在［0，Fs/2］上等间隔采样n点，采样点频率及相应频响值分别记录在f和h中。由用户指定Fs(以Hz为单位)的值。　　h＝freqz(b，a，w)；用于对H(ejw)在［0，2p］上进行采样，采样频率点由矢量w指定。　　h＝freqz(b，a，f，Fs)；用于对H(ejw)在［0，Fs］上采样，采样频率点由矢量f指定。　　freqz(b，a，n)；用于在当前图形窗口中绘制幅频和相频特性曲线。 　　2.angle　　功能：求相角。　　调用格式：　　p＝angle(h)；用于求取复矢量或复矩阵H的相角(以弧度为单位)，相角介于－p和＋p之间。 　　3.grid　　功能：在指定的图形坐标上绘制分格线。　　调用格式：　　grid 紧跟在要绘制分格线的绘图指令后面。例如：plot(t，y)；grid。　　grid on 绘制分格线。　　grid off 不绘制分格线。 　　4.hold　　功能：在当前轴或图形上多次叠绘多条曲线。　　调用格式：　　hold 使当前图形具备刷新性质的双向开关。　　hold on 使当前轴或图形保持而不被刷新，准备接受此后将绘制的新曲线。　　hold off 使当前轴或图形不再具备不被刷新的性质。 　　5.text　　功能：在图形上标注文字说明。　　调用格式：　　Text(xt，yt，￠string￠)；在图面上(xt，yt)坐标处书写文字说明。其中文字说明字符串 必须使用单引号标注。 三、实验原理　　1.离散系统频率响应的基本概念　　已知稳定系统传递函数的零-极点增益(zpk)模型为则系统的频响函数为 其中，系统的幅度频响特性为　　系统的相位频响特性为 　　由公式可见，系统函数与频率响应有着密切的联系。适当地控制系统函数极点、零点的分布，可以改变离散系统的频率响应特性：　　(1)在原点(z＝0)处的零点或极点至单位圆的距离始终保持不变，其值|ejw|＝1，所以对幅度响应不起作用。　　(2)单位圆附近的零点对系统幅度响应的凹谷的位置及深度有明显的影响。　　(3)单位圆内且靠近单位圆附近的极点对系统幅度响应的凸峰的位置及峰度有明显的影响。 　　2.系统的频率响应特性　　MATLAB为求解离散系统的频率响应和连续系统的频率响应，分别提供了freqz和freqs两个函数，使用方法类似。本实验主要讨论离散系统的频率响应。　　例10-1 已知离散时间系统的系统函数为　　求该系统在0～p频率范围内的相对幅度频率响应与相位频率响应。 　　MATLAB程序如下：　b＝［0.1321，0，－0.3963，0，0.3963，0，－0.1321］；　a＝［1，0，0.34319，0，0.60439，0，0.20407］；　freqz(b，a)；　　以上程序采用了freqz不带输出向量的形式，直接出图。执行结果如图10-1所示。　　由图10-1可见，该系统是一个IIR数字带通滤波器。其中幅频特性采用归一化的相对幅度值，以分贝( dB)为单位。 　　 图10-1 例10-1系统的幅度频率响应与相位频率响应 　　例10-2 已知离散时间系统的系统函数，求该系统在0～p频率范围内归一化的绝对幅度频率响应与相位频率响应。　　解 MATLAB程序如下：　　b＝［0.2，0.1，0.3，0.1，0.2］；　　a＝［1，－1.1，1.5，－0.7，0.3］；　　n＝(0：500)*pi/500； %在pi的范围内取501个采样点 　　［h，w］＝freqz(b，a，n)；%求系统的频率响应　　subplot(2，1，1)，plot(n/pi，abs(h))；grid%作系统的幅 度频响图　　axis(［0，1，1.1*min(abs(h))，1.1*max(abs(h))］)；　　ylabel(￠幅度￠)；　　subplot(2，1，2)，plot(n/pi，angle(h))；gr