gpsplot.txt

An Example Simulation of Differential GPS

% GPSPlot.m
% Author: Jim Connor ([e-mail])
% The scatter plot is then plotted and the graph is held so the CEP circle can be drawn


plot(east,north,'x'),title('Scatter Plot')
axis equal
hold;

% Plots Start and End points
plot(east(1),north(1),'m*')
text(east(1),north(1),'START')
plot(east(size(east)),north(size(north)),'m*')
text(east(length(time)),north(length(time)),'END')

% This code plots a red circle of radius r, previosly calculated, on top of the scatter plot
% There has to be a more efficient way to draw transparent circles in MATLAB! Let me know.

n = 50;
th = linspace(0,2*pi,n);
rr = r * ones(1,n);
[xx,yy] = pol2cart(th,rr);
plot(xx,yy,'r')

if R==1
    figure(2),...
    subplot(311),plot(east),title('East Coordinates'),...
    subplot(312),plot(north),title('North Coordinates'),...
    subplot(313),plot(height),title('Height')
elseif R==2
    figure(2),...
    subplot(411),plot(sat),title('Satellite Switching'),...
    subplot(412),plot(east),title('East Coordinates'),...
    subplot(413),plot(north),title('North Coordinates'),...
    subplot(414),plot(height),title('Height'),...

    figure(3),...
    subplot(211),plot(rt20stat),title('Correction Status - 0=best  8=worst'),...
    subplot(212),plot(stat),title('Differential Corrections - 2=DGPS  1=Normal')
else
   
   figure(2),...
    subplot(411),plot(sat),title('Satellite Switching'),...
    subplot(412),plot(east),title('East Coordinates'),...
    subplot(413),plot(north),title('North Coordinates'),...
    subplot(414),plot(height),title('Height'),...
    
   figure(3),...
    subplot(211),plot(rt20stat),title('Correction Status - 0=best  8=worst'),...
    subplot(212),plot(stat),title('Differential Corrections - 2=DGPS  1=Normal')

   
   figure(4)
	 subplot(311),plot(sats),title('# Satellites Used')
	 subplot(312),plot(gdop),title('Geometric Dilution of Precision')
	 subplot(313),plot(pdop),title('Position Dilution of Precision')

   figure(5)
    subplot(311),plot(htdop),title('Horizontal Position and Time Dilution of Precision')
    subplot(312),plot(hdop),title('Horizontal Dilution of Precision')
    subplot(313),plot(tdop),title('Time Dilution of Precision')
   
   
end