[Hb1,w] = freqz (b1, 1); % 分母多项式只有常数项1，没指定返回点数，缺省512点
[Hb2] = freqz (b2); % 这里把分母项省略了，返回的点数和上行一样，所以不再需要 w Hb3 = freqz (b3); if useplot
% 用一句 plot 画3条曲线，第2条用绿色，虚线，第3条红色；grid 画出网格线 subplot (2, 1, 1)
plot (w/pi, abs(Hb1), w/pi, abs(Hb2), 'g--', w/pi, abs(Hb3), 'r'), grid subplot (2, 1, 2)
%plot(w/pi, angle(Hb), w/pi, angle(Hbb), 'r--'), grid
plot(w/pi, unwrap(angle(Hb1)), w/pi, unwrap(angle(Hb2)), 'g--', w/pi, unwrap(angle(Hb3)), 'r') grid
else % 也可以用下面的来显示
% Hb, Hbb 是和 w 同维数的列向量(从Wokspace里可以看到) % H 是 2 列合成的矩阵 H = [Hb1, Hb2, Hb3];
% 用不同的方式显示，如果三种都要同时看，你还得加几句 % 例如用 pause 在时间上隔开，或是用3个窗口 figure (1) freqzplot (H, w) figure (2)
freqzplot (H, w, 'linear') figure (3)
freqzplot (H, w, 'squared') end
- 13 -