如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计

网友回答

技术指标 Wp=0.2*pi,Ws=0.4*pi,Ap=0.25dB,As=50dB
方法一选择海明窗clear all;
Wp=0.2*pi;
Ws=0.4*pi;
tr_wide=Ws-Wp; %过渡带宽度
N=ceil(6.6*pi/tr_wide)+1; %滤波器长度
n=0:1:N-1;
Wc=(Wp+Ws)/2; %理想低通滤波器的截止频率
hd=ideal_lp1(Wc,N); %理想滤波器的单位冲击响应
w_ham=(hamming(N))'; %海明窗
h=hd.*w_ham; %实际海明窗的响应
[db,mag,pha,w]=freqz_m2(h,[1]); %计算实际滤波器的幅度响应
delta_w=2*pi/1000;
Ap=-(min(db(1:1:Wp/delta_w+1))) %实际通带纹波
As=-round(max(db(Ws/delta_w+1:1:501))) %实际阻带纹波
subplot(221)
stem(n,hd)
title('理想单位脉冲响应hd(n)')
subplot(222)
stem(n,w_ham)
title('海明窗')
subplot(223)
stem(n,h)
title('实际单位脉冲响应hd(n)')
subplot(224)
plot(wi/pi,db)
title('幅度响应(dB)')
axis([0,1,-100,10])
方法二Window=blackman(16);
b=fir1(15,0.3*pi ,'low',Window);
freqz(b,128)
注：b=fir1(N,Wn,window)只能设计低通和带通滤波器,并且滤波器的阶数必须为N+1
窗函数你可以换合适的上面只是举个例子