gps算法考试重点整理2 .docxVIP

时间转换主程序： clc; %檫除工作窗中的所有显示内容 %% % *BOLD TEXT* clear; %清除内存中的变量和函数 close all; format long %输出数据为长格式 fprintf(格里历时(Gregorian Calendar)转换为儒略时(Julian Day))% 输出字 符串 jd = julday(2001,9,4,9.6667)%调用julday函数，格里历时转换为儒略时 fprintf(儒略时转换为格里历时) [y,m,d,h]=GreCal(jd)%儒略时转换为格里高利时 fprintf(儒略时转换为GPS周和周内秒) [week,sow] = gps_time(jd)%儒略历转换为 GPS 时 fprintf(GPS时转换为儒略时) JD=week*7+sow/86400+2444244.5 %GPS 时转换为儒略时 % 2444244.5为GPS时间零点1980年1月6日子夜对应的儒略日 %%%%%%%%%%%%%%%%%%%%% end easy1.m %%%%%%%%%%%%%%% %计算2011年7月4日0时0秒的年积日 fprintf(年积日与格里历时间的相互转换)； %计算出当年1月1日的儒略日 jd1=julday(2011,1,1,0); %计算出日历时间的儒略日 jd2=julday(2011,7,4,0); %两个儒略日求差加1，得出年积日 fprintf(两个儒略日求差加1，得出年积日’)； DOY=jd2-jd1+1 %年积日加JD1减1得儒略日 fprintf(年积日加JD1减1得儒略日’)； JD2=jd1+DOY-1 运行结果： 格里历时(Gregorian Calendar)转换为儒略时(Julian Day) jd = 2.452156902779167e+006 儒略时转换为格里历时 y = 2001 m = 9 d = 4 h = 9.66669999808073 儒略时转换为GPS周和周内秒 week = 1130 sow = 2.076001199930906e+005 GPS时转换为儒略时 JD = 2.452156902779167e+006 年积日与格里历时间的相互转换两个儒略日求差加1，得出年积日 DOY = 185 年积日加JD1减1得儒略日 JD2 = 2.455746500000000e+006 大地坐标转直角坐标: function xyz = llh2xyz(llh) %大地坐标转三维直角坐标 phi = llh(1);%大地纬度，弧度表示 lambda = llh(2);%大地经度，弧度表示 h = llh(3);% 大地高 %设置椭球参数 a = 6378137.0000; % earth semimajor axis in meters b = 6356752.3142; % earth semiminor axis in meters e = sqrt (1-(b/a).”2);% 第一偏心率 %以下为求三角函数 sinphi = sin(phi); cosphi = cos(phi); coslam = cos(lambda); sinlam = sin(lambda); tan2phi = (tan(phi))”2; tmp = 1 - e*e; tmpden = sqrt( 1 + tmp*tan2phi ); x = (a*coslam)/tmpden + h*coslam*cosphi; %三维大地直角X坐标 y = (a*sinlam)/tmpden + h*sinlam*cosphi; %三维大地直角Y坐标 tmp2 = sqrt(1 - e*e*sinphi*sinphi); z = (a*tmp*sinphi)/tmp2 + h*sinphi; %三维大地直角Z坐标 %将直角坐标赋给数组xyz xyz(1) = x; xyz(2) = y; xyz(3) = z; 站心直角坐标转空间直角坐标： function xyz = enu2xyz(enu,orgxyz) %站心直角坐标转换空间直角坐标 if nargin2,error(insufficient number of input arguments),end [m,n]=size(enu);if mn,tmpenu=enu;else,tmpenu=enu;end %如果enu矩阵的行数小于列数，将矩阵转置 [m,n]=size(orgxyz);if mn,tmpxyz=orgxyz;else,tmpxyz=orgxyz;end %如果orgxyz矩阵的行数小于列数，将矩阵转置 orgllh = xyz2llh(tmpxyz); %将站心直