% subspace method
[thetas,residues]=sm(R,A,B,2);
arrowb(thetas,residues); hold on;
Ac=compan(eye(1,6));
Bc=eye(5,1);
That=dlsim_complex(Ac,Bc,y');
[th_esprit,r_esprit]=sm(That,Ac,Bc,2); % ESPRIT spectral lines
arrowg(th_esprit,r_esprit);
Subspace methodsWe start with the CarathodoryFejrPisarenko result for Toeplitz matrices. Given a positive definite Toeplitz matrix , let be the smallest number such that is singular and has rank , then has the formE where Accordingly, the power spectrum deomposes as where is the Dirac function. The decomposition corresponds to a set of white noise. See MA decomposition for a decomposition corresponds movingaverage (MA) noise. The above result can be generalized to the case of state covariances [1]. More specifically, let be the unique solution to the Lyapunov equation The matrix is the state covariance when the input is pure white noise. Let now be an arbitray state covariance, let be the smallest eigenvalues of the matrix pencil , and assme that the rank of is , then where is generalization of . The subspace spectral analysis methods rely on the singular value decomposition where is a unitary matrix and , . Partition where and are the first and the last columns of , respectively. Based on this decomposition, there are two ways we can proceed generalizing the MUSIC and ESPRIT methods, respectively [P. Stoica, R.L. Moses, 1997]. Noise subspace analysisThe columns of span the null space of , while the signal , is in the span of the columns of . So the nonnegative trigonometric polynomial has roots at . Given a sample state covariance matrix and an estimate on the number of signals , we let be the matrix of singular vectors of the smallest singular values, and be the corresponding trigonometric polynomial for . Two possible generalization of the MUSIC method are as follows.
Signal subspace analysisA signal subspace method relies on the fact that for the pair , and given above, there is a unique solution to the following matrix equation where is row vector and a matrix. The eigenvalues of are precisely for . If is a Toeplitz matrix and , given as in with a companion matrix, then we recover the ESPRIT result [P. Stoica, R.L. Moses, 1997]. The signal subspace estimation computed using sm.m, whereas music.m and esprit.m implement the MUSIC and ESPRIT methods, respectively. For the example discussed above, the estimated spectral lines are shown in the following figure. For an extensive comparison of the high resolution method sm.m with MUSIC and ESPRIT methods, see the example. % subspace method
[thetas,residues]=sm(R,A,B,2);
arrowb(thetas,residues); hold on;
Ac=compan(eye(1,6));
Bc=eye(5,1);
That=dlsim_complex(Ac,Bc,y');
[th_esprit,r_esprit]=sm(That,Ac,Bc,2); % ESPRIT spectral lines
arrowg(th_esprit,r_esprit);
