A Convergence issue in this Least Squares calculation

AI Thread Summary
The discussion centers on a Least Squares calculation for computing the trajectory of a moving body using DTOA measurements from five stations. The user experiences convergence issues depending on the initialization of the starting point; convergence is successful when initialized within the polygon formed by the stations, but fails when initialized outside. This raises questions about the influence of the initial guess on the iterative process and the potential need for additional data, such as iterative RMS error data, to diagnose the problem. Participants are encouraged to share insights or suggestions for improving the initialization process. The conversation highlights the challenges of achieving convergence in iterative algorithms for trajectory estimation.
ChiPi
Messages
2
Reaction score
0
TL;DR Summary
I am not able to reach convergence if I initialize the starting point with coordinates outside the stations net. The observations are Difference in Time of Arrival.
I'm computing the trajectory of a moving body and my net is composed by 5 stations.
My observations are DTOA: difference in time of Arrival (they have been linearized).
I am trying to use Least Squares with a linear model: Y = Ax + b, where Y are the observed measurements (DTOA), A the design matrix, b is the known terms vector and X is a vector of estimates. The algorithm processes data according to the epoch considered and iterate the process up to a value of 20 times, unless it reaches before a 1mm convergence.

Since it is an iterative process, the system requires an initialization at the starting point with approximated values for the unknowns, and here comes my problem: if I initialize the starting point with coordinates within the polygon formed by the five stations convergence is reached and the solution is successful, but if I initialize the point with coordinates outside the area formed by the stations the convergence is not reached and I’m not able to determine the trajectory.

Does anyone of you have an explanation for this?
 
Physics news on Phys.org
I have a pretty good idea but you should figure this out. Do you have access to the iterative rms error data? Does is say anything to you?
 
hutchphd said:
I have a pretty good idea but you should figure this out. Do you have access to the iterative rms error data? Does is say anything to you?
I don't have access to it but I can add it at the code if it can help in figuring out the problem. What was your idea?
 
Back
Top