# plotting of graph/TDOA localization

20 views (last 30 days)
ali hassan on 10 Oct 2020
HELLO EVERYONEONE
hope everyone is fine and doing great.just need some help guys.
this is my code based on TDOA(Time Difference Of Arrival) and it localizes an RF transmitter which is the output in 3D coordinates.now i want to plot the accuracy of my algorithm on matlab.i want to plot a graph of accuracy vs an array of increasing values from current values and increasing 150km.
CODE:
x_p = 4000; y_p = 5000; z_p = 20;
x_1 = 8000; y_1 = 9000; z_1 = 50;
x_2 = 2000; y_2 = 5000; z_2 = 10;
x_3 = 6000; y_3 = 1000; z_3 = 30;
c=3.0*10^8;
t1 = 2.923281245*10^-6;
t2 = -1.284077874*10^-6;
t3 = -1.40649205*10^-5;
syms xs ys zs %our unknowns
eqn1 = sqrt((xs-x_p)^2+(ys-y_p)^2+(zs-z_p)^2)-sqrt((xs-x_1)^2+(ys-y_1)^2+(zs-z_2)^2)-(c*t1);
eqn2 = sqrt((xs-x_p)^2+(ys-y_p)^2+(zs-z_p)^2)-sqrt((xs-x_2)^2+(ys-y_2)^2+(zs-z_2)^2)-(c*t2);
eqn3 = sqrt((xs-x_p)^2+(ys-y_p)^2+(zs-z_p)^2)-sqrt((xs-x_3)^2+(ys-y_3)^2+(zs-z_3)^2)-(c*t3);
sol = solve([eqn1, eqn2, eqn3], [xs ys zs]);
m = 1;
for n = 1:length(sol.xs)
possibleSol(1,m) = double(sol.xs(n));
possibleSol(2,m) = double(sol.ys(n));
possibleSol(3,m) = double(sol.zs(n));
m=m+1
end
idx = any(possibleSol < 0) | any(imag(possibleSol) ~=0);
possibleSol(:, idx) = []
OUTPUT:
BEST REGARDS