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

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

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

Autoregressive power spectral density estimate — modified covariance method

`pxx = pmcov(x,order)`

`pxx = pmcov(x,order,nfft)`

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

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

`[pxx,w] = pmcov(x,order,w)`

`[pxx,f] = pmcov(x,order,f,fs)`

`[___] = pmcov(x,order,___,freqrange)`

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

`pmcov(___)`

returns
the power spectral density estimate, `pxx`

= pmcov(`x`

,`order`

)`pxx`

, of
a discrete-time signal, `x`

, found using the modified
covariance method. 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`

. `pxx`

is
the distribution of power per unit frequency. The frequency is expressed
in units of rad/sample. `order`

is the order of
the autoregressive (AR) model used to produce the PSD estimate.

uses `pxx`

= pmcov(`x`

,`order`

,`nfft`

)`nfft`

points
in the discrete Fourier transform (DFT). For real `x`

, `pxx`

has
length (`nfft`

/2+1) if `nfft`

is
even, and (`nfft`

+1)/2 if `nfft`

is
odd. For complex–valued `x`

, `pxx`

always
has length `nfft`

. If you omit `nfft`

,
or specify it as empty, then `pmcov`

uses a default
DFT length of 256.

`[`

returns
a frequency vector, `pxx`

,`f`

] = pmcov(___,`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/second (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`

).

`[`

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

,`f`

] = pmcov(`x`

,`order`

,`f`

,`fs`

)`f`

. The vector, `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).

`[___,`

returns
the `pxxc`

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

)`probability`

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

.

`pmcov(___)`

with no output arguments
plots the AR PSD estimate in dB per unit frequency in the current
figure window.

Was this topic helpful?