MATLAB: Plotting in matlab using while loop

debugplotting

HELLO EVERYONE! i have a simple question regarding plotting in MATLAB using while loop.plzz help me out. my problem is that when i plot using while loop, my green color is coming first and then it is changing and then my blue color is displayed.i want to know what command is making it happen??
CODE: 
clear clc close
%% Receivers coordinates 
iter = 0;
lat_tgt = 34.0151;
long_tgt = 71.5249;
z_s = 0.01/375;
while iter < 100 
    iter=iter+1;
    %lat_rand = rand/10;
    long_rand = rand/10;
    lat_rand = (-0.5 + (0.5+0.5)*rand)/10;
    long_rand = (-0.5 + (0.5+0.5)*rand)/10;
    lat_tgt = lat_tgt+lat_rand; long_tgt = long_tgt+long_rand;
    lats = [34.1989 34.0105 34.067894 34.1166 lat_tgt]; longs = [72.0231 71.9876 71.992783 72.0216 long_tgt];
    [x,y] = grn2eqa(lats,longs,[34.1166, 72.0216]);
    % x=[8 0 -8 0 ].*100; %%x=[x_1 x_2 x_3 x_p]
    %x(5) = -x(5); y(5) = -y(5);
    % y=[-4 8 -4 5].*100;  %%y=[y_1 y_2 y_3 y_p]
    z=[0 0 0 0];   %%z=[z_1 z_2 z_3 z_p]
    c=2.997924580*10^8;
    
    %%Source TDOA calculation
    z_s=z_s+(-0.5 + (0.5+0.5)*rand)/37500;
    t1 = (sqrt((x(5)-x(4))^2+(y(5)-y(4))^2+(z_s-z(4))^2)-sqrt((x(5)-x(1))^2+(y(5)-y(1))^2+(z_s-z(1))^2))/c;
    t2 = (sqrt((x(5)-x(4))^2+(y(5)-y(4))^2+(z_s-z(4))^2)-sqrt((x(5)-x(2))^2+(y(5)-y(2))^2+(z_s-z(2))^2))/c;
    t3 = (sqrt((x(5)-x(4))^2+(y(5)-y(4))^2+(z_s-z(4))^2)-sqrt((x(5)-x(3))^2+(y(5)-y(3))^2+(z_s-z(3))^2))/c;
    
    %%Source localization
    syms xs ys zs     %our unknowns
    eqn1 = sqrt((xs-x(4))^2+(ys-y(4))^2+(zs-z(4))^2)-sqrt((xs-x(1))^2+(ys-y(1))^2+(zs-z(1))^2)-(c*t1)==0;
    eqn2 = sqrt((xs-x(4))^2+(ys-y(4))^2+(zs-z(4))^2)-sqrt((xs-x(2))^2+(ys-y(2))^2+(zs-z(2))^2)-(c*t2)==0;
    eqn3 = sqrt((xs-x(4))^2+(ys-y(4))^2+(zs-z(4))^2)-sqrt((xs-x(3))^2+(ys-y(3))^2+(zs-z(3))^2)-(c*t3)==0;
    sol = solve([eqn1, eqn2, eqn3], [xs ys zs]);
    
    %%Debug
    % figure(1)
    % fimplicit3(eqn1,[-0.008 0.001 -0.0020 0.002 -0.01 0.02],'FaceAlpha',0.5,'LineStyle','none')
    % hold on


    % fimplicit3(eqn2,[-0.008 0.001 -0.0020 0.002 -0.01 0.02],'FaceAlpha',0.5,'LineStyle','none')
    % hold on
    % fimplicit3(eqn3,[-0.008 0.001 -0.0020 0.002 -0.01 0.02],'FaceAlpha',0.5,'LineStyle','none')
    %%Debug End
    
    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
    %%Filtering Results
    
    %idx = any(possibleSol < 0,1) | any(imag(possibleSol) ~=0) %possibleSol(3,:) < 0 | any(imag(possibleSol) ~=0)
    idx = possibleSol(3,:) < 0 | any(imag(possibleSol) ~=0);
    possibleSol(:, idx) = [];
    [lat,long] = eqa2grn(possibleSol(1),possibleSol(2),[34.1166, 72.0216]);
    
    %Plotting on 3D coordinates
    % fprintf('Xs=%g\n', possibleSol(1,1));
    % fprintf('Ys=%g\n', possibleSol(2,1));
    % fprintf('Zs=%g\n', possibleSol(3,1));
    % scatter3(x, y, z, 'o');
    % hold on
    % scatter3(possibleSol(1,1), possibleSol(2,1), possibleSol(3,1), '+');
    % hold off
    % legend({'Receivers', 'Source'})
    posSolMoving(:,iter) = possibleSol;
    x_s=x(length(x)); y_s=y(length(y)); x = x(1:length(x)-1); y = y(1:length(y)-1);
    %fprintf('Xs=%.f\nYs=%.f\nZs=%.f\n', possibleSol);
    figure(1)
    hold on
    grid on
    view(3);
    plot3(x,y,z, 'ro', 'LineWidth', 2, 'MarkerSize', 10);
    plot3(posSolMoving(1,:),posSolMoving(2,:),posSolMoving(3,:), 'b+', 'LineWidth', 4, 'MarkerSize', 4)
    plot3(x_s,y_s,z_s,'g+', 'LineWidth', 2, 'MarkerSize', 4)
    if iter > 1
        q = quiver3(posSolMoving(1,iter-1),posSolMoving(2,iter-1),posSolMoving(3,iter-1),posSolMoving(1,iter)-posSolMoving(1,iter-1),posSolMoving(2,iter)-posSolMoving(2,iter-1),posSolMoving(3,iter)-posSolMoving(3,iter-1),0);
        q.LineWidth=0.75;
    end
    legend({'Receivers', 'Source','exactvalue'})
    view(-14.3,33.2)
    pause(1)
end
output:

Best Answer

The answer ought to be, step throw the code and observe what each line of code makes to your diagram. However, the problem is that the figure window constantly disappears under other windows. There are to ways to keep the figure window visible.
I've modified your code somewhat to better see what's going on. 
        fh = figure(1);
        WindowAPI( fh,'Topmost')    % <<<<<<<<<<<<<<<<<
        hold on
        grid on
        view(3);
        plot3(x,y,z, 'ro', 'LineWidth', 2, 'MarkerSize', 10);
        plot3(posSolMoving(1,:),posSolMoving(2,:),posSolMoving(3,:), 'b+'   ...
            , 'LineWidth', 2, 'MarkerSize', 6)
        plot3(x_s,y_s,z_s,'go', 'LineWidth', 2, 'MarkerSize', 8)
In the first iteration (iter==1)
  • The first plot3 draws four red circles
  • The second plot3 draws a blue cross and changes the axis limits.
  • The third plot3 adds (in my case) a green circle on top of the blue cross with . "my green color is coming first and then it is changing" no it's hidden by the green marker that is added on top.
Now it's your turn to step throw the remaining 99 iterations.
Related Question