利用MATLAB编制的GPS单点定位程序.docVIP

function t=TimetoJD(Y,M,D,h,f,s) if(Y=80) Y=Y+1900; else Y=Y+2000; end if M=2 Y=Y-1; M=M+12; end JD=fix(365.25*Y)+fix(30.6001*(M+1))+D+h/24+f/1440+s/24/3600+1720981.5; t=JD-2444244.5; function [head,obs]=ReadObsData %读接收机观测数据文件 %HeadODat :a structure stores header information if o-file % .ApproXYZ[3]; //approximate coordinate % .interval; //intervals of two adjacent epochs % .SiteName; //point name % .Ant_H; //antenna height % .Ant_E; //east offset of the antenna center %