《语音信号处理 matlab》由会员分享，可在线阅读，更多相关《语音信号处理 matlab（49页珍藏版）》请在金锄头文库上搜索。
1、目录摘要-3 关键词-3 1课程设计的目的、意义-3 1.1 目的-31.2 意义-32 设计任务及技术指标-3 2.1 设计任务-32.2技术指标-43 设计方案论证-43.1 信号的抽样与恢复-4 3.2 滤波器设计-4 3.2.1模拟滤波器设计原理-4 3.2.2 模拟滤波器数字化原理-53.2.3数字高通、带通、带阻滤波器的设计-54 设计内容 -6 4.1必选部分-64.2 拓展部分-65 实验结果与分析-66 总结-29 参考文献-30题目语音信号分析系统设计摘要:语音信号分析系统设计涵盖了基于matlab对语音信号进行处理的原理、方法及源程序,包括语音信号的采集与时域频域分析、对
2、语音信号进行加噪滤波处理、噪声滤除前后的语音进行时频域分析、语音信号的重采样和回放、语音信号部分时域参数的提取及GUI用户友好界面的实现,实现了对简单的单音节语音的综合处理与分析。本设计程序界面简练,操作简便,有一定的实际意义。最后对语音信号的处理发展方向提出了自己的看法。关键词:语音,matlab,GUI,滤波,噪声,采样1课程设计的目的、意义1.1 目的设计一个简单的语音信号分析系统,实现对语音信号时域波形显示、进行频谱分析,利用滤波器滤除噪声、对语音信号的参数进行提取分析等功能。采用Matlab设计语言信号分析相关程序,并且利用GUI设计图形用户界面。(1采集语音信号。(2对原始语音信号
3、加入干扰噪声,对原始语音信号及带噪语音信号进行时频域分析。(3针对语音信号频谱及噪声频率,设计合适的数字滤波器滤除噪声。(4对噪声滤除前后的语音进行时频域分析。(5对语音信号进行重采样,回放并与原始信号进行比较。(6对语音信号部分时域参数进行提取。(7设计图形用户界面(实现以上功能。1.2 意义本次设计基本完成了语音信号的分析与处理,在下面方面有实际意义:1.从语音信号中得到某些语音特征参数,如:短时平均幅度,短时平均能量,短时平均过零率,短时自相关函数。2.对原信号,带噪信号进行简单的时域频域分析。3.分析不同采样频率下的语音信号的不同。2 设计任务及技术指标2.1设计任务设计一个简单的语音
4、信号分析系统,实现对语音信号时域波形显示、进行频谱分析,利用滤波器滤除噪声、对语音信号的参数进行提取分析等功能。采用Matlab设计语言信号分析相关程序,并且利用GUI设计图形用户界面。具体任务是:(1采集语音信号。(2对原始语音信号加入干扰噪声,对原始语音信号及带噪语音信号进行时频域分析。(3针对语音信号频谱及噪声频率,设计合适的数字滤波器滤除噪声。(4对噪声滤除前后的语音进行时频域分析。(5对语音信号进行重采样,回放并与原始信号进行比较。(6对语音信号部分时域参数进行提取。(7设计图形用户界面(实现以上功能。2.2技术指标低通滤波器的技术指标WP,WS,RP,AS,采样频率的技术指标FS,
5、帧长,帧移,帧重叠采样点数3 设计方案论证3.1 信号的抽样与恢复利用抽样脉冲把一个连续信号变为离散时间样值的过程称为抽样,在满足抽样定理条件下,抽样信号保留了原信号的全部信息,并且从抽样信号中可以无失真的恢复出原始信号。抽样定理在通信系统、信息传输理论方面占有十分重要的地位。数字通信系统是以此定理作为理论基础。抽样过程关系到通信设备整个系统的性能指标。3.2 滤波器设计3.2.1模拟滤波器设计原理(1模拟巴特沃思滤波器原理巴特沃斯滤波器具有单调下降的幅频特性:在小于截止频率的范围内,具有最平幅度c的响应,而在c 后,幅频响应迅速下降。 巴特沃思低通滤波器幅度平方函数为:221(1(a NcH
6、 j =+ (3-1式中N 为滤波器阶数,c 为3dB 截止角频率。将幅度平方函数写成s 的函数:21(1(a a NcH s H s s j -=+ (3-2该幅度平方函数有2N 个等间隔分布在半径为c 的圆上的极点121(22 k j Nk c s e +=,0,1,.21k N =- 为了形成稳定的滤波器,取左半平面的N 个极点构成(a H s ,即:1(N N a ckk H s s s -=- (3-3为使设计统一,将频率归一化,得到归一化极点121(22k j Nk p e+=,相应的归一化系统函数为:1(1(N a k k H p p p -=- (3-4多项式形式为: 01(1
7、(.N a H p b b p p =+ (3-5(2模拟切比雪夫滤波器原理切比雪夫滤波器的幅频特性具有等波纹特性,有两种形式,在通带内等波纹、阻带单调的是I 型滤波器,在通带内单调、在阻带内等波纹的是II 滤波器。以I 型滤波器为例。 切比雪夫滤波器的幅度平方函数为:22221(1(aN pA H j C =+ (3-6为小于1的正数,表示通带内幅度波动的程度。p 称为通带截止频率。令=/p ,称为对p 的归一化频率。C N (x为N 阶切比雪夫多项式。幅度平方函数的极点是分布在b p 为长半轴,a p 为短半轴的椭圆上的点。同样取s 平面左半平面的极点构成(a H s :11(2(NN N
8、 a pi i H s s s -=- (3-7进行归一化,得到:11(12(NN a i i H p p p -=- (3-8其中(21(21sincos22k k p ch jch i NN-=-+ ,11(Arsh N=3.2.2 模拟滤波器数字化原理将模拟滤波器转化为数字滤波器在工程上常用的有脉冲响应不变法和双线性变换法。脉冲响应不变法时一种时域上的转换方法,它是数字滤波器的单位取样响应在抽样点上等于模拟滤波器的单位冲激响应,即:(h n h nT a = (3-9设模拟滤波器只有单阶极点,其系统函数为:1(Ni a i iA H s s s =-(3-10对(a H s 进行拉氏反变
9、换得到(a h t ,对(a h t 进行等间隔抽样,得到(h n h nT a =,对(h n进行Z 变换,得到数字滤波器系统函数:11(1i Nis Ti A H z ez-=- (3-11这种方法s 和z 的关系是:sT z e =。该方法的优点是频率坐标变换时线性的切数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应,时域特性逼近好;缺点是会产生频谱混叠现象,适合低通、带通滤波器的设计,不适合高通、带阻滤波器的设计。双线性变换法为了克服频谱混叠现象,采用非线性频率压缩方法,将整个频率轴上的频率范围压缩到/T 之间,再用sT z e =转换到Z 平面上。这种方法s 和z 的关系是:
10、11(2/(1/1s T z z -=-+。该方法克服了频谱混叠现象,但带来了频率坐标变换的非线性:(2/tan(/2T =,由模拟滤波器系统函数转换为数字滤波器系统函数公式为:11211(|a z z T z H z H s -=+= (3-123.2.3数字高通、带通、带阻滤波器的设计这些滤波器可以借助于模拟滤波器的频率变换设计一个所需类型的模拟滤波器, 再通过双线性变换法将其转换成所需类型的数字滤波器。首先确定所需类型数字滤波器的技术指标;然后将数字滤波器技术指标按照公式(2/tan(/2T =转换成所需类型滤波器的模拟域技术指标;将所需类型滤波器的模拟域技术指标转换成低通滤波器技术指标
11、;设计归一化模拟低通滤波器;去归一化得到模拟低通滤波器的系统函数;将模拟低通滤波器转换为所需类型的模拟滤波器;最后通过双线性变换法转换成所需类型的数字滤波器。4 设计内容 4.1必选部分利用MATLAB 软件对语音信号进行频谱分析;并对语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声。(1根据设计要求分析系统功能,掌握设计中所需理论、阐明设计原理(抽样频率、量化位数的概念,抽样定理;信号的FFT 分析;数字滤波器设计原理和方法,各种不同类型滤波器的性能比较。(2对语音信号以不同采样频率进行重采样,回放并与原始信号进行比较(如原语音信号采样频率为16KHz ,对其进行
12、8KHz ,4KHz 的重采样。(3对语音信号做FFT ,进行频谱分析,画出信号的时域波形图和频谱图。(4对语音信号加入干扰噪声,对语音信号进行回放(利用sound 函数,感觉加噪前后声音的变化。对其做FFT ,进行频谱分析,比较加噪前后语音信号的波形及频谱,对所得结果进行分析。(5根据带噪语音信号的特点,设计合适的数字滤波器,绘制所设计滤波器的幅频和相频特性。(6用所设计的滤波器对带噪语音信号进行滤波。对滤波后的语音信号进行FFT 频谱分析。记录处理过程中所得各种波形及频谱图。(7对语音信号进行回放,感觉滤波前后声音的变化。比较滤波前后语音信号的波形及频谱,对所得结果和滤波器性能进行频谱分析。4.2 拓展部
