基于matlab德语音信号处理

基于Matlab的数字滤波器设计及其对语音信号的应用摘要MATLAB 语言是一种数据分析和处理功能十分强大的计算机应用软件 ,它可以将声音文件变换为离散的数据文件 , 然后利用其强大的矩阵运算能力处理数据 ,如数字滤波、傅里叶变换、时域和频域分析、声音回放以及各种图的呈现等 , 信号处理是MATLAB 重要应用的领域之一。

本课程设计介绍了基于Matlab的对语音信号采集、处理及滤波器的设计，并使之实现的过程。

理解与掌握课程中的基本概念、基本原理、基本分析方法，用Matlab进行数字语音信号处理，并阐述了课程设计的具体方法、步骤和内容。

综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。

关键词： MATLAB 工具信号语音采集滤波器一课程设计的目的：1．掌握数字信号的频谱分析及数字滤波器的基本概念；2．掌握基于MATLAB 的FIR和IIR数字滤波器的设计方法；3．掌握在MATLAB下语音信号的采集、时域波形分析、频谱分析及滤波设计；4．培养学生分析、解决问题的能力；5. 初步培养学生撰写论文的能力，为下一步毕业设计做准备。

二滤波器设计方法综述：1)．直接设计数字滤波器；数字滤波器是数字信号处理的一个重要组成部分。

数字滤波实际上是一种运算过程，其功能是将一组输入的数字序列通过一定的运算后转变为另一组输出的数字序列。

数字滤波器一般可以用两种方法来实现：一种是根据描述数字滤波器的数学模型或信号流图，用数字硬件装配成一台专门的设备，构成专用的信号处理机；另一种方法就是直接利用通用计算机，将所需要的运算编成程序让计算机来执行，这也是用软件来实现数字滤波器。

数字滤波器是离散时间系统，所处理的信号是离散时间信号。

一般时域离散系统或网络可以用差分方程、单位脉冲响应以及系统函数进行描述。

如果输入、输出服从N阶差分方程接近于线性关系，Ω进一步增加时，ω增长变得缓慢， (ω终止于折叠频率处)，所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。

双线性变换法的缺点：Ω与ω的非线性关系，导致数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。

双线性变换比脉冲响应法的设计计算更直接和简单。

由于s与z之间的简单代数关系，所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。

4)．使用MATLAB下的Filter Designed Tool3。

在Command Window窗口中输入fdatool，出现滤波器设计工具界面图1 滤波器设计界面设计思路：首先将数字高通滤波器的技术指标转化为模拟滤波器的技术指标，再采用Buttetworth逼近的方法求得其系统函数H(S)，再利用冲激响应不变法将模拟滤波器的系统函数H(s)变为相应的数字滤波器的系统函数H(Z)，然后用MATALB进行仿真，然后确定数字滤波器的结构。

5)．IIR滤波器设计的基本思想将数字滤波器的设计为模拟滤波器的设计。

设计满足技术指标的模拟滤波器。

将模拟滤波器转换为数字滤波器。

无限长单位脉冲响应(IIR)数字滤波器的单位脉冲响应h(n)是一个无限长序列，器滤波器的结构采用递归结构，即在信号流图中包含反馈支路。

基本网络结构有直接型、级联型和并联型三种。

其中，直接型有可以分为直接Ⅰ型和直接Ⅱ型两种。

直接Ⅰ型()()()M Niii i y n b x n i a y n i ===-+-∑∑直接Ⅱ型121()()()1Miii Nii i b zH z H z H z a z-=-===-∑∑级联型11111(1)()1(1)MMi iii i N Nii ii i c zb zH z Aa z dz --==--==-==--∑∏∑∏并联型012()()()()N H z A H z H z H z =++++0111Ni i iA A dz -==+-∑6)．设计FIR 滤波器用窗函数法设计上面要求的3种滤波器。

可以利用函数fir1设计FIR 滤波器。

滤波器的单位脉冲响应h(n)长度为N ，那么1()()N nn H z h n z--==∑H(z)是1z -的N-1次多项式(FIR 滤波器的阶数是N-1)，该系统在z 平面上有N-1个零点，在z=0处有一个N-1重极点。

系统永远稳定。

与IIR 滤波器比较，设计方法差别很大，由于h(n)长度为有限的，因此FIR 滤波器设计通常是找出这有限个值。

直接型10()()N nn H z h n z--==∑级联型11201201()()()MN nii i n i H z h n zaa za z----====++∑∏频率采样型'1()()()c k k H z H z H z N==∑三课程设计的总体内容：给定滤波器的性能指标，在MATLAB平台下设计各种FIR和IIR数字滤波器，画出滤波器的频率响应；录制一段语音信号，画出语音信号的时域波形和频谱图，然后给原始语音信号加入噪声，然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱,并对滤波前后的信号频谱图进行对比，分析信号的变化，回放语音信号。

四课程设计的准备工作：1. 准备相关参考书，查阅相关文献，复习数字信号处理中各种滤波器的性能指标及设计方法；2. 准备基于MATLAB的语音信号的录音、读取、时域波形、频谱分析等基本知识。

3. 根据你对滤波器原理及其应用的理解，制定自己的设计方案，使用Matlab进行设计，写出源代码，得出仿真结果。

五具体步骤与要求：1.设计数字滤波器，画出其频率响应曲线各滤波器的性能指标：(1)低通滤波器性能指标fp＝1000Hz, fs＝1200Hz, As＝40dB, Ap＝1dB。

(2)高通滤波器性能指标fs＝4000Hz, fp＝4300Hz, As＝40dB, Ap＝1dB(3)带通滤波器性能指标fp1＝1200 Hz, fp2＝3000 Hz, fs1＝1000 Hz, fs2＝3200 Hz,As＝40dB, Ap＝1dB。

要求：(1)设计IIR滤波器：利用函数butter设计低通、利用函数cheby1设计高通，利用函数ellip设计带通。

(2)(选做)设计IIR滤波器：用双线性变换法设计上面要求的3种滤波器。

(3)(选做)设计FIR滤波器：用窗函数法设计上面要求的3种滤波器。

可以利用函数fir1设计FIR滤波器。

(4)函数freqz画出各滤波器的频率响应。

①低通滤波器图1 cheby1滤波器设计流程图程序为：[y,fs,bit]=wavread('xm.wav');fp=1000;fr=1200;ap=1;as=40;[n,fn]= buttord(fp/(fs/2),fr/(fs/2),ap,as,'z');[b,a]= buttap(n,ap,fn);freqz(b,a);[h,w]=freqz(b,a);subplot(221);plot(w*fs/(2*pi),abs(h));title(' buttap低通滤波');y1=filter(b,a,y);subplot(222);plot(y);title('原始语音信号');subplot(223);plot(y1);title('IIR滤波后语音信号');X=fft(y1,4096);subplot(224);plot(abs(X));title('滤波后信号频谱');显示为：020004000600000.511.5butter 低通滤波02000400060008000-1-0.500.51原始语音信号2000400060008000-0.50.51IIR 滤波后语音信号02000400060008000-4-2024-8FIR 滤波后语音信号图1 butter 高通滤波(a) 原始语音信号(b) IIR 滤波后语音信号(c) 滤波后信号频谱(d)②高通滤波器图2 cheby1滤波器设计流程图程序为：[y,fs,bit]=wavread('xm.wav');fp=4300;fr=4000;ap=1;as=40;[n,fn]=cheb1ord(fp/(fs/2),fr/(fs/2),ap,as,'z');[b,a]=cheby1(n,ap,fn,'high');freqz(b,a);[h,w]=freqz(b,a);subplot(221);plot(w*fs/(2*pi),abs(h));title('cheby1高通滤波');y1=filter(b,a,y);subplot(222);plot(y);title('原始语音信号');subplot(223);plot(y1);title('IIR滤波后语音信号');X=fft(y1,4096);subplot(224);plot(abs(X));title('滤波后信号频谱');显示为：图2 cheby1高通滤波(a) 原始语音信号(b)IIR滤波后语音信号(c) 滤波后信号频谱(d) ③带通滤波器程序为：[y,fs,bit]=wavread('xm.wav');fp1=1200;fr1=1000;fp2=3000;fr2=3200;ap=1;as=40;fp=[fp1,fp2];fr=[fr1,fr2];[n,fn]=cheb1ord(fp/(fs/2),fr/(fs/2),ap,as,'z');[b,a]=cheby1(n,ap,fn);[h,w]=freqz(b,a);subplot(221);plot(w*fs/(2*pi),abs(h));title('cheby1带通滤波');y1=filter(b,a,y);subplot(222);plot(y);title('原始语音信号');subplot(223);plot(y1);title('IIR滤波后语音信号');X=fft(y1,4096);subplot(224);plot(abs(X));title('滤波后信号频谱');显示为：图3 cheby1带通滤波(a) 原始语音信号(b)IIR滤波后语音信号(c) 滤波后信号频谱(d)2.语音信号的采集录制一段自己的话音,或利用老师给的语音，在Matlab软件平台下,利用函数wavread对语音信号进行采样,记住采样频率和采样点数。