MATLAB语音预处理主体程序

%语音读入global X T xx,fs,bits=wavread(4.wav); sound(x,fs,bits); T=1/fs;X=fft(x,4000);axes(handles.axes1),plot(x);title(原始信号时域图);xlabel(时间)ylabel(幅值)axes(handles.axes2),plot(abs(X); title(原始信号频谱图);xlabel(频率/Hz);ylabel(幅值);%预加重global x xx xx=double(x);xx=filter(1 -0.9375,1,xx);%通过一个一阶高通滤波器进行滤波XX=fft(xx,4000);%figure(2);axes(handles.axes3),plot(xx);title(预加重后信号波形);xlabel(时间)ylabel(幅值)axes(handles.axes4),plot(abs(XX);title(预加重后信号频谱);xlabel(频率/Hz);ylabel(幅值)分帧% 设定帧移和帧长step=st; framelength =fr;% 分帧后如果有结余，则补足最后一帧数据lengthOfTailPad = framelength - rem( length(xx),framelength); %求余数tailPad=xx(length(xx)-lengthOfTailPad+1:length(xx) );nf=fix(length(xx)-framelength+step)/step); %截尾取整 计算帧数sn=zeros(nf,framelength); %219*256indf=step*(0:(nf-1).;inds=(1:framelength); % 语音分帧u=indf(:,ones(1,framelength)+inds(ones(nf,1),:);%219*256fx=x(u);sn(:)= x(u);t=sn(:);%figure(3);axes(handles.axes5),plot(fx);title(分帧处理后的波形);xlabel(帧数);ylabel(幅值/n);f=fft(fx);axes(handles.axes6),plot(abs(f);title(分帧处理后的频谱图);xlabel(帧数);ylabel(幅值/n);%过零检测global xx fx fr st x nf%figure(5)nt=length(xx);fxx=abs(fx);avap=sum(fxx,2);axes(handles.axes7),plot(avap);title(短时能量波形);xlabel(帧数)ylabel(能量)N=fr;for i=1:nt-1 if xx(i)=0 b(i)=1; else b(i)=-1; end if xx(i+1)=0 b(i+1)=1; else b(i+1)=-1; end w(i)=abs(b(i+1)-b(i);endk=1;j=0;while (k+N-1)nt Zn(k)=0; for i=0:N-1 Zn(k)=Zn(k)+w(k+i); end j=j+1; k=k+N/2;%每次移动半个窗口end for w=1:j Q(w)=Zn(st*(w-1)+1)/(2*N);%短时