Figure 1. Principle of four-satellite TDOA.
, and .
Ideally for satellites location, , Equation (4) can be obtained by using pseudo inverse to solve Equation (3):
The accuracy analysis of four-satellite TDOA localization algorithm is as follows.
, , ,
, , .
We can obtain
Let . Since the main satellite has an error in measuring the arrival time of the radiation source signal, and the error exists in each time difference measurement. The measurement errors of each are related. we assumed that the mean of the measurement error is 0 after the systematic calibration. Moreover, each element of satellite position error is uncorrelated and each satellite position is uncorrelated. The covariance matrix of the target radiation source localization error can be obtained as follows.
where , is the variance
of the measurement error of the distance between the satellite and the main satellite. is the correlation coefficient between and . We assume that the variance of satellite position error on each component is the same, that is . We can obtain the covariance matrix of the satellite position error as Equation (8).
The covariance matrix of the target radiation source position error is as Equation (9).
According to the definition of the GDOP, the GDOP of TDOA localization in three dimensional space can be calculated as Equation (10).
3. Principle of Four-Satellite TDOA/FDOA
On the basis of the four-satellite TDOA localization algorithm, the four-satellite TDOA/FDOA localization algorithm adds frequency difference information. The coordinate of the target is recorded as . The coordinate of the main satellite is recorded as . The speed of the main satellite is recorded as . The coordinate of the subsatellites is recorded as . The speed of the subsatellite is recorded as , where . According to the time difference and frequency difference of the target radiation source arriving at the main satellite and subsatellites, we can obtain that
where . And is the time difference of arrival, is the frequency difference of arrival, is the radiation frequency of radiation source, is the light speed. By solving the least square solution of the equation, the target position can be obtained. The purpose of this paper is to analyze the accuracy of the four-Satellite localization, so the solution process of the equation is no longer exhausted.
Using first order linear approximation expansion for , , , and , we can obtain the Equation (13).
Take the upper equation as
Divide the measurement accuracy of the subsatellite position into measurement accuracy of absolute position and measurement accuracy of relative position. That is
The right is the measurement accuracy of relative position between main satellite and subsatellites. is the measurement accuracy of absolute position of main satellite. We can obtain the Equation (15).
Let , , , , , , , . And we assume that the measurement accuracy of position, the measurement accuracy of speed, the measurement accuracy of TDOA, the measurement accuracy of FDOA are uncorrelated. We can obtain Equation (16) by matrix operation.
where is the correlation matrix of measurement accuracy between TDOA and FDOA. is the correlation matrix of the measurement accuracy of position and the measurement accuracy of speed .
The positioning accuracy can be expressed by the measuring accuracy of 3 orthogonal directions as Equation (17)
From the above analysis results, we can see that the positioning accuracy is calculated by and .
We assume that the mean of the TDOA signal and the FDOA signal measurement accuracy is 0, and they are uncorrelated. We can obtain that
where is a unit matrix， is the measurement accuracy of time difference, is the measurement accuracy of frequency difference.
The absolute position measurement accuracy component are considered to be independent of the random variables. Therefore, the correlation matrix of four-satellite absolute position measurement accuracy is
where is the measurement accuracy of absolute position.
Comprehensively considering the measurement accuracy of relative position and the ranging accuracy, the correlation matrix of relative position measurement accuracy can be represented as
where , and is the measurement accuracy of relative position, is the ranging accuracy.
The correlation matrix of speed measurement accuracy can be represented as where is measurement accuracy of speed.
4. Simulation Results
We assume that the satellite orbit altitude is 200 km and the length of baseline is 60 km. The Y configuration, T configuration and diamond configuration are adopted respectively. The measurement accuracy of time difference is 50 ns. The measurement accuracy of frequency difference is 10 Hz. The measurement accuracy of speed is 1 m/s. The absolute position measurement accuracy is 150 m. The relative position measurement accuracy is 50 m. The GDOP distributions of TDOA localization and TDOA/FDOA localization are obtained by simulation, as shown in Figures 2-4.
Analysis and comparison of Figures 2-4 show that under the same parameters, the accuracy of TDOA/FDOA localization is better than that of TDOA localization. And for the diamond configuration, the improvement of TDOA/FDOA localization accuracy is especially obvious. The GDOP of TDOA/FDOA localization regards subastral point as center and uniformly distributed. But the TDOA localization of GDOP is unevenly distributed.
Figure 2. Y configuration, the GDOP distributions of TDOA localization (left) and TDOA/FDOA localization (right).
Figure 3. T configuration, the GDOP distributions of TDOA localization (left) and TDOA/FDOA localization (right).
Figure 4. Diamond configuration, the GDOP distributions of TDOA localization (left) and TDOA/FDOA localization (right).
Let measurement accuracy of absolute position in continuous 0 - 300 m and continuous value. Select the Y configuration, and the remaining parameters remain unchanged. Taking a point on the 100 km circumference of the subastral point, the relationship of TDOA localization accuracy and TDOA/FDOA localization accuracy with the absolute position measurement accuracy are obtained respectively, as shown in Figure 5 and Figure 6.
Comparing Figure 5 and Figure 6, the localization accuracy of the two algorithms will be improved with the increase of absolute position measurement
Figure 5. Relationship of TDOA localization accuracy with the absolute position measurement accuracy.
Figure 6. Relationship of TDOA/FDOA localization accuracy with the absolute position measurement accuracy.
accuracy. But the localization accuracy of TDOA is greatly affected by absolute position measurement accuracy. The localization accuracy of TDOA/FDOA is affected by absolute position measurement accuracy very little. Compared with TDOA localization, the effect is almost negligible.
The principles of four-satellite TDOA and TDOA/FDOA localization algorithms are analyzed, and the GDOP expressions of the two algorithms are derived respectively. Through simulation analysis, some useful conclusions can be drawn as follows.
With the same parameters, TDOA/FDOA localization algorithm shows better accuracy than TDOA localization algorithm.
When using the diamond configuration, as compared with the TDOA localization algorithm, the TDOA/FDOA localization algorithm not only has better localization accuracy, but also improves the GDOP distribution uniformity. It solves the problem of uneven distribution of diamond configuration and increases the engineering feasibility.
As compared with the TDOA localization, the absolute position measurement accuracy of satellite has almost no effect on the localization accuracy of TDOA/ FDOA, which reduces the localization error caused by unavoidable absolute position measuring factors.
This work is supported by the National Natural Science Foundation of China under Grant (No. 61201277, No. 61371184, No. 61671137, No. 61771114) and the Fundamental Research Funds for the Central Universities (No. ZYGX2016J028).
 Guo, X., Chu, L. and Sun, X. (2016) Accurate Localization of Multiple Sources Using Semidefinite Programming Based on Incomplete Range Matrix. IEEE Sensors Journal, 16, 5319-5324. https://doi.org/10.1109/JSEN.2016.2558184