return min(fitted_variance, naive_variance)
@splittableKwargsFunction()
-def fit_psd(freq_axis, psd_data, minFreq=500, maxFreq=7000) :
+def fit_psd(freq_axis, psd_data, minFreq=500, maxFreq=25000) :
"""
freq_axis : array of frequencies in Hz
psd_data : array of PSD amplitudes for the frequencies in freq_axis