3D TDOA/AOA Localization in MIMO Passive Radar with Transmitter and Receiver Position Errors

This paper deals with the problem of determining the position of a single target from time diﬀerence of arrival (TDOA) and angle of arrival (AOA) measurements using multitransmitter multireceiver passive radar system with widely separated antennas. A practically motivating scenario where the transmitter and receiver positions are contaminated by errors is addressed. First, the reduction in localization accuracy due to the presence of transmitter and receiver position errors is derived through the Cram´er-Rao lower bound (CRLB) analysis. Then, a novel algebraic localization algorithm based on weighted least squares minimization is proposed that takes the transmitter and receiver position errors into consideration to reduce the estimation error. The proposed solution is shown theoretically to reach the CRLB even when the transmitter and receiver positions have errors. Simulation results also verify the theoretical developments and the performance improvement of the proposed solution over existing algorithms.


Introduction
Passive radar detects and tracks potential targets by exploiting noncooperative transmitters as their sources of radar transmission [1][2][3][4]. Dispensing with the need for a dedicated transmitter makes passive radar inherently low cost and hence attractive for a broad range of applications. Recently, by taking the advantage of the spatial diversity from widely separated antenna configuration, passive radar with multiple separated receivers and noncooperative transmitters, also known as multiple-input multiple-output (MIMO) passive radar, has received growing attention due to its enormous potential in improving the detection and localization performances [5][6][7][8].
Target localization is one of the salient issues in the MIMO passive radar field. e time difference of arrival (TDOA) measurement, which usually comes from the crosscorrelation (CC) processing between the reference signal and the reflected target echo [9], is a common measurement used to determine the target position. And over the years, many localization methods have been developed based on TDOA measurements [10][11][12][13][14]. However, the TDOA-based localization with terrestrial transmitters and receivers suffers from poor accuracy in estimating the target height [11] and is overly dependent on the TDOA measurement accuracy. To overcome the fundamental flaw of TDOA-based localization, as suggested in [11], angle of arrival (AOA) measurement of the reflected target echo, which can be determined by subspace-based estimators [15], can be jointly utilized [16]. However, unlike TD-based localization which has been extensively studied [10][11][12][13][14], the hybrid TDOA/AOA localization is potentially more challenging due to the higher nonlinearity between the desired target position and TDOA/AOA measurements, and less effort has been devoted to hybrid TDOA/AOA localization.
Recently, borrowing the two-stage weighted least squares (TSWLS) idea originally proposed for radiation source localization [17], A. Noroozi1 et al. [18] developed a TSWLS algebraic solution for target localization with a MIMO passive radar using TDOA and AOA measurements. e performance analysis indicates that it can achieve the Cramér-Rao lower bound (CRLB) under mild noise conditions. By using a different way to linearize the TDOA and AOA measurement equations, R. Amiri et al. [19] explored a different algebraic localization algorithm based on weighted least squares (WLS) minimization, which is also shown theoretically and numerically to achieve the CRLB. Unlike Noroozi1's method which requires two WLS stages, Amiri's method determines the target position in only one WLS stage. Nevertheless, the above studies are based on the ideal assumption that the transmitter and receiver positions are exactly known, which is certainly not practical. Actually, the transmitter and receiver positions need to be estimated before the localization of an unknown target can be achieved, and the transmitter and receiver positions cannot be precisely known, especially when the antennas are mounted on moving platforms [20][21][22], GPS signal is sheltered, or the transmitters are extremely noncooperative (like the hostile radar radiation whose position could usually only roughly determined by electronic reconnaissance techniques [23]). e performance degradation created by transmitter and receiver position errors has been known for a while in the TDOA-based localization with MIMO passive radar [24]. And it was shown that the transmitter and receiver position errors can significantly degrade the localization performance of MIMO passive radars. Consequently, the errors in transmitter and receiver positions need to be taken into consideration in practical applications during the design of localization algorithms in MIMO passive radars.
In this paper, we address the target localization problem from TDOA and AOA measurements in the presence of transmitter and receiver position errors. We evaluate how much degradation the target localization accuracy is with respect to the amount of transmitter and receiver position errors by deriving the CRLB in the presence of transmitter and receiver position errors and examining the increase in CRLB due to the transmitter and receiver position errors. en, a novel closed-form solution is proposed for the localization problem to reduce the performance degradation created by the transmitter and receiver position errors. e proposed solution is shown analytically, under some mild approximations, to reach the CRLB, even in the presence of transmitter and receiver position errors. Some numerical simulations will be conducted to support the theoretical development of the proposed solution.

Notations.
is paper involves numerous symbols. By convention, uppercase and lowercase bold fonts denote matrices and vectors, respectively. e operations (·) T and (·) − 1 represent transpose and inverse, respectively. I k is the k-by-k identity matrix, and 0 k is a k-by-1 vector with all elements equal to zero. diag a { } is a diagonal matrix with the elements of a on the main diagonal. E(·) − 1 , tr(·) − 1 , and || · || stand for the statistical expectation, trace operation, and 2norm, respectively. e superscript (·) o is the true value of the noisy vector.  e paper is organized as follows. Section 2 presents the localization scenario and introduces the symbols involved. Section 3 evaluates the CRLB in the presence of transmitter and receiver position errors. Section 4 presents a novel proposed algebraic solution as well as a theoretical accuracy analysis. Section 5 contains the simulation results to verify the localization performance of the proposed solution, and Section 6 is the conclusion.

Problem Formulation
Consider a MIMO passive radar like the one illustrated in Figure 1. is MIMO passive radar system is comprised of N geographically separated receivers with positions be unknown position coordinates of the target. e transmitters send out a set of waveforms, which are then reflected by potential targets. e receivers sense the direct path signals from the transmitters as well as the reflected signals from the target. After some processing, each receiver extracts TDOAs and one AOA pair (an elevation angle and an azimuth angle), which are then sent to the fusion center to determine the target position.
In a practical localization scenario, the actual positions of all transmitters and receivers are not known, and the available positions are as follows: where s t,m and s r,n represent the available (noisy) positions of transmitter m and receiver n, and Δs t,m and Δs r,n represent the corresponding position errors. For notation simplicity, we put the positions of the transmitters and receivers into a 3(M + N)-by-1 column vector as where the covariance matrix Q β can be determined from theoretical analysis and actual measurement [25].

Mathematical Problems in Engineering
Based on the above geometry, the range between transmitter m and the target is R o t,m � ‖u o − s o t,m ‖, and the baseline distance between transmitter m and receiver n is R o t,m,r,n � ‖s o t,m − s o r,n ‖. us, the time difference between the direct path signal from transmitter m and the corresponding reflected signal arriving at receiver n can be given by e true AOA pair for receiver n, that is, the elevation angle denoted by θ o n and the azimuth angle denoted by φ o n , is, respectively, given by Considering the unavoidable measurement noises in reality, the erroneous version of the TDOA and AOA measurements can be expressed as where τ m,n , θ n , φ n represent the noisy TDOA and AOA measurement, and Δτ m,n , Δθ n , Δφ n are the corresponding measurement noises. For easier manipulation, we put these TDOA and AOA measurements into an (MN+2N)-by-1 column vector form as where ] T stands for the noisy TDOA and AOA measurement vector, T stands for the true TDOA and AOA vector, and Δφ N ] T stands for the TDOA and AOA measurement noise vector assumed to be zero-mean Gaussian with covariance matrix.
where the covariance matrix Q α can be determined from the specific signal conditions [26].
In this work, we are interested in identifying the unknown target position as accurately as possible, using the noisy transmitter/receiver positions and the TDOA/AOA measurements, as well as their error statistical characteristics. Nevertheless, this is a potentially challenging task since the desired target position is highly nonlinear with respect to the TDOA and AOA measurements.

CRLB Analysis
CRLB traces out the lowest possible variance of unbiased estimators and is often used as a benchmark for performance evaluation. In this section, we shall characterize the influence of transmitter and receiver position errors on the localization accuracy by deriving a novel CRLB for target position estimation with transmitter/receiver position errors and comparing it with the one without transmitter/receiver position error.

CRLB with Transmitter and Receiver Position Errors.
In addition to TDOA/AOA measurement noises, the presence of transmitter and receiver position errors is included. erefore, the unknown parameter vector for the CRLB evaluation is η o � [(u o ) T , (β o ) T ] T , and the observation vector is z � [α T , β T ] T . As described in Section 2, the TDOA/AOA measurement noise vector Δα and transmitter/ receiver position error vector Δβ are independent of each other, and zero-mean Gaussian with covariance matrices Q α and Q β , respectively. Hence, the logarithm of the joint probability density function of z under η o is where κ is a constant independent of η o . According to the structure of η o , the Fisher information matrix (FIM) will be a (3+3M+3 N)-by-(3+3M+3 N) matrix given by and the blocks X, Y, and Z are given by where zα o /zu o and zα o /zβ o are the partial derivatives of the parametric form of α o with respect to u o and β o , respectively, with their elements given by Mathematical Problems in Engineering for m � 1, 2, . . . , M, n � 1, 2, . . . , N, and zeros elsewhere. By definition, the CRLB of η o , denoted by CRLB(η o ), is equal to the inverse of the FIM, that is, Substituting (14) in (16) and then invoking the partitioned matrix inversion formula and the matrix inversion lemma [27,28], we have Comparing (17) with the CRLB in [18,19] indicates that X − 1 is the CRLB when there is no transmitter and receiver position errors, and is the change in CRLB brought by the transmitter and receiver position errors.

Performance Degradation from Transmitter and Receiver
Position Errors. We shall characterize analytically the performance degradation in the presence of transmitter and receiver position errors by establishing the positive definiteness of the second term in (17), that is, has a full column rank of 3. Moreover, it can be inferred from the expression of (zα o /zβ o ) that (zα o /zβ o ) T has a full column rank of (MN+2N). erefore, utilizing the expression of Y in (11), Y T has a full column rank.
Next, we proceed to prove Z − Y T X − 1 Y is positive definite. From the expressions of Y and Z in (11), we obtain (15) after some mathematical manipulations where Q −1 β is an obvious positive definite matrix. By using the expression of X in (11) and performing the Cholesky decomposition on Q −1 α , that is, Q −1 α � L α L T α , we can rewrite the first term on the right side of (15) as (15) is positive definite as long as Synthesizing the above results, that is, Y T has a full column rank, and at is to say, we can obtain (17) from (14) that where us, it can be concluded that the presence of transmitter and receiver position errors indeed deteriorates the target localization accuracy, at least at the CRLB level.

Algorithm Development and Analysis
e degradation of localization accuracy in the presence of transmitter and receiver position errors has been shown in Section 3 through the derivation and analysis of the CRLB. In what follows, to minimize the influence of transmitter and receiver position errors on target localization accuracy, we will proceed to design a novel algebraic localization algorithm for the aforementioned practical localization scenario. After that, a theoretical analysis will be performed to show that the proposed solution achieves the CRLB when satisfying some mild conditions.

Localization Algorithm.
e proposed solution is derived based on transforming the nonlinear TDOA and AOA equations into linear ones, from which the target position can be estimated using a simple WLS minimization.
To achieve this, we first rearrange the TDOA equation in (5) as Squaring both sides of (18), and then rearranging it as Since the nuisance parameter R o r,n are nonlinearly related to the target position u o , (19) is pseudo-linear with respect to u o . In order to thoroughly linearize (19), rewrite (6) as Putting (20) into (19) yields To linearize the AOA equations, rewrite (6) and (7)
for m � 1, 2, . . . , M, n � 1, 2, . . . , N. e error vector Δh can be expressed as R t,m,r,n − 2 cτ m,n + R t,m,r,n z r,n − z t,m R t,m,r,n + 2z r,n − 2 cτ m,n + R t,m,r,n sin φ n , for m � 1, 2, . . . , M, n � 1, 2, . . . , N, and zeros elsewhere. Note that (25) is a linear set of equations with respect to the target position u o , from which the target position u o can be estimated using the WLS minimization. e WLS solution of (25), which minimizes the cost function E(Δh T WΔh) with respect to the target position u o , is given by where W is the weighting matrix, and its optimal choice that gives the minimum parameter variance is given by However, as suggested in (30) Ignoring the second-and higher-order error terms in (31) and taking expectation result in that E(Δu) � 0 3 , that is, the target position estimate u in (29) is asymptotically unbiased. By multiplying (34) by its transpose and taking expectation, we have the covariance matrix of u as (32)  (14). By inserting (33) into (35), we have after some mathematical manipulations where with Comparing (36) with the CRLB in (17), we observe that CRLB(u o ) and cov(u) have the same structural form. Assuming the transmitter/receiver position errors and the TDOA/AOA measurement noises are sufficiently small, we can derive after tedious algebraic manipulations that Based on this, it can be inferred that In other words, the covariance matrix of the proposed solution accomplishes the CRLB given sufficiently small transmitter/receiver position errors and TDOA/AOA measurement noises.

Numerical Examples
is section contains some numerical simulations to evaluate the performance of the proposed solution. e localization scenario is set as shown in Figure 2 Table 1. MN � 24 TDOAs and N � 6 AOA pairs are extracted to determine the target position.
In order to simulate a practical localization scenario, zero-mean Gaussian noises with known covariance matrices Q α � di ag σ 2 t I MN , σ 2 a I 2N and Q β � σ 2 s I 3(M+N) are added to actual true TDOAs/AOAs and transmitter/receiver positions, respectively. σ t represents TDOA measurement noise level, σ a represents AOA measurement noise level, and σ s represents the transmitter/receiver position error level. e localization accuracy is assessed using root mean squares error (RMSE) defined as where u l is the estimation of u o at the lth Monte Carlo trial, and the number of Monte Carlo trials for each simulation is set as L � 5000.

CRLB Comparison.
In this section, in order to evaluate how sensitive the target localization accuracy is with respect to the transmitter and receiver position errors, we compare the CRLBs with and without transmitter/receiver position errors, under different TDOA/AOA measurement noise and transmitter/receiver position error levels. e comparison results are exhibited in Figure 3. Figure 3(a) presents the CRLBs with and without transmitter/receiver position error as TDOA measurement noise level σ t changes from 0.01us to 10us while σ a � 0.1°and σ s � 50m. It is readily seen that the CRLB with transmitter and receiver position errors is generally above the one without and with the former larger than the latter about one order of magnitude. Figure 3(b) compares the traces of the two CRLBs as AOA measurement noise level σ a changes from 0.01°to 10°while σ t � 0.1us and σ s � 50m. When the AOA measurement noise is small, the localization accuracy with transmitter and receiver position errors level is apparently lower than that without. As the AOA measurement noise level increases, the gap between the two CRLBs is becoming smaller and smaller. When the AOA measurement noise level is sufficiently large, the two CRLBs are nearly the same. e CRLB comparison versus the transmitter/receiver position error level σ s is plotted in Figure 3(c) where the TDOA/AOA measurement noise level is fixed at σ t � 0.1us and σ a � 0.1°. It can be seen that the difference between the two CRLBs becomes greater and greater with the increase of the transmitter/receiver position error level σ s . At a transmitter/receiver location error level of σ s � 100m which is not rare in practical applications, the localization CRLB with transmitter and receiver position errors is ten times as the localization CRLB without.
is reconfirms the significance of including the transmitter and receiver position errors in the design of the localization algorithm for MIMO passive radar.

Performance
Comparison. Now, we proceed to evaluate the localization RMSE of the proposed solution by comparing with two typical algorithms, that is, Noroozi1's method in [18] and Amiri's method in [ [19]], under different measurement noise and transmitter/receiver position error levels. In order to achieve a more comprehensive insight on the performance of the proposed solution, we consider two cases, that is, an ideal case where the transmitter and receiver position errors are negligibly small and a nonideal case where the transmitter and receiver position errors are significant. We first address the ideal case. e comparison results are presented in Figure 4.     Figure 4 depicts the localization RMSE curves of the algorithms for different TDOA/AOA measurement noise levels in the ideal case. As expected, when there are no transmitter and receiver position errors, three algorithms generally perform comparably on the whole. At small TDOA/AOA measurement noise levels, the RMSE curves of the three algorithms match the CRLB very well. With the increase of TDOA/AOA measurement noise levels, the RMSE curves rise correspondingly and deviate gradually from the CRLB. And the deviation from the CRLB, known as the thresholding phenomenon, is owing to the discarded second and higher-order error terms in the derivation of the algorithms, which is invalid for large measurement noise levels. However, it is worth noting that, after deviating from CRLB, the localization RMSE of Noroozi's method is slightly higher than that of Amiri's method and the proposed solution. e reason may be that Noroozi's method requires two WLS stages while Amiri's method and the proposed solution require only one WLS stage, which means Noroozi's method has to discard more second and higher-order error terms.
Next, we consider a more practical case, where the transmitter and receiver position errors cannot be ignored. e results are given in Figure 5. is is because both Noroozi's method and Amiri's method assume transmitter and receiver positions are accurate but in fact have errors. By contrast, the proposed solution provides better localization performance, with its RMSE reaching the CRLB at small TDOA/AOA measurement noise levels. Although the RMSE of the proposed solution starts to deviate from the CRLB when the TDOA/DOA measurement noises are large, it is still much smaller than that of Noroozi's method and Amiri's method. Figure 5(c) presents the RMSE comparison of the algorithms as transmitter and receiver position errors level σ s varies from 1m to 1000m, and the TDOA/AOA measurement noise level is fixed at σ t � 0.1us and σ a � 0.1°. Within expectation, when the transmitter and receiver position errors are sufficiently small, the RMSE curves of the three algorithms are very close. However, with the increase CRLB without transmitter and receiver position error CRLB with transmitter and receiver position error  of transmitter and receiver position errors level, the performance deterioration brought by the presence of transmitter and receiver position errors becomes more and more obvious for Noroozi's method and Amiri's method. us, the proposed solution performs much better than Noroozi's method and Amiri's method over large transmitter and receiver position errors regions. is verifies again that the presence of transmitter and receiver location errors can deteriorate the localization accuracy and the necessity of taking the transmitter and receiver position errors into consideration in the design of localization algorithm for MIMO passive radar systems.

GDOP Analysis.
In order to assess the effect of varying the target position on the proposed solution, we use the contour plots of the geometrical dilution of precision (GDOP) values defined as GDOP(u) � ���������� σ 2 x + σ 2 y + σ 2 z where σ x , σ y , and σ z are the variances of localization on the x, y, and z axes, respectively. In all GDOP simulations, the TDOA/AOA measurement noise level and transmitter/receiver position error level are set to be σ t � 0.1us, σ a � 0.1°, σ s � 50m, respectively. It should be pointed out that the setting of σ t � 0.1us, σ a � 0.1°, and σ s � 50m is of practical value under the circumstance of the target echoes being contaminated by moderate noises and GPS signal being sheltered or jammed. e contour plots of the GDOP values in an 4km × 4km area around the MIMO passive radar are provided in Figure 6. Figure 6(a) presents the GDOP contours at z o � 50m, from which we observe that the GDOP plots have a relatively symmetric form corresponding to the transmitter/receiver layout. As the target position moves towards the receivers, the GDOP values decrease, so that the minimum level occurs near the receiver positions. Figure 6(b) presents the GDOP contours at z o � 500m. Comparing Figure 6(a) with Figure 6(b), we see that by increasing the z-coordinate of the target, the GDOP contours become more symmetric. Furthermore, when the target x-coordinate and y-coordinate are the same, the larger the target z-coordinate, the larger the GDOP values, and the smaller the influence of the transmitter/receiver position on target localization accuracy. is simulation shows that when the target z-coordinate becomes larger, the GDOP also becomes larger and the influence of transmitter/receiver position errors decreases. is can be useful in our application.

Computation Complexity Comparison.
e computational complexity is also an important index for performance evaluation. In what follows, to evaluate the proposed solution in terms of computational complexity, we count the average running time of the algorithms from 5000  Table 2.
As presented in Table 2, the time cost of Noroozi1's method is almost twice higher than that of Amiri's method.
is is because Noroozi1's method needs two WLS stages while Amiri's method identifies the target position in only one WLS stage. e proposed solution incurs the highest time cost among the algorithms, approximately 3 times higher than Amiri's method. is is not surprising since the proposed solution took the transmitter and receiver position errors into account while Noroozi1 and Amiri's methods do not. at is to say, it is at the expense of the higher computation cost that the proposed solution achieves higher localization accuracy. However, in view of the significant performance enhancement, the increased computation cost is worthy and acceptable.

Conclusions
Target localization from time difference of arrival (TDOA) and angle of arrival (AOA) measurements using multitransmitter multireceiver passive radar system requires very precise knowledge of the transmitter and receiver positions. A small error in the transmitter and receiver positions may result in a significant degradation in target localization accuracy. Hence, this paper addresses a practically motivated scenario, where the transmitter and receiver positions are not known perfectly and only the nominal values are available for processing. To minimize the influence of transmitter and receiver position errors on target localization accuracy, we proposed a novel algebraic solution to improve the target position estimate. By taking the transmitter and receiver position errors into consideration in the measurement model, the proposed solution can achieve the CRLB, no matter whether there exist transmitter and receiver position errors or not. Both theoretical performance analysis and numerical simulations are performed to demonstrate the superiority of the proposed solution over existing algorithms.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.