Note: This page has been translated by MathWorks. Please click here

To view all translated materials including this page, select Japan from the country navigator on the bottom of this page.

To view all translated materials including this page, select Japan from the country navigator on the bottom of this page.

Multitaper power spectral density estimate

`pxx = pmtm(x)`

`pxx = pmtm(x,nw)`

`pxx = pmtm(x,nw,nfft)`

`[pxx,w] = pmtm(___)`

`[pxx,f] = pmtm(___,fs)`

`[pxx,w] = pmtm(x,nw,w)`

`[pxx,f] = pmtm(x,nw,f,fs)`

`[___] = pmtm(___,method)`

`[___] = pmtm(x,e,v)`

`[___] = pmtm(x,dpss_params)`

`[___] = pmtm(___,'DropLastTaper',dropflag)`

`[___] = pmtm(___,freqrange)`

`[___,pxxc] = pmtm(___,'ConfidenceLevel',probability)`

`pmtm(___)`

returns
Thomson’s multitaper power spectral density (PSD) estimate, `pxx`

= pmtm(`x`

)`pxx`

,
of the input signal, `x`

. When `x`

is
a vector, it is treated as a single channel. When `x`

is
a matrix, the PSD is computed independently for each column and stored
in the corresponding column of `pxx`

. The tapers
are the discrete prolate spheroidal (DPSS), or Slepian, sequences.
The time-halfbandwidth, `nw`

, product is 4. By
default, `pmtm`

uses the first 2`nw`

– 1 DPSS sequences. If `x`

is
real-valued, `pxx`

is a one-sided PSD estimate.
If `x`

is complex-valued, `pxx`

is
a two-sided PSD estimate. The number of points, `nfft`

,
in the discrete Fourier transform (DFT) is the maximum of 256 or the
next power of two greater than the signal length.

`[`

returns a frequency vector, `pxx`

,`f`

] = pmtm(___,`fs`

)`f`

, in cycles per unit time. The sampling
frequency, `fs`

, is the number of samples per unit time. If the unit of time
is seconds, then `f`

is in cycles/sec (Hz). For real–valued signals,
`f`

spans the interval [0,`fs`

/2] when
`nfft`

is even and [0,`fs`

/2) when
`nfft`

is odd. For complex-valued signals, `f`

spans the
interval [0,`fs`

). `fs`

must be the fourth input to
`pmtm`

. To input a sample rate and still use the default values of the
preceding optional arguments, specify these arguments as empty, `[]`

.

`[`

returns the two-sided multitaper PSD estimates at the frequencies specified in the vector,
`pxx`

,`f`

] = pmtm(`x`

,`nw`

,`f`

,`fs`

)`f`

. `f`

must contain at least two elements. The
frequencies in `f`

are in cycles per unit time. The sampling frequency,
`fs`

, is the number of samples per unit time. If the unit of time is
seconds, then `f`

is in cycles/second (Hz).

`[___] = pmtm(`

uses
the cell array, `x`

,`dpss_params`

)`dpss_params`

, to pass input arguments
to `dpss`

except the number of elements in the
sequences. The number of elements in the sequences is the first input
argument to `dpss`

and is not included in `dpss_params`

.
An example of this usage is `pxx = pmtm(randn(1000,1),{2.5,3})`

.

`[___] = pmtm(___,'DropLastTaper',`

specifies
whether `dropflag`

)`pmtm`

drops the last taper in the computation
of the multitaper PSD estimate. `dropflag`

is a
logical. The default value of `dropflag`

is `true`

and
the last taper is not used in the PSD estimate.

`[___] = pmtm(___,`

returns
the multitaper PSD estimate over the frequency range specified by `freqrange`

)`freqrange`

.
Valid options for `freqrange`

are `'onesided'`

, `'twosided'`

,
and `'centered'`

.

`[___,`

returns
the `pxxc`

] = pmtm(___,'ConfidenceLevel',`probability`

)`probability`

× 100%
confidence intervals for the PSD estimate in `pxxc`

.

`pmtm(___)`

with no output arguments
plots the multitaper PSD estimate in the current figure window.

[1] Percival, D. B., and A. T. Walden, *Spectral
Analysis for Physical Applications: Multitaper and Conventional Univariate
Techniques*. Cambridge, UK: Cambridge University Press,
1993.

[2] Thomson, D. J., “Spectrum estimation
and harmonic analysis.” *Proceedings of the IEEE ^{®}*.
Vol. 70, 1982, pp. 1055–1096.

`dpss`

| `periodogram`

| `pwelch`

Was this topic helpful?