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,
because otherwise the function interprets it as `nfft`

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