Encontrando los coeficientes cepstrales.

cepstrum.m

function [c] = cepstrum(segs)

%Encontrar los coeficientes cepstrales de cada segmento
[M,N] = size(segs); 
% M filas, N columas, 
% Definir Mo=output frame size(# de %frecuencias a las que se aplica FFT)
[c] = zeros(M,10); 
% Inicialiando el vector para guardar los coeficientes cepstrales. 
% P=10

% Calculando los coeficientes cepstrales reales
for i=1:M
    r = rceps(segs(i,:));
    c(i,:)= r(1:10);
end

% Normalizacion de los coeficientes cepstrales
avg_c = sum(c)/M;

% Sustrayendo minimo, solo las que estan arriba del promedio de los fonemas
for i=1:M
    c(i,:)= c(i,:)-avg_c; % smoothened spectrum
end