% Title: ESC model
% Author: Lilei
% Data: 2008.6.2
% Description: Default parameters is fm=150,f0=500,k=0;
% Change: 
% Addition: 

function    [ ESCSpectrum , ESCf ] = ESCModel( f0 , fm , k , n )
f = 1 : n;
ESCSpectrum = (( fm + k .* ( f + f0 )) ./ ( fm .^ 2 + ( f + f0 ) .^ 2) + ( fm - k .* ( f - f0 )) ./ ( fm .^ 2 + ( f - f0 ) .^ 2 )) ./ pi;
ESCSpectrum = (ESCSpectrum - min( ESCSpectrum )) ./ (max( ESCSpectrum ) - min( ESCSpectrum ));
ESCf = f;
figure;
plot( ESCf , ESCSpectrum );