语音信号处理实验报告 精选 .pdfVIP

语音信号处理实验报告

实验一

1用Matlab读取一段话音（自己录制一段，最好其中含有汉语四种声调变化，

该段话音作为本课实验原始材料），绘制原始语音波形图。

2.用Matlab计算这段语音的短时平均过零率、短时平均能量和短时平均幅度，

并将多个波形同步显示绘图。

3.观察各波形在不同音情况下的参数特点，并归纳总结其中的规律。

clcclc;

[x,fs]=wavread(benpao.wav);

figure

plot(x);

axis([0length(x)min(x)max(x)]);

title(原始语音波形)

xlabel(时间)

f=enframe(x,300,100);

[m,n]=size(f);

fori=1:m

energy(i)=sum(f(i,1:n).^2);

mn(i)=sum(abs(f(i,1:n)));

end

figure

plot(energy);

axis([0length(energy)min(energy)max(energy)]);

title(短时能量)

figure

plot(mn);

axis([0length(mn)min(mn)max(mn)]);

title(短时幅度)

lingd=zeros(m);

forx=1:m

fory=1:n-1

temp=f(x,y)*f(x,y+1);

iftemp=0

lingd(x)=lingd(x)+1;

end

end

%temp1=num(x,1)/300;

%count(x)=temp1;

end

figure

plot(lingd);

%axis([0length(lingd)min(lingd)max(lingd)]);

title(短时过零率)

子函数：

functionf=enframe(x,win,inc)%定义函数。

nx=length(x(:));%x(:)的作用是把x给弄成一个向量，x为一行，则变成一列，如果为

矩阵，则按每一列的顺序排成一列。得出的nx为序列的数据个数。

nwin=length(win);

if(nwin==1)

len=win;%如果win中就一个数，则len就=该数，此例中为256个点。即每帧长

else

len=nwin;%如果有多个数，则len=个数。

end

if(nargin3)%nargin返回的是函数输入的个数，如果中间有变量，返回的是负值。

inc=len;%也就是说，如果函数enframe的输入只有两个的话，系统就自动赋inc

end

nf=fix((nx-len+inc)/inc);

%这个比较关键，nf为分帧的组数

f=zeros(nf,len);%构成以组数为行，帧长为列的矩阵。

indf=inc*(0:(nf-1)).;%indf为一列nf个数据，即0到nf-1的inc倍，即分好了每

幀起点。

inds=(1:len);%构成了长度为len的一行。

f(:)=x(indf(:,ones(1,len))+inds(ones(nf,1),:));

if(nwin1)%nwin大于1的情况就不说了。

w=win(:);

f=f.*w(ones(nf,1),:);

end

原始语音波形

0.02

0.015

0.01

0.005

0

-0.005

-0.01

-0.015

-0.02

00.511.52

时间5

x10

短时能量

0.035

0.03

0.025

0.02