//exec('pspect.1.ex')
//test modified periodogram method
//generate white data
rand('normal');
rand('seed',0);
x=rand(1:1024-33+1);
//make low-pass filter with eqfir
nf=33;
bedge=[0 .1;.125 .5];
des=[1 0];
wate=[1 1];
hn=eqfir(nf,bedge,des,wate);
//filter white data to obtain colored data 
h1=[hn 0*ones(1:maxi(size(x))-1)];
x1=[x 0*ones(1:maxi(size(hn))-1)];
hf=fft(h1,-1);
xf=fft(x1,-1);
yf=hf.*xf;
y=real(fft(yf,1));
