基于Matlab的FIR带通滤波器设计及DSP实现
本文采用第二种方法。选择菜单File→Export打开Export对话框，选择Coefficient File(ASCII)，指定输出文件名，生成*fcf文件，就可以直接把生成的*fcf文件中的系数添加到DSP程序的变量中了。
在Matlab软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数。采样值放在向量y中，fS表示采样频率(单位：Hz)，bits表示采样位数。根据语音信号的特点给出数字带通FIR滤波器的性能指标：设置通带频率fp1=1 200 Hz，fp2=3 000 Hz，阻带频率fs1=1 000 Hz，fs2=3 200 Hz，抽样频率fs=8 000 Hz，阻带最小衰减As=50 dB，通带最大衰减Ap=1 dB。这里采用窗函数设计法设计FIR带通滤波器。首先根据阻带衰减选择窗形状，海明窗和布莱克曼窗等窗函数均可提供大于50 dB的衰减。由于海明窗可提供较小的过渡带，所以选择海明窗。根据ωp1=2fp1／fs；ωs1=2fs1／fs，B(带宽)=ωp1-ωs1，N(窗口长度)=ceil(6．6／B)，n(阶数)=N-1；算出阶数为132。在FDATo ol设计界面，按要求填写以上参数，得到FIR数字带通滤波器的仿真图像，如图1所示。本文引用地址：//m.amcfsurvey.com/article/186064.htm
然后再用设计的FIR数字带通滤波器对加噪的语音信号进行滤波。在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波。编写Matlab程序如下，实现上述功能。
程序中“quzao”为录制的波形音频文件(．wav)，对加噪的语音信号进行滤波仿真处理后，输出结果如图2所示。
评论