MATLAB语音预处理主体程序(共3页)

1、精选优质文档-倾情为你奉上%语音读入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('幅值')

2、;%预加重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('幅值')分帧%

3、 设定帧移和帧长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:framelen

4、gth); % 语音分帧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('帧数');ylab

5、el('幅值/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);%短时平均过零检测 endaxes(handles.a