语音实验三信号的特征提取 .pdfVIP

语音实验三信号的特征提取

实验三语音信号的特征提取

一、实验目的

1、熟练运用MATLAB软件进行语音信号实验。

2、熟悉短时分析原理、LPC的原理。

3、学习运用MATLAB编程进行LPC的提取。

4、学会利用短时分析原理提取LPC特征序列。二、实验仪器设

备及软件MATLAB三、实验原理

由于频率响应H(ω

je

)反映声道的频率响应和被分析信号的谱包络，因此用log)

(ω

je

H做反傅里叶变换求出的LPC倒谱系数。

通过线性预测分析得到的合成滤波器的系统函数为H(z)=1

-∑=-p

iiiza11/,其冲激响应为h(n)。h(n)的倒谱为∧

h(n)，∧

H(z)=

)(1

nhn∑+∞=∧

n

z

-就是说∧H(z)的逆变换∧

h(n)是存在的。设

∧h(0)=0，将式∧

H(z)=)(1

nhn∑+∞=∧

nz-两边同时对1-z求导，得

-??∑=--piiza111

11

logz=1-??z)(1

nhn∑+∞=∧

n

z

-得到

1

1

)(+-+∞

=∧

∑nnz

nhn=

∑∑=-=+--p

ii

iiii

zaz

ia1

p

1

1

1，于是有

∑=-p1i1ia-1z1

1

)(+-+∞

=∧

∑nnz

nhn令其左右两边z的各次幂前系数分别相等，得到

h(n)和

间的递推关系

按其可直接从预测系数｛a｝求得倒谱h(n).这个倒谱是根据线性

预测模型得到的，又称为LPC倒谱。LPC倒谱由于利用线性预测中声

道系统函数H（z）的最小相位特性，因此避免了一般同态处理中求

复对数的麻烦。

四、实验步骤及程序LPC(1)、实验步骤

1、输入原始语音

2、对样本语音进行加窗处理

3、计算LPC系数

4、建立语音正则方程

5、输出原始样本语音、预测语音波形和预测误差

6、输出LPC谱

7、求出预测误差的倒谱

8、输出原始语音和预测语音的语谱图

（2）实验流程

（3）实验源程序

[y,Fs]=wavread(1.wav);

e_x=(frame(y,lpc_spectrum,Fs));

%plot(e_x(2,:));

plot(e_x(:,89))；

holdon;

c=melcepst(y,Fs);

plot(c(89,:),k)

frame定义

%function[y]=frame(x,func,SAMP_FREQ,l,step)

%whereyisoutputonaframebyframebasis,xisinput

speech,

%andlisthewindowsize.landstepareoptionalparameters,

%bydefaultSAMP_FREQis8000,lis200,andstepis100.

%funcisastringe.g.pitchthatdeterminesafunctionthat

youwant%toapplytoxonashort-timebasis.

%

%Writtenby:LeventArslanApr.11,1994

%

function[yy]=frame(x,func,SAMP_FREQ,l,step)

[m,n]=size(x);

ifmn

n=