语音信号处理.docVIP

PAGE PAGE 3 学号： 2010 - 2011学年 第 2 学期 《语音信号处理》 课 程 设 计 报 告 题 目： 语音信号的分析与处理 专 业： 电子信息工程 班 级： 姓 名： 指导教师： 冯 锁 成 绩： 电气工程系 2011年6月7日 目 录 TOC \o 1-3 \h \z \u 一.课程设计目的 3 二.设计内容 3 三．设计意义 3 四.设计原理 3 4.1采样定理 3 4.2 DFT和FFT原理 4 4.3 MATLAB中计算DFT(FFT)的函数 4 4.4 窗函数法 4 4.5 窗函数法设计线性相位FIR滤波器步骤 5 五.详细设计过程 6 5.1 绘出采样后语音信号的时域波形和频谱 6 5.2对语音信号加噪声，绘出加噪声后语音信号的时域波形和频谱图 6 5.3通过窗函数法设计FIR滤波器 7 六 设计总结 9 七 参考文献 10 一.课程设计目的 综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。 二.设计内容 基本要求：录制一段自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标，采用窗函数法设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化。 三．设计意义 语言是人类持有的功能．声音是人类常用的工具，是相互传递信息的最主要的手段。语音信号是人们构成思想疏通和感情交流的最主要的途径，通过语音传递倍息是人类最重要、最有效、最常用和最方便的交换信息的形式。现在随着时代的发展，人类已开始进入了信息化时代，用现代手段研究语音信号，使人们能更加有效地产生、传输、存储、获取和应用语音信息，这对于促进社会的发展具有十分重要的意义。现在社会衍生了很多现代的语音通信方式，手机语音、语音聊天软件以及语音小说等，涉及语音的方面几乎遍及了我们的大部分生活，和我们生活息息相关，因此语音信号处理是具有现实意义的研究课程。 四.设计原理 4.1采样定理 在进行模拟/采样信号的转换过程中，当采样频率大于模拟最高频率的2倍时，则采样后的采样信号完整地保留了原始信号中的信息，采样定理又称奈奎斯特定理。 4.2 DFT和FFT原理 长度为N的序列x(n)的离散傅立叶变换为X(k)： 首先按n的奇偶把时间序列x(n)分解为两个长为N/2点的序列 r=0,1,...,N/2-1 r=0,1,...,N/2-1 其中X1(k) 和X2(k)分别为x1(n) 和x2(n)的N/2点DFT。因为X1(k) 和X2(k)均是以N/2为周期的，且。因此可将N点DFT X(k)分解为下面的形式 k=0,1,...,N/2-1 k=0,1,...,N/2-1 通过上面的推导可以看出，N点的DFT可以分解为两个N/2点的DFT，每个N/2点的DFT又可以分解为两个N/4点的DFT。依此类推，当N为2的整数次幂时()，由于每分解一次降低一阶幂次，所以通过M次的分解，最后全部成为一系列2点DFT运算。以上就是按时间抽取的快速傅立叶变换(FFT)算法。 4.3 MATLAB中计算DFT(FFT)的函数 函数fft用来求序列的DFT,调用格式为：[Xk]=fft(x,N) 其中x 为有限长序列，N为序列x的长度，Xk为序列xn的DFT. 函数ifft用来求IDFT,调用格式为：[x]=ifft(Xk,N) 其中，Xk为有限长序列，N为序列Xk的长度，x为序列Xk的IDFT 4.4 窗函数法 设欲设计的滤波器的理想频率响应为，单位脉冲响应为，与是一对傅式变换，因此有 根据给定的求得的一般是无限长的且是非因果的。为了得到一个因果的有限长的滤波器，最直接的方法是截断，或者说用一个窗口函数对进行加窗处理 成为实际设计FIR滤波器的单位脉冲响应,其频率响应为为： 其中N为窗口的长度。窗口函数的形状和窗口长度N决定了窗函数法设计出的FIR滤波器的性能。 4.5 窗函数法设计线性相位FIR滤波器步骤 数字滤波器的性能要求：临界频率{ωk}、滤波器单位脉冲响应长度N； 2）根据性能要求，合理选择