pshdsi
wmpim得到理想低通滤波器ki
put所有取值boxcar1hammi
g2ha
i
g3blackma
4
kifk1BboxcarN矩形窗stri
gBoxcarN
um2strNelseifk2Bhammi
gNstri
gHammi
gN
um2strNelseifk3Bha
i
gNstri
gHa
i
gN
um2strNelseifk4Bblackma
Nstri
gBlackma
N
um2strNe
dhhdB得到FIR数字滤波器Hmfreqzh11024whole求其频率响应magabsH得到幅值db20log10magepsmaxmagdb计算公式phaa
gleH得到相位figure1subplot221
0N1
3
f用窗函数发设计4FIR数字滤波器
stem
haxis0N10103holdo
0N1xzerosNplot
xxlabel
ylabelh
title实际低通滤波器的h
text03N027stri
gholdoffsubplot222plotmpidbaxis011000xlabelwpiylabeldBtitle衰减特性dBgridsubplot223plotmphaholdo
07xzeros8plot
xtitle相频特性xlabel频率radylabel相位radaxis031544subplot224subplot224plotmmagtitle频率特性xlabel频率Wradylabel幅值axis0315005text0912stri
gN15时
4
f用窗函数发设计5FIR数字滤波器
实际低通滤波器的h
03
Ha
i
gN15
02
衰减特性dB
0
Ha
i
gN15
dB
h
相位rad
01
50
0
010
4
5
10


相频特性
2
0
2
4
0
1
2
3
频率rad
N33时
实际低通滤波器的h
03
Ha
i
gN33
02
幅值
1000
04030201
00
0
05
1
wpi
频率特性
1
2
3
频率Wrad
衰减特性dBHa
i
gN33
dB
h
相位rad
01
50
0
01
100
0
10
20
30
0
05
1


wpi
相频特性
频率特性
4
2
04
幅值
03
002
2
01
4
0
1
2
3
频率rad
0
0
1
2
3
频率Wrad
分析：通过对上述两幅图的分析，我们可以看出，当N增大时过渡带会相应地减少，也就是过渡带宽会变窄。N越大，阻带衰减越严重，越符合理想低通滤波器。随着N的增大，阻带波动越来越快。
5
f用窗函数发设计6FIR数字滤波器
b、N33，WcΠrad，用四种窗函数设计线性相位低通滤波器。绘制相4
应的幅频特性曲线，观察3dB带宽和20dB带宽以及阻带最小衰减，比较四种窗函数对滤波特性的影响。矩形窗
h
相位rad
实际低通滤波器的h
03
BoxcarN33
02
01
0
01
0
10
20
30


相频特性
4
2
0
2
4
0
1
2
3
频率rad
汉宁窗
实际低通滤波器的h
03
Ha
i
gN33
02
h
相位rad
01
0
010
4
10
20
30


相频特性
2
0
2
4
0
1
2
3
频率rad
幅值
dB
幅值
dB
0
50
1000
04030201
00
0
50
1000
04030201
00
衰减特性dBBoxcarN33
05
1
wpi
频率特性
1
2
3
r