fig5_16.m

来自「数字信号处理Matlab演示文件,其中各个文件加放置了不同的matlab子文件」· M 代码 · 共 21 行

echo on;
% Fig5_16.m
% Comparison of low-pass,band-pass,high-pass
% designs using Chebychev prototype filter
%
echo off;
[B1,A1] = cheby1(2,0.5,1/3);
[B2,A2] = cheby1(2,0.5,2/3,'high');
[B3,A3] = cheby1(2,0.5,[1/3,2/3]);
[H1,W1] = freqz(B1,A1,512);
[H2,W2] = freqz(B2,A2,512);
[H3,W3] = freqz(B3,A3,512);
H1(512) = 1e-100;
H2(1) = 1e-100;
H3(1) = 1e-100;
H3(512) = 1e-100;
W = 0:9000/512:9000*(511/512);
plot(W,20*log10(abs(H1)),W,20*log10(abs(H2)),W,20*log10(abs(H3)));
axis([0 9000 -40 10]);
pause;