《语音信号滤波去噪使用脉冲响应不变法设计的巴特沃斯滤波器》由会员分享，可在线阅读，更多相关《语音信号滤波去噪使用脉冲响应不变法设计的巴特沃斯滤波器（26页珍藏版）》请在金锄头文库上搜索。
1、王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页语音信号滤波去噪使用脉冲响应不变法设计的巴特沃斯滤波器学生姓名：王雄杰指导老师：喻小虎摘 要 本课程设计主要内容是设计利用脉冲响应不变法设计的巴特沃斯滤波器，对一段含噪语音信号进行滤波去噪处理并根据滤波前后的波形和频谱分析滤波性能。本课程设计仿真平台为MATLAB7.0，开发工具是M语言编程。录制一段语音信号，并人为加 入一单频噪声，然后对信号进行频谱分析以确定所加噪声频率，即fn=3000Hz,并设计滤波器进行滤波去噪处理，最后比较滤波前后的波形和频谱，得出结果为，滤波器后的 语音信号与原始信号基本一致，即设计的巴
2、特沃斯滤波器能够去除信号中所加单频噪 声，达到了设计目的。关键词 课程设计；滤波去噪；巴特沃斯滤波器；脉冲响应不变法；MATLAB本课程设计需要录制一段频率为 8000Hz，8位的单声道语音信号，并绘制波形观察 其频谱，再对其进行加噪处理，同样要绘制加噪后的频谱图，再用 MATLAB利用脉冲 响应不变法设计一个巴特沃斯滤波器，将该语音信号进行滤波去噪处理。1.1课程设计目的课程设计是教学的最后一个步骤，课程设计有利于基础知识的理解，我们掌握了基 础知识和基本技能，但是要真正接触才能真正理解课程的深入部分；还有利于逻辑思维 的锻炼，在许多常规学科的日常教学中，我们不难发现这样一个现象，不少学生的
3、思维 常常处于混乱的状态，写起作文来前言不搭后语，解起数学题来步骤混乱，这些都是缺 乏思维训练的结果，所以我们可以通过实践来分析问题、解决问题、预测目标等目的； 同时也有利于与其他学科的整合，例如我们这次的课程设计就要运用MATLAB软件的 帮助才能实现；最重要的有利于治学态度的培养，在课程设计中，我们可能经常犯很多 小错误，可能要通过好几次的反复修改、调试才能成功，但这种现象会随着学校的深入 而慢慢改观。这当中就有一个严谨治学、一丝不苟的科学精神的培养，又有一个不怕失 败、百折不饶品格的锻炼。数字信号处理课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的
4、目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设 课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等 步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力， 提高实际应用水平。1.2课程设计的要求（1）滤波器指标必须符合工程设计。（2）设计完后应检查其频率响应曲线是否满足指标。（3）处理结果和分析结论应该一致，而且应符合理论。（4）独立完成课程设计并按要求编写课程设计报告。1.3设计平台MATLAB是美国MathWorks公司生产的一个为科学和工程计算专门设计的交互式 大
5、型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。 它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医学、金融和其他 需要进行复杂数值计算的领域得到了广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的数学 工具，在世界各地的高和大型计算机上运行， 适用于Windows、UNIX等多种系统平台。MATLAB作为一种科学计算的高级语言之所以受欢迎，就是因为它有丰富的函数资 源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再去编写大量繁琐的 程序代码，从而减轻了编程人员的工作负担， 被称为
6、第四代编程语言的MATLAB最大的 特点就是简洁开放的程序代码和直观实用的开发环境 。王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页2设计原理本课程设计需要录制的一段频率为 8000Hz, 8位的单声道语音信号，并绘制波形观 察其频谱，再对其进行加噪处理，同样要绘制加噪后的频谱图，再用 MATLAB禾U用脉 冲响应不变法设计一个巴特沃斯滤波器，将该语音信号进行滤波去噪处理。2.1 IIR滤波器IIR滤波器具有无限长脉冲响应，因此能够与模拟滤波器相匹敌；一般来说，所有 的模拟滤波器都有无限长脉冲响应。因此，IIR滤波器设计的基本方法是利用复值映射 将大家熟知的模拟
7、滤波器变换为数字滤波器。这一方法的优势在于各种模拟滤波器设计 (AFD )表格和映射在文献中普遍都能获得。这种 IIR滤波器设计的基本方法存在两种 途径：途径1:设计模拟低通滤波器实行频带变换S T S实行滤波器 变换STZ匚=期望的IIR滤波器途径2:设计模拟低| 通滤波器=实行滤波器变 换STZ实行频带变 换ZtZ匚二期望的IIR滤波器我们将研究途径二的设计基本方法，在这种IIR滤波器设计方法将按下列步骤进行:(1) 设计模拟低通滤波器。(2) 研究并实行滤波器变换以得到数字低通滤波器。(3) 研究并实行频带变换以便从数字低通滤波器得到其他数字滤波器o22巴特沃斯滤波器巴特沃斯滤波器是用上
8、述IIR滤波器的性质表征的：它的幅度响应在通带和阻带都是平坦的。一个N阶低通滤波器的幅度平方响应给出为：卫囚汕(円)卜1调(2-1)式中，N为整数，是滤波器的阶，ic(rad/s)是截止频率。称为滤波器的阶数，N越大,王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页通带和阻带的近似性越好，过渡带也越陡。幅度平方响应的图如下所示。、N=N=1002 IN=1RkN=200JLJ图2.1幅度平方响应图从这张图可看出下面几个性质:(1) 0 =0, Ha(j0)2*=1，对全部 N。2 1(2) c, Ha(Zc)，对全部N，这意味着在门c有3dB衰减。2(3) |Ha
9、(j0)|是0的单调下降的函数。(4) |Ha(j0)|随Nt珀向一个理想低通滤波器趋近。(5) |Ha(j0)|2在0 =0世最大平坦，因为在这里所有阶的导数存在且等于零。为了确定系统函数Ha(s)现将得到Ha(S)Ha(-S)二一1S丿(rc)2N(rc)2N2N2N 2Ns + jc)卄)k=1由(2-2)式分母多项式的根(或Ha(s)Ha(-s)的极点)给出为(2-2)Pk =(T)叫厂 c)*页(2宀 1)J誌 k = 0,1,M2N_1(2-3)(2-3)式的解释是:王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页(1) 2N个极点均匀分布在半径为 Q
10、c的圆上，相隔pi/N(2) N 为奇数时，pk= Qcej2pik/N(3) N 为偶数时，pk= Qcej(pi/2N+kpi/N)(4) 极点对虚轴(j Q对称分布(5) 极点永远不会落在虚轴上，且仅当N为奇数时才会落在实轴上(6) 选取极点全部落在左半平面可得到因果稳定的Ha(s)作为一个例子，三阶和四阶巴特沃斯滤波器的极点分布如图所示：JOijQ(a)(A)图2.2巴特沃斯滤波器的极点图通过选取在左半面的极点就能给出一个稳定和因果的Ha(s)，并且能将Ha(s)写成:Ha(S)(2-4)2.3脉冲响应不变法 (s-久)kT脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲
11、响应序h( n)模仿模拟滤波去的冲击响应ha(t)，使h(n)正好等于ha(t)的采样值，即h( n)=ha( nT)T为采样周期。如以Ha(s)及H(z)分别表示ha(t)的拉式变换及h(n)的z变换 NR变换前：H a ( S )kk =1 s 一 p k(2-5)(2-6)王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页变换后:NZk =1z -1(2-7)由于z=ejw是在单位圆上，s=j jx是在虚轴上，所以有下面从s平面到z平面的变换z=esT王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页(2-8)系统函数H (z)和H
12、a( s)是经过频域混叠公式联系的:H(z)| z严1 f-Z Ha s+T m_oc (2-9)王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页在(2-8)式的映射关系下，复平面的变换如图2.3王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页图2.3脉冲响应不变法中的复平面映射映射特征：(1) 记c = Re(s): (T 映射到|z|0央射到|z|1 (单位圆外)。(2) 多个
13、s映射到一个z :多点对一点的映射，每个宽度为2n /T的全部半无限 带都映射到单位圆内。(3) 因果稳定的AF仍然映射为因果稳定的DF。(4) 如果滤波器不是真正带限会产生混叠误差 (Aliasi ng )3.设计步骤3.1设计流程图设计流程图如语音信号滤波去噪一一使用脉冲响应不变法设计的巴特沃斯滤波器,下3.1所示：图3.1脉冲响应不变法巴特沃斯滤波器对语音信号去噪流程图3.2语言信号的采集人工录制一段为.wav的语音信号，时间在8-9s左右，要求为8000Hz, 8位单声道的王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第 8页共20页 音频格式。然后在Matlab软件
14、平台下，利用函数wavread对语音信号进行采样，函数为x,fs,bits=wavread(music.wav),记住采样频率 fs 和采样点数。所下载音乐的属性设置如图3.2：图3.2语音信号参数设置3.3语音信号的频谱分析首先使用MATLAB分别画出原始语音信号的时域波形，再对原始语音信号进行快 速傅里叶变换，画出原始信号的频率 一一幅度谱，然后再原始语音信号的基础上加入单 频噪声后，噪声频率为fn=3000Hz,再所得到的音乐做时间一一幅度谱，同样对它做傅里 叶变换，得到此信号的频率 一一频谱特性。可得到图形3.3和图3.4。其中傅里叶变换表示为：X=abs(fft(x); Y=abs(fft(y)王雄杰语音信号滤波去噪 -使用脉冲响应不变法设计的巴特沃斯滤波器第#页共20页9 Edit or - d; IATLAB7vorkUnt itlcd Eile Edit Iect Cell Tgls Debug Qeskt gp Window Help勺
