m0721.m

来自「MATLAB仿真技术与实例应用教程」· M 代码 · 共 16 行

%使用Kaiser窗设计低通滤波器
Wp=0.2*pi;
Ws=0.3*pi;
Rp=0.25;
As=50;
w=[Wp/pi Ws/pi];
tr_width=Ws-Wp;
M=ceil(6.6*pi/tr_width)+1 %计算滤波器的长度
beta=0.1102*(As-8.7);
kai_w=(kaiser(M,beta));
b=fir1(M-1,w,'stop',kai_w);
[h,w]=freqz(b,1,512,2);
plot(w,20*log10(abs(h)));
grid on
title('低通滤波器幅频特性')