卫星导航定位算法与程序设计——实验报告_2.doc

2013 级测绘工程专业 卫星导航定位算法与程序设计 实 验 报 告 实验名称： 卫星导航基本程序设计 班 级： 学 号： 姓 名： 实验时间： 2016年6月28日~2016年6月30 中 国 矿 业 大 学 目录 TOC \o 1-3 \h \z \u 实验一 时空基准转换 2 一、实验目的 2 二、实验内容 2 三、实验过程 2 四、实验感想 6 实验二 RINEX文件读写 7 一、实验目的 7 二、实验内容 7 三、实验过程 7 实验三 卫星轨道计算 12 一、实验目的 12 二、实验内容 12 三、实验过程 12 四、实验感想 15 实验一 时空基准转换 一、实验目的 1、加深对时空系统及其之间转换关系的理解 2、掌握常用时空基准之间的转换模型与软件实现 3、每人独立完成实验规定的内容 二、实验内容 本实验内容包括： 内容一：编程实现GPS起点1980年1月6日0时对应的儒略日 内容二：编程实现2011年11月27日对应的GPS周数与一周内的秒数 内容三：在WGS84椭球的条件下，编程实现当中央子午线为117度时，计算高斯坐标x = , y = 对应的经纬度坐标？ 内容四：WGS84椭球下，表面x=-2408000; y=4698000;z= 3566000处的地平坐标系坐标为: e=;n=;u=的点对应的直角坐标为多少？ 三、实验过程 1.针对第一、二部分内容： 解决思路：先建立” ”的头文件，将格里高利历、GPS时间结构、儒略日时间结构共结构体的方式放在里面；在建立“TimeTr”的头文件，建立类“CTimeTr”，创建变量“GPSTime”、“Time”、”JulDay”，并且申明函数“TIME2JUL”、“TIME2GTIME”等，用这些函数分别实现所需要的转换。 具体的实现函数： “TIME2JUL”函数： double CTimeTr::TIME2JUL()XT”文件即可。 、主要函数代码: “get_eph”函数： function eph = get_eph(ephemeridesfile) fide = fopen(ephemeridesfile); [eph, count] = fread(fide, Inf, double); noeph = count/22; eph = reshape(eph, 22, noeph); “rinexe”函数（部分）： function rinexe(ephemerisfile, outputfile) fide = fopen(ephemerisfile); head_lines = 0; while 1 head_lines = head_lines+1; line = fgetl(fide); answer = findstr(line,END OF HEADER); if ~isempty(answer), break; end; end; head_lines 主函数中输出结果得函数部分： af0=data(19);%卫星中差 M0=data(3); roota=data(4); deltan=data(5); ecc=data(6); omega=data(7); cuc=data(8); cus=data(9); crc=data(10); crs=data(11); i0=data(12); idot=data(13); toe=data(18); af1=data(20);%对所要输出的参数赋值 fprintf(fid,\n卫星编号：%d\n卫星钟差：%d\n平近点角距：%d\n轨道长半轴的平方根：%d\n平均运动修正量：%d\n轨道偏心率：%d\n近地点角距：%d\n纬度幅角的余弦调和项改正的振幅,prn,af0,M0,roota,deltan,ecc,omega,cuc); fprintf(fid,纬度幅角的正弦调和项改正的振幅：%d\n轨道半径的余弦调和项改正的振幅：%d\n轨道半径的正弦调和项改正的振幅：%d\n轨道倾角：%d\n轨道倾角变化率：%d\n星历参考时刻：%d\n,cus,crc,crs,i0,idot,toe) fclose(fid); 、输出结果 2、针对第二部分内容： 编程实现O文件读入，并采用中文标注出主要