课程设计说明书语音信号的数字滤波-应用MATLAB仿真.docVIP

课程设计说明书 NO.1 语音信号的数字滤波 ------应用MATLAB仿真 一、课程设计的目的 课程设计的目的，是使学生融会贯通本课程所学专业理论知识，完成一个较完整的设计过程，以加深学生对所学理论的理解与应用，培养学生综合运用基础理论知识和专业知识去解决实际工程设计问题的能力。课程设计是工科学不可少的一个综合性实践环节。 课程设计是课程学习后的一个综合性实践环节，是对课程理论和课程实验的综合和补充，是学完一门课后应用本课知识及以前的知识积累而进行的综合性、开放性的训练，是培养学生工程意识和创新能力的重要环节，培养学生综合运用已学过的理论和技能去分析和解决问题的能力，对加深课程理论的理解和应用具有重要意义。 二、课程设计方案与论证 线性相位FIR滤波器通常采用窗函数法设计。窗函数法设计FIR滤波器的基本思想是：根据给定的滤波器技术指标，选择滤波器长度N和窗函数ω（n），使其具有最窄宽度的主瓣和最小的旁瓣。其核心是从给定的频率特性，通过加窗确定有限长单位脉冲响应序列h(n)。工程中常用的窗函数共有6种，即矩形窗、巴特利特（Bartlett）窗、汉宁（Hanning）窗、汉明（Hamming）窗、布莱克曼（Blackman）窗和凯塞（Kaiser）窗。设计一个线性相位带通FIR滤波器指标如下： fn=[1000,1375,3625,4000]; a=[0,1,0]; dev=[0.0005,0.05,0.0005];最小衰减来选择凯塞（Kaiser）窗，利用MATLAB的求阶函数和FIR滤波器的设计函数，可以快速地设计出所需的数字滤波器。这两个设计函数如下：[N,Wn,beta,ftype]=kaiserord(fn,a,dev,fs) b=fir1(n,Wn,ftype,Kaiser(n+1,beta)； 最后，利用的滤波器分析函数freqz分析所设计出的滤波器的幅频特性和相频特性，并用图形显示函数plot将它们显示出来，如图1所示。由图1可见，设计结果满足指标要求。p0几何对称，因此ws1=wp0- (ws2-wp0)=0.3π 通带截止频率wc1=0.4π,wc2=0.6π;阻带截止频率wr1=0.3π，wr2=0.7π；阻带最小衰减αs=3dB和通带最大衰减αp=15dB; (2)用Ω=2/T*tan(w/2)对带通数字滤波器H(z)的数字边界频率预畸变，得到带通模拟滤波器H(s)的边界频率主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1，ωs2的转换。 为了计算简便，对双线性变换法一般T=2s 沈 阳 大 学 课程设计说明书 NO.3 通带截止频率wc1=(2/T)*tan(wp1/2)=tan(0.4π/2)=0.7265 wc2=(2/T)*tan(wp2/2)=tan(0.6π/2)=1.3764 阻带截止频率wr1=(2/T)*tan(ws1/2)=tan(0.3π/2)=0.5095 wr2=(2/T)*tan(ws2/2)=tan(0.7π/2)=1.9626 阻带最小衰减αs=3dB和通带最大衰减αp=15dB; (3)运用低通到带通频率变换公式λ=(((Ω^2)-(Ω0^2))/(B*Ω))将模拟带通滤波器指标转换为模拟低通滤波器指标。 B=wc2-wc1=0.6499 normwr1=(((wr1^2)-(w0^2))/(B*wr1))=2.236 normwr2=(((wr2^2)-(w0^2))/(B*wr2))=2.236 normwc1=(((wc1^2)-(w0^2))/(B*wc1))=1 normwc2=(((wc2^2)-(w0^2))/(B*wc2))=1 得出，normwc=1，normwr=2.236 模拟低通滤波器指标：normwc=1，normwr=2.236，αp=3dB，αs=15dB (4)设计模拟低通原型滤波器。用模拟低通滤波器设计方法得到模拟低通滤波器的传输函数Ha(s);借助巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等。 (5)调用lp2bp函数将模拟低通滤波器转化为模拟带通滤波器。 (6)利用窗函数变换法将模拟带通滤波器Ha(s)转换成数字带通滤波器H(z). 语音信号的频谱分析： 语音文件“ZJH.wav”频谱分析程序： z=wavread(Z