倒谱计算与分析

1、视频语音处理技术倒谱计算与分析学院名称：计算机与信息工程学院专业名称：计算机科学与技术年级班级：姓名：学号：计算机与信息技术学院综合性、设计性实验报告专业：计算机科学技术年级/班级：2011级20122013学年第一学期课程名称视频语音处理技术指导教师张新明本组成员学号姓名1108114153王彦秋实验地点计科楼324实验时间项目名称倒谱计算与分析实验类型设计性一、实验目的：对语音信号进行同态分析可得到语音信号的倒谱参数。 语音的倒谱是将语音的短时谱取对数后再进行 IDFTIDFT 得到的，所以浊音信号的激励反映在倒谱上是同样周期的冲激，借此，可从倒谱波形中估计出基音周期。对倒谱进行低时窗选，

2、通过语音倒谱分析的最后一级，进彳 fDFTfDFT 后的输出即为平滑后的对数模函数，这个平滑的对数谱显示了特定输入语音段的谐振结构，即谱的峰值基本上对应于共振峰频率，对于平滑过的对数谱中的峰值进行定位，即可估计共振峰。对于倒谱计算与分析的设计实验可作如下训练：1 1、复倒谱的几种计算方法：2 2、最小相位信号法和递归法；3 3、基音检测；4 4、共振峰检测。二、实验仪器或设备：windowsXPwindowsXP 下的 MatlabMatlab 编程环境三、总体设计(设计原理、设计方案及流程等)1.1.复倒谱的几种计算方法：在复倒谱分析中，z z 变换后得到的是复数，所以取对数时要进行复对数运

3、算。这时存在相位的多值性问题，称为“相位卷绕”。设信号为x(n)=x(n)X2(n)则其傅里叶变换为X(ej)=X1(ej)X2(ej)对上式取复对数为lnX(ej)=lnX1(ej).lnXzlej)则其幅度和相位分别为：lnX(e，|=lnX4ej制+lnX2(e)()CD)()=l()：2()上式中，虽然9),%()的范围均在(冗，冗)内，但中心)的值可能超过(%江)范围。计算机处理时总相位值只能用其主值表示，然后把这个相位主值“展开”，得到连续相位。所以存在下面的情况：甲(0)二+2kn(K为整数)此时即产生了相位卷绕。下面介绍几种避免相位卷绕求复倒谱的方法。最小相位信号法这是解决相位

4、卷绕的一种较好的方法。但它有一个限制条件：被处理的信号想 x(n)x(n)必须是最小相位信号。实际上许多信号就是最小相位信号，或可以看作是最小相位信号。语音信号的模型就是极点都在 z z 平面单位圆内的全极点模型，或者极零点都在 z z 平面单位圆内的极零点模型。设信号x( (n) )的z变换为X( (z)=)=N( (z)/)/D(z),则有?N(z)4(z)=lnX(z)=ln-D(z)根据z变换的微分特性有D(z)N(z)-N(z)D(z)l二-zN(z)D(z)若x( (n) )是最小相位信号，则X(必然是稳定的因果序列。由 HilbertHilbert 变换的性质可知，任一因果复倒谱

5、序列都可分解为偶对称分量和奇对称分量之和：X(n)=Xe(n)X,(n)其中xe(n)-除n)x(n)1/2xo(n)=b?(n)一夕(-n)1/2这两个分量的傅里叶变换分别为？(的傅里叶变换的实部和虚部。父(ej)=&n)e二力武)j)?I(ej)n二二：从而可得0n：二0 x(n)=xe(n)n=02)?e(n)n0此即复倒谱的性质 3,3,山就是说一个因果序列可由其偶对称分量来恢复。如果引入一个oO二=-z义(z)=-zdz一喘辅助因子g( (n),上式可写作0g(n)=12可知x(0)=zInzx(0)U-Inx(0)、(n)=Inx(0)如果x( (n) )是最大相位序列，则

6、变为0 0n n0 0g(n)=g(n)= 1n=01n=02n02n0其中n n：0 0n n=0=0dd-zX(z)-X(z)=-z-X(z)dzdz对上式求逆z变换，根据z变换的微分特性，有n?(n)Lx(n)=nx(n)所以oOx(n)二三k=一二kiR(k)x(n-k)n)设x( (n) )是最小相位序列，而最小相位信号序列一定为因果序列，所以有nknx(n)=()*k)x(n-k)二二k=0nk=0n(k)x(n-k)及(n)x(0)由于5?(k)=0(kn)可得递推公式x(n)x(n)x(0)k田)x(n-k)x(0)递归运算后由复倒谱定义X(n)=z44nzk(n)。=z。ro

7、o、，lnZx(n)zL.n=oOJ最小相位信号法求复倒谱原理框图如下递归法这种方法仅限于是最小相位信号的情况。根据 z z 变换的微分特性得犯”(与?(k)nn0.08&indLF;a=fs/(LF+ind);elsea=0;endpitch=a%!图figure(1);subplot(3,1,1);plot(time1,y);title(语音波形);%axistightylim=get(gca,ylim);line(time1(startlndex),time1(startlndex),ylim,color,r);line(time1(endlndex),time1(endlnde

8、x),ylim,color,r);xlabel(样点数);ylabel(幅度);subplot(3,1,2);plot(frame);%axis(0,400,-0.5,0.5)title(一帧语音);xlabel(样点数);ylabel(幅度)subplot(3,1,3);time2=-199:1:-1,0:1:200;plot(time2,cepstrum1);%axis(-200,200,-0.5,0.5)title(,一帧语音的倒谱,);加汉明窗求倒谱基因周期的范围是70500Hzxlabel（样点数）;ylabel（幅度）;运行结果如下图：倒谱法求浊音的基音周期样点数样点数一帧语音的倒

9、谱0硼业郴际的如岫曲啊的即H泮w蚓M垢N脑出晒晒I1嶙口511 1 1,aso10a15D20a250aoa35040a语音波形jog0511.522.5335样点数/十一帧语音曾曾0.40.20-0.2-0.4语音波形语音波形0.511522.533.E样点数一帧语音样点数一帧语音0.51111111样点数样点数清音的倒谱3 3. .共振峰检测程序waveFile=qinghua.wav;y,fs,nbits=wavread(waveFile);time=(1:length(y)/fs;frameSize=floor(40*fs/1000);%帧长startIndex=round(15000

10、);%起始序号endIndex=startIndex+frameSize-1;%结束序号frame=y(startIndex:endIndex);%取出该帧frameSize=length(frame);frame2=frame.*hamming(length(frame);%力口汉明窗rwy=rceps(frame2);%倒谱求ylen=length(rwy);cepstrum=rwy(1:ylen/2);%基首检测LF=floor(fs/500);HF=floor(fs/70);cn=cepstrum(LF:HF);mx_cepind=max(cn);%找到最大的突起的位置%共振峰检测核心

11、代码NN=ind+LF;ham=hamming(NN);cep=cepstrum(1:NN);ceps=cep.*ham;%汉明窗formant1=20*log(abs(fft(ceps);formant(1:2)=formant1(1:2);fort=3:NN%dosomemedianfilteringz=formant1(t-2:t);md=median(z);formant2(t)=md;endfort=1:NN-1ift=2formant(t)=formant1(t);elseformant(t)=formant2(t-1)*0.25+formant2(t)*0.5+formant2(

12、t+1)*0.25;endendsubplot(3,1,1);plot(cepstrum);title(倒谱);xlabel(样点数);ylabel(幅度)%axis(0,220,-0.5,0.5)spectral=20*log(abs(fft(frame2);subplot(3,1,2);xj=(1:length(spectral)/2)*fs/length(spectral);plot(xj,spectral(1:length(spectral)/2);title(频谱);xlabel(频率/Hz);ylabel(幅度/dB)%axis(0,500,-100,50)subplot(3,1,3);xi=(1:NN/2)*fs/NN;plot(xi,formant(1:NN/2);title(平滑对数幅度谱);xlabel(频率/Hz);ylabel(幅度/dB);%axis(0,5500,-80,0)运行结果如图所示：0.50-0.5020406080100120140160180200220样点数0500100015