数字信号处理课程设计-语音信号的数字滤波——fir数字滤波器的汉宁窗函数设计(编辑修改稿)

数字信号处理课程设计-语音信号的数字滤波——fir数字滤波器的汉宁窗函数设计(编辑修改稿)内容摘要：

下图 5 所示为 hanning 窗 及系统的 hd和 hn 函数特性图。 图 5 汉宁窗的函数特性图 ⑸ 滤波器进行滤波 如下图 6 所示，为滤波器滤波后的效果图。 第一个是滤波后的原始语音信号，第二个和第三个是滤波后的频谱图。 由于参数中已经设定 3dB 截止频率为 1100Hz,所以频率在 1100Hz的信号均被滤除 ，达到滤波效果。 课程设计说明书 沈 阳 大 学 图 6 滤波后的波形图 ⑹ 设计过程中完整的代码 clear close all %设参数。 ***************************** fc=1100 N=42 [s,fs]=wavread(39。 39。 )。 %*************************************** jieduan=s(1000:1000+fs1)。 wavwrite(jieduan,fs,39。 39。 ) jiazao=awgn(jieduan,20)。 叠加噪声 wavwrite(jiazao,fs,39。 39。 ) JIEDUAN=fft(jieduan) JIAZAO=fft(jiazao) 课程设计说明书 沈 阳 大 学 sound(jieduan) sound(jiazao) figure(1) plot(s) figure(2) subplot(311) plot(jieduan) title(39。 截短后时域波形 39。 ) subplot(312) plot(abs(JIEDUAN)) title(39。 截短波形 FFT 幅频特性 39。 )。 subplot(313)。 %在图形窗口第三个分区中画图 plot(abs(JIEDUAN))。 %绘制被处理的语音信号的实际幅频特性（正频率段） axis([0,4000,0,700]) title(39。 实际幅频特性 39。 )。 figure(3) subplot(311)。 %把一个图形窗口中分三行一列，在第一个分区中画图 plot(jiazao)。 %绘制被处理的语音信号的时间域波形 title(39。 滤波前信号波形 39。 )。 subplot(312)。 %在图形窗口第二个分区中画图 plot(abs(JIAZAO))。 %绘制被处理的语音信号的 FFT 幅频特性 grid。 title(39。 FFT 幅频特性 39。 )。 subplot(313)。 %在图形窗口第三个分区中画图 plot(abs(JIAZAO))。 %绘制被处理的语音信号的实际幅频特性（正频率段） grid。 axis([0,4000,0,700]) title(39。 实际幅频特性 39。 )。 figure(4) 课程设计说明书 沈 阳 大 学 N=42。 wn=hanning(N)。 subplot(311) plot(wn) grid。 title(39。 汉宁窗函数特性 39。