简易讲解：使用MATLAB设计低通滤波器（LPF）教程，附带Matlab源代码及第323期刊登

%设计FIR数字低通滤波器，截止频率为pi/4，在不同窗口长度（N=15，N=33）下，分别求出
%h(n),通过幅频特性和相频特性，观察3dB带宽和20dB带宽，总结窗口长度N对滤波特性的影响。

close all;
clc;
clear;
N=15;
N1=N+mod(N+1,2);
wcc=pi/4;
wc=wcc/pi;
hn=fir1(N1-1,wc,hanning(N1)); %矩形窗boxcar(N) 三角窗bartlett(N) 汉宁窗 hanning(N) 哈明窗hamming(N)
% 布莱克曼窗函数blackman(N) 凯瑟-贝赛尔窗函数 kaiser(N,beta)
figure(1)
subplot(3,1,1)
freqz(hn,1)
ylim([-150,50]);
hold on;
plot([-1 1],[-3 -3],'r')
plot([-1 1],[-20 -20],'y')
title('N=15')
legend('N=15','3dB','20dB')

clc;
clear;
N=33;
N1=N+mod(N+1,2);
wcc=pi/4;
wc=wcc/pi;
hn=fir1(N1-1,wc,hanning(N1)); %矩形窗boxcar(N) 三角窗bartlett(N) 汉宁窗 hanning(N) 哈明窗hamming(N)
% 布莱克曼窗函数blackman(N) 凯瑟-贝赛尔窗函数 kaiser(N,beta)
%以N=33，wc=pi/4为技术指标，用三种窗函数设计线性相位低通滤波器，通过对相应的幅频特性分析，
%观察3dB和20dB带宽以及阻带最小衰减，比较三种窗函数对滤波特性的影响。
%假设：通带边界频率wp为pi/5,阻带截止频率ws为pi*3/10,此时通带截止频率wc为pi/4.
%结论：观察波形和数值，发现窗函数会影响滤波器带内波动和阻带最小衰减，
% 其中三角窗的幅频特性为单调衰减，相频特性在整个频域均为线性。矩形窗的阻带内最小衰减最大，衰减最快。

close all;
clc;
clear;
N=33;
N1=N+mod(N+1,2)
wcc=pi/4;
wc=wcc/pi;
hn=fir1(N1-1,wc,hanning(N1)); %矩形窗boxcar(N) 三角窗bartlett(N) 汉宁窗 hanning(N) 哈明窗hamming(N)
%布莱克曼窗函数blackman(N) 凯瑟-贝赛尔窗函数 kaiser(N,beta)
figure(1)
freqz(hn,1); %汉宁窗3dB带宽为0.2246*pi,20dB带宽为0.3074*pi,阻带内最小衰减为17.73
hold on;
plot([-1 1],[-3 -3],'r')
plot([-1 1],[-20 -20],'y')
title('N=33 汉宁窗')

hn1=fir1(N1-1,wc,boxcar(N1)); %矩形窗3dB带宽为0.2385*pi,20dB带宽为0.2776*pi,阻带内最小衰减为21.41
figure(2)
freqz(hn1,1);
hold on;
plot([-1 1],[-3 -3],'r')
plot([-1 1],[-20 -20],'y')
title('N=33 矩形窗')

hn2=fir1(N1-1,wc,bartlett(N1)); %三角窗3dB带宽为0.2261*pi,20dB带宽为0.315*pi,阻带内最小衰减为16.09
figure(3)
freqz(hn2,1);
hold

上一篇： 理解dB和dBm的基本概念及其计算方法

下一篇： 详解:dBi、dBd、dB、dBm和dBc的区别与含义