Anonymous

% sig = loadbin('csvowels\BLOCK113\V11302\V11302_A.BIN');
fs=48000; %hz
fmin = 0; %hz
fmax = fs/2; %hz
tlen = 0.015; %s
tstep=0.08; %s
wlen=floor(fs*tlen);
wstep=floor(fs*tstep);
n=5;
nfft=2^n;
M=22;
Wmel=melbf(M,fs,nfft,fmin,fmax);
Bmel =[];
   norm_specgram=spectrogram(sig,wlen,wlen-wstep,nfft,fs,'yaxis');
Bmel = Wmel*abs(norm_specgram);
[av,Epv,tv]=speechlpc(sig,fs,16,tlen*1000,tstep*1000); % vybereme prvnich 10 LPC koeficientu ( tlen a tstep fce vyzaduje v milisekundach )
meanLPC = mean(av(:,2:11));% prumer delame podle zadani z prvnich deseti koeficientu
mel_pwr_res.a =[];
mel_pwr_res.e =[];
mel_pwr_res.i =[];
mel_pwr_res.o =[];
mel_pwr_res.u =[];
mel_pwr_res.aa =[];
lpc_coeff_res.a =[];
lpc_coeff_res.e =[];
lpc_coeff_res.i =[];
lpc_coeff_res.o =[];
lpc_coeff_res.u =[];
lpc_coeff_res.aa =[];
final = [];
for v = 1:30
 mel_pwr_res.a = [mel_pwr_res.a (mel_pwr.a(:,v) - mean(Bmel,2));]; 
 mel_pwr_res.e = [mel_pwr_res.e (mel_pwr.e(:,v) - mean(Bmel,2));];  
 mel_pwr_res.i = [mel_pwr_res.i (mel_pwr.i(:,v) - mean(Bmel,2));];
 mel_pwr_res.o = [mel_pwr_res.o (mel_pwr.o(:,v) - mean(Bmel,2));];  
%  mel_pwr_res.aa = [mel_pwr_res.aa (mel_pwr.aa(:,v) - mean(Bmel,2));];  
%  mel_pwr_res.i = [mel_pwr_res.i (mel_pwr.i(:,v) - mean(Bmel,2));];
lpc_coeff_res.a  = [lpc_coeff_res.a; (lpc_coeff.a(v,:) - meanLPC); ];
lpc_coeff_res.e  = [lpc_coeff_res.e; (lpc_coeff.e(v,:) - meanLPC); ];
lpc_coeff_res.i  = [lpc_coeff_res.i; (lpc_coeff.i(v,:) - meanLPC); ];
lpc_coeff_res.o  = [lpc_coeff_res.o; (lpc_coeff.o(v,:) - meanLPC); ];
end
[h,l] = size(mel_pwr.u);
for v = 1:h
     mel_pwr_res.aa = [mel_pwr_res.aa (mel_pwr.aa(:,v) - mean(Bmel,2));]; 
 mel_pwr_res.u = [mel_pwr_res.u (mel_pwr.u(:,v) - mean(Bmel,2));]; 
 lpc_coeff_res.aa  = [lpc_coeff_res.aa; (lpc_coeff.aa(v,:) - meanLPC); ];
 lpc_coeff_res.u  = [lpc_coeff_res.u; (lpc_coeff.o(v,:) - meanLPC); ];
end
er_res_a = sum(sum(mel_pwr_res.a,2));
er_res_e = sum(sum(mel_pwr_res.e,2));
er_res_i = sum(sum(mel_pwr_res.i,2));
er_res_o = sum(sum(mel_pwr_res.o,2));
er_res_u = sum(sum(mel_pwr_res.u,2));
er_res_aa = sum(sum(mel_pwr_res.aa,2));
lpc_res_a = sum(sum(lpc_coeff_res.a,2));
lpc_res_e = sum(sum(lpc_coeff_res.e,2));
lpc_res_i = sum(sum(lpc_coeff_res.i,2));
lpc_res_o = sum(sum(lpc_coeff_res.o,2));
lpc_res_u = sum(sum(lpc_coeff_res.u,2));
lpc_res_aa = sum(sum(lpc_coeff_res.aa,2));
final = [er_res_a er_res_e er_res_i er_res_o er_res_u er_res_aa;lpc_res_a lpc_res_e lpc_res_i lpc_res_o lpc_res_u lpc_res_aa;];