数字语音滤波系统设计..docxVIP

信息与通信工程学院 数字信号处理实验报告系：通信工程系专 业：通信工程班 级：通信133实验题目：语音信号分析与处理学生姓名：xxx指导教师：xx 学 号xxxxxxxx 完成日期：2015年 12月 18 日设计目的 综合运用数字信号处理和MATLAB设计语音数字滤波系统。熟悉MATLAB的使用，掌握MATLAB的程序设计和GUI用户界面设计的方法。学习并掌握在MATLAB环境下语音信号采集的方法。理解并运用数字信号处理的基本概念、基本理论和基本方法解决问题。掌握MATLAB设计数字滤波器的方法。学会用MATLAB对信号进行分析和处理。设计原理1.采样定理：在进行模拟/数字信号的转换过程中，当采样频率大于等于信号最高截止频率的2倍时，即：，则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高截止频率的5～10倍；采样定理又称奈奎斯特定理。l2.时域信号的FFT分析:FFT即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。?在MATLAB的信号处理工具箱中函数FFT和IFFT用于快速傅立叶变换和逆变换。函数FFT用于序列快速傅立叶变换，其调用格式为y=fft(x)，其中，x是序列，y是序列的FFT，x可以为一向量或矩阵，若x为一向量，y是x的FFT且和x相同长度；若x为一矩阵，则y是对矩阵的每一列向量进行FFT。如果x长度是2的幂次方，函数fft执行高速基－2FFT算法，否则fft执行一种混合基的离散傅立叶变换算法，计算速度较慢。函数FFT的另一种调用格式为y=fft(x,N)，式中，x，y意义同前，N为正整数。函数执行N点的FFT，若x为向量且长度小于N，则函数将x补零至长度N；若向量x的长度大于N，则函数截短x使之长度为N；若x?为矩阵，按相同方法对x进行处理3.数字滤波器设计原理和方法3.1窗函数法如果所希望的滤波器的理想频率响应函数为 ，则其对应的单位脉冲响应为 = （2-1）窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近。由于往往是无限长序列，且是非因果的，所以用窗函数将截断，并进行加权处理，得到： ＝ (2-2)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为：＝ (2-3)式中，N为所选窗函数的长度。 由第七章可知，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。 这样选定窗函数类型和长度N后，求出单位脉冲响应＝·，并按式（2-3）求出。是否满足要求，要进行验算。一般在尾部加零使长度满足于2的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果不满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。如果要求线性相位特性，则还必须满足 (2-4) 根据上式中的正负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据设计的滤波特性正确选择其中一类。例如，要设计线性低通特征，可选择一类，而不能选一类。3.2脉冲不变响应法如果从模拟到数字滤波器我们想要保留脉冲响应的形状，那么就得到一种方法称为脉冲不变响应法的变换方法。?脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序h(n)模仿模拟滤波去的冲击响应ha(t)，使h(n)正好等于ha(t)的采样值，即?h(n)=ha(nT)T为采样周期。如以Ha(s)及H(z)分别表示ha(t)的拉式变换及h(n)的z变换，即?Ha(s)=L[ha(t)]H(z)=Z[h(n)]则根据采样序列z变换与模拟信号拉式变换的关系，得：?上式表明，采样脉冲响应不变法将模拟滤波器变换为数字滤波器时，它所完成的S平面到Z平面的变换，正是以前讨论的拉式变换到Z变换的标准变换，即首先对Ha(s)作周期严拓，然后再经过的映射关系映射到Z平面上。3.3双线性变换法脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Z平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-π/T～π/T一条横带里；第二步再通过标准变换关系T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。 双线性变换法的S域与Z域间的变换关系为： （1）由式(1)可以看出,z和s之间可以直接代换,由于这是一非线性变换,需考察一下是否能把jΩ映射成单