Note: This page has been translated by MathWorks. Click here to see

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

To view all translated materials including this page, select Country 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?