第一步:读入已经录好的 WAV 声音文件;
[k,fs]=wavread('zmkm');%读入一个声音文件
_____________
第二步:预处理
%下面的这段程序用于提取.WAV 文件中的有声部分
[k,fs]=wavread('zmkm');
h=length(k);
[M,N]=size(k); %[m,n]是大小为 K 的矩阵
x=linspace(0,M,M); %linspace 线性间矩矢量
subplot(2,1,1);
plot(x,k);
p1=abs(k);
for i=1:M,
B(i)=k(M-i+1); %取数组 k 的倒序给 B
end
p2=abs(B);
th=0.035;
for i=1:M,
if p1(i)>th
break;
end
end
for j=1:M,
if p2(j)>th
break;
end
end
new=k(i:(M-j));
L=length(new);
subplot(2,1,2);
y=linspace(0,L,L);
plot(y,new);
________________
%本段程序用于对语音文件进行预加重
%输入.WAV 语音文件,输出预加重后的文件 B
[k,fs]=wavread('zmkm');
% [B,no]=readwavdata('zmkm');
wavplay(k,fs ); %播放语音文件
h=length(k); %计算语音文件的
长度