F-test



In [2]:

    
import sys
sys.path.append("..")
import splitwavepy as sw

import numpy as np
import matplotlib.pyplot as plt

The energy in the noise is chi-squared distributed.



In [28]:

    
nsamps = 1001
sigma = 1
nptsmooth = 1
energy = [ np.sum(sw.core.noise(nsamps,sigma,nptsmooth)**2) for _ in range(1000) ]
# plot a histogram of the energy on randomly generated smoothed white noise
plt.hist(energy)
plt.show()

The amplitudes in the noise are Gaussian distributed.



In [29]:

    
nsamps = 10001
sigma = 1
nptsmooth = 1
noise = sw.core.noise(nsamps,sigma,nptsmooth)
# plot histogram
plt.hist(noise)
plt.show()



In [ ]:



In [ ]: