FREQCOR_functions.fit_func

FREQCOR_functions.fit_func(freq_binned, a, b, d, PS_original)[source]

@author: Aslan et al. 2021 Translated by A. Faurès from Matlab to Python on 21/10/2024

Function to fit power spectra. The function is defined as the sum of an attenuated turbulent signal component and an independent white noise compontent.

Parameters:

  • freq_binned (numpy array) – Natural frequency, binned.

  • a (float64) – Low-pass filtering time constant (s).

  • b (float64) – Transfer function normalisation factor.

  • d (float64) – Intercept of white noise line with the y-axis.

  • PS_original (numpy array) – Binned power spectra of sonic temperature (unattenuated, noise-free). Is in the form: f*sp(T)/var(T)

Returns:

Attenuated and noise-contaminated power spectra.

Return type:

Numpy array