function D=spec_dist_ceps(c1,c2);
% Syntax: D=spec_dist_ceps(c1,c2)
%
% Calcule la distortion spectrale en dB defini par
% D=sqrt(1/2pi * integral de -pi a pi (10 log10 A1 - 10 log10 A2)^2 domega)
% ou A1, A2 sont les spectres de puissance LPC: P(omega)=1/|A(omega)|^2
% a partir du cepstre complexe en utilisant le theoreme de Parceval.
% c1, c2 sont les vecteurs cepstraux de type
% c=[0 c(1) c(2) ... c(N)] calcules par la fonction a_to_cepst

D2=100/((log(10))^2)*2*sum((c1-c2).^2); % constante pour la conversion en DB

% constante 2 - les coeffs. ont une partie
% symetrique qu'il faut aussi prendre en compte !
D=sqrt(D2);