samp5_17.m

MATLAB7.x数字信号处理 光盘内容

%Samp5_17
ws=[1000 2000]*2*pi;wp=[500 2500]*2*pi;Rp=1;Rs=50;  %滤波器设计参数
[N,Wn]=cheb2ord(wp,ws,Rp,Rs,'s');   %求滤波器的最小阶数和截止频率
[b,a]=cheby2(N,Rs,Wn,'stop','s');  %设计带阻滤波器
w=linspace(1,3000,1000)*2*pi;   %设置计算复数频率响应的频率点
H=freqs(b,a,w);   %计算给定频率点的复数频率响应
magH=abs(H);phaH=unwrap(angle(H));  %计算幅频响应和相频响应
plot(w/(2*pi),20*log10(magH));  %以频率为横轴绘制幅频响应
xlabel('频率/Hz');ylabel('振幅/dB');
title('Chebyshev II 模拟带阻滤波器')
hold on; plot([1000 1000],ylim, 'r');plot([2000 2000],ylim, 'r');%阻带边界
grid on