Question

如何用MATLAB设计fir低通滤波器

设计具有下面指标的低能FIR滤波器
ωp=0.2rad ,Ap=0.25dB,ωs=0.3rad ,As=50dB
分别选择矩形窗,汉宁窗,汉明窗来实现这个滤波器,画出其幅频特性曲线和相频特性曲线,找出不同
求高手解答,小弟感激不尽!

Answer 1

2-4用海明窗设计一个FIR滤波器，其中Wp=0.2*pi，Ws=0.3*pi,通带衰减不大于0.25dB,阻带衰减不小于50dB。
代码：
wp=0.2*pi;
ws=0.3*pi;
tr_width=ws-wp;
M=ceil(6.6*pi/tr_width)+1;
n=[0:M-1];
wc=(ws+wp)/2;%ideal LPF cutoff frequency
hd = ideal_lp(wc,M);
w_ham=(hamming(M))';
h=hd.*w_ham;% .*
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
Rp=-(min(db(1:1:wp/delta_w+1:1:501))) %Min stopband attenuation
%Plots
subplot(2,2,1);
stem(n,hd);
title('Ideal Impulse Response')
axis([0 M-1 -0.1 0.3]);
xlabel('n');ylabel('hd(n)');
%%%%%%
subplot(2,2,2);
stem(n,w_ham);
title('Hamming Window')
axis([0 M-1 0 1.1]);
xlabel('n');ylabel('w(n)');
%%%%%%
subplot(2,2,3);
stem(n,h);
title('Actual Impulse Response')
axis([0 M-1 -0.1 0.3]);
xlabel('n');ylabel('h(n)');
%%%%%%
subplot(2,2,4);
plot(w/pi,db)
title('Magnitude Response In DB');grid
axis([0 1 -100 10]);
xlabel('Frequency In pi Units');ylabel('Decibels');

Answer 2

这个直接用MATLAB中的FDAtool就可以的。