GPS测量与数据处理实验报告

《GPS测量与数据处理实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《GPS测量与数据处理（N文件）》实验（上机）报告班级·学号姓名实验日期6月10号任课教师实验名称实践型一、实验目的及要求1.通过编程N文件，深入理解matlab编程方式与运行原理。2..通过老师给定的N文件，看懂并在理解的基础上，加入自己的内容。3.学习高级语言软件MATLAB.二、实验平台1.微机一台2.matlab软件一套三、实验内容请写出实验内容的操作步骤（本表不够填写可加页）：（请写出实验内容各环节的详细操作步骤（可附屏幕截图）和常用功能的实现代码清单）三、实验内容1.打开MATLAB程序及面如下图所示2.打开已有的RENDNRINEX文件，并点击运行，运行成果如下（具体程序

2、段见附录）。3.通过本组编辑后的程序及运行结果如下（具体程序段见附录）。◆已知程序functionEphDat=Read_N_Renix(filename)%ReadtheN-filewiththerenixformat%Input:%filename:thenameoffileincludingthefilepath%Output:%EphDat:astructure,storingtheephemeris%.SatPRN;//theSatcode%.toc;//gpsweeksec%.a0;//differentofthesatclock(second)%.a1;//satcloc

3、kdrift(s/s)%.a2;//therateoftheSatclockdrift(s/s^2)%.IODE;//星历数据有效期IODE=Toe-Tl;%.Crs;//轨道半径正弦调和改正项振幅(m)%.DetaN;//平均运动修正量%.M0;//Toe时的平近点角%.Cuc;//纬度幅角余弦调和改正项振幅%.e;//卫星轨道扁心率%.Cus;//纬度幅角正弦调和改正项振幅%.sqrtA;//卫星轨道长半径方根%.Toe;//星历参考时间%.Cic;//轨道倾角余弦调和改正项振幅%.OMG0;//升交点赤经%.Cis;//轨道倾角正弦调和改正项振幅%.I0;//轨道倾角%.Crc;

4、//轨道半径余弦调和改正项振幅(m)%.omg;//近地点角距%.OMG0dot;//升交点赤经变化率%.I0dot;//轨道倾角变化率%.ISL2;//L2数据标志%.GpsWn;//GPSweeknumber%.ISL2P;//L2P数据标志%.SatAccu;//Theaccuracyofthesatellite(m)%.SatHth;//ThehealthoftheSatellite(MSB)%.Tgd;//单频接收机延迟改正数%.IODC;//时钟数据有效期%Beginprogram[filename,pathname]=uigetfile('*.**N','读取GPS广播星

5、历文件');fp=fopen(strcat(pathname,filename),'rt');if(fp==-1)error('errortoopenthe'+filename);end%readthefileheaderwhile(1)strTemp=fgets(fp);if(strTemp==-1)error('errorintheheaderfileof'+filename);endif(length(strTemp)<73)strTemp(length(strTemp)+1:73)='X';endif(strTemp(61:73)=='ENDOFHEADER')break;en

6、dend%readtheNfiledata%readthefirstepochstrTemp=fgets(fp);EphDat(1).SatPRN=str2num(strTemp(1:2));Year=str2num(strTemp(4:5));Mon=str2num(strTemp(7:8));Day=str2num(strTemp(10:11));THour=str2num(strTemp(13:14));TMin=str2num(strTemp(16:17));TSec=str2num(strTemp(19:22));if(Year>80)%GPS从20C80年代上商用Year=

7、Year+1900;elseYear=Year+2000;endEphDat(1).toc=ConvertGpsTime(Year,Mon,Day,THour,TMin,TSec);%转换成GPS周秒EphDat(1).a0=str2num(strTemp(23:41));%D,E需要转化。。。。。。。。。。EphDat(1).a1=str2num(strTemp(42:60));EphDat(1).a2=str2num(strTemp(61: