海明窗课程设计.docVIP

课 程 设 计 任 务 书 专 业 通信工程 班 级 14通信4班 姓 名 张永余 设 计 起 止 日 期 设计题目： 语音信号的数字滤波——FIR数字带通滤波器的（三角）窗函数法设计 设计任务（主要技术参数）： 语音信号的采样录音（、N=60001）； 语音信号的频谱分析（） 语音信号的加噪（SNR1=8dB）； 数字滤波器的设计（，滤波器的阶数N=27）； 滤波结果输出及效果分析（SNR2=11dB） 指导教师评语： 成绩： 签字： 年 月 日 设计题目 语音信号的数字滤波——FIR数字带通滤波器的（三角）窗函数法设计 一、课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；掌握利用数字滤波器对信号进行滤波的方法。。 1.语音信号的采集 利用Windows系统采集：按“开始－＞程序－＞附件－＞娱乐－＞录音机”的顺序操作 打开Windows系统中的录音机软件 图1录音机 采集语音信号信息： （1）语音格式：*.wav 语音名称： （3）语音长度： （4）语音属性：kHz 窗函数设计法的原理 （1）窗函数的基本思想：先选取一个理想滤波器（它的单位抽样响应是非因果、无限长的），再截取（或加窗）它的单位抽样响应得到线性相位因果FIR滤波器。这种方法的重点是选择一个合适的窗函数和理想滤波器。 设x(n)是一个长序列，是长度为N的窗函数，用截断，得到N点序列，即 在频域上则有 由此可见，窗函数不仅仅会影响原信号在时域上的波形，而且也会影响 到频域内的形状。 （2）海明窗函数简介：海明窗（Hamming）和汉明窗是一个意思。简单的说汉明窗就是个函数，它的形状像窗，所以类似的函数都叫做窗函数。海明窗也是余弦窗的一种，又称改进的升余弦窗。海明窗与汉宁窗都是余弦窗，只是加权系数不同。海明窗加权的系数能使旁瓣达到更小。 分析表明，海明窗的第一旁瓣衰减为一42dB．海明窗的频谱也是由3个矩形时窗的频谱合成，但其旁瓣衰减速度为20dB/（10oct），这比汉宁窗衰减速度慢。海明窗与汉宁窗都是很有用的窗函数。公式如下： 将语音信号进行截取 编写MATLAB编码实现对信号的截取处理，将原始语音信号和截短语音信号在同一 图上画出来 图2原始信号频谱图 图3 截短信号频谱图 代码如下所示： [x,fs]=wavread(250.wav); x1=x(:,1); plot (x1, DisplayName, x1, YDataSource, x1); figure(gcf) x2=x(9500:67201); plot (x2, DisplayName, x2, YDataSource, x2); figure(gcf) 3.语音信号的加噪处理 (1)产生一个长度与你所截取的语音信号等长的噪声信号，为“250.wav”，所用指令为wavewrite(),所加噪声如图6所示。 (2)叠加噪声可以用awgn()语句，在叠加噪声时，应选取适当的值，尽量能体现出能量主要集中的频段。 (3)将语音与噪声合成的信号“250.wav”送入所设计的基于凯泽窗的FIR数字滤波器中进行滤波。 (4)将滤波结果导出保存成文件“250.wav”. 得出加噪后信号如下： 图4 加噪后的语音信号 Y = AWGN(x2,8,measured); 一个理想的带通滤波器应该有一个完全平坦的通带，在通带内没有放大或者衰减，有源带通滤波器电路并且在通带之外所有频率都被完全衰减掉，另外，通带外的转换在极小的频率范围完成。 4.原语音信号的频谱图 图5 语音信号频谱图1 图6 语音信号频谱图2[x,fs]=wavread(250.wav); x1=x(:,1); plot (x1, DisplayName, x1, YDataSource, x1); figure(gcf) x2=x(9500:67201); plot (x2, DisplayName, x2, YDataSource, x2); figure(gcf) X2=fft(x2); AX2=abs(X2); plot (AX2, DisplayName, AX2, YDataSource, AX2); figure(gcf) AX2=fftshift(AX2); N=length(AX2); plot([-N/2