matlab zplane多画出零点,上机实验2 零极点分析与幅频特性.ppt

上机实验2 零极点分析与幅频特性

L o g o Digital Signal Processing, * 实验2 离散事件系统频域分析 数字信号处理实验 Digital Signal Processing * 频率响应的 Matlab 实现 在 Matlab 中，提供了一个 zplane函数 、freqz 函数来计算幅度和相位响应，该函数有如下几种调用方式。 1)[H,w] = freqz(b,a,N) b 和 a 分别表示分子和分母的系数向量，与 filter(b,a,x) 函数中的相同。此函数在上半个单位园 [0，π) 上等间距的 N 点计算频率响应，返回该系统的 N 点频率矢量 ω 和对应的 N 点复数频率响应矢量 H。 2)[H,w] = freqz(b,a,N,’whole’) 在整个单位园上等间距的 N 点计算频率响应。 3)H=freqz(b,a,w) 它返回矢量 ω 指定的那些频率点上的频率响应，通常在 0 到 π 之间。 系统函数、稳定性、零极点 * 例2.14 已知因果系统 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 画出它的零极点图。 2.8 系统函数、稳定性、零极点 * % 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 函数假设标量是零点或极点。) 2.8 系统函数、稳定性、零极点 * 利用 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(b,a,N)，沿单位园的上半圆取 100 个点。 2.8 系统函数、稳定性、零极点 * 这种情况下计算区间是 0≤w≤0.99π，而在 w=π 点不确定。这是由于在 Matlab 中，单位园的下半圆是从 w=π 开始的。 为了克服这个问题，我们使用 freqz 函数的第二种形式 freqz(b,a,N,’whole’)，程序如下： ? % Using the second usage of freqz: freqz(b,a,'whole') [H,w]=freqz(b,a,200,'whole'); magH = abs(H(1:101)); phaH = angle(H(1:101)); ? w = w(1:101); figure(3); subplot(2,1,1);plot(w/pi,magH);grid subplot(2,1,2);plot(w/pi,phaH/pi);grid 2.8 系统函数、稳定性、零极点 * 使用 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 xla