基于窗函数法的FIR数字低通滤波器设计
3.3 采用布莱克曼窗设计FIR数字低通滤波器
0-50-100-15000.511.522.5布莱克曼窗滤波器幅頻响应33.50-20-40-60-8000.10.20.30.40.50.60.7布莱克曼窗相频响应0.80.91
图3-3 采用布莱克曼窗设计FIR数字低通滤波器的仿真图
9
基于窗函数法的FIR数字低通滤波器设计
0.0150.010.005005101520冲激响应序列2530350.20.150.10.0500510冲激响应傅里叶变换15
图3-4采用布莱克曼窗冲激响应仿真图
1510x 104
50
0246810含噪语音信号121416x 104
3002001000
0246810滤波后语言信号121416x 104
图3-5采用布莱克曼窗对含噪语音信号及滤波后的语音信号仿真图
10
基于窗函数法的FIR数字低通滤波器设计
参考文献
[1] 从玉良.数字信号处理原理及其MATLAB实现[M].北京：电子工业出版社.2009.7
[2] 胡广书.数字信号处理理论、算法与实现[M].北京：清华大学出版社.2003,8 [3] 万永革.数字信号处理的MATLAB实现[M].北京: 科学出版社.2007. [4] 薛山.MATLAB基础教程[M].北京: 清华大学出版社.2011.
[5] 陈怀琛.数字信号处理教程—MATLAB释义与实现[M].北京: 电子工业出版社.2002
[6] 缪家鼎，徐文娟． 光电技术[M]． 杭州：浙江大学出版社，1994：112-156．
[7] 李晓东，张庆红，叶瑾琳．气候学研究的若干理论问题[J]．北京大学学报：
自然科学版，1999，35(1)：101-106．
[8] 郑开青．通讯系统模拟及软件[D]．北京：清华大学无线电系，1987． [9]
Online
Computer
Library
Center,
Inc.
History
of
OCLC[EB/OL]．[2000-01-08]． http: //www. oclc. org/about/ history/default. htm.
11
基于窗函数法的FIR数字低通滤波器设计
附 录
矩形窗设计程序 passrad=0.4*pi; w2=boxcar(16);
n=1:1:16;
hd=sin(passrad*(n-8))./(pi*(n-8)); hd(8)=passrad/pi; h2=hd.*rot90(w2);
title('designed by Hanning window'); [mag2,rad]=freqz(h2); subplot(2,1,1);
plot(rad,20*log10(abs(mag2))); grid on;
[h2,w2]=freqz(h2,1,100,2); subplot(2,1,2);
plot(w2,unwrap(angle(h2))); grid on;
汉明窗设计程序 passrad=0.4*pi; w2=hamming(32); n=1:1:32;
hd=sin(passrad*(n-16))./(pi*(n-16)); hd(16)=passrad/pi; h2=hd.*rot90(w2);
title('designed by Hanning window'); [mag2,rad]=freqz(h2); subplot(2,1,1);
12
基于窗函数法的FIR数字低通滤波器设计
plot(rad,20*log10(abs(mag2))); grid on;
[h2,w2]=freqz(h2,1,100,2); subplot(2,1,2);
plot(w2,unwrap(angle(h2))); grid on;
布莱克曼窗函数设计程序 passrad=0.4*pi; w2=blackman(64); n=1:1:64;
hd=sin(passrad*(n-32))./(pi*(n-32)); hd(32)=passrad/pi; h2=hd.*rot90(w2);
title('designed by Hanning window'); [mag2,rad]=freqz(h2); subplot(2,1,1);
plot(rad,20*log10(abs(mag2))); grid on;
[h2,w2]=freqz(h2,1,100,2); subplot(2,1,2);
plot(w2,unwrap(angle(h2))); grid on;
13