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

Authors: Donghua Huang; Yongsheng Zhao; Dexiu Hu; Yongjun Zhao
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

This paper deals with the problem of determining the position of a single target from time difference 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ér-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.


---

## Body

## 1. Introduction

Passive radar detects and tracks potential targets by exploiting noncooperative transmitters as their sources of radar transmission [1–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–8].Target localization is one of the salient issues in the MIMO passive radar field. The 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–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–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. The 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–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]). The 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. Then, a novel closed-form solution is proposed for the localization problem to reduce the performance degradation created by the transmitter and receiver position errors. The 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.


### 1.1. Notations

This paper involves numerous symbols. By convention, uppercase and lowercase bold fonts denote matrices and vectors, respectively. The operations·T and ·−1 represent transpose and inverse, respectively. Ik is the k-by-k identity matrix, and 0k is a k-by-1 vector with all elements equal to zero. diaga 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 2-norm, respectively. The superscript ·o is the true value of the noisy vector. ak and aj:k are indexing expressions for a column vector a, with ak representing the kth element of a and aj:k representing the jth to the kth elements of a. Ak,:, A:,k, Aj:k,i:l, Aj:k,:, A:,j:k are indexing expressions for a matrix A, with Ak,: representing the kth column of matrix A, A:,k representing the kth row of matrix A, Aj:k,i:l representing a submatrix including the elements of A indexed by j:k in the first dimension and i:l in the second dimension, Aj:k,: and A:,j:k include all subscripts in one of the dimensions but use the vector j:k to index in the other dimension.The paper is organized as follows. Section2 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.


## 1.1. Notations

This paper involves numerous symbols. By convention, uppercase and lowercase bold fonts denote matrices and vectors, respectively. The operations·T and ·−1 represent transpose and inverse, respectively. Ik is the k-by-k identity matrix, and 0k is a k-by-1 vector with all elements equal to zero. diaga 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 2-norm, respectively. The superscript ·o is the true value of the noisy vector. ak and aj:k are indexing expressions for a column vector a, with ak representing the kth element of a and aj:k representing the jth to the kth elements of a. Ak,:, A:,k, Aj:k,i:l, Aj:k,:, A:,j:k are indexing expressions for a matrix A, with Ak,: representing the kth column of matrix A, A:,k representing the kth row of matrix A, Aj:k,i:l representing a submatrix including the elements of A indexed by j:k in the first dimension and i:l in the second dimension, Aj:k,: and A:,j:k include all subscripts in one of the dimensions but use the vector j:k to index in the other dimension.The paper is organized as follows. Section2 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.


## 2. Problem Formulation

Consider a MIMO passive radar like the one illustrated in Figure1. This MIMO passive radar system is comprised of N geographically separated receivers with positions sr,no=xr,no,yr,no,zr,noT, n=1,2,…,N, and M noncooperative transmitters with positions st,mo=xt,mo,yt,mo,zt,moT, m=1,2,…,M. Let uo=xo,yo,zoT be unknown position coordinates of the target. The transmitters send out a set of waveforms, which are then reflected by potential targets. The 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.Figure 1
A practical localization scenario with MIMO passive radar.In a practical localization scenario, the actual positions of all transmitters and receivers are not known, and the available positions are as follows:(1)st,m=st,mo+Δst,m,m=1,2,…,M,sr,n=sr,no+Δsr,n,n=1,2,…,N,where st,m and sr,n represent the available (noisy) positions of transmitter m and receiver n, and Δst,m and Δsr,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(2)β=βo+Δβ,where β=stT,srTT with st=st,1T,st,2T,…,st,MTT and Δsr=Δsr,1T,Δsr,2T,…,Δsr,NTT denotes the noisy transmitter and receiver position vector, βo=stoT,sroTT with sto=st,1oT,st,2oT,…,st,MoTT and sro=sr,1oT,sr,2oT,…,sr,NoTT stands for the actual transmitter and receiver position vector, and Δβ=ΔstT,ΔsrTT with Δst=Δst,1T,Δst,2T,…,Δst,MTT and Δsr=Δsr,1T,Δsr,2T,…,Δsr,NTT stands for the transmitter and receiver position errors vector assumed to be zero-mean Gaussian with covariance matrix(3)EΔβΔβT=Qβ,where the covariance matrix Qβ can be determined from theoretical analysis and actual measurement [25].Based on the above geometry, the range between transmitterm and the target is Rt,mo=uo−st,mo, and the baseline distance between transmitter m and receiver n is Rt,m,r,no=st,mo−sr,no. Thus, the time difference between the direct path signal from transmitter m and the corresponding reflected signal arriving at receiver n can be given by(4)τm,no=1cRt,mo+Rr,no−Rt,m,r,no.The true AOA pair for receivern, that is, the elevation angle denoted by θno and the azimuth angle denoted by φno, is, respectively, given by(5)θno=arc tanyo−yr,noxo−xr,no,φno=arc tanzo−zr,noxo−xr,no2+yo−yr,no2.Considering the unavoidable measurement noises in reality, the erroneous version of the TDOA and AOA measurements can be expressed as(6)τm,n=τm,no+Δτm,n,θn=θno+Δθn,φn=φno+Δφn,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(7)α=αo+Δα,where α=θT,φT,τTT with τ=τ1T,τ2T,…,τMTT, τm=τm,1,τm,2,…,τm,NT, θ=θ1,θ2,…,θNT, and φ=φ1,φ2,…,φNT stands for the noisy TDOA and AOA measurement vector, αo=θoT,φoT,τoTT with τo=τ1oT,τ2oT,…,τMoTT, τmo=τm,1o,τm,2o,…,τm,NoT, θo=θ1o,θ2o,…,θNoT, and φo=φ1o,φ2o,…,φNoT stands for the true TDOA and AOA vector, and Δα=ΔθT,ΔφT,ΔτTT with Δτ=Δτ1T,Δτ2T,…,ΔτMTT, Δτm=Δτm,1,Δτm,2,…,Δτm,NT, Δθ=Δθ1,Δθ2,…,ΔθNT, and Δφ=Δφ1,Δφ2,…,ΔφNT stands for the TDOA and AOA measurement noise vector assumed to be zero-mean Gaussian with covariance matrix.(8)EΔαΔαT=Qα,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.


## 3. 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.


### 3.1. CRLB with Transmitter and Receiver Position Errors

In addition to TDOA/AOA measurement noises, the presence of transmitter and receiver position errors is included. Therefore, the unknown parameter vector for the CRLB evaluation isηo=uoT,βoTT, and the observation vector is z=αT,βTT. 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(9)lnpz|ηo=lnpα|ηo+lnpβ|ηo=κ−12α−αoTQα−1α−αo−12β−βoTQβ−1β−βo,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(10)FIMηo=E∂lnpz|ηo∂ηo∂lnpz|ηo∂ηoT,=XYYTZ,and the blocks X, Y, and Z are given by(11)X=∂αo∂uoTQα−1∂αo∂uo,Y=∂αo∂uoTQα−1∂αo∂βo,Z=∂αo∂βoTQα−1∂αo∂βo+Qβ−1,where ∂αo/∂uo and ∂αo/∂βo are the partial derivatives of the parametric form of αo with respect to uo and βo, respectively, with their elements given by(12)∂αo∂uon,1=−yo−yr,noxo−xr,no2+yo−yr,no2∂αo∂uon,2=xo−xr,noxo−xr,no2+yo−yr,no2∂αo∂uoN+n,1=−xo−xr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂uoN+n,2=−yo−yr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂uoN+n,3=xo−xr,no2+yo−yr,no2Rr,no2∂αo∂uo2N+m−1N+n,1:3=uo−st,moTcRt,mo+uo−sr,noTcRr,no∂αo∂βon,3M+3n−2=−yr,no−yoxo−xr,no2+yo−yr,no2∂αo∂βon,3M+3n−1=xr,no−xoxo−xr,no2+yo−yr,no2∂αo∂βoN+n,3M+3n−2=−xo−xr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂βoN+n,3M+3n−1=−yo−yr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂βoN+n,3M+3n=xo−xr,no2+yo−yr,no2Rr,no2∂αo∂βo2N+m−1N+n,3m−2:3m=st,mo−uoTcRt,mo−st,mo−sr,noTcRt,m,r,no∂αo∂βo2N+m−1N+n,3M+3n−2:3M+3n=sr,no−uoTcRr,no−sr,no−st,moTcRt,m,r,nofor 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,(13)CRLBηo=FIMηo−1,CRLBηo is a (3+3M+3 N)-by-(3+3M+3 N) matrix, of which the upper left 3-by-3 block gives the CRLB of uo. Substituting (14) in (16) and then invoking the partitioned matrix inversion formula and the matrix inversion lemma [27, 28], we have(14)CRLBuo=X−YZ−1YT−1=X−1+X−1YZ−YTX−1Y−1YTX−1.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 X−1YZ−YTX−1Y−1YTX−1 is the change in CRLB brought by the transmitter and receiver position errors.


### 3.2. 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, X−1YZ−YTX−1Y−1YTX−1. In order to prove X−1YZ−YTX−1Y−1YTX−1 is positive definite, we only need to prove YZ−YTX−1Y−1YT is positive definite since X−1 is full rank.SinceX=∂αo/∂uoTQα−1∂αo/∂uo is invertible, we can deduce Qα−1∂αo/∂uo has a full column rank of 3. Moreover, it can be inferred from the expression of ∂αo/∂βo that ∂αo/∂βoT has a full column rank of (MN+2N). Therefore, utilizing the expression of Y in (11), YT has a full column rank.Next, we proceed to proveZ−YTX−1Y is positive definite. From the expressions of Y and Z in (11), we obtain (15) after some mathematical manipulations(15)Z−YTX−1Y=∂αo∂βoTQα−1−Qα−1∂αo∂uoX−1∂αo∂uoTQα−1∂αo∂βo+Qβ−1,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(16)Qα−1−Qα−1∂αo∂uoX−1∂αo∂uoTQα−1=LαI−ΦΦTΦ−1ΦTLαT,where Φ=LαT∂αo/∂uo. It is readily to see from (16) that Qα−1−Qα−1∂αo/∂uoX−1∂αo/∂uoTQα−1 is positive semidefinite since I−ΦΦTΦ−1ΦT is a projection matrix. Therefore, Z−YTX−1Y in (15) is positive definite as long as Qβ−1≠O.Synthesizing the above results, that is,YT has a full column rank, and X−1 and Z−YTX−1Y are positive definite, we can infer that X−1YZ−YTX−1Y−1YTX−1 is a positive definite matrix. That is to say, we can obtain (17) from (14) that(17)CRLBuo≻X−1,where CRLBuo≻X−1 represents CRLBuo−X−1 is positive definite. This condition implies trCRLBuo>trX−1. trCRLBuo and trX−1 have the physical meaning, that is, the minimum possible mean-square position errors for the target localization with and without transmitter/receiver position errors, respectively. Thus, 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.


## 3.1. CRLB with Transmitter and Receiver Position Errors

In addition to TDOA/AOA measurement noises, the presence of transmitter and receiver position errors is included. Therefore, the unknown parameter vector for the CRLB evaluation isηo=uoT,βoTT, and the observation vector is z=αT,βTT. 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(9)lnpz|ηo=lnpα|ηo+lnpβ|ηo=κ−12α−αoTQα−1α−αo−12β−βoTQβ−1β−βo,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(10)FIMηo=E∂lnpz|ηo∂ηo∂lnpz|ηo∂ηoT,=XYYTZ,and the blocks X, Y, and Z are given by(11)X=∂αo∂uoTQα−1∂αo∂uo,Y=∂αo∂uoTQα−1∂αo∂βo,Z=∂αo∂βoTQα−1∂αo∂βo+Qβ−1,where ∂αo/∂uo and ∂αo/∂βo are the partial derivatives of the parametric form of αo with respect to uo and βo, respectively, with their elements given by(12)∂αo∂uon,1=−yo−yr,noxo−xr,no2+yo−yr,no2∂αo∂uon,2=xo−xr,noxo−xr,no2+yo−yr,no2∂αo∂uoN+n,1=−xo−xr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂uoN+n,2=−yo−yr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂uoN+n,3=xo−xr,no2+yo−yr,no2Rr,no2∂αo∂uo2N+m−1N+n,1:3=uo−st,moTcRt,mo+uo−sr,noTcRr,no∂αo∂βon,3M+3n−2=−yr,no−yoxo−xr,no2+yo−yr,no2∂αo∂βon,3M+3n−1=xr,no−xoxo−xr,no2+yo−yr,no2∂αo∂βoN+n,3M+3n−2=−xo−xr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂βoN+n,3M+3n−1=−yo−yr,nozo−zr,noRr,no2xo−xr,no2+yo−yr,no2∂αo∂βoN+n,3M+3n=xo−xr,no2+yo−yr,no2Rr,no2∂αo∂βo2N+m−1N+n,3m−2:3m=st,mo−uoTcRt,mo−st,mo−sr,noTcRt,m,r,no∂αo∂βo2N+m−1N+n,3M+3n−2:3M+3n=sr,no−uoTcRr,no−sr,no−st,moTcRt,m,r,nofor 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,(13)CRLBηo=FIMηo−1,CRLBηo is a (3+3M+3 N)-by-(3+3M+3 N) matrix, of which the upper left 3-by-3 block gives the CRLB of uo. Substituting (14) in (16) and then invoking the partitioned matrix inversion formula and the matrix inversion lemma [27, 28], we have(14)CRLBuo=X−YZ−1YT−1=X−1+X−1YZ−YTX−1Y−1YTX−1.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 X−1YZ−YTX−1Y−1YTX−1 is the change in CRLB brought by the transmitter and receiver position errors.


## 3.2. 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, X−1YZ−YTX−1Y−1YTX−1. In order to prove X−1YZ−YTX−1Y−1YTX−1 is positive definite, we only need to prove YZ−YTX−1Y−1YT is positive definite since X−1 is full rank.SinceX=∂αo/∂uoTQα−1∂αo/∂uo is invertible, we can deduce Qα−1∂αo/∂uo has a full column rank of 3. Moreover, it can be inferred from the expression of ∂αo/∂βo that ∂αo/∂βoT has a full column rank of (MN+2N). Therefore, utilizing the expression of Y in (11), YT has a full column rank.Next, we proceed to proveZ−YTX−1Y is positive definite. From the expressions of Y and Z in (11), we obtain (15) after some mathematical manipulations(15)Z−YTX−1Y=∂αo∂βoTQα−1−Qα−1∂αo∂uoX−1∂αo∂uoTQα−1∂αo∂βo+Qβ−1,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(16)Qα−1−Qα−1∂αo∂uoX−1∂αo∂uoTQα−1=LαI−ΦΦTΦ−1ΦTLαT,where Φ=LαT∂αo/∂uo. It is readily to see from (16) that Qα−1−Qα−1∂αo/∂uoX−1∂αo/∂uoTQα−1 is positive semidefinite since I−ΦΦTΦ−1ΦT is a projection matrix. Therefore, Z−YTX−1Y in (15) is positive definite as long as Qβ−1≠O.Synthesizing the above results, that is,YT has a full column rank, and X−1 and Z−YTX−1Y are positive definite, we can infer that X−1YZ−YTX−1Y−1YTX−1 is a positive definite matrix. That is to say, we can obtain (17) from (14) that(17)CRLBuo≻X−1,where CRLBuo≻X−1 represents CRLBuo−X−1 is positive definite. This condition implies trCRLBuo>trX−1. trCRLBuo and trX−1 have the physical meaning, that is, the minimum possible mean-square position errors for the target localization with and without transmitter/receiver position errors, respectively. Thus, 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.


## 4. Algorithm Development and Analysis

The degradation of localization accuracy in the presence of transmitter and receiver position errors has been shown in Section3 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.


### 4.1. Localization Algorithm

The 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(18)cτm,no+Rt,m,r,no−Rr,no=Rt,mo.Squaring both sides of (18), and then rearranging it as(19)2sr,no−st,moTuo+2cτm,no+Rt,m,r,noRr,no=cτm,no+Rt,m,r,no2+sr,noTsr,no−st,moTst,mo.Since the nuisance parameterRr,no are nonlinearly related to the target position uo, (19) is pseudo-linear with respect to uo. In order to thoroughly linearize (19), rewrite (6) as(20)sinφnoRr,no=zo−zr,no.Putting (20) into (19) yields(21)2xr,no−xt,moxo+2yr,no−yt,moyo+2zr,no−zt,mo+2cτm,no+Rt,m,r,nosinφnozo=2cτm,no+Rt,m,r,nozr,nosinφno+cτm,no+Rt,m,r,no2+sr,noTsr,no−st,moTst,mo.Since only the noisy values ofsr,no, τm,no, and φno are available, we replace them by st,mo=st,m−Δst,m, sr,no=sr,n−Δsr,n, τm,no=τm,n−Δτm,n, and φno=φn−Δφn and expand them around the noisy values and then ignore the second- and higher-order error terms. Then, we have(22)2xr,n−xt,mxo+2yr,n−yt,myo+2zr,n−zt,m+2cτm,n+Rt,m,r,nsinφnzo=2cτm,n+Rt,m,r,nzr,nsinφn+sr,nTsr,n−st,mTst,m+cτm,n+Rt,m,r,n2+−2uo+zo−zr,nsinφnst,m−sr,nRt,m,r,n−2cτm,n+Rt,m,r,nst,m−sr,nRt,m,r,n+2st,mTΔst,m+2uo+zo−zr,nsinφnsr,n−st,mRt,m,r,n−2cτm,n+Rt,m,r,nsr,n−st,mRt,m,r,n+2sr,nTΔsr,n−2cτm,n+Rt,m,r,nsinφnΔzr,n−2cτm,n+Rt,m,r,nzo−zr,nsinφn2Δφn−2Rt,mocΔτm,n.To linearize the AOA equations, rewrite (6) and (7) as(23)xotanθno−yo=xr,notanθno−yr,no,xotanφno1+tanθno2−zo=xr,notanφno1+tanθno2−zr,no.Substituteθno=θn−Δθn, φno=φn−Δφn, and sr,no=sr,n−Δsr,n into (25) and (23), expand them in Taylor’s series keeping only terms below second order:(24)tanθnxo−yo=tanθnxr,n−yr,n+xo−xr,ncosθn2Δθn−tanθnΔxr,n+Δyr,n,xotanφn1+tanθn2−zo=tanφn1+tanθn2xr,n−zr,n−tanφn1+tanθn2Δxr,n+Δzr,n+xo−xr,ntanφntanθn1+tanθn2cosθn2Δθn+xo−xr,n1+tanθn2cosφn2Δφn.Now, collecting (25), (28), and (29) with respect to theM transmitters and N receivers, we can rearrange them in matrix form as(25)Guo=h+Δh,where G is a (MN+2N)-by-3 matrix and h is a (MN+2N)-by-1 matrix, respectively, with their elements given by(26)Gn,:=tanθn,−1,0,GN+n,:=tanφn1+tanθn2,0,−1,G2N+m−1N+n,1:2=2xr,n−xt,m,2yr,n−yt,m,G2N+m−1N+n,1:2=2cτm,n+Rt,m,r,nsinφn,hn=tanθnxr,n−yr,nhN+n=xr,ntanφn1+tanθn2−zr,nh2N+m−1N+n=2cτm,n+Rt,m,r,nzr,nsinφn+cτm,n+Rt,m,r,n2+sr,nTsr,n−st,mTst,m.form=1,2,…,M, n=1,2,…,N. The error vector Δh can be expressed as(27)Δh=BΔα+DΔβ,where B and D are (MN+2N)-by-(MN+2N) matrix and (MN+2N)×(3M+3 N) matrix, respectively, with their elements given by(28)Bn,n=xo−xr,ncosθn2BN+n,n=xo−xr,ntanφntanθn1+tanθn2cosθn2BN+n,N+n=xo−xr,n1+tanθn2cosφn2B2N+m−1N+n,N+n=−2cτm+Rt,m,r,nzo−zr,nsinφn2B2N+m−1N+n,2N+m−1N+n=−2cRt,moDn,3M+3n−2=−tanθnDn,3M+3n−1=1DN+n,3M+3n−2=−tanφn1+tanθn2DN+n,3M+3n=1D2N+m−1N+n,3m−2:3m=−2uo+zo−zr,nsinφnst,m−sr,nRt,m,r,n−2cτm+Rt,m,r,nst,m−sr,nRt,m,r,n+2st,mTD2N+m−1N+n,3M+3n−2=2xo+zo−zr,nsinφnxr,n−xt,mRt,m,r,n−2cτm,n+Rt,m,r,nxr,n−xt,mRt,m,r,n+2xr,nD2N+m−1N+n,3M+3n−1=2yo+zo−zr,nsinφnyr,n−yt,mRt,m,r,n−2cτm,n+Rt,m,r,nyr,n−yt,mRt,m,r,n+2yr,nD2N+m−1N+n,3M+3n=2zo+zo−zr,nsinφnzr,n−zt,mRt,m,r,n−2cτm,n+Rt,m,r,nzr,n−zt,mRt,m,r,n+2zr,n−2cτm,n+Rt,m,r,nsinφ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 uo, from which the target position uo can be estimated using the WLS minimization. The WLS solution of (25), which minimizes the cost function EΔhTWΔh with respect to the target position uo, is given by(29)u=GTWG−1GTWh,where W is the weighting matrix, and its optimal choice that gives the minimum parameter variance is given by(30)W=EΔhΔhT−1=BQαBT+DQβDT−1.However, as suggested in (30), the computation of W relies on the unknown target position uo via B and D. Thus, it is unrealistic to evaluate W directly. To resolve this difficulty, a realistic and workable way is to form a least squares (LS) solution of target position uo using (29) with W=IMN+2N. Using this LS solution, we update weighting matrix W via (33) and employ (32) once again to form a more accurate estimate of the target position uo.Rewritinguo as the identical equation uo=GTWG−1GTWGuo and subtracting it from sides of (29) yield(31)Δu=GTWG−1GTWΔh.Ignoring the second- and higher-order error terms in (31) and taking expectation result in that EΔu=03, 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)covu=GTWG−1.


### 4.2. Performance Analysis

We shall evaluate the efficiency of the proposed solution by comparing its covariance matrix with the CRLB in (14). By inserting (33) into (35), we have after some mathematical manipulations(33)covu=X⌢−Y⌢Z⌢−1Y⌢T−1,where(34)X⌢=GαTQα−1Gα,Y⌢=GαTQα−1Gβ,Z⌢=GβTQα−1Gβ+Qβ−1,with(35)Gα=B−1G,Gβ=B−1D.Comparing (36) with the CRLB in (17), we observe that CRLBuo and covu 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(36)Gα≃∂αo∂uo,Gβ≃∂αo∂βo.Based on this, it can be inferred that(37)covu≃CRLBuo.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.


## 4.1. Localization Algorithm

The 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(18)cτm,no+Rt,m,r,no−Rr,no=Rt,mo.Squaring both sides of (18), and then rearranging it as(19)2sr,no−st,moTuo+2cτm,no+Rt,m,r,noRr,no=cτm,no+Rt,m,r,no2+sr,noTsr,no−st,moTst,mo.Since the nuisance parameterRr,no are nonlinearly related to the target position uo, (19) is pseudo-linear with respect to uo. In order to thoroughly linearize (19), rewrite (6) as(20)sinφnoRr,no=zo−zr,no.Putting (20) into (19) yields(21)2xr,no−xt,moxo+2yr,no−yt,moyo+2zr,no−zt,mo+2cτm,no+Rt,m,r,nosinφnozo=2cτm,no+Rt,m,r,nozr,nosinφno+cτm,no+Rt,m,r,no2+sr,noTsr,no−st,moTst,mo.Since only the noisy values ofsr,no, τm,no, and φno are available, we replace them by st,mo=st,m−Δst,m, sr,no=sr,n−Δsr,n, τm,no=τm,n−Δτm,n, and φno=φn−Δφn and expand them around the noisy values and then ignore the second- and higher-order error terms. Then, we have(22)2xr,n−xt,mxo+2yr,n−yt,myo+2zr,n−zt,m+2cτm,n+Rt,m,r,nsinφnzo=2cτm,n+Rt,m,r,nzr,nsinφn+sr,nTsr,n−st,mTst,m+cτm,n+Rt,m,r,n2+−2uo+zo−zr,nsinφnst,m−sr,nRt,m,r,n−2cτm,n+Rt,m,r,nst,m−sr,nRt,m,r,n+2st,mTΔst,m+2uo+zo−zr,nsinφnsr,n−st,mRt,m,r,n−2cτm,n+Rt,m,r,nsr,n−st,mRt,m,r,n+2sr,nTΔsr,n−2cτm,n+Rt,m,r,nsinφnΔzr,n−2cτm,n+Rt,m,r,nzo−zr,nsinφn2Δφn−2Rt,mocΔτm,n.To linearize the AOA equations, rewrite (6) and (7) as(23)xotanθno−yo=xr,notanθno−yr,no,xotanφno1+tanθno2−zo=xr,notanφno1+tanθno2−zr,no.Substituteθno=θn−Δθn, φno=φn−Δφn, and sr,no=sr,n−Δsr,n into (25) and (23), expand them in Taylor’s series keeping only terms below second order:(24)tanθnxo−yo=tanθnxr,n−yr,n+xo−xr,ncosθn2Δθn−tanθnΔxr,n+Δyr,n,xotanφn1+tanθn2−zo=tanφn1+tanθn2xr,n−zr,n−tanφn1+tanθn2Δxr,n+Δzr,n+xo−xr,ntanφntanθn1+tanθn2cosθn2Δθn+xo−xr,n1+tanθn2cosφn2Δφn.Now, collecting (25), (28), and (29) with respect to theM transmitters and N receivers, we can rearrange them in matrix form as(25)Guo=h+Δh,where G is a (MN+2N)-by-3 matrix and h is a (MN+2N)-by-1 matrix, respectively, with their elements given by(26)Gn,:=tanθn,−1,0,GN+n,:=tanφn1+tanθn2,0,−1,G2N+m−1N+n,1:2=2xr,n−xt,m,2yr,n−yt,m,G2N+m−1N+n,1:2=2cτm,n+Rt,m,r,nsinφn,hn=tanθnxr,n−yr,nhN+n=xr,ntanφn1+tanθn2−zr,nh2N+m−1N+n=2cτm,n+Rt,m,r,nzr,nsinφn+cτm,n+Rt,m,r,n2+sr,nTsr,n−st,mTst,m.form=1,2,…,M, n=1,2,…,N. The error vector Δh can be expressed as(27)Δh=BΔα+DΔβ,where B and D are (MN+2N)-by-(MN+2N) matrix and (MN+2N)×(3M+3 N) matrix, respectively, with their elements given by(28)Bn,n=xo−xr,ncosθn2BN+n,n=xo−xr,ntanφntanθn1+tanθn2cosθn2BN+n,N+n=xo−xr,n1+tanθn2cosφn2B2N+m−1N+n,N+n=−2cτm+Rt,m,r,nzo−zr,nsinφn2B2N+m−1N+n,2N+m−1N+n=−2cRt,moDn,3M+3n−2=−tanθnDn,3M+3n−1=1DN+n,3M+3n−2=−tanφn1+tanθn2DN+n,3M+3n=1D2N+m−1N+n,3m−2:3m=−2uo+zo−zr,nsinφnst,m−sr,nRt,m,r,n−2cτm+Rt,m,r,nst,m−sr,nRt,m,r,n+2st,mTD2N+m−1N+n,3M+3n−2=2xo+zo−zr,nsinφnxr,n−xt,mRt,m,r,n−2cτm,n+Rt,m,r,nxr,n−xt,mRt,m,r,n+2xr,nD2N+m−1N+n,3M+3n−1=2yo+zo−zr,nsinφnyr,n−yt,mRt,m,r,n−2cτm,n+Rt,m,r,nyr,n−yt,mRt,m,r,n+2yr,nD2N+m−1N+n,3M+3n=2zo+zo−zr,nsinφnzr,n−zt,mRt,m,r,n−2cτm,n+Rt,m,r,nzr,n−zt,mRt,m,r,n+2zr,n−2cτm,n+Rt,m,r,nsinφ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 uo, from which the target position uo can be estimated using the WLS minimization. The WLS solution of (25), which minimizes the cost function EΔhTWΔh with respect to the target position uo, is given by(29)u=GTWG−1GTWh,where W is the weighting matrix, and its optimal choice that gives the minimum parameter variance is given by(30)W=EΔhΔhT−1=BQαBT+DQβDT−1.However, as suggested in (30), the computation of W relies on the unknown target position uo via B and D. Thus, it is unrealistic to evaluate W directly. To resolve this difficulty, a realistic and workable way is to form a least squares (LS) solution of target position uo using (29) with W=IMN+2N. Using this LS solution, we update weighting matrix W via (33) and employ (32) once again to form a more accurate estimate of the target position uo.Rewritinguo as the identical equation uo=GTWG−1GTWGuo and subtracting it from sides of (29) yield(31)Δu=GTWG−1GTWΔh.Ignoring the second- and higher-order error terms in (31) and taking expectation result in that EΔu=03, 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)covu=GTWG−1.


## 4.2. Performance Analysis

We shall evaluate the efficiency of the proposed solution by comparing its covariance matrix with the CRLB in (14). By inserting (33) into (35), we have after some mathematical manipulations(33)covu=X⌢−Y⌢Z⌢−1Y⌢T−1,where(34)X⌢=GαTQα−1Gα,Y⌢=GαTQα−1Gβ,Z⌢=GβTQα−1Gβ+Qβ−1,with(35)Gα=B−1G,Gβ=B−1D.Comparing (36) with the CRLB in (17), we observe that CRLBuo and covu 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(36)Gα≃∂αo∂uo,Gβ≃∂αo∂βo.Based on this, it can be inferred that(37)covu≃CRLBuo.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.


## 5. Numerical Examples

This section contains some numerical simulations to evaluate the performance of the proposed solution. The localization scenario is set as shown in Figure2, where a MIMO passive radar system comprised of N = 6 geographically separated receivers and M = 4 transmitters is deployed to locate a target at position uo=100,200,100Tm. The positions of the transmitters and receivers are listed in Table 1. MN = 24 TDOAs and N = 6 AOA pairs are extracted to determine the target position.Figure 2
Localization geometry for simulations.Table 1
Positions of the transmitters and receivers.
TransmitterTX1TX2TX3TX4ReceiverRX1RX2RX3RX4RX5RX6xt,mom–200–200200200xr,nom–4504500600–6000yt,mom–300300300–300yr,nom–45045060000–600zt,mom25010080120zr,nom200100200150150100In order to simulate a practical localization scenario, zero-mean Gaussian noises with known covariance matricesQα=diagσt2IMN,σa2I2N and Qβ=σs2I3M+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. The localization accuracy is assessed using root mean squares error (RMSE) defined as(38)RMSEu=1L∑l=1Lul−uo2,where ul is the estimation of uo at the lth Monte Carlo trial, and the number of Monte Carlo trials for each simulation is set as L=5000.


### 5.1. 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. The comparison results are exhibited in Figure3.Figure 3
Comparisons of the CRLBs with and without transmitter/receiver position error. (a) CRLB comparison as TDOA measurement noise levelσt increases and σa=0.1°, σs=50m. (b) CRLB comparison as AOA measurement noise level σa increases and σt=0.1us, σs=50m. (c) CRLB comparison as transmitter/receiver position error level σs increases and σt=0.1us, σa=0.1°.
(a)(b)(c)Figure3(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. The 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. This reconfirms the significance of including the transmitter and receiver position errors in the design of the localization algorithm for MIMO passive radar.


### 5.2. 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. The comparison results are presented in Figure 4.Figure 4
RMSE comparison of different algorithms when there are no transmitter and receiver position errors. (a) RMSE comparison as TDOA measurement noise levelσt increases and σa=0.1°. (b) RMSE comparison as AOA measurement noise level σa increases and σt=0.1us.
(a)(b)Figure4 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. The 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. The results are given in Figure5.Figure 5
RMSE comparison of different algorithms when there are transmitter and receiver position errors. (a) RMSE comparison as TDOA measurement noise levelσt increases and σa=0.1°, σs=50m. (b) RMSE comparison as AOA measurement noise level σa increases and σt=0.1us, σs=50m. (c) RMSE comparison as transmitter/receiver position error level σs increases and σt=0.1us, σa=0.1°.
(a)(b)(c)Figures5(a) and 5(b) give the RMSE curves of the algorithms under different TDOA measurement noise levels and AOA measurement noises, respectively. The transmitter and receiver position errors level is set to σs=50m. As can be seen from Figures 5(a) and 5(b), when there are non-negligible transmitter and receiver position errors, the localization RMSEs of Noroozi’s method and Amiri’s method cannot reach the CRLB, even at small TDOA/AOA measurement noise regions. This 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 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. Thus, the proposed solution performs much better than Noroozi’s method and Amiri’s method over large transmitter and receiver position errors regions. This 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.


### 5.3. 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 asGDOPu=σx2+σy2+σz2 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. The contour plots of the GDOP values in an 4km×4km area around the MIMO passive radar are provided in Figure 6.Figure 6
GDOP contours withσt=0.1us, σa=0.1°, and σs=50m. (a) zo=50m. (b) zo=500m.
(a)(b)Figure6(a) presents the GDOP contours at zo=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 zo=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 inﬂuence of the transmitter/receiver position on target localization accuracy. This simulation shows that when the target z-coordinate becomes larger, the GDOP also becomes larger and the inﬂuence of transmitter/receiver position errors decreases. This can be useful in our application.


### 5.4. Computation Complexity Comparison

The 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 independent Monte Carlo runs. The main configuration of the computer is shown as follows: Intel(R) Core(TM) CPU i5-7200U@2.50GHz; 8.00 G RAM; Windows 10 64 bit Operating System; Matlab 2018a Software. The comparison results are given in Table2.Table 2
Time cost of the algorithms.
AlgorithmsAverage run time (ms)Noroozi1’s method1.76Amiri’s method0.98Proposed solution2.66As presented in Table2, the time cost of Noroozi1’s method is almost twice higher than that of Amiri’s method. This is because Noroozi1’s method needs two WLS stages while Amiri’s method identifies the target position in only one WLS stage. The proposed solution incurs the highest time cost among the algorithms, approximately 3 times higher than Amiri’s method. This is not surprising since the proposed solution took the transmitter and receiver position errors into account while Noroozi1 and Amiri’s methods do not. That 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.


## 5.1. 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. The comparison results are exhibited in Figure3.Figure 3
Comparisons of the CRLBs with and without transmitter/receiver position error. (a) CRLB comparison as TDOA measurement noise levelσt increases and σa=0.1°, σs=50m. (b) CRLB comparison as AOA measurement noise level σa increases and σt=0.1us, σs=50m. (c) CRLB comparison as transmitter/receiver position error level σs increases and σt=0.1us, σa=0.1°.
(a)(b)(c)Figure3(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. The 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. This reconfirms the significance of including the transmitter and receiver position errors in the design of the localization algorithm for MIMO passive radar.


## 5.2. 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. The comparison results are presented in Figure 4.Figure 4
RMSE comparison of different algorithms when there are no transmitter and receiver position errors. (a) RMSE comparison as TDOA measurement noise levelσt increases and σa=0.1°. (b) RMSE comparison as AOA measurement noise level σa increases and σt=0.1us.
(a)(b)Figure4 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. The 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. The results are given in Figure5.Figure 5
RMSE comparison of different algorithms when there are transmitter and receiver position errors. (a) RMSE comparison as TDOA measurement noise levelσt increases and σa=0.1°, σs=50m. (b) RMSE comparison as AOA measurement noise level σa increases and σt=0.1us, σs=50m. (c) RMSE comparison as transmitter/receiver position error level σs increases and σt=0.1us, σa=0.1°.
(a)(b)(c)Figures5(a) and 5(b) give the RMSE curves of the algorithms under different TDOA measurement noise levels and AOA measurement noises, respectively. The transmitter and receiver position errors level is set to σs=50m. As can be seen from Figures 5(a) and 5(b), when there are non-negligible transmitter and receiver position errors, the localization RMSEs of Noroozi’s method and Amiri’s method cannot reach the CRLB, even at small TDOA/AOA measurement noise regions. This 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 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. Thus, the proposed solution performs much better than Noroozi’s method and Amiri’s method over large transmitter and receiver position errors regions. This 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.


## 5.3. 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 asGDOPu=σx2+σy2+σz2 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. The contour plots of the GDOP values in an 4km×4km area around the MIMO passive radar are provided in Figure 6.Figure 6
GDOP contours withσt=0.1us, σa=0.1°, and σs=50m. (a) zo=50m. (b) zo=500m.
(a)(b)Figure6(a) presents the GDOP contours at zo=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 zo=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 inﬂuence of the transmitter/receiver position on target localization accuracy. This simulation shows that when the target z-coordinate becomes larger, the GDOP also becomes larger and the inﬂuence of transmitter/receiver position errors decreases. This can be useful in our application.


## 5.4. Computation Complexity Comparison

The 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 independent Monte Carlo runs. The main configuration of the computer is shown as follows: Intel(R) Core(TM) CPU i5-7200U@2.50GHz; 8.00 G RAM; Windows 10 64 bit Operating System; Matlab 2018a Software. The comparison results are given in Table2.Table 2
Time cost of the algorithms.
AlgorithmsAverage run time (ms)Noroozi1’s method1.76Amiri’s method0.98Proposed solution2.66As presented in Table2, the time cost of Noroozi1’s method is almost twice higher than that of Amiri’s method. This is because Noroozi1’s method needs two WLS stages while Amiri’s method identifies the target position in only one WLS stage. The proposed solution incurs the highest time cost among the algorithms, approximately 3 times higher than Amiri’s method. This is not surprising since the proposed solution took the transmitter and receiver position errors into account while Noroozi1 and Amiri’s methods do not. That 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.


## 6. 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.



---

*Source: 9964332-2022-03-15.xml*