Implementation Issues with Subspace Based Estimation Algorithms
MUSIC and ESPRIT are considered the two main subspace based estimation algorithms. They both require the theoretical or known correlation matrices. However, in reality, the problem starts often from a sequence of measured data. Therefore, three issues must be resolved from the given data sequence for subspace based estimation algorithms.
Correlation matrix and singular value decomposition
It is sufficient to simply use the unnormalized estimate of the correlation matrix,
Define as the with the first column removed. The principal component can be obtained by the following,
Where is the first column vector of , and is the rank pseudoinverse of the data matrix, or,
The above procedure simplifies the singular value decomposition greatly.
Number of Signals
In most common situations the number of signals present is not known. The second issue is how to decide the number of signals from the given sequence.
When the signal power is much larger than the noise power, estimating the number of signals is usually not a problem. Akaike information criterion (AIC) can often help.
AIC is defined in terms of the likelihood function for the observations,
Where is the number of free parameters represented by . If the observation vector is Gaussian with zero mean, AIC can be further simplified to compute the ratio of geometric over arithmetic means.
Spectrum is more easily computed utilizing FFT, then pseudospectrum should also be computed in similar ways. For MUSIC, the spectrum can be written
the denominator can also be simplified with FFT.