基于语音信号去噪处理的FIR低通滤波器设计.
语音信号滤波去噪基于布莱克曼窗的FIR的滤波器

语音信号滤波去噪——基于布莱克曼的FIR滤波器学生姓名:指导老师:摘要本课程设计主要内容是设计利用窗口设计法选择布莱克曼窗设计一个FIR滤波器,希望对一段含噪语音信号进行滤波去噪处理并根据滤波前后的波形和频谱分析滤波性能。
本课程设计仿真平台为MATLAB7.0,开发工具是M语言编程。
用麦克风采集一段语音信号,绘制波形并观察其频谱特点,加入一个带外单频噪声,用布拉克曼窗设计一个满足指标的FIR滤波器,对该含噪语音信号进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析,根据结果和学过的理论得出合理的结论。
最终在通过程序的调试得出最后的分析图,且基本符合设计要求。
关键词语音信号去噪;FIR滤波器;布莱克曼窗;MATLAB1 引言本课程设计主要解决在含噪情况下对语音信号的滤波去噪处理,处理时采用的是利用窗口设计法选择布莱克曼窗设计的FIR滤波器[1]。
在我们日常生活中,不可避免的要通过电话的形式与他人沟通、交流,在语音通信的过程中将不可避免的会受到来自周围环境的干扰,例如在信息传输中,传输媒介引入的噪声,通信设备内部的热噪声,或者是其他人在进行语音信号传输交流的信号与其混叠等。
在这些干扰噪声的存在下,接收者接收到的语音已经不是原始的纯净语音信号,而是受噪声干扰的带噪声语音信号。
而本课程设计就是利用MATLAB集成环境用布莱克曼窗的方法设计一个FIR滤波器,对语音信号进行滤波去噪处理,并将虑噪前后的频谱图进行对比。
1.1 课程设计目的数字信号处理是一门理论性和实践性都很强的学科,通过课程设计可以加深理解掌握数字信号处理的基本理论,培养我们的分析问题以及解决问题的能力,同时也为将来我们走向工作岗位奠定坚实的基础。
在这次课程设计综合运用本课程的理论知识进行频谱分析以及滤波器设计,通过理论推导得出相应结论,并利用MATLAB作为编程工具进行计算机实现,从而复习巩固了课堂所学的理论知识,提高了对所学知识的综合应用能力,并从实践上初步实现了对数字信号的处理。
基于语音信号去噪处理的FIR滤波器设计

基于语音信号去噪处理的FIR滤波器设计概述:在许多语音处理应用中,如通信系统、语音识别和语音合成等,语音信号往往受到各种噪声的污染影响。
为了减少或去除这些噪声的影响,需要设计有效的滤波器。
本文将介绍基于语音信号去噪处理的FIR滤波器的设计方法。
1.FIR滤波器的原理FIR滤波器是一种线性时不变系统,由一组滤波器系数和延迟单元组成。
其输出信号为输入信号与滤波器系数之间的卷积运算。
根据语音信号的频谱特性,可以设计一组合适的滤波器系数,用于去除特定频段的噪声。
2.FIR滤波器的设计方法2.1频率采样法频率采样法是一种直接设计FIR滤波器的方法,它通过在频域中指定所需的频率响应来设计滤波器。
首先,根据信号的采样率和带宽要求,确定滤波器的阶数和过渡带宽。
然后,在频域中绘制所需的频率响应曲线,并进行插值得到滤波器的频率响应。
最后,通过傅立叶逆变换得到滤波器的时域响应,即滤波器系数。
2.2窗函数法窗函数法是一种常用的FIR滤波器设计方法,它通过在时域中选择适当的窗函数来设计滤波器的频率响应。
常用的窗函数有矩形窗、汉宁窗、汉明窗等。
首先,根据信号的采样率和带宽要求,确定滤波器的阶数和过渡带宽。
然后,选择合适的窗函数,并在时域中将滤波器的频率响应与窗函数相乘。
最后,通过傅立叶逆变换得到滤波器的时域响应,即滤波器系数。
3.基于语音信号的去噪处理基于语音信号的去噪处理可以通过FIR滤波器来实现。
首先,通过对语音信号进行时频分析,分析出语音信号的频谱特性和噪声的频谱特性。
然后,根据噪声的频谱特性设计一个合适的FIR滤波器,使其在噪声频段上具有较高的增益,而在语音频段上具有较低的增益。
最后,将噪声信号输入到设计好的FIR滤波器中,并将滤波器的输出与原始语音信号相减,得到去噪后的语音信号。
4.总结本文介绍了基于语音信号去噪处理的FIR滤波器设计方法。
通过选择合适的设计方法和滤波器系数,可以有效地去除语音信号中的噪声。
在实际应用中,可以根据具体的需求和信号特性,选择合适的设计方法和优化算法,以获得更好的去噪效果。
FIR数字低通滤波器的(汉宁)窗函数法设计

)(9cos 15.0)(12cos 15.0)(1919n R n n R N n n w ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛--=ππ2.3进行语音信号的采集(1)按“开始”-“程序”-“附件”-“娱乐”-“录音机”的顺序操作打开Window s系统中的录音机软件。
如图1所示。
图1 wi ndows 录音机(2)用麦克风录入自己的声音信号并保存成wav 文件。
如图2所示。
图2 保存文件保存的文件按照要求如下:① 音信号文件保存的文件名为“y uxueji ao.wav ”。
②语音信号的属性为“8.000KH z,8位,单声道 7KB /秒” ,其它选项为默认。
2.4语音信号的分析将“yu xuejiao .wav ”语音文件复制到计算机装有Ma tlab 软件的磁盘中相应Mat lab目录中的“work ”文件夹中。
打开Matlab 软件,在菜单栏中选择“File ”-图3语音信号的截取处理图在图3中,其中第一个图为原始语音信号;第二个图是截短后的信号图。
图4频谱分析图其中第二个图是信号的FFT 结果,其横坐标的具体值是X (k)中的序号k;第三个图是确定滤波频率范围的参考图,其横坐标的具体值应当是遵循DFT 定义式和频率分辨率求得的:∑-===10)()]([)(N n k N W n x n x DFT k X π当k等于0时, 020j kn Njk knNe eW ==⋅-=π,从数字角频率上看,对应的正好是0=ω即直流的位置,也就是说,在取滤波频段时,当将主要能量(即红色框的部分)保留,其余频段部分的信号滤除。
)]([)(n x DFT k X =相当于是信号)(n x 的实际频谱)]([)(n x DFT ej X w =采样,而)(n x 又是连续时间语音信号)(t x 的采样。
)(k X 的每两个相邻取值之间的频率间隔大小对应到语音信号)(n x 的频谱中去,其频率间隔大小正好是采样结果的长度采样速率===∆L f f f s det f ∆称频率分辨率,其中Hz f s 8000=,10000=L ,p2=sum(s2.^2)-sum(s1.^2);SNR1=10*log10(p1/p2);p3=sum(s4.^2)/8000;p4=sum(s3.^2)/8000-sum(s4.^2)/8000;SNR2=10*log10(p3/p4);2.6 噪声叠加图5 语音信号与加噪声后语音信号对比图五为语音信号与加噪声后语音信号对。
基于语音信号去噪处理的FIR低通滤波器设计要点

摘要本次课程设计分析了FIR数字滤波器的基本原理,在MATLAB环境下利用窗函数设计FIR低通滤波器,实现了FIR低通滤波器的设计仿真。
本文根据滤波后的时域图和原始语音信号时域图的比较,以及滤波后信号的频谱图和原始语音信号频谱图的比较,最后回放滤波后语音信号,滤波后的语音信号与原始语音信号一样清晰,仿真结果表明,设计的FIR滤波器的各项性能指标均达到了指定要求,设计过程简便易行。
该方法为快速、高效地设计FIR滤波器提供了一个可靠而有效的途径。
关键词:DSP ;FIR;低通滤波器;语音信号;MATLAB目录第一章引言 (1)1.1 设计目的及意义 (1)1.2 设计任务及要求 (2)1.3 课程设计平台 (2)第二章基本原理 (3)2.1 FIR滤波器的基本概念 (3)2.2 FIR滤波器的特点 (3)2.3 FIR滤波器的种类 (4)第三章FIR数字低通滤波器的设计 (5)3.1 FIR低通滤波器设计原理 (5)3.2 FIR低通滤波器的设计方法 (5)3.2.1 频率采样法 (5)3.2.2 最优化设计 (6)3.2.3 窗函数法 (6)3.3 窗函数法设计步骤 (8)第四章详细设计 (9)4.1 语音信号的采集 (9)4.2 语音信号的读入与打开 (10)4.3 语音信号的FFT变换 (11)4.4 含噪信号的合成 (12)4.5 利用FIR滤波器滤波 (13)4.6 结果分析 (14)总结 (15)参考文献 (16)附录 (17)致谢 (21)第一章引言随着信息科学和计算机技术的不断发展,数字信号处理(DSP,Digital Signal Processing)的理论和技术也得到了飞速的发展,并逐渐成为一门重要的学科,它的重要性在日常通信、图像处理、遥感、声纳、生物医学、地震、消费电子、国防军事、医疗方面等显得尤为突出。
在我们面临的信息革命中,数字信号处理几乎涉及了所有的工程技术领域。
数字信号处理是一种将信号以数字形式进行处理的一种理论和技术,它的目的是将真实世界中的一些信号进行分析并滤波,最后得出其中的有用的信号。
FIR滤波器的设计与性能评估

FIR滤波器的设计与性能评估一、引言滤波器在信号处理中起到了至关重要的作用。
滤波器可以根据信号的频率特性对信号进行处理,使我们可以去除噪声、增强感兴趣的频段等操作。
本文将介绍FIR(Finite Impulse Response)滤波器的设计原理和性能评估方法。
二、FIR滤波器的设计方法FIR滤波器是一种经典的数字滤波器,它利用有限的输入响应对输入信号进行滤波处理。
FIR滤波器具有线性相位和稳定性等优点,因此在许多应用中得到广泛应用。
1. 理想低通滤波器设计首先,我们需要确定FIR滤波器的设计参数,其中最基本的是滤波器的类型。
假设我们需要设计一个低通滤波器,即只保留低于一定频率的信号分量。
可以采用理想低通滤波器的方法进行设计。
2. 频率响应的离散化接下来,我们需要将理想低通滤波器的频率响应离散化,得到滤波器的系数。
常用的方法有频率采样法和窗函数法。
频率采样法通过在频域上均匀采样理想滤波器的频率响应得到系数,而窗函数法则需要选择一个窗函数来对离散化后的频率响应进行加窗。
3. 系数计算与滤波器实现根据离散化后的频率响应,可以通过逆变换得到滤波器的系数。
然后,我们可以将这些系数用于实现FIR滤波器。
常见的实现方式包括直接形式(Direct Form)、级联形式(Cascade Form)和线性相位形式(Linear Phase Form)等。
三、FIR滤波器的性能评估方法设计完成后,我们需要对FIR滤波器进行性能评估,以确保其能够满足我们的需求。
1. 幅频响应和相频响应在性能评估中,我们通常关注滤波器的幅频响应和相频响应。
幅频响应可以反映滤波器对不同频率分量的衰减或增益情况,而相频响应则描述了信号在滤波器中的相位变化。
2. 截止频率和过渡带宽对于低通滤波器而言,截止频率和过渡带宽是评估性能的重要指标。
截止频率是指滤波器开始起作用的频率,而过渡带宽则是指截止频率和衰减区域之间的频率范围。
3. 线性相位特性FIR滤波器具有线性相位的特点,这意味着不同频率分量的信号在滤波器中的延迟是相同的。
语音信号的FIR滤波器处理.

(语音信号的FIR滤波器处理)课程名称:DSP原理及应用报告提交日期2011年12月13日项目答辩日期2011年12月13日1 导论以及设计要求1.1引言随着信息时代和数字世界的到来,数字信号处理已成为如今一门极其重要的学科和技术领域。
数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。
数字信号处理器,也称DSP芯片,是针对数字信号处理需要而设计的一种具有特殊结构的微处理器,它是现代电子技术、相结合的产物。
一门主流技术,随着信息处理技术的飞速发展,计算机技术和数字信号处理技术逐渐发展成为使它在电子信息、通信、软件无线电、自动控制、仪表技术、信息家电等高科技领域得到了越来越广泛的应用。
数字滤波是语音处理、图像处理、频谱分析等应用中的基本处理算法。
DSP 是一种处理数字信号的专用微处理器, 主要应用于实时快速地实现各种信号的数字处理算法。
用 DSP 芯片实现数字滤波具有稳定性好、精确度高、不受环境影响等优点。
数字滤波器分为有限冲激响应滤波器(FIR滤波器)和无限冲激响应滤波器(IIR滤波器)。
FIR 滤波器属于经典滤波器,优点就是由于不存在系统极点,FIR 滤波器是绝对稳定的系统,FIR 滤波器还确保了线性相位,在信号处理中占有极其重要的地位。
数字滤波器一直以来就是数字信号处理器(DSP)最广为人知的应用,FIR 滤波器的单位冲激响应 b(n)为有限长序列,若 b(n)为实数, 且满足偶对称:b(n)= b(N- 1- n)的条件, 称为系数对称FIR 滤波器。
系数对称 FIR 滤波器在数字信号处理中应用十分广泛。
1.2设计的目的通过课程设计,加深对DSP芯片TMS320C54x的结构、工作原理的理解,获得DSP应用技术的实际训练,掌握设计较复杂DSP系统的基本方法。
通过使用汇编语言编写具有完整功能的图形处理程序或信息系统,使学生加深对所学知识的理解,进一步巩固汇编语言讲法规则。
数字信号处理课程设计-语音信号的数字滤波——fir数字滤波器的(汉宁)窗函数设计[管理资料]
![数字信号处理课程设计-语音信号的数字滤波——fir数字滤波器的(汉宁)窗函数设计[管理资料]](https://img.taocdn.com/s3/m/5ad56004551810a6f42486cc.png)
语音信号的数字滤波——FIR数字滤波器的(汉宁)窗函数设计设计题目:语音信号的数字滤波——FIR数字滤波器的(汉宁)窗函数设计一、课程设计的目的通过对常用数字滤波器的设计和实现,掌握数字信号处理的工作原理及设计方法;掌握利用数字滤波器对信号进行滤波的方法。
并能够对设计结果加以分析。
二、设计步骤⑴按“开始->程序->附件->娱乐->录音机”的顺序操作打开Windows系统中的录音机软件⑵用麦克风录入自己的声音信号并保存成文件(语音信号的长度不得少于1秒)⑶记录以下内容:语音信号文件保存的文件名为“”、格式PCM,8位,单声道,如图1和图2所示。
语音信号的采样速率为8000Hz/s。
图1语音信号的采集图2 原始语音信号⑴将上一步骤中保存下来的语音信号文件“*.wav”复制到计算机装有Matlab 软件的磁盘中相应Matlab目录中的“work”文件夹中⑵双击桌面上Matlab软件的快捷图标,打开Matlab软件⑶在菜单栏中选择“File->new->M-File”或是点击快捷按钮,打开m 文件编辑器⑷在m文件编辑器中输入相应的指令将自己的语音信号导入Matlab工作台。
程序部分首先用语音文件将自己的录音导入,指令为wavread(),本设计中为waveread(’’),然后将处理后的语音信号导出,指令为:wavwrite(‘’);本设计中录入的是单声道语音。
一般情况下录入的双声道语音信号中()右导入交保存为变量后,其变量应当是一个二列的二维数组,其中每一列对应一个声道,数组的行数等于采样速率与时间的乘积(即单声道的采样点数);本课程设计过程中的语音原始信号存为“”;截短后的输出语音为:“”;叠加噪声后的语音为“”。
具体程序段见小标题⑹,频谱分析如下图3。
图3截短后语音信号的时域和频域波形如上图3所示,上面的图是原始声音截去大部分空白后的截短语音,这样有利于频谱分析;中间的图是截短后的声音在频域的分析,首先对语音进行采样,采样频率大于信号最高频率的2倍即可。
基于FIR滤波器的语音信号处理设计_1

毕业设计选题:基于FIR滤波器的语音信号处理设计设计思路:1.介绍数字滤波器的概念及分类主要分数字滤波器的概念、发展现状、按照不同性能的分类,突出IIR和FIR滤波器。
2.重点介绍FIR滤波器的设计原理主要从FIR的设计原理(理论)、性能指标、分类、设计方法(重点,一共三种,分别介绍其理论和实现方法)以及优缺点入手分析3.介绍语音信号的处理过程及原理简单介绍一下采样原理,再介绍语音信号的采集流程(录音、转码、采样、打开、回放)4.简单介绍MATLAB的功能及使用方法简单介绍MATLAB在信号处理方面的功能(信号的采集、频谱分析、滤波器的设计)以及它的使用方法5.设计一个能处理声音信号的FIR滤波器全文的核心,先录制一段语音,通过MATLAB采集显示,分析时域波形及频谱显示;再设置滤波后的性能指标,用MATLAB分别用三种方法设计出FIR滤波器;将设计好的滤波器加入原语音信号中,观察处理后的时域波形及频谱;最后分析比较三种方法以及与无滤波器的显示比较,得出结论。
(如果可以,首先可以在原音中加入噪音,观察加入滤波器前后波形的差异,说明滤波器在语音信号处理方面的作用;只通过编程来完成滤波器的设计,还不够形象,最好可以通过Simulink将整个系统仿真出来,说明结果的可靠性。
)6.简单分析一下所设计滤波器的优缺点由于FIR滤波器设计原理本身存在一些误差,导致仿真出来的结果必然存在一些误差,分析这些误差的来源,再说明FIR 滤波器在语音信号处理方面的优点,得出结论。
7. 总结全文的思路和研究方法总结在毕业设计过程中遇到的问题,以及解决的方法;最后感谢老师的指导,通过毕业设计学到了哪些东西。
8. 参考文献主要是图书馆借阅的书籍、翻看的论文和报刊,还有网上资料的文献来源。
(每引用一些知识点,就注明其出处,但不能完全照抄,按照自己的话来重新组织。
)9. 附录将MATLAB 的编程代码全部附上,如果有其他的图表也附上研究对象:FIR 滤波器辅助对象:语音信号工具软件:MATLAB FilterMATLAB频谱显示+语音信号噪 声语音信号的滤波流程图摘要本文设计介绍了基于Matlab的对语音信号采集、处理及FIR滤波器的设计,并使之实现的过程。
(完整版)基于MATLAB的FIR滤波器语音信号去噪

*****************实践教学******************兰州理工大学计算机与通信学院2013年春季学期《信号处理》课程设计题目:基于MATLAB的FIR滤波器语音信号去噪专业班级:姓名:学号:指导教师:成绩:摘要本次课程设计是基于MATLAB的FIR滤波器语音信号去噪,在设计过程中,首先录制一段不少于10秒的语音信号,并对录制的信号进行采样;其次使用MATLAB会出采样后的语音信号的时域波形和频谱图;然后在给原始的语音信号叠加上噪声,并绘出叠加噪前后的时域图及频谱图;再次设计FIR滤波器,针对语音信号的性质选取一种适合的窗函数设计滤波器进行滤波;最后对仿真结果进行分析。
设计出的滤波器可以满足要求。
关键词: FIR滤波器;语音信号;MATLAB仿真目录一 FIR滤波器设计的基本原理 (1)1.1滤波器的相关介绍 (1)1.1.1数字滤波器的概念 (1)1.1.2 IIR和FIR滤波器 (1)1.2利用窗函数法设计FIR滤波器 (1)1.2.1窗函数法设计FIR滤波器的基本思想 (1)1.2.2窗函数法设计FIR滤波器的步骤 (2)1.2.2窗函数法设计FIR滤波器的要求 (2)1.2.3常用窗函数的性质和特点 (3)1.2.4 语音处理中的采样原理 (3)二语音信号去噪实现框图 (5)三详细设计 (7)3.1 信号的采集 (7)3.2 语音信号的读入与打开 (7)3.3 语音信号的FFT变换 (8)3.4含噪信号的合成 (9)3.5 FIR滤波器的设计 (10)3.6 利用FIR滤波器滤波 (11)3.7 结果分析 (14)总结 (15)参考文献 (16)附录 (17)致谢 (21)一 FIR滤波器设计的基本原理1.1滤波器的相关介绍1.1.1数字滤波器的概念数字滤波器(Digital Filter,简称为DF)是指用来对输入信号进行滤波的硬件和软件。
所谓数字滤波器,是指输入、输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。
语音信号的FIR滤波器处理课程设计

DSP课程设计报告语音信号的FIR滤波器处理姓名:班级: 11电信4 学号:指导老师:日期: 2014.06.03~2014.06.13摘要DSP( Digital Signal Processing) 也就是我们常说的数字信号处理, 它是利用计算机或专用处理设备, 以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理, 以得到符合人们需要的信号形式。
DSP 芯片在信号处理、通信、雷达等许多领域得到广泛的应用。
Matlab是一款强大的软件,它将数值分析、矩阵计算、科学数据可视化等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。
CSS集成开发环境使用CCS内置的软件仿真simulator对程序进行编译,调试和运行,主要用于检测目标程序运行的正确性和连贯性,并能通过仿真器与目标板连接,在目标板上实时观察效果。
在本次设计中,我们选择的课题是基于DSP的语言信号的FIR滤波处理。
首先利用MATLAB进行了仿真,得到滤波前后的时域波形和频谱。
然后通过调用MATLAB的分析工具FDATOOL,根据仿真结果导出了滤波器的相关参数,将原始信号数据和滤波器参数输入CCS进行DSP编程。
最后在DSP中实现了FIR低通滤波,并通过CCS的频谱分析功能查看了最终DSP的滤波效果。
关键词: 语音信号 DSP FIR滤波 MATLAB CCS目录1.前言 (4)2.设计原理 (4)2.1数字信号处理器 (4)2.2滤波器简介 (5)2.3 fir滤波器原理 (6)2.4窗函数简介 (7)2.5滤波器的特点 (7)3.FIR滤波器的MATLAB设计 (8)3.1总体方案的设计 (8)3.2语音信号的采集 (9)3.3MATLAP的具体处理 (9)3.4FIR滤波器的系统参数 (11)4.FIR滤波器的CCS设计与仿真 (11)4.1CCS程序流程图 (12)4.2CCS仿真结果 (15)5.问题及解决办法 (18)6.设计感想 (19)参考文献 (19)7.附录 (20)1 前言随着信息与数字技术的发展,数字信号处理已经成为当今极其重要而学科与技术领域之一。
FIR带通滤波器设计

1.2利用窗函数法设计FIR滤波器
1.2.1窗函数法设计FIR滤波器的基本思想
窗函数法的基本思想是要选取某一种合适的理想频率选择性滤波器,这种滤波器总是有一个非因果,无限长的脉冲响应,然后将它的脉冲响应截断(或加窗)以得到一个线性相位和因果的FIR滤波器,因为必须设计一个因果可实现的FIR滤波器[1]。
(5)布莱克曼窗
布莱克曼窗的幅度函数主要由五部分组成,他们的位移都不同,其幅度也是不同的W (w)使旁瓣再进一步抵消。旁瓣峰值幅度进一步增加,其幅度谱主瓣宽度是矩形窗的3倍。设计程序时用backman函数调用。
(6)凯塞窗
以上几种窗函数是各以一定主瓣加宽为代价,来换取某种程度的旁瓣抑制,窗函数的主瓣宽度和旁瓣峰值衰耗是矛盾的,一项指标的提高总是以另一项指标的下降为代价,窗口选择实际上是对两项指标作权衡。而两项指标是跳变的,于是有人提出可调整窗,适当修改参数,可在这两项指标间作连续的选择。常用的可调整窗是凯塞(Kaiser)窗。而凯窗则是全面地反映主瓣与旁瓣衰减之间的交换关系,可以在它们两者之间自由地选择它们的比重。
频带为F的连续信号 f(t)可用一系列离散的采样值f(t1),f(t1±Δt),f(t1±2Δt),...来表示,只要这些采样点的时间间隔Δt≤1/2F,便可根据各采样值完全恢复原来的信号f(t)。这是时域采样定理的一种表述方式。时域采样定理的另一种表述方式是:当时间信号函数f(t)的最高频率分量为fm时,f(t)的值可由一系列采样间隔小于或等于1/2fm的采样值来确定,即采样点的重复频率f≥2fm。
课程设计--语音信号滤波去噪使用BOHMAN窗设计的频率采样型FIR滤波器

课程设计--语音信号滤波去噪使用BOHMAN窗设计的频率采样型FIR滤波器第2页共23页语音信号滤波去噪——使用BOHMAN窗设计的频率采样型FIR滤波器学生姓名:指导老师:摘要本课程设计主要内容是,利用窗口设计法设计频率采样型FIR滤波器并对加噪语音信号进行滤波去噪处理。
仿真平台为MATLAB7.0,开发工具是M语言编程。
在课程设计中,首先利用Windows下的录音机工具录制一段格式为.wav的语音信号,然后在MATLAB中对语音信号进行加噪,并绘制原始语音信号和加噪语音信号的时域和频域波形,进行频谱分析以确定所加噪声频率,再利用Bohman窗设计FIR滤波器,并检测是否达到指标,最后使用滤波器对信号进行滤波去噪处理,并通过对比原始信号,加噪信号,滤波去噪信号的时域和频域波形,或回放语音信号,检测是否设计成功。
通过程序调试及完善,本课程设计滤波后的语音信号与原始语音信号基本一致,即设计的滤波器能够从含噪信号中滤除单频噪声,还原原始信号,达到了设计目的。
关键词课程设计;滤波去噪;FIR滤波器;BOHMAN窗;MATLAB1 引言本课程设计主要解决在含噪情况下如何对语音信号进行滤波去噪处理。
设计采用BOHMAN窗设计的频率采样型FIR滤波器对一段含噪语音信号进行滤波去噪处理。
数字滤波器是一种用来过滤时间离散信号的数字系统,它是通过对抽样数据进行数字处理来达到频域滤波的目的。
数字信号处理与模拟信号处理是信号处理的子集。
数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。
随着现代通信的数字化,数字滤波器变得更加重要。
数字滤波器的种类很多,但总的来说可以分成两大类,一类是经典滤波器,另一类可称为现代滤波器。
从滤波特性方面考虑,数字滤波器可分为数字高通、数字低通、第3页共23页数字带通和数字帯阻等滤波器。
从实现方法上考虑,将滤波器分成两种,一种称为无限脉冲响应滤波器,简称IIR滤波器,另一种为FIR滤波器。
基于MATLAB的FIR滤波器语音信号去噪_课程设计

*****************实践教学******************兰州理工大学计算机与通信学院《信号处理》课程设计题目:基于MATLAB的FIR滤波器语音信号去噪摘要本次课程设计是基于MATLAB的FIR滤波器语音信号去噪,在设计过程中,首先录制一段不少于10秒的语音信号,并对录制的信号进行采样;其次使用MATLAB会出采样后的语音信号的时域波形和频谱图;然后在给原始的语音信号叠加上噪声,并绘出叠加噪前后的时域图及频谱图;再次设计FIR滤波器,针对语音信号的性质选取一种适合的窗函数设计滤波器进行滤波;最后对仿真结果进行分析。
设计出的滤波器可以满足要求。
关键词: FIR滤波器;语音信号;MATLAB仿真目录一 FIR滤波器设计的基本原理 (1)1.1滤波器的相关介绍 (1)1.1.1数字滤波器的概念 (1)1.1.2 IIR和FIR滤波器 (1)1.2利用窗函数法设计FIR滤波器 (1)1.2.1窗函数法设计FIR滤波器的基本思想 (1)1.2.2窗函数法设计FIR滤波器的步骤 (2)1.2.2窗函数法设计FIR滤波器的要求 (2)1.2.3常用窗函数的性质和特点 (3)1.2.4 语音处理中的采样原理 (3)二语音信号去噪实现框图 (5)三详细设计 (7)3.1 信号的采集 (7)3.2 语音信号的读入与打开 (7)3.3 语音信号的FFT变换 (8)3.4含噪信号的合成 (9)3.5 FIR滤波器的设计 (10)3.6 利用FIR滤波器滤波 (11)3.7 结果分析 (14)总结 (15)参考文献 (16)附录 (17)致谢 (21)一 FIR滤波器设计的基本原理1.1滤波器的相关介绍1.1.1数字滤波器的概念数字滤波器(Digital Filter,简称为DF)是指用来对输入信号进行滤波的硬件和软件。
所谓数字滤波器,是指输入、输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。
用FIR低通滤波器对语音信号进行处理

一、实验目的对一段已知的音频进行信号处理,将音频中2500Hz的噪声用低通滤波器滤除。
二、实验思路1.用MATLAB设计滤波器系数,确定滤波器阶数;2.利用CCS设计FIR滤波器,将滤波器系数添加进来;3.用MATLAB将音频wav转换成CCS可以读的输入dat文件;4.将输入dat文件作为输入,经过卷积运算,输出输出dat文件;5.再将输出dat文件用MATLAB还原成语音wav格式。
三、实验步骤1.用MATLAB设计FIR滤波器系数(1)键入fdatool(2)选择lowpass(3)选择FIR滤波器(4)Specify ordor(设计的是65阶)(5)Window选择hamming窗(6)采样频率fs:8000Hz(7)截止频率fc:(2000~2200)设计为2100Hz(8)点击design(9)File、 export、 Mat-file、 export、保存2.用CCS设计FIR滤波器,程序如下:.title "FIR4.asm".mmregs.def _c_int00.bss y,1 ;给y分配一个存储单元xn .usect "xn",66 ;定义xn的缓冲区域a0 .usect "a0",66PA0 .set 0 ;I/O口地址赋值PA1 .set 1.datatable .word -5 ;定义滤波器系数.word 27.word 1.word -35.word 5.word 49.word -17.word -70.word 38.word 94.word -72.word -120.word 123.word 144.word -195.word -160.word 292.word 161.word -416.word -139.word 79.word -777.word 35.word 1039.word -240.word -1408.word 614.word 2014.word -1417.word -3421.word 4282.word 15302.word 15302.word 4282.word -3421.word -1417.word 2014.word 614.word -1408.word -240.word 1039.word 35.word -777.word 79.word 575.word -139.word -416 .word 161 .word 292.word -160.word -195.word 144.word 123.word -120.word -72.word 94.word 38.word -70.word -17.word 49.word 5.word -35.word 1.word -5.text_c_int00: SSBX FRCTSTM #a0,AR1 ;AR1指向a0的首地址RPT #65 ;重复执行下条语句65次MVPD table,*AR1+ ;从table首址重复传递66个数据STM #xn+65,AR3 ;AR3指向x(n-65)STM #a0+65,AR4 ;AR4指向a65STM #66,BK ;循环缓冲区长度BK=5STM #-1,AR0 ;AR0=-1,双操作数减量LD #xn,DPPORTR PA1,@xn ;输入xnSTM #1000h,AR6 ;输出数据缓冲区首址为#1000hSTM #9fffh,AR7 ;循环计算40959个样本点FIR: RPTZ A,#65 ;A清0,共迭代65次MAC *AR3+0%,*AR4+0%,A ;双操作数乘法累加STH A,@y ;保存y(n)PORTW @y,PA0 ;输出y(n);BD FIRPORTR PA1,*AR3+0% ;输入新数据STH A,*AR6+ ;保存y(n)BANZ FIR,*AR7- ;循环40959次end: B end.END3.用MATLAB将音频wav转换成CCS可以读的input.dat文件,程序如下:y=wavread('E:\SpeechAndNosie.wav');%读入wav数据fid=fopen('E:\input.dat','w');%打开input.dat文件(没有就自动创建),获取文件ID fprintf(fid,'1651?1?0?1?0?\n');%写入dat文件的头几个参数fprintf(fid,'0x%x \n',round(y*100)+(y<0)*2^16);%写入音频信号的数据,以16进制显示,负数用补码显示4.将input.dat文件作为输入,经过卷积运算,输出output.dat文件;5.再将output.dat文件用MATLAB还原成语音wav格式,程序如下:fid=fopen('E:\DSPprogram\FIR4\output.dat','rt');%打开output.dat文件,获取文件ID(fid)y=fscanf(fid,'%x \n');%将数据以16进制写入到变量y中y1=(y-(y>2^15)*(2^16))/100;%16进制数转换为小数wavwrite(y1,8000,'E:\DSP program\FIR4\output.wav');%输出output.wav文件四、实验结果1.FIR低通滤波器对语音信号的响应及频谱:可以看到在2500Hz处,滤波器对语音信号的衰减很大2.通过处理后的语音与原始语音进行比较,可以明显的发现声音清晰,没有噪声的干扰。
FIR滤波器设计与实现实验报告

FIR滤波器设计与实现实验报告目录一、实验概述 (2)1. 实验目的 (3)2. 实验原理 (3)3. 实验设备与工具 (4)4. 实验内容与步骤 (6)5. 实验数据与结果分析 (7)二、FIR滤波器设计 (8)1. 滤波器设计基本概念 (9)2. 系数求解方法 (10)频谱采样法 (11)最小均方误差法 (14)3. 常用FIR滤波器类型 (15)线性相位FIR滤波器 (16)非线性相位FIR滤波器 (18)4. 设计实例与比较 (19)三、FIR滤波器实现 (20)1. 硬件实现基础 (21)2. 软件实现方法 (22)3. 实现过程中的关键问题与解决方案 (23)4. 滤波器性能评估指标 (25)四、实验结果与分析 (26)1. 实验数据记录与处理 (27)2. 滤波器性能测试与分析 (29)通带波动 (30)虚部衰减 (31)相位失真 (32)3. 与其他设计方案的对比与讨论 (33)五、总结与展望 (34)1. 实验成果总结 (35)2. 存在问题与不足 (36)3. 未来发展方向与改进措施 (37)一、实验概述本次实验的主要目标是设计并实现一个有限脉冲响应(Finite Impulse Response,简称FIR)滤波器。
FIR滤波器是数字信号处理中常用的一种滤波器,具有线性相位响应和易于设计的优点。
本次实验旨在通过实践加深我们对FIR滤波器设计和实现过程的理解,提升我们的实践能力和问题解决能力。
在实验过程中,我们将首先理解FIR滤波器的基本原理和特性,包括其工作原理、设计方法和性能指标。
我们将选择合适的实验工具和环境,例如MATLAB或Python等编程环境,进行FIR滤波器的设计。
我们还将关注滤波器的实现过程,包括代码编写、性能测试和结果分析等步骤。
通过这次实验,我们期望能够深入理解FIR滤波器的设计和实现过程,并能够将理论知识应用到实践中,提高我们的工程实践能力。
本次实验报告将按照“设计原理设计方法实现过程实验结果与分析”的逻辑结构进行组织,让读者能够清晰地了解我们实验的全过程,以及我们从中获得的收获和启示。
fir滤波器的设计实验报告

fir滤波器的设计实验报告fir滤波器的设计实验报告引言:滤波器是信号处理中常用的工具,它可以对信号进行去噪、降噪、频率调整等操作。
在本次实验中,我们将设计一种fir滤波器,通过对信号进行滤波处理,实现对特定频率成分的增强或抑制。
本报告将详细介绍fir滤波器的设计原理、实验步骤和结果分析。
一、设计原理:fir滤波器是一种无限冲激响应滤波器,其特点是具有线性相位和稳定性。
其基本原理是通过对输入信号和滤波器的冲激响应进行线性卷积运算,得到输出信号。
fir滤波器的冲激响应由一组有限长的系数决定,这些系数可以通过不同的设计方法得到,如窗函数法、最小二乘法等。
二、实验步骤:1. 确定滤波器的频率响应需求:根据实际应用需求,确定滤波器需要增强或抑制的频率范围。
2. 选择滤波器的设计方法:根据频率响应需求和系统要求,选择合适的fir滤波器设计方法。
3. 设计滤波器的冲激响应:根据所选设计方法,计算得到fir滤波器的冲激响应系数。
4. 实现滤波器的数字滤波器:将fir滤波器的冲激响应系数转换为差分方程,得到数字滤波器的差分方程表示。
5. 实现滤波器的数字滤波器:将fir滤波器的冲激响应系数转换为差分方程,得到数字滤波器的差分方程表示。
6. 通过编程实现滤波器:使用编程语言(如MATLAB)编写代码,实现fir滤波器的数字滤波器。
7. 信号滤波处理:将待滤波的信号输入到fir滤波器中,通过数字滤波器进行滤波处理,得到输出信号。
8. 结果分析:对滤波后的信号进行分析,评估滤波器的性能和效果。
三、实验结果分析:在本次实验中,我们设计了一个fir滤波器,并对一段音频信号进行滤波处理。
通过实验结果分析,我们发现滤波器能够有效地增强或抑制指定频率范围内的信号成分。
滤波后的音频信号听起来更加清晰,噪音得到了有效的抑制。
同时,我们还对滤波器的性能进行了评估。
通过计算滤波器的幅频响应曲线和相频响应曲线,我们发现滤波器在指定频率范围内的增益和相位变化符合预期。
基于FIR滤波的语音信号处理系统设计

信号与系统课程设计报告课题五基于FIR滤波的语音信号处理系统设计课题五基于FIR滤波的语音信号处理系统设计摘要:本课题研究数字语音信号的初步分析方法、FIR数字滤波器的设计及应用。
Matlab中实现对语音信号进行采集(读取),播放以及重采样的功能。
对于加入噪声的语音信号设计合适的FIR数字滤波器进行滤波处理,并通过FFT 变换对语音信号的频谱进行分析。
还通过enframe函数实现对语音信号的分帧,计算出短时能量和短时平均过零率等参数,通过时域波形和自相关的处理可以看到清音帧和浊音帧的区别。
本次设计中还采用LabVIEW进行仿真设计,实现了包括语音文件的读取及播放、滤波前后语音信号时域波形及频谱的显示、滤波器设计、滤波器幅频特性的显示等功能。
最后,通过Express声音采集VI实现采集语音信号的功能,并对采集的语音信号进行后续分析。
本设计很好的完成了设计课题的要求,并在一定程度上进行了延伸。
关键词:语音信号采样 FFT FIR滤波器频谱分析 LabVIEW 仿真Abstract:This paper researches the initial analysis of the digital voice signal and the design and application of the digital filter.We use Matlab to achieve the function,such as the acquisition of voice,play the voice and the resample of voice signal.We design the proper FIR digital filter to obtain the pure signal and analyze the voice spectrum by FFT transformation.We also use “enframe” to achieve the framing of voice signal and calculate the short-time energy.We can see the distinguishbetween unvoiced speech frame and the voice frame through the voice transaction.Also,the paper introduce the LabVIEW simulation design.The design realize the functions such as read and play the audio files,display of the time domain waveform and the spectrum, filter design,display of the filter amplitude-frequency characteristic.Finally,we use VI named sound acquisition to collect sound signal and analyze.The research completes the requirement of the task well and make more further study.Keywords: Voice signal sampling FFT FIR filter spectral analysis LabVIEW simulation1课程设计介绍1.1 课程设计目的本设计课题主要研究数字语音信号的初步分析方法、FIR数字滤波器的设计及应用。
fir数字滤波器的设计与实现

fir数字滤波器的设计与实现一、引言数字滤波器是数字信号处理中的重要组成部分,它可以用于去除信号中的噪声,平滑信号等。
其中,fir数字滤波器是一种常见的数字滤波器。
本文将介绍fir数字滤波器的设计与实现。
二、fir数字滤波器概述fir数字滤波器是一种线性相位、有限脉冲响应(FIR)的数字滤波器。
它通过一系列加权系数对输入信号进行卷积运算,从而实现对信号的过滤。
fir数字滤波器具有以下特点:1. 稳定性好:由于其有限脉冲响应特性,使得其稳定性优于IIR(无限脉冲响应)数字滤波器。
2. 线性相位:fir数字滤波器在频域上具有线性相位特性,因此可以保持输入信号中各频率分量之间的相对时延不变。
3. 设计灵活:fir数字滤波器可以通过改变加权系数来实现不同的频率响应和截止频率。
三、fir数字滤波器设计步骤1. 确定需求:首先需要确定所需的频率响应和截止频率等参数。
2. 选择窗函数:根据需求选择合适的窗函数,常用的有矩形窗、汉明窗、布莱克曼窗等。
3. 计算滤波器系数:利用所选窗函数计算出fir数字滤波器的加权系数。
常见的计算方法有频率采样法、最小二乘法等。
4. 实现滤波器:将计算得到的加权系数应用于fir数字滤波器中,实现对信号的过滤。
四、fir数字滤波器实现方法1. 直接形式:直接将计算得到的加权系数应用于fir数字滤波器中,实现对信号的过滤。
该方法简单易懂,但是需要大量运算,不适合处理较长的信号序列。
2. 快速卷积形式:利用快速傅里叶变换(FFT)来加速卷积运算。
该方法可以大大减少计算量,适合处理较长的信号序列。
五、fir数字滤波器应用案例1. 语音处理:fir数字滤波器可以用于去除语音信号中的噪声和杂音,提高语音质量。
2. 图像处理:fir数字滤波器可以用于图像去噪和平滑处理,提高图像质量。
3. 生物医学信号处理:fir数字滤波器可以用于生物医学信号的滤波和特征提取,如心电信号、脑电信号等。
六、总结fir数字滤波器是一种常见的数字滤波器,具有稳定性好、线性相位和设计灵活等优点。
语音信号滤波去噪——使用凯塞窗设计的FIR滤波器

摘要本课程设计的主要目的是在MATLAB环境下,使用窗口设计法设计一个滤波器,并对语音信号进行滤波去噪。
开发平台为MATLAB7.0,设计方法为窗口设计法。
用麦克风采集一段语音信号,绘制波形并观察其频谱,给定相应技术指标,用凯塞窗设计一个满足指标的FIR滤波器,对该语音信号进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析。
在滤波前后语音信号波形的变化中,由于我们录制的语音信号噪声不大,所以观察并不明显,但在频域波形中,我们可以明显的看到设计的滤波器对语音信号进行了滤波处理,将噪声进行了滤除。
此次滤波基本达到了要求,完成了设计指标。
关键词滤波去噪;FIR滤波器;凯塞窗;MATLAB1 引言本课程设计是采用凯塞窗设计的FIR滤波器对语音信号滤波去噪。
通过课程设计了解FIR滤波器的原理及使用方法,了解使用MATLAB语言设计FIR滤波器的方法,了解DSP对FIR滤波器的设计及编程方法。
通过观察滤波前后的时域图形,加深对FIR 滤波器作用的理解。
通过对比滤波器前后的波形及回放滤波前后的语音信号,可以看出滤波器对有用信号的无失真放大具有重要作用。
语音信号是基于时间轴上的一维数字信号,在这里主要是对语音信号进行频域上的分析。
在信号分析中,频域往往包含了更多的信息。
对于频域来说,大概有8种波形可以让我们分析:矩形方波,锯齿波,梯形波,临界阻尼指数脉冲波形,三角波,余旋波,余旋平方波,高斯波。
对于各种波形,我们都可以用一种方法来分析,就是傅立叶变换:将时域的波形转化到频域来分析。
于是,本课题就从频域的角度对一段频率为22050Hz,16k的语音信号进行分析,并通过分析频谱来设计出合适的滤波器。
由于MATLAB软件在数字信号处理上有巨大的优势,本课程设计除了录制语音信号外基本都在MATLAB上实现[1]。
本课程设计主要解决在一个综合型超市中,进行某种商品(如某一品牌的服装)的采购、领料、库存与成本等管理的程序设计。
基于语音信号去噪处理的FIR滤波器设计

摘要本次课程设计的题目是基于语音信号去噪处理的FIR滤波器设计,主要是利用MATLAB设计一个FIR数字带通滤波器,对一段噪声环境下的语音信号进行滤波。
在设计过程中,首先要录制一段语音信号,并对录制好的信号进行时域和频域分析;然后对原始的语音信号进行加噪处理,对加噪后的信号进行分析,绘出时域和频域的分析图;最后利用设计出的FIR带通滤波器,针对语音信号的性质选取一种适合的窗函数设计滤波器进行滤波,最后对仿真结果进行分析。
关键词: FIR带通滤波器;语音信号;MATLAB仿真;加噪;滤波目录前言 1第一章基本原理 21.1 语音处理中的采样原理 21.2数字滤波器的设计 21.3窗函数法 4第二章 FIR带通滤波器设计 62.1 FIR滤波器简介 62.2 FIR带通滤波器设计要求 72.3 设计方法 72.4 设计步骤 9第三章 FIR带通滤波器的软件仿真 103.1 程序流程图 103.2 仿真结果及分析 11总结 15参考文献 16附录 17致谢 22前言在信号处理过程中,所处理的信号往往混有噪音,从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。
根据有用信号和噪音信号的不同特性,提取有用信号的过程称为滤波。
实现滤波功能的系统被称为滤波器。
在近代电信设备和各类控制系统中,数字滤波器应用极为广泛。
语音处理是最早应用数字滤波器的领域之一,也是最早推动数字信号处理理论发展的领域之一。
该领域主要包括5个方面的内容:第一,语音信号分析。
即对语音信号的波形特性,统计特性,模型参数等进行分析计算。
第二,语音合成。
即利用专用硬件或在通用计算机上运行软件来产生语音。
第三,语音识别。
即利用专用硬件或计算机识别人的讲话,或者识别说话的人。
第四,语音增强。
即从噪音或者干扰中提取被掩盖的语音信号。
第五,语音编码。
主要用于语音数据的压缩,目前已经建立了一系列语音编码国际标准,大量用于通信和音频处理。
MATLAB是一种面向科学和工程计算的语言,它集数值分析、矩阵运算、信号处理和图形显示于一体,具有编程效率高、调试手段丰富、扩充能力强等特点。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
摘要本次课程设计分析了FIR数字滤波器的基本原理,在MATLAB环境下利用窗函数设计FIR低通滤波器,实现了FIR低通滤波器的设计仿真。
本文根据滤波后的时域图和原始语音信号时域图的比较,以及滤波后信号的频谱图和原始语音信号频谱图的比较,最后回放滤波后语音信号,滤波后的语音信号与原始语音信号一样清晰,仿真结果表明,设计的FIR滤波器的各项性能指标均达到了指定要求,设计过程简便易行。
该方法为快速、高效地设计FIR滤波器提供了一个可靠而有效的途径。
关键词:DSP ;FIR;低通滤波器;语音信号;MATLAB目录第一章引言 (1)1.1 设计目的及意义 (1)1.2 设计任务及要求 (2)1.3 课程设计平台 (2)第二章基本原理 (3)2.1 FIR滤波器的基本概念 (3)2.2 FIR滤波器的特点 (3)2.3 FIR滤波器的种类 (4)第三章FIR数字低通滤波器的设计 (5)3.1 FIR低通滤波器设计原理 (5)3.2 FIR低通滤波器的设计方法 (5)3.2.1 频率采样法 (5)3.2.2 最优化设计 (6)3.2.3 窗函数法 (6)3.3 窗函数法设计步骤 (8)第四章详细设计 (9)4.1 语音信号的采集 (9)4.2 语音信号的读入与打开 (10)4.3 语音信号的FFT变换 (11)4.4 含噪信号的合成 (12)4.5 利用FIR滤波器滤波 (13)4.6 结果分析 (14)总结 (15)参考文献 (16)附录 (17)致谢 (21)第一章引言随着信息科学和计算机技术的不断发展,数字信号处理(DSP,Digital Signal Processing)的理论和技术也得到了飞速的发展,并逐渐成为一门重要的学科,它的重要性在日常通信、图像处理、遥感、声纳、生物医学、地震、消费电子、国防军事、医疗方面等显得尤为突出。
在我们面临的信息革命中,数字信号处理几乎涉及了所有的工程技术领域。
数字信号处理是一种将信号以数字形式进行处理的一种理论和技术,它的目的是将真实世界中的一些信号进行分析并滤波,最后得出其中的有用的信号。
数字滤波器是数字信号处理的一种,一般根据单位脉冲响应h(n)分为无限脉冲响应(IIR)和有限脉冲响应(FIR)系统。
IIR数字滤波器的设计方法简单,特别是采用双线性变换法来设计的数字滤波器不存在频域混叠的现象,但是IIR滤波器存在一个较为明显的缺憾,就是它的相位响应一般都是非线性的,而在传输频带内的相位响应如果不是线性的,就会造成有用信号的传输失真,而FIR数字滤波器不仅可以设计成任意的幅度响应,而且可以设计成在通频带内具有良好的线性相位响应。
FIR数字滤波器的单位脉冲响应h(n)有限长,所以FIR数字滤波器是稳定的,不存在稳定性的问题,且可以通过快速傅里叶变换(FFT)的算法来实现信号滤波,大大的提高的运算效率。
[1]因此,FIR数字滤波器日益引起了人们的关注。
本课程设计是采用kaiser窗设计的FIR滤波器对语音信号进行滤波去噪。
通过课程设计了解FIR滤波器设计的原理和步骤,掌握用Matlab语言设计滤波器的方法,了解DSP对FIR滤波器的设计及编程方法。
通过观察语音信号滤波前后的时域波形的比较,加深对滤波器作用的理解。
通过对比滤波前后波形图的比较和放滤波前后语音信号的对比,可以看出滤波器对有用信号无失真放大具有重大意义。
1.1 设计目的及意义《信号处理》课程设计是现代信号处理技术课程的有效补充部分,通过课程设计,使得学生在设计过程中了解完整的现代信号处理技术的工程实现方法和流程,从而对现代信号处理技术的理论有更深入的认识。
本课程设计的目的是通过学生使用MATLAB等工具,采用窗函数法设计符合一定参数要求的FIR滤波器,并用所设计的滤波器对加噪语音信号进行滤波去噪处理。
1.2 设计任务及要求利用MATLAB语言及其工具箱来完成如下的工作:1、学会MATLAB的使用,掌握MATLAB的程序设计方法;2、录制一段语音信号,完成对信号的采样,画出信号的时域波形和频谱图,确定信号的频谱范围;3、给信号叠加噪声(噪声类型分为如下几种:a白噪声;b单频噪色(正弦干扰);c多频噪声(多正弦干扰);d其它干扰。
),画出受噪声干扰的信号时域波形和频谱图;4、采用窗函数法设计FIR低通滤波器,低通滤波器性能指标fb=1000 Hz,fc=1200 Hz,as=100 dB,ap=1 dB。
画出滤波器的频响特性图;5、用所设计的滤波器对受噪声影响的信号进行滤波,画出滤波后语音信号的时域波形图和频谱图;6、对滤波前后的信号进行对比,分析信号的变化;回放语音信号,并与原始语音信号对比;7、掌握窗函数法设计FIR数字滤波器的方法。
1.3 课程设计平台20世纪70年代后期,时任美国新墨西哥大学计算机科学系主任的Cleve Moler教授出于减轻学生编程负担的动机,为学生设计了一组调用LINPACK和EISPACK库程序的“通俗易用”的接口,此即用FORTRAN编写的萌芽状态的MATLAB。
MATLAB软件包括五大通用功能:数值计算功能(Nemeric);符号运算功能(Symbolic);数据可视化功能(Graphic);数据图形文字统一处理功能(Notebook)和建模仿真可视化功能(Simulink)。
其中,符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLAB命令窗口。
该软件有三大特点:一是功能强大;二是界面友善、语言自然;三是开放性强。
目前,Mathworks 公司已推出30多个应用工具箱。
MATLAB在线性代数、矩阵分析、数理统计和随机信号分析、电路与系统、系统动力学、信号和图像处理、建模和仿真、通信系统、以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。
第二章基本原理2.1 FIR滤波器的基本概念数字滤波器从实现的网络结构或者从单位脉冲响应可分为无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。
FIR滤波器是有限长单位冲激响应滤波器,是数字信号处理系统中最基本的元件,它可以在保证任意幅频特性的同时具有严格的线性相频特性,同时其单位抽样响应是有限长的,因而滤波器是稳定的系统。
[2]因此,FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。
从性能上说,IIR滤波器以非线性相位为代价以较低的阶数获得较高的选择性。
而FIR滤波器想要获得相同的选择性阶数是IIR滤波器的5-10倍,结果成本较高、信号时延也较大:从结构上说,IIR采用递归结构,FIR采用非递归结构;从设计工具上说;IIR可以借助于模拟滤波器的成果,FIR滤波器一般采用没有封闭形式的设计公式;从使用场合上来看,在对相位要求不敏感的场合,如语音通讯等,选用IIR较为合适,可以充分发挥经济高效的特点。
对图像处理、数据传输等以波形携带信息的系统,使用FIR较好。
2.2 FIR滤波器的特点有限长单位冲激响应(FIR)滤波器有以下特点:(1) 系统的单位冲激响应在有限个n值处不为零;(2) 系统函数H(z)在|z|>0处收敛,极点全部在z = 0处(因果系统);(3)结构上主要是非递归结构,没有输出到输入的反馈,但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。
设FIR滤波器的单位冲激响应为一个N点序列,,则滤波器的系统函数为=∑^))(式(2-1)*(H-nzzhn就是说,它有(N—1)阶极点在z = 0处,有(N—1)个零点位于有限z平面的任何位置。
优点:(1)很容易获得严格的线性相位,避免被处理的信号产生相位失真,这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要;(2)可得到多带幅频特性;(3)极点全部在原点(永远稳定),无稳定性问题;(4)任何一个非因果的有限长序列,总可以通过一定的延时,转变为因果序列,所以因果性总是满足;(5)无反馈运算,运算误差小。
缺点:(6)因为无极点,要获得好的过渡带特性,需以较高的阶数为代价;(7)无法利用模拟滤波器的设计结果,一般无解析设计公式,要借助计算机辅助设计程序完成。
[3]2.3 FIR滤波器的种类(1)数字集成电路FIR滤波器[4]一种是使用单片通用数字滤波器集成电路,这种电路使用简单,但是由于字长和阶数的规格较少,不易完全满足实际需要。
虽然可采用多片扩展来满足要求,但会增加体积和功耗,因而在实际应用中受到限制。
设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统,使其系统函数H(z)具有指定的频率特性。
(2)DSP芯片FIR滤波器另一种是使用DSP芯片。
DSP芯片有专用的数字信号处理函数可调用,实现FIR滤波器相对简单,但是由于程序顺序执行,速度受到限制。
而且,就是同一公司的不同系统的DSP芯片,其编程指令也会有所不同,开发周期较长。
(3)可编程FIR滤波器还有一种是使用可编程逻辑器件,FPGA/CPLD。
FPGA有着规整的内部逻辑块整列和丰富的连线资源,特别适合用于细粒度和高并行度结构的FIR滤波器的实现,相对于串行运算主导的通用DSP芯片来说,并行性和可扩展性都更好。
第三章 FIR 数字低通滤波器的设计3.1 FIR 低通滤波器设计原理由于FIR 滤波器冲击响应h[n]是有限长序列,因此这种结构可用非递归结构来实现[1]。
FIR 数字滤波器系统函数一般形式为如式(3-1)所示:∑=-=0)()(n n z k h z H 式(3-1)FIR 滤波器数学表达式可用差分方程(3-2)来表示:∑=-=0)()()(k k n x k h n y 式(3-2)式中:y (n )输出序列;h (k )滤波器系数;n 滤波器阶数;x (k )输入序列。
应用Matlab 设计FIR 滤波器的主要任务就是根据给定的性能指标,设计一个H (z ) ,使其逼近这一指标,进而计算并确定滤波器的系数b(n),再将所设计滤波器的幅频响应、相频响应曲线作为输出,与设计要求进行比较,对设计的滤波器进行优化。
3.2 FIR 低通滤波器的设计方法3.2.1 频率采样法频率采样法是从频域出发,根据频域采样定理,对给定的理想滤波器的频率响应()jw H e 加以等间隔的抽样 ,得到()d h k :N k w jw d d e H k H /)2(|)()(π==,k=0,1,…,N -1 式(3-3) 再利用()d H k 可求得FIR 滤波器的系统函数()H Z 及频率响应()jw H e 。
而在各采样点间的频率响应则是其加权内插函数延伸叠加的结果。
但对于一个无限长的序列,用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的形状,理想频响特性变换越平缓,则内插函数值越接近理想值,误差越小。