Main Content

corrmtx

Data matrix for autocorrelation matrix estimation

collapse all in page

Syntax

H = corrmtx(x,m)
H = corrmtx(x,m,method)
[H,r] = corrmtx(___)

Description

H = corrmtx(x,m) returns an (n+m)-by-(m+1) rectangular Toeplitz matrix H such that H†H is a biased estimate of the autocorrelation matrix for the input vector x. n is the length of x, m is the prediction model order, and H† is the conjugate transpose of H.

H = corrmtx(x,m,method) computes the matrix H according to the method specified by method.

example

[H,r] = corrmtx(___) also returns the (m + 1)-by-(m + 1) autocorrelation matrix estimate r, computed as H†H, for any of the previous syntaxes.

Examples

collapse all

Open Live Script

Generate a signal composed of three complex exponentials embedded in white Gaussian noise. Compute the data and autocorrelation matrices using the 'modified' method.

n = 0:99;
s = exp(i*pi/2*n)+2*exp(i*pi/4*n)+exp(i*pi/3*n)+randn(1,100);
m = 12;
[X,R] = corrmtx(s,m,'modified');

Plot the real and imaginary parts of the autocorrelation matrix.

[A,B] = ndgrid(1:m+1);
subplot(2,1,1)
plot3(A,B,real(R))
title('Re(R)')
subplot(2,1,2)
plot3(A,B,imag(R))
title('Im(R)')

Figure contains 2 axes objects. Axes object 1 with title Re(R) contains 13 objects of type line. Axes object 2 with title Im(R) contains 13 objects of type line.

Input Arguments

collapse all

Input data, specified as a vector.

Prediction model order, specified as a positive real integer.

Matrix computation method, specified as 'autocorrelation', 'prewindowed', 'postwindowed', 'covariance' or 'modified'.

  • 'autocorrelation': (default) H is the (n + m)-by-(m + 1) rectangular Toeplitz matrix that generates an autocorrelation estimate for the length-n data vector x, derived using prewindowed and postwindowed data, based on an mth-order prediction model. The matrix can be used to perform autoregressive parameter estimation using the Yule-Walker method. For more details, see aryule.

  • 'prewindowed': H is the n-by-(m + 1) rectangular Toeplitz matrix that generates an autocorrelation estimate for the length-n data vector x, derived using prewindowed data, based on an mth-order prediction model.

  • 'postwindowed': H is the n-by-(m + 1) rectangular Toeplitz matrix that generates an autocorrelation estimate for the length-n data vector x, derived using postwindowed data, based on an mth-order prediction model.

  • 'covariance': H is the (n – m)-by-(m + 1) rectangular Toeplitz matrix that generates an autocorrelation estimate for the length-n data vector x, derived using nonwindowed data, based on an mth-order prediction model. The matrix can be used to perform autoregressive parameter estimation using the covariance method. For more details, see arcov.

  • 'modified': H is the 2(n – m)-by-(m + 1) modified rectangular Toeplitz matrix that generates an autocorrelation estimate for the length-n data vector x, derived using forward and backward prediction error estimates, based on an mth-order prediction model. The matrix can be used to perform autoregressive parameter estimation using the modified covariance method. For more details, see armcov.

Output Arguments

collapse all

Data matrix, returned for autocorrelation matrix estimation. The size of H depends on the matrix computation method specified in method.

Biased autocorrelation matrix, returned as a (m + 1)-by-(m + 1) rectangular Toeplitz matrix.

Algorithms

The Toeplitz data matrix computed by corrmtx depends on the method you select. The matrix determined by the autocorrelation (default) method is:

H=1n[x(1)0⋯00x(2)x(1)⋯00x(3)x(2)⋯00⋮⋮⋱⋮⋮x(m)x(m−1)⋯x(1)0x(m+1)x(m)⋯x(2)x(1)x(m+2)x(m+1)⋯x(3)x(2)⋮⋮⋱⋮⋮x(n−1)x(n−2)⋯x(n−m)x(n−m−1)x(n)x(n−1)⋯x(n−m+1)x(n−m)0x(n)⋯x(n−m+2)x(n−m+1)⋮⋮⋰⋮⋮00⋯x(n−1)x(n−2)00⋯x(n)x(n−1)00⋯0x(n)].

In the matrix, m is the same as the input argument m to corrmtx and n is length(x). Variations of this matrix are used to return the output H of corrmtx for each method:

  • 'autocorrelation' — (default) H = H.

  • 'prewindowed' — H is the n-by-(m + 1) submatrix of H whose first row is [x(1) … 0] and whose last row is [x(n) … x(n – m)].

  • 'postwindowed' — H is the n-by-(m + 1) submatrix of H whose first row is [x(m + 1) … x(1)] and whose last row is [0 … x(n)].

  • 'covariance' — H is the (n – m)-by-(m + 1) submatrix of H whose first row is [x(m + 1) … x(1)] and whose last row is [x(n) … x(n – m)].

  • 'modified' — H is the 2(n – m)-by-(m + 1) matrix Hmod defined by

    Hmod=12(n−m)[x(m+1)⋯x(1)⋮⋰⋮x(n)⋯x(n−m)x∗(1)⋯x∗(m+1)⋮⋱⋮x∗(n−m)⋯x∗(n)].

References

[1] Marple, S. Lawrence. Digital Spectral Analysis: With Applications. Prentice-Hall Signal Processing Series. Englewood Cliffs, N.J: Prentice-Hall, 1987.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using MATLAB® Coder™.

Version History

Introduced before R2006a

See Also

| | | |