xcon.m

GPS TOOLBOX包含以下内容： 1、GPS相关常量和转换因子； 2、角度变换； 3、坐标系转换： &#61656 点变换； &#61656 矩阵变换； &#61656 向量变换

%                                xcon.m
%  Scope:   This MATLAB program displays the constants specified by the 
%           macros convcon, gpscon, wgs72con, wgs84con, sgs85con, and pz90con.
%           Each macros set the constants as global variables.
%  Usage:   xcon
%  Inputs:   - selection of the set of constants (see the displayed menu)
%  Outputs:  - displayed constants
%  External Matlab macros used: convcon, gpscon, pz90con, sgs85con, wgs72con,
%                               wgs84con
%  Last update: 04/06/00
%  Copyright (C) 1996-00 by LL Consulting. All Rights Reserved.

clear

%  Initialization

mtitle  = 'Make a selection';
sel1 = 'Display the conversion constants';
sel2 = 'Display the GPS constants';
sel3 = 'Display the WGS-72 constants';
sel4 = 'Display the WGS-84 constants';
sel5 = 'Display the SGS-85 constants';
sel6 = 'Display the PZ-90 constants';
sel7 = 'Exit program';

disp(' ');
answer = 'y';
yes = 'y';
while (strcmp(answer,yes) == 1)
   comp = menu(mtitle,sel1,sel2,sel3,sel4,sel5,sel6,sel7);

   if comp == 1

%  Display the conversion constants from macro  convcon

      convcon;

      disp('Conversion constants from macro convcon')
      disp(' ');
      fprintf('sem_rad   =  %16.10e \n',sem_rad  );  
      fprintf('rad_sem   =  %16.10e \n',rad_sem  );  
      fprintf('deg_rad   =  %16.10e \n',deg_rad  );  
      fprintf('rad_deg   =  %16.10e \n',rad_deg  ); 
      fprintf('deg_mrad  =  %16.10e \n',deg_mrad );  
      fprintf('mrad_deg  =  %16.10e \n',mrad_deg ); 
      fprintf('ussf_m    =  %16.10e \n',ussf_m   );  
      fprintf('m_ussf    =  %16.10e \n',m_ussf   );  
      fprintf('if_m      =  %16.10e \n',if_m     );  
      fprintf('m_if      =  %16.10e \n',m_if     );  
      fprintf('nmile_m   =  %16.10e \n',nmile_m  );  
      fprintf('m_nmile   =  %16.10e \n',m_nmile  );  
      fprintf('smile_m   =  %16.10e \n',smile_m  );  
      fprintf('m_smile   =  %16.10e \n',m_smile  );  
      fprintf('smile_if  =  %16.10e \n',smile_if ); 
      fprintf('if_smile  =  %16.10e \n',if_smile ); 
        
   elseif comp == 2

%  Display the gps constants from macro  gpscon

      gpscon;
      
      disp('GPS constants from macro gpscon')
      disp(' ');
      fprintf('f0           =  %16.10e \n',f0          );    
      fprintf('fL1          =  %16.10e \n',fL1         );    
      fprintf('fL2          =  %16.10e \n',fL2         );    
      fprintf('lambdaL1     =  %16.10e \n',lambdaL1    );    
      fprintf('lambdaL2     =  %16.10e \n',lambdaL2    );    
      fprintf('fca_code     =  %16.10e \n',fca_code    );    
      fprintf('fp_code      =  %16.10e \n',fp_code     );  
      fprintf('fL1pL2       =  %16.10e \n',fL1pL2      );    
      fprintf('fL2mL2       =  %16.10e \n',fL1mL2      );    
      fprintf('lambdaL1pL2  =  %16.10e \n',lambdaL1pL2 );    
      fprintf('lambdaL1mL2  =  %16.10e \n',lambdaL1mL2 );    
      fprintf('secpmin      =  %16.10e \n',secpmin     );    
      fprintf('secph        =  %16.10e \n',secph       );    
      fprintf('secpday      =  %16.10e \n',secpday     );    
      fprintf('secpweek     =  %16.10e \n',secpweek    );  
      fprintf('secphday     =  %16.10e \n',secphday    );    
      fprintf('secphweek    =  %16.10e \n',secphweek   );  

   elseif comp == 3

%  Display the wgs72 constants from macro  wgs72con

      wgs72con;

      disp('WGS-72 constants from macro wgs72con')
      disp(' ');
      fprintf('a_smaxis  =  %16.10e \n',a_smaxis );  
      fprintf('b_smaxis  =  %16.10e \n',b_smaxis );  
      fprintf('eccentr   =  %16.10e \n',eccentr  );  
      fprintf('eccentr2  =  %16.10e \n',eccentr2 );  
      fprintf('flatness  =  %16.10e \n',flatness );  
      fprintf('gravpar   =  %16.10e \n',gravpar  );  
      fprintf('rot_rate  =  %16.10e \n',rot_rate );  
      fprintf('c_speed   =  %16.10e \n',c_speed  );  
      fprintf('ucgrav    =  %16.10e \n',ucgrav   );  
      fprintf('mearth    =  %16.10e \n',mearth   );  
      fprintf('g0        =  %16.10e \n',g0       );  

   elseif comp == 4

%  Display the wgs84 constants from macro  wgs84con

      wgs84con;

      disp('WGS-84 constants from macro wgs84con')
      disp(' ');
      fprintf('a_smaxis  =  %16.10e \n',a_smaxis ); 
      fprintf('b_smaxis  =  %16.10e \n',b_smaxis ); 
      fprintf('eccentr   =  %16.10e \n',eccentr  ); 
      fprintf('eccentr2  =  %16.10e \n',eccentr2 ); 
      fprintf('onemecc2  =  %16.10e \n',onemecc2 ); 
      fprintf('eprime    =  %16.10e \n',eprime   ); 
      fprintf('eprime2   =  %16.10e \n',eprime2  ); 
      fprintf('flatness  =  %16.10e \n',flatness ); 
      fprintf('gravpar   =  %16.10e \n',gravpar  ); 
      fprintf('rot_rate  =  %16.10e \n',rot_rate ); 
      fprintf('c_speed   =  %16.10e \n',c_speed  ); 
      fprintf('ucgrav    =  %16.10e \n',ucgrav   ); 
      fprintf('mearth    =  %16.10e \n',mearth   ); 
      fprintf('g0        =  %16.10e \n',g0       ); 
      fprintf('c20       =  %16.10e \n',c20      ); 

elseif comp == 5

%  Display the SGS-85 constants from macro  sgs85con

      sgs85con;

      disp('SGS-85 constants from macro sgs85con')
      disp(' ');
      fprintf('a_smaxis  =  %16.10e \n',a_smaxis ); 
      fprintf('b_smaxis  =  %16.10e \n',b_smaxis ); 
      fprintf('eccentr   =  %16.10e \n',eccentr  ); 
      fprintf('eccentr2  =  %16.10e \n',eccentr2 ); 
      fprintf('flatness  =  %16.10e \n',flatness ); 
      fprintf('gravpar   =  %16.10e \n',gravpar  ); 
      fprintf('gravatm   =  %16.10e \n',gravatm  ); 
      fprintf('rot_rate  =  %16.10e \n',rot_rate ); 
      fprintf('c_speed   =  %16.10e \n',c_speed  ); 
      fprintf('g0        =  %16.10e \n',g0       ); 
      fprintf('gcoratm   =  %16.10e \n',gcoratm  ); 
      fprintf('j2        =  %16.10e \n',j2       ); 

elseif comp == 6

%  Display the PZ-90 constants from macro  pz90con

      pz90con;

      disp('PZ-90 constants from macro pz90con')
      disp(' ');
      fprintf('a_smaxis  =  %16.10e \n',a_smaxis ); 
      fprintf('b_smaxis  =  %16.10e \n',b_smaxis ); 
      fprintf('eccentr   =  %16.10e \n',eccentr  ); 
      fprintf('eccentr2  =  %16.10e \n',eccentr2 ); 
      fprintf('flatness  =  %16.10e \n',flatness ); 
      fprintf('gravpar   =  %16.10e \n',gravpar  );
      fprintf('gravatm   =  %16.10e \n',gravatm  ); 
      fprintf('rot_rate  =  %16.10e \n',rot_rate ); 
      fprintf('c_speed   =  %16.10e \n',c_speed  ); 
      fprintf('g0        =  %16.10e \n',g0       ); 
      fprintf('gcoratm   =  %16.10e \n',gcoratm  ); 
      fprintf('j2        =  %16.10e \n',j2       ); 

   elseif comp == 7

%  Exit the program

      return

   else
      error('Error XCON ;  check the index of the selection');
   end

%  Check if another selection is desired

   disp('  ');
   answer = input('Do you want to make another selection? (y/n) [y] ','s');
   if  isempty(answer)
      answer = 'y';
   end
   disp('  ');

end

disp('End of the program  XCON');
disp('  ');