matlab代码 模板匹配法 语音信号 数字识别（0到1）

clear all close clc; disp('正在计算参考模板的参数...') for i=1:10 fname=sprintf('%d0.mp3',i-1); [x,fs]=audioread(fname); sound(x,fs); [x1 x2]=vad(x); % x1=abs(round(x1)); % x2=abs(round(x2)); m=mfcc(x); m=m(x1-2:x2-4,:); ref(i).mfcc=m; end disp('正在分析语音信号...') for i=1:10 fname=sprintf('%d0.mp3',i-1); %[x,fs,bit]=audioread(fname,[2000,2512]); %采样% [x,fs]=audioread(fname,[2000,2512]); %sound(x,fs); %播放语音信号 figure(i); subplot(3,3,1); plot(x(1:256)); %原始语音信号的时域图形% title('原始信号') subplot(3,3,2) [h,w]=freqz(x,fs); %原始语音信号的频率响应图 hr=abs(h); %求系统幅频响应 plot(w,hr); title('幅频图'); xlabel('Frequency in rad mple') ylabel('Magnitude in dB') subplot(3,3,3) hphase=angle(h); hphase=unwrap(hphase); %求系统相频响应 plot(w,hphase); title('相频图'); xlabel('Frequency in rad mple') ylabel('Phase in degrees') y=fft(x,512); %傅立叶变换% mag=abs(y); mag1=10*log10(mag); f=fs*(0:255)/512; subplot(3,3,4) plot(f,mag(1:256)); %FFT频谱图 % title('fft变换后信号') iff=ifft(y,512); %反傅立叶变换% ifm=abs(iff); subplot(3,3,5) plot(f,ifm(1:256)) title('ifft后信号') % 短时傅里叶变换 Ts=1/fs; %N=T/Ts; N=512; Nw=20; %窗函数长 L=Nw/2; %窗函数每次移动的样点数 Tn=(N-Nw)/L+1; %计算把数据x共分成多少段 nfft=32;%FFT的长度 Tn=round(Tn); TF=zeros(Tn,nfft); %将存放三维谱图，先清零 for i=1:Tn xw=x((i-1)*10+1:i*10+10); %取一段数据 temp=fft(xw,nfft); %FFT变换 temp=fftshift(temp); %频谱以0频为中心 for j=1:nfft TF(i,j)=temp(j); %把谱图存放在TF中 end end subplot(3,3,6) fnew=((1:nfft)-nfft/2)*fs/nfft; tnew=(1:Tn)*L*Ts; [F,T]=meshgrid(fnew,tnew); mesh(F,T,abs(TF)) title('短时傅立叶变换时频图') subplot(3,3,7) contour(F,T,abs(TF)) title('等高线表示') end disp('正在计算测试模板的参数...') for i=1:10 fname=sprintf('%d2.mp3',i-1); x=audioread(fname); [x1 x2]=vad(x); m=mfcc(x); m=m(x1-2:x2-4,:); test(i).mfcc=m; end disp('正在进行模板匹配...') dist=zeros(10,10); for i=1:10 for j=1:10 dist(i,j)=dtw(test(i).mfcc,ref(j).mfcc); end end disp('正在计算匹配结果...') for i=1:10 [d,j]=min(dist(i,:)); fprintf('测试信号%d的识别结果为:%d\n',i-1,j-1); end