第2章离散时间信号与系统的z域分析4.pptVIP

第2章 离散时间信号与系统的频域分析 系统函数、稳定性、零极点 例 求下列系统函数的零极点 解： 用一元二次方程求分母多项式的根，即极点： b=[1,0]; a=[0.8,0.23,0.15]; p=roots(a) zplane(b,a) 例2.12 判断滤波器的形状， 滤波器的系统函数为： 解：滤波器在 z=0.45处有一极点，无零点，其零极点如下图所示，并标出了 ejω 的位置。 从图中可以看到，当 ω=0 时，ejω 到极点的距离最小，而 ω 接近π时，距离最大。因为幅频响应的幅度是这个距离的倒数。所以，幅度在接近 ω=0 时最大，靠近 π 弧度时最小。因此，该滤波器具有低通特性。 例 已知因果系统 y(n) = 0.9y(n-1) + x(n) (1) 求 H(z) ，并画出零极点示意图 (2) 画出 |H(ejw)| 和 ∠H(ejw) 解：差分方程可以变形为： y(n) - 0.9y(n-1) = x(n) ? ??(1) 由差分方程可得: 系统有一个位于 0.9 的极点和一个位于原点的零点。下面用 Matlab 中的函数 zplane 画出它的零极点图。 % Chapter 2: Example 13: % zplane function: b = [1,0]; a = [1, -0.9]; zplane(b,a); title(Pole-Zero Plot); text(0.85,-0.1,0.9);text(0.01,-0.1,0); （注意：令 b=[1,0]，而不是 b=1，是因为 zplane 函数假设标量是零点或极点。） 利用 freqz 函数画出幅度响应和频率响应。 %Drawing the Frequency Response of this system %Using the first usage of freqz: freqz(b,a,N) ? [H,w]= freqz(b,a,100); magH = abs(H); phaH = angle(H); figure(2); subplot(2,1,1); plot(w/pi,magH); grid subplot(2,1,2); plot(w/pi,phaH/pi);grid 使用 freqz 函数的第三种形式的程序如下，运行结果和上图一样： % Using the third usage of freqz: freqz(b,a,w) w = [0:1:100]*pi/100; [H,w] = freqz(b,a,w); magH = abs(H);phaH = angle(H); figure(4) subplot(2,1,1);plot(w/pi,magH); grid xlabel(frequency Unit:pi); ylabel(Magnitude); title(Magnitude Response); subplot(2,1,2);plot(w/pi,phaH/pi); grid xlabel(frequency Unit:pi); ylabel(Phase Unit:pi); title(Phase Response); 信号流图 零极点以单位圆为镜对称 实系数二阶全通滤波器 两个零点（极点）共轭对称 零极点以单位圆镜像对称 图2.6.6　全通滤波器一组零极点示意 全通滤波器的系统函数一般形式如下式： （2.6.12） 二阶滤波器级联形式： （2.6.13） 上面两式中的系数均为实数。容易看出，全通滤波器系统函数H(z)的构成特点是其分子、分母多项式的系数相同，但排列顺序相反。 例 求取 的 Z 反变换。 % Check of residues in the Example b = [0,1]; a = [3,-4,1]; [R,p,k] = residuez(b,a) [b,a] = residuez(R,p,k) 运行结果如下（留数、极点、直接项以及分子分母多项式的系数）： R = 0.5000 -0.5000 p = 1.0000 0.3333 k = [] b = -0.0000 0.3333 a = 1.0000 -1.3333 0.3333 Z 反变换：Matlab 实现 根据收敛域情况，x(n) 表示为右边序列、左边序列或者双边序列。 已知离散系统的差分方程为 y(n)+0.7