viewep.m

MFCC特征提取算法以及语音端点检测源码

clear all;  
fclose('all');

% Fp_Speech=fopen('E:\\已完成的工作\\2001年暑假的工作\\语音\\16k\\M_陈一宁_特定人数码采集.dat','r');%打开语音文件 .dat
Fp_Speech=fopen('E:\PiPi_Demo\皮皮\m_牟耕谷_皮皮_wj.seg','r');%打开语音文件 .dat
Fp_EP1=fopen('D:\语音库\合并的语音库及程序\file1\F_100_NoName_epd.txt','rt');%打开端点文件 .txt
% Fp_EP1=fopen('E:\\已完成的工作\\2001年暑假的工作\\语音\\16k\\M_陈一宁_特定人数码采集_epd.txt','rt');%打开端点文件 .txt
%Fp_EP2=fopen('E:\\已完成的工作\\2001年暑假的工作\\语音\\16k\\M_陈一宁_特定人数码采集_epd.txt','rt');

%读入语音文件总的句子数
fseek(Fp_Speech,198,'bof');
Sentence_Num=fread(Fp_Speech,1,'int16');
%读端点文件
tmpS = fscanf(Fp_EP1,'%s',1);
Sentence_Num1 = fscanf(Fp_EP1,'%s',1);

%tmpS = fscanf(Fp_EP2,'%s',1);
%Sentence_Num2 = fscanf(Fp_EP2,'%s',1);

%循环每一句话
offset_sentence=2; %从第几句话开始读语音

%读入一句话起始端点
for i=0:(offset_sentence-1),
    tmppoint = fscanf(Fp_EP1,'%ld',2);
end
    
for iSentence_Num=offset_sentence:Sentence_Num,

    %读入一句语音数据的首地址
    fseek(Fp_Speech,200,'bof');
    fseek(Fp_Speech,2*4*iSentence_Num,'cof');
    Off_Set=fread(Fp_Speech,1,'int32');
    %读入一句语音数据的总字节数
    SentenceLength_Byte=fread(Fp_Speech,1,'int32');
    fseek(Fp_Speech,Off_Set,'bof');
    %读入语音数据
    VoiceData_Short=fread(Fp_Speech,SentenceLength_Byte/2,'int16');
    
    Start1 = fscanf(Fp_EP1,'%ld',1);
    %读入一句话的终止端点
    End1 = fscanf(Fp_EP1,'%ld',1);
    %读入一句话起始端点
    %Start2 = fscanf(Fp_EP2,'%d',1);
    %读入一句话的终止端点
    %End2 = fscanf(Fp_EP2,'%d',1);

    x=1:SentenceLength_Byte/2;
    x = x %/ 8000;
    Start1 = Start1 / 8000;
    End1 = End1 / 8000;
    y=(-5000):20:(5000);
    %画图
    plot(x,VoiceData_Short,'g');
    hold on;
    plot(Start1,y,'r-',End1,y,'r-');
    axis('tight');
    hold off;
    string=int2str(iSentence_Num);
    title(string);
    %VoiceData_Short=VoiceData_Short/32768;
    %wavwrite(VoiceData_Short,8000,string);     
    pause;
end
fclose('all');