chebyshev.m

来自「Signal and Sytem Using Matlab 书的源代码 注释详尽」· M 代码 · 共 23 行

%切比雪夫滤波器chebyshev ，page351

[z,p,k]=cheb1ap(2,3); %通带内允许有3dB波动的2几点切比雪夫滤波器
%注意cheb1lap里cheb后是一不是L
[b,a]=zp2tf(z,p,k);
w=0:0.01:4;
H=tf(b,a);
% bode(H,w);

%书上图的画法
[mag,phase] = bode(H,w);
mag = squeeze(mag);
phase = squeeze(phase);
clf
subplot(211)
plot(w,mag)
ylabel('|H|')
xlabel('Frequency (rad/sec)')
title('Example 8.27 Bandstop Filter Design')
subplot(212)
plot(w,phase)
ylabel('Angle(H), degrees')
xlabel('Frequency (rad/sec)')