滤波器设计（处理加噪音频，窗函数法）基于matlab的FIR滤波器设计（处理加噪音频，窗函数法）
一、录制音频
二、利用窗函数法设计的低通滤波器处理加有噪声的音频波形f37000所加噪声正弦函数的频率YfsbitswavreadDyi
pi
ihao01wav利用wavread产生音频的函数及采样频率Lle
gthYt01fsL1fs定义时间的范围及步长y0005si
2pif3t
1floorL2所加噪声f10
1fsLYY1sou
dYfs输出加噪前音频Y1yY给音频加噪声FY1absfftY1LFYabsfftYLsou
dY1fs输出加噪后的音频figure1subplot211plott11000Y11000grido
加噪前音频的时域图xlabel时间tylabel幅度Ytitle加噪前录音波形的时域图subplot212plotf1FY1
11grido
加噪前音频的频域图xlabel频率fylabel幅度FY
ftitle加噪前录音波形的频域图figure2subplot211plott11000Y111000grido
加噪后音频的时域图xlabel时间tylabel幅度Y1title加噪声后录音波形的时域图subplot212plotf1FY11
11grido
加噪后音频的频域图xlabel频率fylabel幅度FY1title加噪声后录音波形的频域图m003Mrou
d8mNM1定义滤波器的阶数bfir1N06figure3hffreqzb1512滤波器的幅频特性图plotffs2pi20log10absh参数分别是频率与幅值xlabel频率赫兹ylabel增益分贝title滤波器的增益响应figure4sffilterb1Y1使用filter函数对信号进行滤波FsfabsfftsfLsubplot211plott11000sf11000grido
滤波后音频的时域图xlabel时间tylabel幅度sftitle滤波后录音波形的时域图axis00100500020002subplot212plotf1Fsf1
11grido
滤波后音频的频域图xlabel频率fylabel幅度Fsftitle滤波后录音波形的频域图sou
dsffs
f三、运行结果
加噪前录音波形的时域图和频域图
加噪后录音波形的时域图和频域图
f窗函数法设计低通滤波器的增益响应
滤波后录音波形的时域图和频域图
fr