mne.time_frequency.multitaper_psd(x, sfreq=6.283185307179586, fmin=0, fmax=inf, bandwidth=None, adaptive=False, low_bias=True, n_jobs=1, normalization='length', verbose=None)

Compute power spectrum density (PSD) using a multi-taper method


x : array, shape=(n_signals, n_times) or (n_times,)

The data to compute PSD from.

sfreq : float

The sampling frequency.

fmin : float

The lower frequency of interest.

fmax : float

The upper frequency of interest.

bandwidth : float

The bandwidth of the multi taper windowing function in Hz.

adaptive : bool

Use adaptive weights to combine the tapered spectra into PSD (slow, use n_jobs >> 1 to speed up computation).

low_bias : bool

Only use tapers with more than 90% spectral concentration within bandwidth.

n_jobs : int

Number of parallel jobs to use (only used if adaptive=True).

normalization : str

Either “full” or “length” (default). If “full”, the PSD will be normalized by the sampling rate as well as the length of the signal (as in nitime).

verbose : bool, str, int, or None

If not None, override default verbose level (see mne.verbose).


psd : array, shape=(n_signals, len(freqs)) or (len(freqs),)

The computed PSD.

freqs : array

The frequency points in Hz of the PSD.

See also, mne.Epochs.plot_psd


New in version 0.9.0.