Matlab - Automatizar funcion

clear all

ruta=pwd;

clc

diract = cd;

directorio = uigetdir(cd,...

    'ELEGIR CARPETA DE AUDIOS QUE SE DESEAN CARGAR');

if isequal(directorio,0)

    msgbox(...

        'No ha elegido ninguna carpeta',...

        'MENSAJE',...

        'help')

    return

end

cd(directorio)

lista = dir('*wav');

[n,~] = size(lista);

if isequal(n,0)

    msgbox(...

        'No hay ninguna imagen en la carpeta especificada',...

        'MENSAJE',...

        'help')

    cd(diract)

    return

end

names = cell(n,1);

for i = 1:n

    names{i} = lista(i).name;

end

cd(ruta)

for i = 1:n

   A{i}= audioread([directorio,'\',names{i}],[500,20799]);

    B=cell2mat(A);

    B = (bsxfun(@rdivide, bsxfun(@minus, B, mean(B)), std(B)));

    tw=25;

ts=10;

alpha=0.97;

R=[300 4000];

N=13;

L=22;

hamming = @(N)(0.54-0.46*cos(2*pi*[0:N-1].'/(N-1)));

M=20;

Fs=20000;

% n=512;              %Number of FFT points

% Tf=0.032;           %Frame duration in seconds

% N=Fs*Tf;            %Number of samples per frame

% fn=24;              %Number of mel filters

% l=length(B);        %total number of samples in speech

% Ts=0.016;            %Frame step in seconds

% FrameStep=Fs*Ts;    %Frame step in samples

% a=1;

% b=[1, -0.97];       %a and b are high pass filter coefficients

 [MFCC]= mfcc(B(:,i),20000,tw,ts,alpha,hamming,R,M,N,L);

% opt= mfccOptSet(20000)

% opt.frameSize=3000;

%  [MFCC(:,i)]= wave2mfcc(B(:,i), 20000, opt)

% [MFCC(:,i)]=kannumfcc(13,B(:,i),Fs);

% MFCC=  bsxfun(@minus, MFCC, mean(MFCC));

    MFILA=MFCC';

    [MPF,f]=lpcar2pf(MFILA,127);

%     MPF=sum(MPF)/128;

%     MPF=pow2db(MPF);

%     MPF2=transpose(MPF);

%     norm = transpose(bsxfun(@rdivide, bsxfun(@minus, MPF2, mean(MPF2)), std(MPF2)));

end

 for i =1:n/2

SR(i,:)=MPF(2*i,:);

CR(i,:)=MPF(2*i-1,:);

CB=pdist2(SR,CR,'cityblock');

CH=pdist2(SR,CR,'chebychev');

dCB=diag(CB);

dCH=diag(CH);

HM=pdist2(SR,CR,'hamming');

dHM=diag(HM);

      fig=figure;

      plot(f,MPF(2*i-1,:),f,MPF(2*i,:));

      xlabel('Frecuencia normalizada')

      ylabel('Espectro de potencia MFCC [dB]')

      legend('Posfiltrada','Limpia')

      print(fig,['miplot',num2str(i)],'-dpng');

      pause(0.5)

       eu(i)=disteusq(MPF(2*i-1,:),MPF(2*i,:),'s');

      deu=eu';

       pm=corrcoef(MFILA(2*i-1,:),MFILA(2*i,:)) ;

       pcorM(i)=pm(1,2);

       pcorM=pcorM';

      rmseM(i,:)=sqrt(sum((MFILA(2*i-1,:)-MFILA(2*i,:)).^2))*(1/13);

       EUCLIDIANA= xlswrite('MFCC2.xlsx',deu,'A2:A26');

       cosh(i)=abs(distchpf(MPF(2*i-1,:),MPF(2*i,:)));

      dcosh=10*cosh';

      COSH= xlswrite('MFCC.xlsx',dcosh,'B2:B26');

      DCB= xlswrite('MFCC2.xlsx',dCB,'b2:b26');

      DCH= xlswrite('MFCC2.xlsx',dCH,'c2:c26');

      DHM= xlswrite('MFCC.xlsx',dHM,'E2:E26');

       p=corrcoef(MPF(2*i-1,:),MPF(2*i,:)); 

       pcor(i)=p(1,2);

      pcor=pcor';

      rmse(i,:)=sqrt(sum((MPF(2*i-1,:)-MPF(2*i,:)).^2)*(1/512));

     end

cd(diract)