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

1、实验三 语音信号进行倒谱分析一、 实验目的、要求1 理解倒谱分析的作用2 掌握倒谱分析求基音周期的方法3 了解LPC 倒谱分析方法二、实验原理1 倒谱分析原理同态信号处理也称为同态滤波，实现将卷积关系变换为求和关系的分离处理，即解卷。如 x(n) = x1(n)* x2(n)进行如下 3 步处理(1)Zx(n) X(z)X1(z) X2(z)(2)lnX(z) lnX1(z) lnX2(z) X?1(z) X?2(z) X?(z)(3)Z1X?(z)Z 1X?1(z) X?2(z)x?1(n) x?2(n) x?(n)对于语音信号进行解卷，可将语音信号的声门激励信息及声道响应信息分离开来，从而

2、求得声道共振特征和基音周期，用于语音编码、合成和识别。同态信号处理的基本原理( 1)第一个子系统D* (特征系统)完成将卷积信号转化为加性信号的运算。xi(n)和？2(n)信号也均是时域序列，但它们所处的离散时域显然不同于x(n)所处的离散时域， 故把它称之为复倒频谱域。 x?(n) 是 x(n) 的复倒频谱， 简称为复倒谱，有时也称为对数复倒谱。复倒谱具体计算公式x?(n) = Z- 1ln Z(x(n)其中倒谱计算公式为：DFTx(n) = X(ejw)X?(ejw) = ln X(ejw)x?( n) = IDFT X?(ejw)2 线性预测原理线性预测分析的基本思想由于语音样点之间存在

3、相关性， 所以可以用过去的样点值来预测现在或未来的样点值。 通过使实际语音抽样和线性预测抽样之间的误差在某个准则下达到最小值来决定唯一的一组预测系数， 而这组系数就能反映语音信号的特性， 可以作为语音信号特征参数来用于语音编码、语音合成和语音识别等应用中去。线性预测分析的基本原理每个采样值由前面的 p 个采样值线性组合所构成。记为 x (n) ，有：px(n)= ? akx(n- k)k = 1Pe(n) = x(n)- x(n) = x(n)- ? akx(n- k)k= 1要提高预测精度，就是要预测系数aj的取值使e(n)最小。理论上通常采用均方误差Ee2(n)最小的准则。根据e(n)均方

4、误差最小的原则来求解ak ,有三种方法：自相关法 (Levinson-Durbin算法)、协方差法和格型合成滤波算法。自相关法，就是先解出 Yule-Walker方程，再计算 G。Levinson-Durbin算法 是自相关算法中的一种，形成递推算法。、使用仪器、材料微机(带声卡)、耳机，话筒。1 .2.D2)3)4)四、实验步骤(1)采集语音(浊音，轻音)信号(2)分帧计算语音信号倒谱(3)倒谱作FFT并加短时窗，取大于 25以上的样值，进行IFFT,得到基音周期的倒(3)运用Levinson-Durbin计算一帧语音信号线形预测系数(4)对该帧语音信号进行逆滤波处理(5)对逆滤波后的信号进

5、行倒谱分析确定基音周期五、实验过程原始记录(数据，图表，计算) 倒谱计算公式为：DFTx(n) = X (ejw)*(ejw)= ln X(ejw)5?(n)= IDFT)?(ejw) 线形预测系数计算公式为：E0= Rn(0) i 1 -i 1 Rn(i)aj Rn(i j)K Ena i k i i i- 1i- 1aj = aj - Rai- j75) En (1 if ip go toki2)En16) a j1 j P1)a pLPC倒谱流程如下原始语音信号进行分帧加窗DFTMl高频分量置零基音检测六：实验结果，及分析倒谱分析LPC倒谱分析波形图逆滤波后波形图分析倒谱分析可以观察倒浊

6、音的基音周期。相邻两个峰值之间的距离为基音周 期。七、实验参考程序N=256;R=zeros(1,N);p=12;a=zeros(p,p);En=zeros(1,p);K=zeros(1,p);语音信号处理试验a1.wav);%11000采样频率 qq=max(sp);sp=sp-qq;%相关函数for j=1:pfor n=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)A2)*En(1);%K(2)=(R

7、(3)-a(1,1)*R(2)/En(2);%a(2,2)=K(2);%a(2,1)=a(1,1)-K(2)*a(1,1);for i=2:ppo=0;for j=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);for j=1:i-1a(i,j)=a(i-1,j)-K(i)*a(i-1,i-j);endEn(i+1)=(1-K(i42)*En(i);endsource=sp(1:256);for i=1:12mm(i)=a(12,13-i);endmm1=zeros(1,300);% 逆滤波for k=1:256if kpfor m=1:kmm1(k)=mm1(k)+mm(p-k+m)*source(m);endelseif k256for m=1:pmm1(k)=mm1(k)+mm(m)*source(m+k-p);endelsefor m=k-p+1:256mm1(k)=mm1(k)+mm(m-k+p)*source(m);endendendendDFTmm=abs(fft(mm1);Ln=zeros(1,300);%取对数，高频风量置零for i=1:60Ln(i)=logm(DFTmm(i);e