fir_lp.m

this a DSP program for FIR

%根据技术指标，设计FIR数字低通滤波器
%选哈明窗
clc
wp=0.2*pi;ws=0.3*pi;
deltaw=ws-wp;       %Δω为过渡带的宽度
N1=ceil(6.6*pi/deltaw);  %计算N;;;;;;;;;;;;;;;;;;;;;;;
N=N1+mod(N1+1,2);
%确保N的长度为奇数,实现第一类线性相位滤波

wdham=(hamming(N))';    %求哈明窗
wc=(wp+ws)/2;
%截止频率取两边频率的平均值
hd=ideallp(wc,N);
h=hd.*wdham;
%检验
% % figure(22)
% % freqz(h)
H=fft(h,128);

HdB=20*log10(abs(H+eps)/max(abs(H)));%计算相对振幅
Ap=-min(HdB(1:wp/(2*pi/128)+1))
   %检验通带波动
plot(HdB);
legend('哈明窗,FIR数字低通滤波幅频特性')