4-26.m

matlab数字信号与应用源代码4

%例程4－26    基于切比雪夫一致逼近法设计FIR数字低通滤波器

clear all;               
f=[0 0.6 0.7 1];             %给定频率轴分点 
A=[1 1 0 0];                %给定在这些频率分点上理想的幅频响应
weigh=[1 10];                %给定在这些频率分点上的加权
b=remez(32,f,A,weigh);       %设计出切比雪夫最佳一致逼近滤波器
[h,w]=freqz(b,1,256,1);
h=abs(h);
h=20*log10(h);
subplot(211)
stem(b,'.');
grid;
title('切比雪夫逼近滤波器的抽样值')
subplot(212)
plot(w,h);
grid;
title('滤波器幅频特性(dB)')