pymultifracs.psd.plot_psd#

pymultifracs.psd.plot_psd(signal, fs, n_fft=4096, seg_size=None, n_moments=2, log='log2', ax=None, wt='db')#

Plot the superposition of Fourier-based Welch estimation and Wavelet-based estimation of PSD on a log-log graphic.

Based on the wavelet_fourier_spectrum function of the MATLAB toolbox mf_bs_toolbox-v0.2 [1]

Parameters:
signalndarray, shape (n_samples,)

Time series of sampled values

fsfloat

Sampling frequency of signal

n_fftint, optional

Length of the FFT desired. If seg_size is greater, n_fft = seg_size.

seg_sizeint | None

Length of Welch segments. Defaults to None, which sets it equal to n_fft

n_momentsint

Number of vanishing moments of the Daubechies wavelet used in the Wavelet decomposition.

logstr

Log function to use on the frequency and power axes

References