erg2code.m

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

% -----------------------------------------------------------------------------
% GPSLab:      ERG2CODE.M
%
% Save results after running of DGPSLAB.M:
% mean values, observation statistics    [lst2code.txt]
%
% -----------------------------------------------------------------------------
% (c) iapg 1998-2001 zeb


if ~exist('dgps_is_calc')
   %warndlg(['Ergebnisdatei nicht erstellbar, solange kein DGPS berechnet wurde: '...
   %      'Bitte f黨ren Sie zuerst "DGPS"  -> "berechnen" aus.'] ,...
   %      'GPS Lab: Anwendungsfehler beim Aufruf von ERG2CODE.M');
   return;
end 


fid=fopen('lst2code.txt','wt+');
a1='#########################################################################';
a2='       Results of the Pseudorange Processing (double-difference-DGPS)';
a3='                        GPSLab (c) iapg 2001 zeb';
a4='################################ END ####################################';
a5=['            This file was created ',datestr(datenum(now),0)];
a6='-------------------------------------------------------------------------';
a7='------------------------------------';
a10='lst2dgps.txt (DGPS from dbl.diff.PR): GPS-Wsec, DX, DY, DZ, Bsllength, Dh';
a13='Units                                   [sec]   [m] [m] [m]    [m]     [m]';

fprintf(fid,'%s\n%s\n%s\n%s\n%s\n%s\n\n',a1,a2,a3,a1,a5,a6);

a8='Reference station                   X [m]        Y [m]        Z [m]';
a9='Rover station                       X [m]        Y [m]        Z [m]';
ag8='Reference station (WGS84 ellipsoid) Lat [癩      Lon [癩      H [m]';
ag9='Rover station     (WGS84 ellipsoid) Lat [癩      Lon [癩      H [m]';

if exist('mit_p2')
if mit_p2==1
fprintf(fid,'%s  %s\n%s\n\n','Processing model   (P2 existing) :',what_modl,a6);
else
fprintf(fid,'%s  %s\n%s\n\n','Processing model                 :',what_modl,a6);
end
end

fprintf(fid,'%s\n%s\n%s%12.3f %12.3f %12.3f\n\n',a8,a6,'approximate coordinates        :',xs1,ys1,zs1);
fprintf(fid,'%s\n%s\n%s%12.3f %12.3f %12.3f\n',a9,a6,'approximate coordinates        :',xs2,ys2,zs2);
fprintf(fid,'%s%12.3f %12.3f %12.3f\n','estimated coordinates    (mean):',xxs2,yys2,zzs2);
fprintf(fid,'%s%12.3f %12.3f %12.3f\n\n\n','standard deviations      (mean):',sxxs2,syys2,szzs2);

fprintf(fid,'%s\n%s\n%s%12.8f %12.8f %12.3f\n\n',ag8,a6,'approximate coordinates        :',app1_b,app1_l,app1_h);
fprintf(fid,'%s\n%s\n%s%12.8f %12.8f %12.3f\n',ag9,a6,'approximate coordinates        :',app2_b,app2_l,app2_h);
fprintf(fid,'%s%12.8f %12.8f %12.3f\n','estimated coordinates    (mean):',bbs2,lls2,hhs2);
fprintf(fid,'%s%12.3f %12.3f %12.3f\n\n\n','standard deviations      (mean):',sbbs2,slls2,shhs2);

fprintf(fid,'%s %10.3f\n\n','error of unit weight                (Mittel) [m]:',sig0);

fprintf(fid,'%s %10.3f\n','height difference from approximate coordin.  [m]:',dapp_h);
fprintf(fid,'%s %10.3f\n\n','height difference from DGPS                  [m]:',m_delta_h);
fprintf(fid,'%s %10.3f\n','baseline length from approximate coordinates [m]:',app_bsl);
fprintf(fid,'%s %10.3f\n\n','baseline length from DGPS                    [m]:',m_bsl);

fprintf(fid,'%s %6.3f %s %6.3f %s %6.3f\n\n','RDOP of the DGPS solution - mean:',rdmit,'   min:',rdmin,'   max:',rdmax);

fprintf(fid,'%s %10.3f\n','antenna height reference stat.(phase center) [m]:',ants1);
fprintf(fid,'%s %10.3f\n\n','antenna height at rover stat. (phase center) [m]:',ants2);

fprintf(fid,'\n%s\n%s\n%s\n%s\n','simultaneous observation epochs: all satellites at rover ',a7,'PRN      #',a7);
fprintf(fid,'%3i  %5i\n',[n2;haufig2]);

t1st=s2epo(1,1);
tlast=s2epo(length(s2epo),1);
tspan=tlast-t1st;
tsec=mod(tspan,60);
tmin=(tspan-tsec)/60;
fprintf(fid,'\n%s %6i\n','first  epoch    [GPS week seconds]:  ',t1st);
fprintf(fid,'%s %6i\n','last epoch      [GPS week seconds]:  ',tlast);
fprintf(fid,'%s %6i %s %3i %s\n','time span (processable obs. time) :  ',tmin,' minutes  ',tsec,' seconds ');


fprintf(fid,'\n%s\n%s\n%s\n%s\n%s\n\n',a1,a10,a13,a6,a4);

fclose(fid);

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