mne.time_frequency.psd_array_multitaper(x, sfreq, fmin=0, fmax=inf, bandwidth=None, adaptive=False, low_bias=True, normalization=’length’, n_jobs=1, verbose=None)[source]

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


x : array, shape=(…, 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.

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).

n_jobs : int

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

verbose : bool, str, int, or None

If not None, override default verbose level (see mne.verbose() and Logging documentation for more).


psds : ndarray, shape (…, n_freqs) or

The power spectral densities. All dimensions up to the last will be the same as input.

freqs : array

The frequency points in Hz of the PSD.


New in version 0.14.0.