demon1.m

根据卫星星历及测者位置

plot3(0,0,0);
xlabel('x');
ylabel('y');
zlabel('z');
axis([-3 3 -3 3 -3 3]);
hold on
Delta_t = 360;
for t=0:(Delta_t):12*3600    
    p=SatellitePos(getGPSSatellite( 1 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3);
    %plot(p1,p2);
    hold on;
    pause(0.01);
end
for t=0:(Delta_t):12*3600    
    p=SatellitePos(getGPSSatellite( 5 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3,'k');  
    hold on;
    pause(0.01);
end
for t=0:(Delta_t):12*3600     
    p=SatellitePos(getGPSSatellite( 9 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3,'g');
    hold on;
    pause(0.01);
end
%
for t=0:(Delta_t):12*3600    
    p=SatellitePos(getGPSSatellite( 13 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3);
    %plot(p1,p2);
    hold on;
    pause(0.01);
end
for t=0:(Delta_t):12*3600    
    p=SatellitePos(getGPSSatellite( 17 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3,'y');  
    hold on;
    pause(0.01);
end
for t=0:(Delta_t):12*3600     
    p=SatellitePos(getGPSSatellite( 21 ),pi/3,t);
    p1=p(1)/1.0e7;
    p2=p(2)/1.0e7;
    p3=p(3)/1.0e7;
    plot3(p1,p2,p3,'r');
    hold on;
    pause(0.01);
end