% Title: Main function
% Author: Lilei
% Data: 2008.6.2
% Description: Default parameters is fm=150,f0=500,k=0;
% Change: 
% Addition: 
clc 
clear all
close all
fm=100;
f0=200;
k=0;
n=4000;
[ ESCSpectrum , ESCf ] = ESCModel( f0 , fm , k , n );
[ a ]= ARFilter( ESCSpectrum , 1 , 20 );
[ h , w ] = freqz( 1 , [ 1 , a ] , 1000 );
figure;
plot( abs( h ));
figure;
plot(angle( h ));

WhiteNoise = randn( 1 , 2000 );
Signal = filter( 1 , [ 1 , a ] , WhiteNoise  );
figure;
subplot(221);
plot(WhiteNoise);
subplot(222);
plot(abs(fft(WhiteNoise)));
subplot(223);
plot(Signal);
subplot(224);
plot(pyulear(Signal,200));