卫星导航定位实验报告
gps实验报告
gps实验报告GPS实验报告。
一、实验目的。
本实验旨在通过对GPS（全球定位系统）的原理和使用进行深入研究，掌握GPS的工作原理、定位原理和精度控制方法，以及GPS在实际应用中的一些特点和限制。
二、实验原理。
GPS是由24颗卫星组成的卫星导航系统，其中包括21颗工作卫星和3颗备用卫星。
这些卫星以6个轨道面，每个面上有4颗卫星的方式分布在大气层之外的轨道上，以提供全球范围的导航服务。
GPS接收机接收来自卫星的信号，并计算信号传播时间来确定自身的位置。
通过同时接收多颗卫星的信号，可以实现三维定位和速度测量。
三、实验内容。
1. GPS接收机的基本使用，打开GPS接收机，等待接收卫星信号并进行定位，观察定位结果的精度和稳定性。
2. GPS定位精度的影响因素，在不同环境条件下进行GPS定位实验，观察信号强度、遮挡物、大气层等因素对定位精度的影响。
3. GPS定位的实际应用，通过实际场景模拟，测试GPS在城市、山区、森林等不同环境下的定位效果，并对比不同场景下的定位精度和稳定性。
四、实验结果与分析。
经过一系列实验，我们得出以下结论：1. GPS定位精度受到环境因素的影响较大，如建筑物、树木等遮挡物会导致信号弱或者反射，从而影响定位精度。
2. 在城市环境中，由于高楼大厦的遮挡和信号反射，GPS定位精度可能会受到较大影响，定位结果可能出现偏移。
3. 在山区和森林等复杂环境中，GPS定位精度也会受到影响，但相对于城市环境，精度可能会更高一些。
五、实验总结。
通过本次实验，我们对GPS的工作原理和定位精度有了更深入的了解。
在实际应用中，我们需要注意环境因素对定位精度的影响，合理选择使用场景，以获得更准确的定位结果。
同时，GPS在城市环境下的定位精度仍然存在一定的局限性，需要结合其他定位技术进行辅助，以提高定位精度和稳定性。
六、参考文献。
[1] 赵云. GPS定位精度分析及影响因素研究[J]. 测绘工程, 2015(2): 15-21.[2] 李明. GPS技术在城市环境下的应用研究[J]. 地理信息科学, 2016, 18(3): 45-52.[3] 王强. GPS定位技术及其在森林环境中的应用[J]. 林业科学, 2017, 29(5): 78-84。
gps实验报告
gps实验报告GPS实验报告。
一、实验目的。
本实验旨在通过GPS定位技术，对GPS接收机进行测试，验证其定位精度和稳定性，以及对GPS信号的接收情况进行分析。
二、实验原理。
GPS（Global Positioning System，全球定位系统）是由美国国防部研制的一种卫星导航系统，它能够提供全球范围内的三维定位、速度和时间信息。
GPS系统由24颗运行于地球轨道上的卫星组成，这些卫星以特定的轨道和时间间隔发送信号，接收机通过计算这些信号的传播时间来确定自身的位置。
三、实验步骤。
1. 准备工作，将GPS接收机放置在开阔的空地上，确保周围没有高楼或其他遮挡物，以保证接收到的信号来自卫星而非其他干扰源。
2. 启动GPS接收机，打开GPS接收机的电源，等待其自动搜索卫星信号并进行定位。
3. 数据采集，在接收机显示屏上记录下当前的经度、纬度、海拔高度等信息，并记录下时间。
4. 移动测试，在不同的位置重复步骤2和3，以验证GPS定位的准确性和稳定性。
四、实验结果分析。
通过对实验数据的分析，我们发现在开阔的空地上，GPS接收机能够快速、准确地定位到当前位置的经纬度和海拔高度信息。
在移动测试中，随着移动位置的变化，GPS接收机能够实时更新定位信息，且定位精度较高。
五、实验总结。
本次实验验证了GPS接收机的定位精度和稳定性，证明了GPS定位技术在开阔空地上的可靠性。
然而，在城市高楼林立或密林深处等遮挡物较多的地方，GPS 信号的接收可能会受到影响，导致定位精度下降。
因此，在实际应用中，需要根据具体情况选择合适的定位方式。
六、参考文献。
1. 徐明，李华. GPS原理与应用[M]. 北京，科学出版社，2009.2. 王强，刘明. GPS技术应用实例分析[J]. 测绘通报，2015，(6)，78-82.七、致谢。
感谢实验室的各位老师和同学们在本次实验中的支持与帮助，让我们能够顺利完成实验并取得了丰富的实践经验。
以上即为本次GPS实验的报告内容，希望能对相关领域的学习和研究有所帮助。
gps定位实验报告
gps定位实验报告GPS定位实验报告引言：GPS（全球定位系统）是一种基于卫星导航的定位技术，它利用地球上的卫星系统来确定特定位置的方法。
本实验旨在探究GPS定位的原理和精度，并通过实际操作来验证其可靠性和准确性。
一、GPS定位原理GPS定位原理是基于三角测量的原理。
GPS接收机接收到来自卫星的信号后，通过测量信号的传播时间来计算出距离。
通过同时接收多颗卫星的信号，GPS接收机可以计算出自身与卫星之间的距离差，并根据这些距离差进行三角测量，从而确定自身的位置。
二、实验设备与方法本实验使用了一台GPS接收机和一台笔记本电脑。
首先，将GPS接收机与笔记本电脑通过USB线连接，确保接收机与电脑之间的通信畅通。
然后，打开接收机的电源，并在电脑上打开相应的GPS定位软件。
接下来，等待接收机与卫星建立连接，并获取到足够的卫星信号。
最后，记录下接收机显示的经纬度信息，并与实际位置进行对比。
三、实验结果与分析在进行实验过程中，我们发现GPS接收机的定位速度相对较快，一般在几秒钟内就能够获取到足够的卫星信号进行定位。
通过与实际位置进行对比，我们发现GPS定位的准确性非常高，误差一般在几米以内。
这证明了GPS定位技术的可靠性和精度。
然而，我们也注意到GPS定位的准确性可能会受到一些因素的影响。
例如，高楼大厦、山脉和树木等物体可能会阻碍卫星信号的传播，从而导致定位的不准确。
此外，天气条件也可能对GPS定位的精度产生影响。
在恶劣的天气条件下，如大雨或大雪，卫星信号的传播可能会受到干扰，从而影响定位的准确性。
四、GPS定位的应用GPS定位技术在现代社会中有着广泛的应用。
首先，GPS定位技术在导航领域被广泛使用。
无论是在汽车导航系统中还是在手机导航应用中，GPS定位都能够帮助人们准确地找到目的地。
其次，GPS定位技术在物流和运输领域也发挥着重要作用。
通过实时监控车辆的位置，物流公司可以更好地管理和调度运输车辆，提高物流效率。
卫星导航定位技术实习报告
标题：卫星导航定位技术实习报告一、实习背景与目的随着全球经济一体化和科技发展的日新月异，卫星导航定位技术在各个领域的应用越来越广泛。
为了更好地了解和学习卫星导航定位技术，提高自己在相关领域的实际操作能力，我参加了为期一个月的卫星导航定位技术实习。
实习期间，我参与了导师的科研项目，学习了卫星导航定位技术的原理、应用以及数据处理方法，并对该技术在我国精准农业领域的应用进行了深入研究。
二、实习内容与过程1. 理论学习在实习的开始阶段，导师为我讲解了卫星导航定位技术的基本原理、发展历程以及各类卫星导航系统。
我了解到，卫星导航定位技术是利用导航卫星发射的信号，通过接收器接收并处理这些信号，从而确定用户位置的一种技术。
目前全球主要的卫星导航系统有美国的GPS、中国的北斗、欧盟的伽利略、俄罗斯的格洛纳斯和日本的准天顶。
此外，我还学习了卫星导航定位技术的应用领域，包括交通运输、精确农业、地形测绘、地质勘探等。
2. 实际操作在理论学习的基础上，我开始参与导师的科研项目。
实习期间，我主要负责使用卫星导航定位设备进行实地测量，收集数据，并利用相关软件进行数据处理和分析。
我学会了如何操作卫星导航定位设备，包括发射器的安装、接收器的设置以及数据的采集。
同时，我还掌握了运用专业软件对采集到的数据进行处理和分析的方法，如GPS数据处理软件、地理信息系统（GIS）等。
3. 项目研究在实际操作的过程中，我深入了解了卫星导航定位技术在精准农业领域的应用。
我参与了导师团队的研究项目，针对农业生产中的实际问题，如农田土壤养分监测、作物病虫害防治、灌溉管理等，研究利用卫星导航定位技术进行解决方案的设计。
通过项目研究，我学会了如何将卫星导航定位技术应用于实际生产，提高农业生产效率。
三、实习收获与体会通过这次实习，我对卫星导航定位技术有了更加深入的了解，从理论到实践都有了很大的提高。
我认识到，卫星导航定位技术不仅具有很高的精确度和可靠性，而且在各个领域的应用潜力巨大。
GPS标准实习报告3篇
GPS标准实习报告GPS标准实习报告精选3篇（一）实习报告实习单位：XYZ科技有限公司实习日期：2022年1月1日至2022年1月31日一、实习背景和目的：我作为一名电子信息工程专业的学生，在大学期间学习了很多关于GPS（全球定位系统）的知识。
我选择在XYZ科技有限公司进行实习，是希望能够将课堂上学到的理论知识与实际工作相结合，提升自己在GPS领域的实践能力。
二、实习任务和工作内容：在实习期间，我主要参与了XYZ科技有限公司的GPS标准化工作。
具体工作内容如下：1.学习和熟悉GPS的基本原理和工作原理。
2.调研和分析国内外GPS标准和规范。
3.参与制定GPS标准的技术方案和规范。
4.与相关部门进行配合，收集并整理GPS标准的资料和文档。
5.参与GPS标准化工作的会议和讨论，提出自己的建议和意见。
三、实习收获和体会：通过这次实习，我对GPS标准化工作有了更深入的了解，并且提升了自己的实践能力。
在参与制定GPS标准的过程中，我发现了自己在理论知识上的不足，并且学会了如何通过调研和分析来解决问题。
同时，在与团队成员的合作中我也学到了团队合作的重要性和沟通的技巧。
四、实习总结和建议：通过这次实习，我收获了很多，并且对自己的职业规划也有了更明确的目标。
我觉得在未来的工作中，我应该继续加强自己的理论知识，并且注重实践能力的培养。
同时，我也建议公司在GPS标准化工作中加强与国内外相关部门的合作，提高GPS标准的制定水平。
以上就是我的实习报告，感谢公司给予我这次宝贵的实习机会，希望我能够在未来的工作中发挥自己的专业知识和能力。
GPS标准实习报告精选3篇（二）根据我在GPS测量实习中的经验，我总结出以下几个关键点：1. 准备工作：在实习开始前，我需要熟悉GPS测量的基本原理和使用方法。
这包括了对GPS仪器的了解，熟悉GPS测量的误差来源和影响因素，以及对测区的了解和准备。
2. 测量准确性：在实习过程中，我学会了如何提高测量的准确性。
卫星导航定位技术实习报告
一、实习背景随着科技的不断发展，卫星导航定位技术在各个领域得到了广泛应用。
为了深入了解这一技术，提高自身实践能力，我参加了为期两周的卫星导航定位技术实习。
本次实习旨在通过理论学习和实践操作，掌握卫星导航定位的基本原理、应用领域及操作方法。
二、实习内容1. 理论学习（1）卫星导航定位技术概述实习期间，我首先学习了卫星导航定位技术的基本概念、发展历程、系统组成及工作原理。
通过学习，我了解到全球定位系统（GPS）、全球导航卫星系统（GNSS）等是我国在卫星导航领域的重要成果。
（2）卫星导航定位技术原理我深入学习了卫星导航定位技术的基本原理，包括伪距测量、多普勒测速、双曲面交会定位等。
这些原理是卫星导航定位技术实现的基础。
（3）卫星导航定位技术应用我了解了卫星导航定位技术在农业、交通、测绘、军事等领域的应用。
这些应用使卫星导航定位技术成为现代社会不可或缺的一部分。
2. 实践操作（1）GPS接收机操作在实习过程中，我学习了GPS接收机的操作方法。
通过实际操作，我掌握了GPS接收机的开机、关机、数据采集、数据传输等功能。
（2）卫星导航定位数据处理我学习了卫星导航定位数据处理的流程，包括数据预处理、坐标转换、误差分析等。
通过实际操作，我掌握了数据处理软件的使用方法。
（3）卫星导航定位技术应用实践在实习期间，我参与了以下项目实践：1. 农业机械控制：利用卫星导航定位技术实现变量施肥播种机、联合收割机、无人驾驶拖拉机等农业机械的精准作业。
2. 精准农业：通过卫星导航定位技术，对农田进行精细化管理，提高农作物产量。
3. 测绘：利用卫星导航定位技术进行大地测量、地形测绘等。
三、实习收获1. 理论知识方面通过本次实习，我对卫星导航定位技术的基本原理、应用领域及操作方法有了更深入的了解。
这为我今后从事相关工作奠定了坚实的理论基础。
2. 实践能力方面在实习过程中，我学会了GPS接收机的操作、数据处理软件的使用，以及卫星导航定位技术在实际项目中的应用。
卫星导航原理实验报告
卫星导航原理实验报告实验目的本实验旨在通过实际操作，加深对卫星导航原理的理解，掌握卫星导航的基本工作原理、信号接收与处理方法。
实验原理卫星导航是利用人造卫星在太空中运行，通过卫星定位系统向用户提供空间位置、速度和时间等信息的导航方式。
其原理是通过接收多颗人造卫星发射的信号，利用信号的时间差异和测量误差，计算出用户的三维空间位置。
卫星导航系统由地面控制站、卫星和用户终端组成。
地面控制站负责发送导航信号和控制卫星运行，卫星接收地面控制信号并通过天线以无线电信号形式发送到用户终端，用户终端接收并解码信号，计算用户位置。
实验步骤1. 连接设备：将接收天线连接到接收设备上，确保连接正常；2. 打开接收设备：根据具体型号，按下相应按钮或转动开关打开接收设备；3. 接收卫星信号：对设备进行信号搜索，确保接收到卫星信号；4. 信号处理：接收设备将信号传输到计算机或显示屏上，进行信号处理；5. 计算用户位置：根据接收到的信号，使用相应的算法计算用户的三维空间位置。
实验结果经过一系列操作，最终成功接收到卫星信号，并通过计算机显示用户位置。
实验结果表明，卫星导航系统具备高精度和广域覆盖的能力。
实验总结本实验通过操作接收设备，将卫星信号传输到计算机上进行处理，实现了卫星导航的基本功能。
在实验过程中，我们对卫星导航原理有了更加深入的了解，掌握了信号搜索和处理的方法。
卫星导航在交通、军事和民用领域具有广泛应用前景。
它可以为车辆导航、航空航天、灾害救援等提供准确的定位和导航服务。
此外，随着技术的不断发展，卫星导航系统的精度和覆盖范围将会进一步提高，为人们的生活带来更多的便利。
通过本次实验，我们不仅学习了卫星导航的原理和操作方法，还了解了其应用领域和发展前景。
相信在今后的学习和工作中，我们将会更好地运用卫星导航技术，为社会发展做出贡献。
卫星导航定位算法与程序的设计——实验报告
2013 级测绘工程专业卫星导航定位算法与程序设计实验报告实验名称：卫星导航基本程序设计班级：学号：姓名：实验时间： 2016年6月28日~2016年6月30中国矿业大学目录实验一时空基准转换 (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 = 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”的头文件，建立类“CTimeT r”，创建变量“GPS Time”、“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;}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代码，绕后在函数“function [B,L] = gauss_fansuan （x,y,L0）”中，将坐标x = 3548910.811290287，y = 179854.6172135982，L0 = 117，带入函数的坐边，即可得到所需要的经纬度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
China University of Mining and Technology 《卫星导航定位算法与程序设计》实验报告学号： 07122825姓名：王亚亚班级：测绘12—1指导老师：王潜心/张秋昭/刘志平中国矿业大学环境与测绘学院2015-07-01实验一编程实现读取下载的星历一、实验要求：读取RINEX N 文件，将所有星历放到一个列表（数组）中。
并输出和自己学号相关的卫星编号的星历文件信息。
读取RINEX O文件，并输出指定时刻的观测信息。
二、实验步骤：1、下载2014年的广播星历文件和观测值文件，下载地址如下：ftp:///gps/data/daily/2014/2、要求每一位同学按照与自己学号后三位一致的年积日的数据文件和星历文件，站点的选择必须选择与姓氏首字母相同的站点的数据，以王小康同学为例，学号：07123077，需下载077那天的数据。
有些同学的学号365<后三位<730，则取学号后三位-365，以姜平同学为例：学号10124455，下载455-365=90 天的数据，有些同学的学号730<后三位<=999，则取学号后三位-730，以万伟同学为例：学号：07122854，则下载854-730 = 124天的数据。
可以选择wnhu0124.14n wnhu0124.14o 根据上述要求我下载了2014年第95天的数据，选择其中的wsrt0950.14n和wsrt0950.14o星历文件。
指定时刻（学号后五位对应在年积日对应的秒最相近时刻）的观测值信息如张良09123881，后五位23881，取23881-3600*6= 2281秒，6点38分01秒，最近的历元应该是6点38分00秒的数据。
根据计算与我最接近的观测时刻为2014年4月5日6点20分30.00秒。
3、编程思路：利用rinex函数读取星历文件中第14颗卫星的星历数据并输出显示。
对数据执行762次循环找到对应的2014年4月5日6点20分30.00秒，并输出观测值。
4、程序运行结果：三、编程主要代码：1、读写N 文件主函数代码： clc;%清除屏幕之前打过的代码clear all ;%从内存中释放所有的内存变量和数组以及所有用户自定义的菜单栏, 菜单, 和窗口的定义。
close all ;%关闭当前及所有工作区中所有打开的数据库, 表, 和索引，并选择 1 号工作区 [filename,filepath]=uigetfile('*.14n','选择计算的输入文件');%打开窗口选取计算的输入文件 file=[filepath filename];%获得文件名 rinexe(file,'eph.dat');%函数rinexe 读取一个标准格式化导航信息文件并重新格式化数据为21行列数等于卫星数的矩阵。
Eph = get_eph('eph.dat');%将存放在‘eph.dat ’中的导航电文数据提取出来，为计算卫星位置做准备，返回值-eph 星历矩阵，将读到的数据存入内存（注意：eph 星历矩阵，每颗卫星存放21个参数） Eph[filename,filepath]=uiputfile('王亚亚-作业1读写N 导航文件.txt','选择一个路径对结果文件保存');%选择路径保存输出的结果文件 file=[filepath filename];%获取文件名 fid = fopen(file,'wt+');%文件有两种格式：二进制文件(b)和文本文件(t)。
在Windows 下，打开文件的默认是二进制格式，如果要以文本方式打开，则必须在打方式中加上字符 't'。
fprintf(fid,'%c',' 班级 姓名 学号 作业序号 主要内容 ');fprintf(fid,'%c\n',' ');fprintf(fid,'%c',' 测绘12-1 王亚亚07122825 上机实验作业一 读写N 导航文件 '); fprintf(fid,'%c\n',' '); for k = 1: size(Eph,1)p=num2str(Eph(k,44));%输出第14颗卫星的星历数据fprintf(fid,'%c',p);%以一个字符输出 if k == 1;fprintf(fid,'%c',' %GPS 卫星PRN 编号'); end if k == 2;fprintf(fid,'%c',' %GPS 卫星钟的漂移速度'); endif k == 3; fprintf(fid,'%c',' %参考时刻的升交点赤经M0'); endif k == 4; fprintf(fid,'%c',' %roota sqrt(A)( m1 / 2(根号下) ) 轨道长半轴的平方根'); endif k == 5; fprintf(fid,'%c',' %deltan Δn(rad /s) 平均运动修正量'); endif k == 6; fprintf(fid,'%c',' %ecc e 轨道偏心率'); end if k == 7;fprintf(fid,'%c',' %omega ω(rad) 近地点的角距'); endif k == 8; fprintf(fid,'%c',' %cuc Cuc(rad) 纬度幅角的余弦调和项改正的振幅'); endif k == 9; fprintf(fid,'%c',' %cus Cus(radians) 纬度幅角的正弦调和项改正的振幅'); endif k == 10; fprintf(fid,'%c',' %crc Crc(m) 轨道半径的余弦调和项改正的振幅');endif k == 11; fprintf(fid,'%c',' %crs Crs（m）轨道半径的正弦调和项改正的振幅'); endif k == 12; fprintf(fid,'%c',' %i0 i0(rad) 轨道倾角'); endif k == 13; fprintf(fid,'%c',' %idot i(rad/s)(IDOT) 轨道倾角变化率'); endif k == 14; fprintf(fid,'%c',' %cic Cic(rad) 轨道倾角的余弦调和项改正的振幅'); endif k == 15; fprintf(fid,'%c',' %cisCis(rad) 轨道倾角的正弦调和项改正的振幅'); endif k == 16;fprintf(fid,'%c',' %Omega0Ω(rad)(OMEGA) 参考时刻的升交点的赤经'); endif k == 17;fprintf(fid,'%c',' %Omegadot Ω(rad/s)(OMEGA DOT) 升交点赤经的变化率'); endif k == 18;fprintf(fid,'%c',' %toe TOE星历的参考时刻(GPS周内的秒数) 星历表参考历元'); endif k == 19; fprintf(fid,'%c',' %af0 卫星钟的偏差（s）卫星钟偏差'); endif k == 20; fprintf(fid,'%c',' %af1 卫星钟的漂移（s/s）'); endif k == 21;fprintf(fid,'%c',' %toe 存储TOE星历的参考时刻(GPS周内的秒数) 电文发送时刻'); endfprintf(fid,'%c\n',' ');endfclose(fid);2、读写O文件主函数代码：clc;%清除屏幕之前打过的代码clear all;%从内存中释放所有的内存变量和数组以及所有用户自定义的菜单栏, 菜单, 和窗口的定义。
close all;%关闭当前及所有工作区中所有打开的数据库, 表, 和索引，并选择1 号工作区% We identify the observation file and open it我们确定观测文件并打开它[filename,filepath]=uigetfile('*.14O','选择计算的输入的O文件');file=[filepath filename];%文件名fid1 = fopen(file,'rt');%r：读出观测文件ofile1，t：以文本的格式（因为fopen函数默认打开方式是二进制）[Obs_types1, ant_delta1, ifound_types1, eof11] = anheader(file);%打开观测文件ofile1，从头文件中截取观测类型Obs_types1、天线参数的信息ant_delta1，以及标明是否寻找到观测类型行ifound_types1、是否到文件尾的判定eof11 NoObs_types1 = size(Obs_types1,2)/2;%返回观测类型Obs_types1列数(1:行2：列)的Obs_types1矩阵是1行14列：[C1P1P2L1L2D1D2] size（Obs_types1,2）值是14除以2得观测值类型个数是7种Pos = [];% There are 20 epochs of data in ofile1SITE247J.01N文件中有20个历元的数据for q = 1:761 %循环762次找到对应的2014年4月5日6点20分30.0000000秒[time1, dt1, sats1,eof1,datee] =fepoch_0(fid1);% [time1, dt1, sats1, eof1] =fepoch_0(fid1);%通过函数fepoch_0调用观测文件SITE247J.01O，返回四个参数：time1：周内秒、dt1：接收机时钟的偏差( 单位为s, 为可选项)、sats1：卫星的PRN号矩阵、eof1：返回是否到文件尾的标示符（1：结束0：未结束）NoSv1 = size(sats1,1);%通过查看行向量sats1的长度，获取卫星个数% We pick the observed P2 pseudoranges 我们获取观测的P2码伪距其中伪随机噪声码（PN）即测距码主要有精测距码（P码）和粗测距码（C/A码）两种。