.PAGE/NUMPAGES课程
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
报告书数字信号处理课程设计报告书基于Matlab对语音信号进行频谱分析与滤波学院：专业：设计者：学号：设计周数：完成日期：目录TOC\o"1-3"\h\z\uHYPERLINK\l"_Toc142881860"目录2HYPERLINK\l"_Toc142881861"一、课程设计目的2HYPERLINK\l"_Toc142881862"二、课程设计
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
目3HYPERLINK\l"_Toc142881863"三、课程设计题目描述和要求3HYPERLINK\l"_Toc142881865"五、课程设计
流程
快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计
5HYPERLINK\l"_Toc142881865"六、基于MATLAB的滤波器设计结果分析18HYPERLINK\l"_Toc142881867"八、设计经验
总结
初级经济法重点总结下载党员个人总结TXt高中句型全总结.doc高中句型全总结.doc理论力学知识点总结pdf
PAGEREF_Toc142881867\h4HYPERLINK\l"_Toc142881867"七、参考资料与网址PAGEREF_Toc142881867\h4HYPERLINK\l"_Toc142881868"八、附录PAGEREF_Toc142881868\h4[作者]钟伟雄[机构]XX工业大学在校本科生摘要：数字信号处理是通信工程专业的一门相当重要的学科,对日后就业和科研有重大意义,通过MATLAB,我们可以清晰地理解数字信号处理中难以理解的一面,对理论的知识加以深化。关键字：MATLAB数字信号处理GUI频谱相位滤波器一、课程设计题目应用Matlab对语音信号进行频谱分析与滤波二、课程设计目的数字信号处理是一门以算法为核心,理论和实践性较强的学科。是电子信息工程、通信工程专业、电子信息科学与技术专业的一门重要的专业技术基础课。数字信号处理课程是在学习完数字信号处理的相关理论后,进行的综合性训练课程,其目的是：使学生进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法；增强学生应用Matlab语言编写数字信号处理的应用程序与分析、解决实际问题的能力；三、课程设计内容描述和要求为了巩固所学的数字信号处理理论知识,使学生对信号的采集、处理、传输、显示和存储等有一个系统的掌握和理解,安排了以下的课程设计的内容：录制一段个人自己的语音信号,并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；给定滤波器的性能指标,采用窗函数法和双线性变换设计滤波器,并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波,画出滤波后信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化；回放语音信号；最后,设计一个信号处理系统界面。下面对各步骤加以具体说明。2．1语音信号的采集利用Windows下的录音机,录制一段自己的话音,时间在1s内。然后基于Matlab软件平台下,利用函数[x,fs,bits]=wavread
对语音信号进行采样,记住采样频率和采样点数。通过wavread函数的使用,理解采样频率、采样位数等概念。2．2语音信号的频谱分析首先播放读入的语音信号,画出语音信号的时域波形；然后调用FFT函数对语音号进行快速傅里叶变换,得到信号的频谱特性,调用plot,从而加深对频谱特性的理解。2．3设计数字滤波器和画出其频率响应给出各滤波器的性能指标：<1>低通滤波器性能指标fb＝1000Hz,fc＝1200Hz,As＝100dB,Ap＝1dB。<2>高通滤波器性能指标fc＝4800Hz,fb＝5000HzAs＝100dB,Ap＝1dB。<3>带通滤波器性能指标fb1＝1200Hz,fb2＝3000Hz,fc1＝1000Hz,fc2＝3200Hz,As＝100dB,Ap＝1dB。要求用窗函数法和双线性变换法设计上面要求的3种滤波器。在Matlab中,可以利用函数fir1设计FIR滤波器,可以利用函数butte,cheby1和ellip设计IIR滤波器；利用Matlab中的函数freqz画出各滤波器的频率响应。2．4用滤波器对信号进行滤波要求用自己设计的各滤波器分别对采集的信号进行滤波,在Matlab中,FIR滤波器利用函数fftfilt对信号进行滤波,IIR滤波器利用函数filter对信号进行滤波。2．5比较滤波前后语音信号的波形与频谱要求在一个窗口同时画出滤波前后的波形与频谱。2．6回放语音信号在Matlab中,调用函数sound对声音进行回放。其调用格式：sound；可以感觉滤波前后的声音有变化。2．7设计系统界面为了使编制的程序操作方便,设计处理系统的用户界面。在所设计的系统界面上可以选择滤波器的类型,滤波方式,输入滤波器的参数,显示滤波器的频率响应,选择音频信号等。这里选用matlab的GUI工具进行界面设计。四、课程设计进度安排序号设计内容所用时间1熟悉Matlab程序设计方法,了解数字信号处理工具箱使用1天2分析题目,设计程序框图,编写程序代码1天3上机调试程序,修改并完善设计,并完成设计报告1天合计3天五、课程设计基本流程5．1熟悉MATLAB基本操作,基本语法,以与GUI界面设计中各控件的调用方式。MATLAB是一款功能超强的数学软件,应用于各个行业。而其基本操作却很大众化,操作起来很人性化。其中所附带的GUI面向对象的友好编程方式更便于学术交流和设计开发。5．2熟悉各种滤波器设计
方案
气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载
以与和收集与滤波器相关的MATLAB函数与其调用格式。一．IIR数字滤波器〔一〕IIR数字滤波器的传递函数与特点设IIR滤波器的输入序列为x,则IIR滤波器的输入序列x与输出序列y之间的关系可以用下面的方程式表示：其中,和是滤波器的系数,其中中至少有一个非零。与之相对应的差分方程为：由传递函数可以发现无限常单位冲激响应滤波器有如下特点：a.单位冲激响应h是无限长的。b.系统传递函数H在有限z平面上有极点存在。c.结构上存在着输出到输入的反馈,也就是结构上是递归型的。<二>IIR数字滤波器的设计与实现IIR数字滤波器的设计有多种方法,如频率变换法、数字域直接设计以与计算辅助设计等。下面只介绍频率变换设计法。首先考虑由模拟低通滤波器到数字低通滤波器的转换,其基本的设计过程如下:A.将数字滤波器的技术指标转换为模拟滤波器的技术指标；B.设计模拟滤波器G；C.将G转换成数字滤波器H；在低通滤波器的设计基础上,可以得到数字高通、带通、带阻滤波器的设计流程如下：A.给定数字滤波器的设计要求〔高通、带阻、带通〕；B.转换为模拟〔高通、带阻、带通〕滤波器的技术指标；C.转换为模拟低通滤波器的指标；D.设计得到满足第三步要求的低通滤波器传递函数；E.通过频率转换得到模拟〔高通、带阻、带通〕滤波器；F.变换为数字〔高通、带阻、带通〕滤波器。<三>在matlab中设计IIR滤波器的方法与其它们所用到的函数如表所示。表1.基于matlab中设计IIR滤波器的方法列表方法描述函数模拟原型法采用经典低通滤波器作为连续域上的设计模型,通过频率变换得到IIR数字滤波器,最后进行离散化处理完整设计函数：Beself,butter,cheby1,cheby2,ellip滤波器的阶估计函数：Buttord,cheb1ord,cheb2ord,ellipord低通模拟滤波器原型函数：beselap,buttap,cheb1ap,cheb2ap,ellipap频域变换函数：Lp2bp,lp2bs,lp2hp,lp2lp其他函数：Bilinear,impinvar直接设计方法直接在离散时域上估计线性的幅度响应yulewalk通用butterworth设计方法使用butterworth设计低通数字滤波器Maxflat参数建模方法寻找接近于所需要设计的滤波器的通用原型时域上的建模函数：Lpc,prony,stmcb频域上的建模函数：Invfreqs,invfreqz表2频率转换函数列表频率转换转换函数低通到低通[numt,dent]=lp2lp[At,Bt,Ct,Dt]=lp2lp低通到高通[numt,dent]=lp2hp[At,Bt,Ct,Dt]=lp2hp低通到带通[numt,dent]=lp2bp[At,Bt,Ct,Dt]=lp2bp低通到带阻[numt,dent]=lp2bs[At,Bt,Ct,Dt]=lp2bs〔四〕标准数字滤波器设计函数Matlab提供了一组标准的数字滤波器设计函数,大大简化了滤波器的设计过程。butter功能：Butterworth模拟/数字滤波器设计格式：[b,a]=butter[b,a]=butter说明：选项中加入‘S’用于设计各种模拟Butterworth滤波器；不加设计各种数字Butterworth滤波器Ftype为缺省,设计低通滤波器Ftype＝hign,设计高通滤波器Ftype＝stop,设计带阻滤波器[实例1]设计一个5阶Butterworth数字高通滤波器,阻带截止频率为250Hz。设采样频率为1000Hz。源代码如下：[b,a]=butter<5,250/500,'high'>[z,p,k]=butter<5,250/500,'high'>freqz程序运行后,产生结果如下所示。b=0.0528-0.26390.5279-0.52790.2639-0.0528a=1.0000-0.00000.6334-0.00000.0557-0.0000z=11111p=0.0000+0.7265i0.0000-0.7265i0.0000+0.3249i0.0000-0.3249i0.0000k=0.0528图15阶Butterworth数字高通滤波器Cheby1、Cheby2功能：chebyshevI、chebyshevII型模拟/数字滤波器设计格式：[b,a]=cheby1[b,a]=cheby2[实例2]设计一个7阶chebyshevII型数字低通滤波器,截止频率为3000Hz,Rs＝30dB。设采样频率为1000Hz。解：源程序如下：[b,a]=cheby2<7,30,300/500'>；[z,p,k]=butter<5,250/500,'high'>；freqz程序运行后,产生如图6－2所示的波形。图27阶chebyshevII型数字低通滤波器〔五〕冲激响应不变法一般来说,在要求时域冲激响应能模仿模拟滤波器的场合,一般使用冲激响应不变法。冲激响应不变法一个重要特点是频率坐标的变换是线性的,因此如果模拟滤波器的频响带限于折叠频率的话,则通过变换后滤波器的频响应可不失真的反映原响应与频率的关系。[实例3]设计一个中心频率为500HZ,带宽为600Hz的数字带通滤波器,采样频率为1000Hz。解：源代码如下：[z,p,k]=buttap<3>;[b,a]=zp2tf;[bt,at]=lp2bp;[bz,az]=impinvar;%将模拟滤波器变换成数字滤波器freqz程序运行后,产生如图6－3所示的波形。图3用冲激响应不变法设计数字低通滤波器〔六〕双线性变换法与冲激响应不变法比较,双线性变换的主要优点是靠频率的非线性关系得到S平面与Z平面的单值一一对应关系,整个值对应于单位圆一周。所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。[实例4]设计一个截止频率为200Hz的数字低通滤波器,采用频率为1000Hz。解：源代码如下：[z,p,k]=buttap<3>;[b,a]=zp2tf;[bt,at]=lp2lp;[bz,az]=bilinear;freqz结果如图6－4所示：bz=0.07530.22590.22590.0753az=1.0000-0.82660.5154-0.0865图4用双线性变换法设计数字低通滤波器[实例5]基于Butterworth模拟滤波器原型,使用双线性转换设计数字滤波器,其中参数指标为：通带截止频率：通带波动值：阻带截止频率：阻带波动值：解：首先确定滤波器的阶数N,同时根据确定＝0.5。接着使用bilinear进行双线性转换,最后绘制在频域上的各种图像,其源代码如下：％数字滤波器指标wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;％将数字滤波器指标反转变化为模拟滤波器的参数T=1;fs=1/T;omegap=<2/T>*tan;omegas=<2/T>*tan;ep=sqrt<10^
浙公网安备 33021202002483