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

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

ID:48210242

大小:190.00 KB

页数:16页

时间:2020-01-22

卫星导航定位算法与程序设计——实验报告.doc_第1页
卫星导航定位算法与程序设计——实验报告.doc_第2页
卫星导航定位算法与程序设计——实验报告.doc_第3页
卫星导航定位算法与程序设计——实验报告.doc_第4页
卫星导航定位算法与程序设计——实验报告.doc_第5页
资源描述:

《卫星导航定位算法与程序设计——实验报告.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、.2013级测绘工程专业卫星导航定位算法与程序设计实验报告实验名称:卫星导航基本程序设计班级:学号:姓名:实验时间:2016年6月28日~2016年6月30中国矿业大学..目录实验一时空基准转换2一、实验目的2二、实验内容2三、实验过程2四、实验感想2实验二RINEX文件读写2一、实验目的2二、实验内容2三、实验过程2实验三卫星轨道计算2一、实验目的2二、实验内容2三、实验过程2四、实验感想2..实验一时空基准转换一、实验目的1、加深对时空系统及其之间转换关系的理解2、掌握常用时空基准之间的转换模型与软件实现3、每人独立完成

2、实验规定的内容二、实验内容本实验内容包括:内容一:编程实现GPS起点1980年1月6日0时对应的儒略日内容二:编程实现2011年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的点对应

3、的直角坐标为多少?三、实验过程1.针对第一、二部分内容:1.1解决思路:先建立”TimeStruct.h”的头文件,将格里高利历、GPS时间结构、儒略日时间结构共结构体的方式放在里面;在建立“TimeTr”的头文件,建立类“CTimeTr”,创建变量“GPSTime”、“Time”、”JulDay”,并且申明函数“TIME2JUL”、“TIME2GTIME”等,用这些函数分别实现所需要的转换。1.2具体的实现函数:“TIME2JUL”函数:doubleCTimeTr::TIME2JUL()//TIMETime,JULIAND

4、AY&JulDay{doublem,y;doubleD;//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.lDa

5、y=int(D);JulDay.lSecond=D-int(JulDay.lDay);return0;}“TIME2GTIME”:voidCTimeTr::TIME2GTIME(){doubleJD;longm,y;intWN;doubleWsecend;//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.b

6、yMonth;}JD=int(365.25*y)+int(30.6001*(m+1))+Time.byDay+Time.byHour/24.0+1720981.5;WN=floor((JD-2444244.5)/7.0);GpsTime.lWeek=WN;Wsecend=(JD-2444244.5-7*WN)*604800;GpsTime.lSecond=Wsecend;}1.3实验结果:..2针对第三部分内容:2.1解决思路:运用实验指导书中提供的matlab高斯反算的代码,进行解算;将高斯反算的公式直接输成matlab代

7、码,绕后在函数“function[B,L]=gauss_fansuan(x,y,L0)”中,将坐标x=3548910.811290287,y=179854.6172135982,L0=117,带入函数的坐边,即可得到所需要的经纬度。2.2主要函数的代码:function[B,L]=gauss_fansuan(x,y,L0)a=6378137;f=1/298.257223563;b=a-a*f;c=a^2/b;e=sqrt(a^2-b^2)/a;e1=sqrt(a^2-b^2)/b;Beta0=1-(3/4)*e1^2+(45

8、/64)*e1^4-(175/256)*e1^6+(11025/16384)*e1^8;Beta2=Beta0-1;Beta4=(15/32)*e1^4-(175/384)*e1^6+(3675/8192)*e1^8;Beta6=-(35/96)*e1^6+(735/2048)*e1^8;

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。