1、1实验三 语音信号进行倒谱分析一、 实验目的、要求1 理解倒谱分析的作用2 掌握倒谱分析求基音周期的方法3 了解 LPC 倒谱分析方法二、实验原理1倒谱分析原理同态信号处理也称为同态滤波，实现将卷积关系变换为求和关系的分离处理，即解卷。如进行如下 3 步处理 )()()()()3(lnlln2)()()1( 21211 21 nxnxzXZzXzXx 对于语音信号进行解卷，可将语音信号的声门激励信息及声道响应信息分离开来，从而求得声道共振特征和基音周期，用于语音编码、合成和识别。同态信号处理的基本原理（1）第一个子系统 D*（特征系统）完成将卷积信号转化为加性信号的运算。和 信号也均是时域序列
2、，但它们所处的离散时域显然不同于 x(n)所)(1nx)(2处的离散时域，故把它称之为复倒频谱域。 是 x(n)的复倒频谱，简称为复)(nx倒谱，有时也称为对数复倒谱。复倒谱具体计算公式其中倒谱计算公式为：2 线性预测原理线性预测分析的基本思想由于语音样点之间存在相关性，所以可以用过去的样点值来预测现在或未来的样点值。通过使实际语音抽样和线性预测抽样之间的误差在某个准则下达到最小值来决定唯一的一组预测系数，而这组系数就能反映语音信号的特性，可以作为语音信号特征参数来用于语音编码、语音合成和语音识别等应用中去。线性预测分析的基本原理12x=*1()ln()xZx-=()l()jwjwjDFTXe
3、exnI=2每个采样值由前面的 p 个采样值线性组合所构成。记为 x(n)，有：要提高预测精度，就是要预测系数 的取值使 e(n)最小。理论上通常采用均ka方误差 Ee2(n)最小的准则。根据 e(n)均方误差最小的原则来求解 ，有三种方法：自相关法k（Levinson-Durbin 算法） 、协方差法和格型合成滤波算法。自相关法，就是先解出 Yule-Walker 方程，再计算 G。Levinson-Durbin 算法是自相关算法中的一种，形成递推算法。三、使用仪器、材料微机（带声卡） 、耳机，话筒。四、 实验步骤（1）采集语音（浊音，轻音）信号（2）分帧计算语音信号倒谱（3）倒谱作 FFT
4、 并加短时窗，取大于 25 以上的样值，进行 IFFT，得到基音周期的倒谱。（3）运用 Levinson-Durbin 计算一帧语音信号线形预测系数（4）对该帧语音信号进行逆滤波处理（5）对逆滤波后的信号进行倒谱分析确定基音周期五、实验过程原始记录（数据，图表，计算）1倒谱计算公式为：2线形预测系数计算公式为：1） 0()nER=2） 1)(iniji jak3) iia1()()pkxnaxn=-1()pkexn=-()ljwjwjDFTxnXeeI34) 5) 12)(iniinEkif ip go to 1)6) pjja1jpLPC 倒谱流程如下六：实验结果，及分析倒谱分析LPC 倒谱
5、分析波形图1iiijjja-=4逆滤波后波形图分析倒谱分析可以观察倒浊音的基音周期。相邻两个峰值之间的距离为基音周期。七、实验参考程序N=256;R=zeros(1,N);p=12;a=zeros(p,p);En=zeros(1,p);K=zeros(1,p);sp=wavread(D:语音信号处理试验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)
6、/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); 5%a(2,1)=a(1,1)-K(2)*a(1,1); for i=2:p po=0; for j=1:i-1 po=po+a(i-1,j)*R(i-j+1); end K(i)=(R(i+1)-po)/En(i); a(i,i)=K(i); for j=1:i-1 a(i,j)=a(i-1,j)-K(i)*a(i-1,i-j); end En(i+1)=(1-K(i)2)*En(i); end source=sp(1
7、:256); for i=1:12 mm(i)=a(12,13-i); end mm1=zeros(1,300);6% 逆滤波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); end endendendDFTmm=abs(fft(mm1);Ln=zeros(1,300);%取对数，高频风量置零for i=1:60Ln(i)=logm(DFTmm(i);end% IDFTifftLn=real(ifft(Ln);mm=ifftLn;for i=1:150daopu1(i)=mm(151-i);endfor i=1:150daopu1(150+i)=mm(301-i);endplot(daopu1)
