(完整版)语音信号进行倒谱分析 .doc

(完整版)语音信号进行倒谱分析 .doc

ID:60937967

大小:160.99 KB

页数:6页

时间:2021-01-04

(完整版)语音信号进行倒谱分析   .doc_第1页
(完整版)语音信号进行倒谱分析   .doc_第2页
(完整版)语音信号进行倒谱分析   .doc_第3页
(完整版)语音信号进行倒谱分析   .doc_第4页
(完整版)语音信号进行倒谱分析   .doc_第5页
资源描述:

《(完整版)语音信号进行倒谱分析 .doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、实验三语音信号进行倒谱分析一、实验目的、要求理解倒谱分析的作用3)掌握倒谱分析求基音周期的方法4)了解LPC倒谱分析方法5)二、实验原理1.倒谱分析原理同态信号处理也称为同态滤波,实现将卷积关系变换为求和关系的分离处理,即解卷。如x(n)=x1(n)*x2(n)进行如下3步处理(1)Z[x(n)]X(z)X1(z)X2(z)?(z)X?(z)X?(z)ilnX(z)lnX1(z)lnX2(z)X1211XjZ[X?(z)]Z[X?1(z)?2(z)]x?1(n)x?2(n)x?(n)对于语音信号进行解

2、卷,可将语音信号的声门激励信息及声道响应信息分离开来,从而求得声道共振特征和基音周期,用于语音编码、合成和识别。同态信号处理的基本原理(1))第一个子系统D*[](特征系统)完成将卷积信号转化为加性信号的运算。x?1(n)和x?2(n)信号也均是时域序列,但它们所处的离散时域显然不同于x(n)所处的离散时域,故把它称之为复倒频谱域。x?(n)是x(n)的复倒频谱,简称为复倒谱,有时也称为对数复倒谱。复倒谱具体计算公式x?(n)=Z-1[lnZ(x(n))]其中倒谱计算公式为:DFT[x(n)]=X(e

3、jw)X?(ejw)=lnX(ejw)x?(n)=IDFT[X?(ejw)]2线性预测原理线性预测分析的基本思想由于语音样点之间存在相关性,所以可以用过去的样点值来预测现在或未来的样点值。通过使实际语音抽样和线性预测抽样之间的误差在某个准则下达到最小值来决定唯一的一组预测系数,而这组系数就能反映语音信号的特性,可以作为语音信号特征参数来用于语音编码、语音合成和语音识别等应用中去。线性预测分析的基本原理每个采样值由前面的p个采样值线性组合所构成。记为x(n),有:1px¢(n)=?akx(n-k)k=1

4、pe(n)=x(n)-x¢(n)=x(n)-?akx(n-k)k=1要提高预测精度,就是要预测系数{ak}的取值使e(n)最小。理论上通常采用均方误差E[e2(n)]最小的准则。根据e(n)均方误差最小的原则来求解ak,有三种方法:自相关法(Levinson-Durbin算法)、协方差法和格型合成滤波算法。自相关法,就是先解出Yule-Walker方程,再计算G。Levinson-Durbin算法是自相关算法中的一种,形成递推算法。三、使用仪器、材料微机(带声卡)、耳机,话筒。四、实验步骤(1)采集语

5、音(浊音,轻音)信号(2)分帧计算语音信号倒谱(3)倒谱作FFT并加短时窗,取大于25以上的样值,进行IFFT,得到基音周期的倒谱。(3)运用Levinson-Durbin计算一帧语音信号线形预测系数(4)对该帧语音信号进行逆滤波处理(5)对逆滤波后的信号进行倒谱分析确定基音周期五、实验过程原始记录(数据,图表,计算)1.倒谱计算公式为:DFT[x(n)]=X(ejw)X?(ejw)=lnX(ejw)x?(n)=IDFT[X?(ejw)]2.线形预测系数计算公式为:0=R(0)1)Enn(2)1Rn(

6、i)i1ajRn(ij)ki(3)1i12)Eni1.aiki2.aii-1i-1=a-kajjii-j2i2i1E(1k)E3.ninifi

7、');%11000采样频率qq=max(sp);sp=sp-qq;%相关函数forj=1:pforn=j:N-1R(j)=R(j)+sp(n)*sp(n-j+1);endend%(Levinson-Durbin算法)预测系数%完成程序En(1)=R(1);K(1)=R(2)/R(1);a(1,1)=K(1);En(2)=(1-K(1)^2)*En(1);%K(2)=(R(3)-a(1,1)*R(2))/En(2);%a(2,2)=K(2);4%a(2,1)=a(1,1)-K(2)*a(1,1);for

8、i=2:ppo=0;forj=1:i-1po=po+a(i-1,j)*R(i-j+1);endK(i)=(R(i+1)-po)/En(i);a(i,i)=K(i);forj=1:i-1a(i,j)=a(i-1,j)-K(i)*a(i-1,i-j);endEn(i+1)=(1-K(i)^2)*En(i);endsource=sp(1:256);fori=1:12mm(i)=a(12,13-i);endmm1=zeros(1,300);5%逆滤波fork=1

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。