dsp课程设计matlab,dsp课程设计基于matlab的语音信号滤波处理.doc

dsp课程设计基于matlab的语音信号滤波处理.doc

数字信号处理课程设计报告基于MATLAB的语音信号滤波处理系(院)专业电子信息科学与技术班级学号姓名摘要滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用MATLAB作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，并利用MATLAB作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。关键词数字滤波器MATLAB窗函数法目录1绪论12课程设计内容23课程设计的具体实现231语音信号的采集232语音信号的时频分析233语音信号加噪与频谱分析334设计FIR数字滤波器435用滤波器对加噪语音信号进行滤波736比较滤波前后语音信号的波形及频谱737回放语音信号10200940620239数字信号处理课程设计DESIGNEDBY杨慧文01绪论数字信号处理是利用计算机或专用处理设备，以数值计算的方法对信号进行采集、抽样、变换、综合、估值与识别等加工处理，借以达到提取信息和便于应用的目的。它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等优点。数字滤波器,是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应FIR，FINITEIMPULSERESPONSE滤波器和无限冲激响应IIR，INFINITEIMPULSERESPONSE滤波器。FIR滤波器结构上主要是非递归结构，没有输出到输入的反馈，系统函数HZ在处收敛，极点全部在Z0处(因果系统)，因而只能用较高的阶数达到高的选择性。FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变，这是很好的性质。FIR数字滤波器是有限的单位响应也有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。FIR滤波器因具有系统稳定，易实现相位控制，允许设计多通带(或多阻带)滤波器等优点收到人们的青睐。滤波器的设计可以通过软件或设计专用的硬件两种方式来实现。随着MATLAB软件及信号处理工具箱的不断完善，MATLAB很快成为应用学科等领域不可或缺的基础软件。它可以快速有效地实现数字滤波器的设计、分析和仿真，极大地减轻了工作量,有利于滤波器设计的最优化。2课程设计内容语音信号滤波处理查阅资料，了解WAVREAD和SOUND函数的使用方法，利用WAVREAD函数对语音信号进行采集。将语音信号转换成计算机能够运算的有限长序列，加长一定的噪声，并通过选择适当的滤波器进行滤波，观察滤波后的效果,200940620239数字信号处理课程设计DESIGNEDBY杨慧文1并通过SOUND函数回放。(学号141,142,143,144,145,,239,240,241,242,243,343,344,345)。3课程设计的具体实现31语音信号的采集利用PC机上的声卡和WINDOWS操作系统可以进行数字信号的采集。将话筒输入计算机的语音输入插口上,启动录音机。按下录音按钮,接着对话筒说话“语音信号处理”,说完后停止录音,屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所录音的重现。以文件名“DSP”保存入D\MATLAB\WORK中。可以看到,文件存储器的后缀默认为WAV,这是WINDOWS操作系统规定的声音文件存的标准。32语音信号的时频分析利用MATLAB中的“WAVREAD”命令来读入(采集)语音信号，将它赋值给某一向量。再对其进行采样，记住采样频率和采样点数。接下来，对语音信号DSPWAV进行采样。其程序如下Y,FS,NBITSWAVERED DSPWAV 把语音信号加载入MATLAB仿真软件平台中然后，画出语音信号的时域波形，再对语音信号进行频谱分析。在本次设计中，我们利用FFT对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。其程序如下Y,FS,NBITSWAVREAD DSPWAV SOUNDY,FS,NBITS回放语音信号NLENGTHY求出语音信号的长度YFFTY,N傅里叶变换SUBPLOT2,1,1PLOTYTITLE 原始信号波形 SUBPLOT2,1,2PLOTABSYTITLE 原始信号频谱 程序结果如下图200940620239数字信号处理课程设计DESIGNEDBY杨慧文2图1原始信号波形及频谱33语音信号加噪与频谱分析在本次课程设计中，我们是利用MATLAB中的随机函数RAND或RANDN产生噪声加入到语音信号中，模仿语音信号被污染，并对其频谱分析。在这里，我们选用RANDN(M,N)函数。语音信号添加噪声及其频谱分析的主要程序如下Y,FS,NBITSWAVREAD DSPWAV NLENGTHY求出语音信号的长度NOISE001RANDNN,2随机函数产生噪声SIYNOISE语音信号加入噪声SOUNDSISUBPLOT2,1,1PLOTSITITLE 加噪语音信号的时域波形 SFFTSI傅里叶变换SUBPLOT2,1,2PLOTABSSTITLE 加噪语音信号的频域波形 程序结果如下图200940620239数字信号处理课程设计DESIGNEDBY杨慧文3图2加噪后的波形及频谱分析34设计FIR数字滤波器用窗函数法设计FIR滤波器的步骤。如下(1)根据对阻带衰减及过渡带的指标要求，选择串窗数类型(矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等)，并估计窗口长度N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。(2)构造希望逼近的频率响应函数。(3)计算HN。(4)加窗得到设计结果。接下来，我们根据语音信号的特点给出有关滤波器的技术指标①低通滤波器的性能指标FP1000HZ，FC1200HZ，AS100DB,AP1DB②高通滤波器的性能指标FP3500HZ，FC4000HZ，AS100DB，AP1DB；③带通滤波器的性能指标FP11200HZ，FP23000HZ，FC11000HZ，FC23200HZ，AS100DB，AP1DB在MATLAB中，可以利用函数FIR1设计FIR滤波器，利用MATLAB中的函数FREQZ画出各步步器的频率响应。下面我们将给出FIR数字滤波器的主要程序。FIR低通滤波器用窗函数设计低通滤波器的程序如下200940620239数字信号处理课程设计DESIGNEDBY杨慧文4