ephplot.m

这是国外关于卫星导航方面一书的源代码

function []=ephplot(x,y,z,tsat,rela,typ,wie)
% -----------------------------------------------------------------------------
%  GPSLab - EPHPLOT.M  Darstellung der Orbits etc. aller im Rinex-NAV-File 
%                      vorhandenen Satelliten im WGS84 (erdfest).
% -----------------------------------------------------------------------------
%
%
% (c) iapg 1998 zeb
% -----------------------------------------------------------------------------

figure(99);
subplot(1,1,1);

switch typ
   
case 1
%--------------- GPS-Satellitenorbits im erdfesten System ----------------------
k=plot3(x,y,z);
% grid;
axis equal;
box;
title('GPS-Orbits in chosen system (all satellites from NAV-File)');
xlabel('meter');
ylabel('meter');
zlabel('meter');
view(3);
% colormap('gray');
% colormap('white');
colormap('summer');
hold on;
[a,b,c]=sphere(18);
surfl(a*6380000,b*6380000,c*6380000,[0,-23]);
hold off;

case 2
%-------------- Sicht auf 膓uatorebene --------------------------
k=plot3(x,y,z);
% grid;
box;
title('X-Y-plane of the GPS orbits (24h) in chosen system');
xlabel('meter');
ylabel('meter');
zlabel('meter');
view([0 0 1]);
% colormap('gray');
% colormap('white');
colormap('summer');
hold on;
[a,b,c]=sphere(18);
surfl(a*6380000,b*6380000,c*6380000,[0,-23]);
hold off;
axis equal;

   
case 3
%-------------- Sicht auf Meridianebene Greenwich --------------------------
k=plot3(x,y,z);
% grid;
box;
title('X-Z-plane of the GPS orbits (24h) in chosen system');
xlabel('meter');
ylabel('meter');
zlabel('meter');
view([0 1 0]);
% colormap('gray');
% colormap('white');
colormap('summer');
hold on;
[a,b,c]=sphere(18);
surfl(a*6380000,b*6380000,c*6380000,[0,-23]);
hold off;
axis equal;


case 4
%--------------- GPS-Satellitenspuren zur Beobachtungszeit ----------------------
if wie==0; 
   warndlg(['The option "groundtracks" is nonsense for ECI representation.' ...
         ' Calculate the orbits in the ECEF System at first.'] ,...
         'GPS Lab: Break at the 4th option of EPHPLOT.M');
   return;
end 
   
load coast0;		% K黶tenlinien-Daten laden
[b,l,h]=xyz2blh(x,y,z);
colormap('default');
plot(lam,phi,'k'); 		% liefert Weltkarte
hold on;		% Plot "黚ermalen"
% hg=plot(l,b,'p');
hg=plot(l,b,'.');
% set(hg,'LineWidth',5,'MarkerSize',5);   % nur f黵 Liniendarst. in '.'
xlabel('L鋘ge');
ylabel('Breite');
title('GPS - groundtracks for 24 hours, every 10 minutes');
text(-180,-90,'all satellites from NAV-File');
hold off;		% "躡ermalen" aus
   
case 5
%--------------- relativistischer Effekt wegen Bahnexzentrizit鋞 ----------------------
plot(rela);
title('relativistic effect due to orbit excentricity (all satellites from NAV-File)');
ylabel('meter');
xlabel('24 hours, every 10 minutes (144 values)');
   
   
case 6
%--------------- Bahnh鰄e der GPS-Satelliten ----------------------
[b,l,h]=xyz2blh(x,y,z);
plot(h);
title('orbit height above WGS84-Ellipsoid (all satellites from NAV-File)');
ylabel('meter');
xlabel('24 hours, every 10 minutes (144 values)');
   
   
case 7
%--------------- Anzahl von Ephemeriden-Sets ----------------------
load s1eph;
satnrs=s1eph(1:24:length(s1eph));
hist(satnrs,[1:32]);
colormap cool;
title('ephemeris sets of GPS satellites in NAV-File');
xlabel('PRN-number of the satellites');
ylabel('number of ephemeris sets depending on satellite');
grid;


case 8
subplot(2,1,1);
[yydummy,ii]=max(max(abs(tsat)));
clear yydummy;
plot(tsat(:,ii).*1000000,'.');
title('maximal occurring satellite clock error (absolutes)');
%xlabel('24 Stunden im 10-Minutentakt (144 Werte)');
ylabel('祍ec');
grid;

subplot(2,1,2);
[yydummy,jj]=min(min(abs(tsat)));
clear yydummy;
plot(tsat(:,jj).*1000000,'.');
title('minimal occurring satellite clock error (absolutes)');
xlabel('for each satellite during 24 hours, every 10 minutes (144 values)');
ylabel('祍ec');
grid;


otherwise
   
   warndlg(['unvalid variable "typ": value of "typ" has to be ' ...
         ' 1,2,3,4,5,6,7 or 8 !'], 'error during call of EPHPLOT.M');
   return;
   
end   % end switch   
   
%---------------------------- Ende --------------------------

% -----------------------------------------------------------------------------
% (c) iapg 1998 zeb
% -----------------------------------------------------------------------------