卫星导航定位算法与程序的设计_实验报告

2024 级测绘工程专业卫星导航定位算法与程序设计实验报告实验名称： 卫星导航基本程序设计 班 级：学 号：姓 名：实验时间： 2024 年 6 月 28 日 ~ 2024 年 6 月 3 0 中 国 矿 业 大 学目录实验一 时空基准转换 2一、实验目的 2二、实验容 2三、实验过程 2四、实验感想 6实验二 RINEX 文件读写 7一、实验目的 7二、实验容 7三、实验过程 7实验三 卫星轨道计算 12一、实验目的 12二、实验容 12三、实验过程 12四、实验感想 15实验一 时空基准转换一、实验目的1、加深对时空系统与其之间转换关系的理解2、掌握常用时空基准之间的转换模型与软件实现3、每人独立完成实验规定的容二、实验容本实验容包括：容一：编程实现 GPS 起点 1980 年 1 月 6 日 0 时对应的儒略日容二：编程实现 2024 年 11 月 27 日对应的 GPS 周数与一周的秒数容三：在 WGS84 椭球的条件下，编程实现当中央子午线为 117 度时，计算高斯坐标 x = 3548910.811290287, y = 179854.6172135982 对应的经纬度坐标？容四：WGS84 椭球下，表面 x=-2408000; y=4698000;z= 3566000 处的地平坐标系坐标为: e=704.8615;n=114.8683;u=751.9771 的点对应的直角坐标为多少？三、实验过程1.针对第一、二部分容：1.1 解决思路：先建立” TimeStruct.h”的头文件，将格里高利历、GPS时间结构、儒略日时间结构共结构体的方式放在里面；在建立“TimeTr”的头文件，建立类“CTimeTr”，创建变量“GPSTime”、“Time”、”JulDay”，并且申明函数“TIME2JUL”、“TIME2GTIME”等，用这些函数分别实现所需要的转换。1.2 具体的实现函数：“TIME2JUL”函数：double CTimeTr::TIME2JUL()//TIME Time,JULIANDAY &JulDay{double m,y;double D;//h =Time.byHour+Time.byMinute/60.0+Time.dSecond/3600.00;if(Time.byMonth<=2){y=Time.wYear-1;m=Time.byMonth+12;}else{y=Time.wYear;m=Time.byMonth;}D=floor(365.25*(y+4716))+floor(30.6001*(m+1))+Time.byDay+Time.byHour/24.0-1537.5;JulDay.lDay = int(D);JulDay.lSecond = D-int(JulDay.lDay);return 0;}“TIME2GTIME”：void CTimeTr::TIME2GTIME(){double JD;long m,y;int WN;double Wsecend;//UT=Time.byHour+Time.byMinute/60.0+Time.dSecond/3600.00;if(Time.byMonth<=2){y=Time.wYear-1;m=Time.byMonth+12;}else{y=Time.wYear;m=Time.byMonth...