I would like to compute a power spectrum in which the frequencies are logarithmically spaced.
In Welch's method there is a trade-off between the frequency resolution of the resulting power spectrum and the number of averages (i.e. error in the result). I would like this trade-off to be dynamic, i.e. do fewer averages for low-frequency points in order to have a finer resolution at low frequency.
Is there a standard way to do this?
I suppose one way would be to initially do
pwelch with a very high resolution (low number of averages), and then rebin the resulting spectrum using logarithmic binning.