Application of the Total Least Square ESPRITMethod to Estimation of Angular Coordinates of Moving Objects

The TLS ESPRIT method is investigated in application to estimation of angular coordinates (angles of arrival) of two moving objects at the presence of an external, relatively strong uncorrelated signal. As a radar antenna system, the 32-element uniform linear array (ULA) is used. Various computer simulations have been carried out in order to demonstrate good accuracy and high spatial resolution of the TLS ESPRIT method in the scenario outlined above. It is also shown that accuracy and angle resolution can be significantly increased by using the proposed preprocessing (beamforming). The most of simulation results, presented in a graphical form, have been compared to the corresponding equivalent results obtained by using the ESPRIT method and conventional amplitude monopulse method aided by the coherent Doppler filtration.


Introduction
One of the important problems of today's radar technology is the detection and estimation of angular coordinates of fast moving objects, the radar echoes of which are received on a background of strong uncorrelated interferences and signals reflected from other objects located nearby.It is obvious that in such situations the conventional monopulse methods are rather noneffective [1,2].This noneffectiveness results from the essence of them, since they enable the correct estimation of angular coordinates merely in situations when just one sufficiently strong useful signal is being received in a specific moment.Unfortunately, in the real electromagnetic environment, besides the useful signal (reflected from a detected object), different spurious signals are also received by the radar antenna unit.For instance, these spurious signals may be reflected from the earth, weather precipitation, or other masking objects flying around.It follows from the literature that problems of this kind can be effectively solved by modern radar systems using the angular superresolution methods.Thus, in this paper, one of them, namely, the Total Least Square ESPRIT (TLS ESPRIT), has been used for this aim [3][4][5][6].The corresponding theory is presented in Section 2. In the next, that is, Section 3, it is shown how this method can be effectively applied to solve different variants of the general problem formulated above.A novel component of the proposed approach is an algorithm of preprocessing the samples of signals received simultaneously by a linear uniform antenna array.This algorithm creates a new "virtual" linear array with smaller number of radiating/receiving elements.It is important that signal-to-noise ratio evaluated at each of its virtual elements is greater, thanks to integration of signals from consecutive radiating/receiving elements of the original array.Consequently, the angular coordinates of the detected object can be evaluated more accurately.The angular resolution of the applied method is also increased.A short discussion on the results of the corresponding computer simulations is presented in Sections 3, 4, 5, 6, and 7.

TLS ESPRIT Method
The TLS ESPRIT method can be treated as an improved version of the ESPRIT method, the theoretical basis of which is described in many publications, for example, [5][6][7].To recall the essence of the TLS ESPRIT method, let us assume that the signals arriving from different directions are received 2 International Journal of Antennas and Propagation by two identical subarrays, X and Y , mutually displaced by distance Δ, Figure 1.
The resultant signals received by the arrays are the sum of the K signals s k (t), k = 1, . . ., K being received from specific directions determined by the angles θ k .The vectors of signals received by the first subarray (X) and second subarray (Y ), as shown in Figure 1, are described by the following equations, respectively: where where ω 0 is the angular frequency of carrier wave of the signal received by the antenna array, λ is the wavelength of the carrier wave of the signal received by the antenna array, a i (θ k ) is the parameter representing the directional gain of the ith radiating element, and τ i (θ k ) is the delay of the signal reaching the ith radiating element, having s(n) as the vector of samples received by the antenna array, s(n) = [s 1 (n), s 2 (n), . . ., s K (n)] T , u(n) as the vector of noise samples in the radiating elements of the first subarray, and v(n) as the vector of noise samples in the radiating elements of the second subarray.The vector of samples z(n), which is composed of samples of signals received from both subarrays, is defined as follows: where For the vector z(n), the spatial correlation matrix is formulated where  1).All generalized eigenvectors v i satisfy the equation where λ i denotes the ith eigenvalue of matrix R Z .On the basis of eigenvectors corresponding to the K biggest eigenvalues, the vector matrix of the signal space is created The vectors included in matrix (8) can be presented as a combination of the directional vectors: since the directional vectors of received signals (written into columns of matrix A) span over the same signal space.Like in the case of the matrix of directional vectors A, the eigenvector matrix V S has the structure that allows us to determine the phase shift between the signals received by subarrays X and Y , see Figure 1 That means that for the matrix fulfilling the following equation, [6,8]: By using the matrix Ψ = −F X F Y −1 , (12) can be written in the form ATΨ = AΦT, which results in Here, it should be noticed that (13) has similar form as equation Ψ = VΛV H , which defines eigenvalues and eigenvectors of matrix Ψ.Therefore, columns of matrix T −1 are equal to eigenvectors of matrix Ψ, and values of diagonal elements of matrix Φ are equal to related eigenvalues of matrix Ψ.When the values of the phase shifts included in the diagonal elements of Φ are known, one can determine the directions represented by the angles θ k , from which the signals s k (t) are received.
Unfortunately, the method described above requires a very accurate estimate of the correlation matrix R Z , for which the following equation should be satisfied: In practice, ( 14) is not satisfied precisely due to an inaccurate estimate of the correlation matrix R Z determined on basis of limited number of measurements of the vector z(n).The estimate of the correlation (covariance) matrix is where N ret denotes the number of pulse returns (realizations) employed to estimate the values of the spatial correlation matrix, g is the index of range gate, and G is the number of range gates obtained after sampling of the received pulse return.In such case, the generalized eigenvectors, obtained as the result of the decomposition of matrices R Z and Σ, are inaccurate.The sought matrix Ψ, for which the error of mapping the matrix V X into matrix V Y is the least, can be found using the criterion of the total least square error.
According to this criterion, matrix F is evaluated, for which the value of square of the Euclidean (Frobenius) norm of the product of the matrices V XY and F, that is is the least.Additionally, a limitation as described by the relationship below is imposed on the value of the matrix F The solution of the aforediscussed problem is the matrix F of dimension 2K × K composed of K columns, which are equal to eigenvectors corresponding to K least eigenvalues of the matrix product V H XY V XY .The eigenvalues of the matrix Ψ correspond to the sought values of the elements lying on the diagonal of the matrix Φ. Knowing the values of diagonal elements of matrix Φ, values of angles θ k can be found.To summarize, the TLS ESPRIT method can be implemented in the following way.
Step 1. Determine the estimate of the spatial correlation matrix R Z , based on the measurement of the signal vector z(n).
Step 2. Decompose the matrices R Z and Σ with respect to their eigenvectors and eigenvalues, see ( 6) and (7), which satisfy the following condition: where V is a matrix of eigenvectors and Λ is a diagonal matrix of eigenvalues.
Step 3. Determine the number of signal sources K.
Step 4. Out of the matrix V, select vectors corresponding to K biggest eigenvalues and create a matrix V s of them.Step 5. Divide the matrix V s into matrices V X and V Y and create a matrix V XY of them Step 6. Decompose the matrix product V H XY V XY with respect to its eigenvectors and eigenvalues, which satisfy the equation Step 7. Split the eigenvector matrix E into four submatrices of dimensions K × K each Step 8. Based on the matrices E 12 and E 22 , the columns of which correspond to the K smallest eigenvalues, determine a matrix Ψ by using the transformation Step 9. Find the eigenvalues φ k of matrix Ψ.
Step 10.Based on complex coefficients φ k , calculate the angles, which determine the directions of arrival of signals s k (t), see Figure 1.

Application of the TLS ESPRIT Method to Evaluate Angular Coordinates of Two Moving Objects
In this section, the TLS ESPRIT method has been used to solve the problem of estimation of angular coordinates of two objects, moving with velocity v = 100 m/s, Figure 2. Let us assume that both objects bear at α = −45 • north and the azimuth angle of the first of them is θ 1 = −14.1 • .The azimuth angle θ 2 of the second one is changing in the range of −14.1 • ≤ θ 2 ≤ −10.1 • , although without changing its bearing and velocity.Gradually the angular distance between these two objects, is changing in the range of Δθ = 0 • ÷ 4 • .In order to estimate angular coordinates of these objects a radar system with electronically controlled antenna (nonrotating mechanically) has been used.This linear array is composed of N = 32 radiating/receiving elements.The distance between two adjacent elements is equal to Δ = 0.7λ, where λ is a length of the received electromagnetic wave.It is also assumed that the radar system performs signal processing using the MTD (moving target detection) algorithm, [2,9].The processing is performed on the signal composed of 24 pulses grouped in three 8-pulse bursts which differ by their pulse repetition interval, Figure 3.The number of pulses emitted by the radar in given direction is smaller than the number of elements of the array.In order to reduce an influence of inaccurate estimate of correlation matrix R Z on estimation of angular coordinates, the appropriate "diagonal loading" has been used ( 22) where σ 2 is the variance of noise in elements of the array and I is the identity matrix.
The signals received by the antenna array are used to form the sum and difference antenna patterns.Thus, the angular coordinates of objects being detected can be evaluated by means of amplitude monopulse method.
The samples s(n) of echoes received by the sum beam Σ are subjected to coherent integration with the use of several sets of weight coefficients.The so-understood integration is called Doppler filtering, because it allows the discrimination of signals, which originate from objects moving with different radial velocities with respect to the radar.The signal obtained as a result of the so-understood integration is described by the general relationship: where s doppl (g) is the sample of signal in a given Doppler filter, s(g) is the sample of the pulse return received by sum antenna beam, g is the index of the range gate which determines the distance between the object and the radar, G is the number of range gates, m is the index of a pulse within the pulse sequence called "burst", and w m is appropriately chosen weight coefficient.After detection, that is, after distinguishing the range gates in which the level of detected signal exceeds a decision threshold, the radar system performs an estimation of angular coordinates.Naturally, it can be performed by using the conventional amplitude monopulse method.Moreover, the samples of signals appearing at sum and difference antenna outputs can be additionally filtered by the same Doppler filter.The more efficient solution uses the ESPRIT and TLS ESPRIT methods which both use the complex samples of the signal received by the 32-element antenna array; the samples form the vector x(n).The samples are obtained at outputs of matched filters MF, that is, before the beamforming network, see Figure 4.An attempt to evaluate angular coordinates of two closely spaced objects by means of the conventional amplitude monopulse method was met with failure.This failure resulted from the fact that monopulse methods are not able to discriminate objects seen in the same antenna beam when these objects are located at a similar distance to the antenna.The angular resolution of angular coordinates of these objects is determined by 3 dB beam width of azimuth and elevation receive characteristics.Usually the second object, for instance, an aircraft causes errors of angular coordinates of the first one when it is seen within the same antenna beam.This effect can be partially suppressed using the additional Doppler filtering.It allows the attenuation of signals, which sustain another frequency shift due to Doppler effect.It is well known that a change of angular position of an object causes a change of the corresponding radial velocity and Doppler frequency shift related to it.Unfortunately, at the higher radial velocities, the normalized Doppler frequencies of echoes of these objects may be similar and cause errors of estimation of angular coordinates.This effect is illustrated in Figure 5. Thus, Doppler filtering can be used to increase signal-to-interference ratio when a radar receives an echo reflected from a single object together with interferences (weather clutter, ground clutter, etc.).However; it is not an appropriate approach to the resolution of two closely spaced objects.When these objects are moving with relatively small velocities (e.g., v = 100 m/s), their Doppler frequencies, normalized with respect to pulse repetition frequency (PRF) f d-norm = f d % F p , are different, but they may appear in 3 dB passband of the same Doppler filter.In the opposite situation, when two objects are moving with high speed (e.g., v = 1000 m/s) and PRF (e.g., F p = 1/T p = 250 Hz) is small in comparison to their Doppler frequencies, the normalized Doppler frequencies of corresponding echoes may be similar or equal.The errors caused by this effect are illustrated in Figure 5.
For comparison, Figure 6 presents the diagrams of the estimates of the angular coordinates of two closely spaced objects found using the TLS ESPRIT method.

The Preprocessing Preceding the TLS ESPRIT Method
The results illustrated in Figure 6 have been obtained using the TLS ESPRIT method for signal samples recorded at outputs of the 24-element equivalent (virtual) antenna array.This "virtual" antenna array has been created in the manner described below, known as preprocessing.The samples of signals x i (n), where 1 ≤ i ≤ 32, received by the real 32element antenna array, are used to compute the following samples: Samples (24) correspond to particular elements of the 24element antenna array being formed.According to Figure 1,  subarray X includes 23 "virtual" elements, at the outputs of which, samples x (e) k (n) (1 ≤ k ≤ 23) are evaluated.Similarly, samples x (e) k (n) 2 ≤ k ≤ 24 are evaluated at the outputs of 23 "virtual" elements of subarray Y .Such preprocessing scheme allows to improve the TLS ESPRIT method applied for weaker signals, at S/N noise < 6 dB.The signal to noise ratio S/N noise has been determined at outputs of the matched filters MF, see Figure 4.The signal to interference (jamming) ratio S/N jamm has been determined at inputs of matched filters (MF).It is justified by fact, that jamming signals may be correlated to emitted chirp modulation pulse.

Estimation of Angular Coordinates in the Presence of Interferences
An important advantage of the TLS ESPRIT method is its high immunity against interference.To illustrate this property, let us assume that a point-like interference source is located in the direction represented by angle θ jamm = 1.4 • .It has been assumed also that the surface power density of the jamming signal is 60 dB higher than the surface power density of both useful signals.In the presence of such interference, conventional radar systems (with the receive patterns nonadapting to the interference) will be jammed.An effective way of suppressing the influence of interference is to use the adaptive antenna with digitally formed receive patterns.Such receive patterns should have deep minima (maximum attenuation) at directions from which the strong interfering signals are received.The corresponding results of simulations are summarized in Table 1.

Reducing the Errors of Angular Coordinates by Using Bursts with Variable Pulse Repetition Intervals
The presence of the second object within the antenna beam causes errors of angular coordinates when both objects have similar or equal normalized Doppler frequencies.Fortunately, the TLS ESPRIT method under discussion is able to discriminate both objects and sources of interference when a variable pulse repetition interval is applied, see Figure 3.This property is confirmed by estimates of angles evaluated by using the TLS ESPRIT method.The corresponding results are illustrated in Figures 8, 10 been obtained by using the initial preprocessing described in Section 4. Simulations have been carried out for S/N noise = 0 dB, S/N noise = 6 dB, and S/N noise = 30 dB, respectively.For comparison, the equivalent results, obtained by means of the ESPRIT method, are illustrated in Figures 7, 9, and 11, [7].
From the analysis of the diagrams obtained by means of the ESPRIT and the TLS ESPRIT method, it results that TLS ESPRIT offers better accuracy and higher resolution.The corresponding values of mean square errors for angular coordinates are given in Table 1.
Results presented in Figures 5-12 and Table 1 have been obtained for bursts with variable pulse repetition intervals.For comparison, Figure 13 presents mean square errors for angular coordinates occurring when bursts with the same pulse repetition interval are used and when the pulse repetition frequency is small (T p = 4 ms) versus the Doppler frequency of signals reflected from objects moving with high velocity (v = 1000 m/s).The maximum of error observed in Figure 13 occurs when normalized Doppler frequencies of both objects are equal.The comparison of results obtained by the TLS ESPRIT for signal samples from the 32-element antenna array to corresponding results obtained for signal samples formulated for equivalent "virtual" 24-element antenna array states that the use of preprocessing increases efficiency of the TLS ESPRIT for weaker signals, for which the S/N noise < 6 dB.The results presented in Figures 14 and 15 justify this conclusion.
From the analysis performed above, results that proposed additional preprocessing significantly improve the effectiveness of the TLS ESPRIT method in application under consideration.

Conclusions
The results of numerous computer simulations illustrated in Figures 6, 8, 10, and 12 confirm the applicability of the TLS ESPRIT method to solve different variants of the problem under consideration.The comparison of these results with their correspondents obtained by using the ESPRIT method, see Figures 7, 9, and 11, indicates that the TLS ESPRIT method offers better angular resolution and more reliable estimation of angular coordinates of two moving objects, especially when weak echoes are being received in the presence of strong jamming signal.Another, also important, result of the performed simulations is the confirmation of the usefulness of the proposed signal preprocessing, the essence of which is expressed by formula (24).From the obtained results, see Figures 8,10,14, and 15, it follows that the proposed preprocessing contributes to increasing the accuracy of angular coordinates and consequently makes the TLS ESPRIT method more effective, naturally in the application under consideration.In other words, the preprocessing increases the effectiveness of the TLS ESPRIT method especially for the weaker echoes of useful signals.In the first approximation, this essential improvement can be explained in the following way: as a result of the addition of signals x i (n), see (24), the uncorrelated noise waveforms are summed.Thereby, in elements of antenna array, increasing of a signal-to-noise ratio is observed.In this paper, an emphasis has been put on the fact that the effective use of the TLS ESPRIT method requires an appropriate method of scanning.When scanning signals with identical pulse repetition intervals in all three bursts used, see Figure 3, significant errors may occur for some angle locations.This effect is illustrated by the simulation results shown in Figure 13.This parasitic effect can be reduced by changing pulse repetition interval in subsequent bursts, see Figure 3. Here, it deserves noting that bursts with variable pulse repetition intervals are widely used in order to avoid errors related to the so-called "blind speeds" [2,9].

Figure 1 :Figure 2 :
Figure 1: A uniform linear array antenna composed of N sensors and their (N − 1)-sensor subarrays X and Y .

Figure 3 :
Figure 3: A structure of a radar signal used in the MTD processing.

Figure 4 :
Figure 4: The functional scheme of an antenna unit including the digital beamformer.A/D: analog to digital converter, QPD: quadrature phase detector, MF: digital matching filter.

Figure 5 :
Figure 5: Histogram of polar coordinates evaluated by means of the amplitude monopulse method aided by the coherent Doppler filtration, v = 1000 m/s, T p = 4 ms, and S/N noise = 30 dB.

Figure 6 :Figure 7 :
Figure 6: Histogram of polar coordinates evaluated by using the TLS ESPRIT method for the plane velocity v = 1000 m/s, T p = 4 ms, and S/N noise = 30 dB.

Figure 8 :
Figure 8: Histogram of polar coordinates evaluated by using the TLS ESPRIT method at S/N noise = 0 dB and S/N jamm = −60 dB.

5 SFigure 9 :Figure 10 :
Figure 9: Histogram of polar coordinates evaluated by using the ESPRIT method at S/N noise = 6 dB and S/N jamm = −60 dB.

Figure 11 :Figure 12 :
Figure 11: Histogram of polar coordinates evaluated by using the ESPRIT method at S/N noise = 30 dB and S/N jamm = −60 dB.

Figure 13 :Figure 14 :
Figure 13: Histogram of mean square error of polar coordinates evaluated by using the TLS ESPRIT method at S/N noise = 10 dB, in absence of an external noise signal and a constant repetition time of radar pulses, T p = 4 ms, v = 1000 m/s.

Figure 15 :
Figure 15: Histogram of polar coordinates evaluated by means of the TLS ESPRIT method at S/N noise = 6 dB, S/N jamm = −60 dB, and without initial preprocessing.
1, 1, . . ., 1} is a unitary matrix of the (2N − 2)-th degree.The generalized eigenvectors v i , i = 1, . . ., 2N − 2 of the matrix R Z create a space, in which the directional vectors of signals s k (n), received by both N − 1 element antenna arrays lie (Figure

Table 1 :
The angular resolution and mean-square error of angular coordinates evaluated by means of the TLS ESPRIT method.