A matlab toolbox for cross-wavelet and wavelet coherence analysis
This project is maintained by grinsted
It is usually an appropriate null hypothesis if the theoretical AR1 spectrum is degreesa good model degrees for the power decay in the observed spectrum. I recommend to simply visually compare the two power spectra:
X=rednoise(200,.8);
[P,freq]=pburg(zscore(X),7,[],1);
aa=ar1(X);
Ptheoretical=(1-aa.^2)./(abs(1-aa.*exp(-2*pi*i*freq))).^2;
semilogy(freq,P/sum(P),freq,Ptheoretical/sum(Ptheoretical),'r');
legend('observed',sprintf('Theoretical AR1=%.2f',aa),'location','best')