p=noiseModel( f, sigma_w, f_knee, alpha) Returns natural log og the noise power spectrum for given frequencies f with noise statistics sigma_w, f_knee and alpha p = ( 1 + (f_knee./x).^alpha ).*(sigma_w^2)
0001 function logp = noiseModel( f, sigma_w, f_knee, alpha ) 0002 % p=noiseModel( f, sigma_w, f_knee, alpha) 0003 % Returns natural log og the noise power spectrum for given frequencies f 0004 % with noise statistics sigma_w, f_knee and alpha 0005 % p = ( 1 + (f_knee./x).^alpha ).*(sigma_w^2) 0006 emConst = 5.7721566490153286060; 0007 logp = real( 2.*log(sigma_w) + log( 1 + (f_knee./f).^alpha ) )+emConst; 0008 %p = ( 1 + (f_knee./x).^alpha ).*(sigma_w^2); 0009 end