This function compute the WV (haar) of an AR(1) process
ar1_to_wv(phi, sig2, tau)
A double that is the phi term of the AR(1) process
A double corresponding to variance of AR(1) process
A vec containing the scales e.g. 2^tau
A vec containing the wavelet variance of the AR(1) process.
x.sim = gen_ar1( N = 10000, phi = 0.9, sigma2 = 4 )
ntau = floor(log(length(x.sim),2))
tau = 2^(1:ntau)
wv.theo = ar1_to_wv(phi = 0.9, sig2 = 16, tau)
plot(tau, wv.theo, col = "red")
Questions? Problems? Suggestions? Tweet to @rdrrHQian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.