GPS绝对定位的程序设计修改02要点

1、目录摘 要IABSTRACTII1 绪论11.1 引言11.2 国内外发展动态11.3 本文主要研究工作和章节安排22 GPS卫星轨道坐标计算数学模型23 GPS单点定位数学模型64 程序设计84.1 VB介绍84.2 绝对定位程序的总体设计94.2.1 系统描述与需求分析94.2.2 软件总体设计94.2.3 设计软件系统结构94.2.4 数据结构与数据库设计104.2.5 程序中的模块104.2.6 程序流程105 计算实例及误差分析126 结论14致 谢15参考文献16GPS绝对定位的程序设计摘 要GPS单点定位是根据计算出的卫星坐标，再以伪距或者载波相位为观测量，通过空间后方交会原理计

2、算出地面测量站点的位置。本文主要介绍了GPS卫星轨道坐标计算数学模型，单点定位数学模型，并根据最小二乘原理，利用VB语言对以上过程进行了实现。本文对此程序的开发过程进行了阐述，包括时间转换程序、利用广播星历计算卫星坐标程序和地面点近似坐标计算程序，并总结了程序开发过程中对误差进行处理的方法和一些关键技术问题的解决思路。最后，选取实例进行计算并进行精度分析。关键词：GPS 单点定位 坐标计算 精度分析PROGRAM DESIGN OF GPS ABSOLUTE POSITIONINGABSTRACTSingle-point positioning of GPS calculates the lo

3、cation of the site ground-based measurements by the principle of space resection according to the calculated satellite coordinates, which measure the pseudorange or carrier phase for the concept. My article describes a mathematical model of the GPS satellite orbit coordinates calculating, a mathemat

4、ical model of single point of positioning, and achieve the above process in accordance with the principle of least squares, using VB language . my article describes the process of the program development, and summarizes the ways to deal with the error and the idea of some key technical issues in the

5、 program development process. Finally, selects the instance to calculate and analyse the accuracy .KEY WORDS GPS Single-point positioning Coordinate calculation Accuracy Analysis1 绪论1.1 引言GPS是美国从20 世纪70 年代开始研制的, 于1994 年全面建成, 具有在海、陆、空进行全方位实时三维导航与定位能力的新一代卫星导航与定位系统。尤其是经过近几年的研究,GPS 更在测绘、航空遥感和气象等方面有了新的应用

6、, 并以全天候、高精度、自动化、高效益等显著特点, 赢得广大用户的信赖。随着对定位精度要求的不断提高, 人们对GPS卫星星历的精度和实时性提出了越来越高的要求。卫星的星历, 是描述有关卫星运动轨道的信息。利用GPS进行定位, 就是根据已知的卫星轨道信息和用户的观测资料, 通过数据处理来确定接收机的位置及其载体的航行速度。所以, 精确的轨道信息是精密定位的基础。GPS 的卫星星历按照精度可分为精密星历和广播星历。精密星历是由国际GPS服务中心( IGS) 通过Internet 发布,它的轨道精度可达到10cm 左右, 足以满足精密定位的需要。但是精密星历只能在卫星观测的11d后获得, 无法为实时

7、定位、导航、气象等实时性要求很强的应用提供有效的服务。广播星历是通过接收机接收卫星发射的含有轨道信息的导航电文, 经过解码获得的卫星星历推算得到卫星位置, 可以实现实时的导航和定位。本程序以2009年11月21日上海跟踪站（SHAO）的RINEX格式广播星历shao3250.09n和观测数据shao3250.09o为例，取了200个连续观测历元，在不同历元求出坐标值，最后求出坐标平差值，对平差值的各分量作比较。1.2 国内外发展动态绝对定位按天线所处状态可以分为动态绝对定位和静态绝对定位；按观测量的性质，可分为测码伪距绝对定位和测码伪距相对定位。尤其是动态绝对定位的方法，由于它的实时、省力、简

8、洁等优点具有广阔的应用天地。近十年GPS理论以趋于完善，特别是GPS静态定位方面，已有许多可靠和成功的软件，但在绝对定位方面，尤其是实时绝对定位方面，由于GPS定位中包含许多误差，利用传统的方法很难将其影响除掉，使得GPS在某些场合的应用收到一定限制。GPS定位的主要误差源包括：1）卫星测量误差。可分为卫星时钟误差、星历误差、电离层的附加延迟误差、对流层的附加延迟误差、多径效应以及接收机本身的噪音。2）卫星几何位置造成的定位误差。GPS正展现出愈来愈巨大的应用潜力，应用不仅从定位范围和定位方式上提出了各种要求，定位精度的要求也有大幅提高，从米级到亚毫米级，而且提出的较高要求指标还包括可靠性和可

9、用性。为了增强可靠性和可用性，在实际应用中还发展了与GLONASS、无线信标、WARS、INS等系统进行联合集成，同时欧空局正在积极推进GALILEO系统，国际上也随之在酝酿GPS和GALILEO定位系统的集成。比较可靠地定位系统还必须有可靠的定位软件与计算机硬件的辅助，定位方面的软件开发也是各个研究机构比较注重的方面之一。有了好的软件，人们才能在生活和工作中更方便、快捷、准确、可靠和充分的利用定位系统给人类带来的方便之处。1.3 本文主要研究工作和章节安排本文主要介绍了GPS卫星轨道坐标计算数学模型，单点定位数学模型，并根据最小二乘原理，利用VB语言对以上过程进行了实现。整个论文内容安排如下

10、：第1章简单介绍了GPS发展，阐述了国内外GPS绝对定位的研究状况，以及本文所做的工作及论文结构第2章介绍了卫星轨道参数与广播星历参数以及GPS卫星轨道坐标计算数学模型。第3章介绍了GPS绝对定位的算法实现以及数学模型第4章介绍了编程语言VB的基本情况，GPS绝对定位的程序的设计目的、思路、流程以及程序中设计模块的应用第5章研究了软件的使用效能，也就是实例的计算与分析第6章为本论文的结论2 GPS卫星轨道坐标计算数学模型广播星历就是卫星GPS将含有轨道信息的导航电文发送给用户接收机，然后经过解码获得的卫星星历。GPS用户通过卫星广播星历，可以获得16个卫星星历参数，其中，1个参考时刻，6个相应

11、参考时刻的开普勒轨道参数和9个摄动力影响的参数。这些参数的定义如下表所示：表1导航电文中的参考参数参数名称参数名称星历参数的参考历元平均运行速度差轨道长半轴的方根升交点赤经变化率轨道偏心率轨道倾角变化率参考时刻的轨道倾角升交距角的调和改正项振幅参考时刻的升交点赤经轨道倾角的调和改正项振幅近地点角距卫星地心距的调和改正项振幅参考时刻的平近点角AODE星历数据的龄期其中，AODE表示从最后一次注入电文起外推星历时 0的外推时间间隔，它反映了外推星历的可靠程度。根据上述数据，便可外推出观测时刻t的轨道参数，从而计算卫星在不同参考系中的相应坐标。用广播星历参数计算卫星位置：在利用GPS信号进行导航定位

12、时，为了解算用户在地心坐标系中的位置，GPS接收机需要测定测站到卫星的距离并且要知道同一卫星在同一时刻的地心坐标2。卫星的地心坐标是从卫星的导航电文中提供的开普勒轨道参数和轨道摄动修正量按一定公式计算的。1)计算卫星运行的平均角速度：卫星的平均角速度用下式计算： (1)式中，GM=398600.5是WGS-84坐标系中地球引力常数。利用导航电文中给出的摄动改正数，用下式求卫星运行的平均角速度： (2)2)对观测时刻做卫星钟差改正： (3)在计算卫星钟差改正时，可近似取。3)观测时刻的平近点角的计算： （4）4)计算偏近点角： （5）（5）式可用迭代法进行计算，即先令代入上式，求出再代入上式计算

13、，由于偏心率e很小（只有0.01），因此收敛很快，只需迭代两次便可求出偏近点角。5)真近点角的计算： （6）6)计算升交角距及轨道摄动改正项：升交角距： 摄动改正项： （7）7)计算经过摄动改正的升交角距、卫星到地心距离、轨道倾角 （8）8)计算卫星轨道平面坐标系中的坐标：卫星在轨道平面坐标系中的坐标为 （9）9)计算观测时刻升交点经度：升交点经度为该时刻升交点赤经与格林尼治恒星时GAST之差，即 （10）观测时刻的升交点赤经为参考历元的升交点赤经加上观测时刻与参考历元之间的升交点的赤经变化，即 （11）另外，卫星电文中提供了一周开始时刻（星期六子夜）以秒计算的格林尼治恒星时GATS。由于地球

14、的自转作用，GAST也不断增加。增加量与地球自转速率有关=7.29211567。所以，观测时刻GAST用下式计算： GAST=GAST+ （12）考虑到（11）式和(12)式，则 (13)因为 考虑到和都是从开始起算，即0，则(13)式为 （14）10)计算卫星在地心坐标系中空间直角坐标： （15）11)如果考虑极移影响，可求在协议地球坐标系中的空间直角坐标： （16）3 GPS单点定位数学模型由于接收机测量的是伪距，在观测值中存在着接收机钟差，加之测量点的三维坐标为待求值，一共有4个未知数。要求解出这4个未知数，必须有4个方程式。为此，要实现单点绝对定位必须同时观测4颗卫星，才能组成定位的基

15、本方程4。设为伪距观测量，为接收机到卫星的真距离，为接收机钟差，则观测方程为 （17）式中，假定伪距观测量已经过星历中的对流层和电离层改正；为卫星的瞬时地心坐标，可由卫星星历电文中求出；为接收机的地心坐标，是待求量。为了求解方便和数据处理的需要，将式（17）进行微分，作线性化处理，并将接收机的概略坐标作为初始值代入，得到 （18）式中，为接收机钟差对应的空间距离，从式（18）中看出，三个坐标分量的系数是接收机到卫星的单位矢径分别向三个坐标轴投影的方向余弦。采用符号 （19）规定上标为卫星号，下标为测站号，则组成伪距定位的基本方程 （20）采用矩阵表示 观测量 状态矩阵 未知数则式（20）变为

16、（21）对式（21）求解，便得到接收机地心坐标的唯一解 （22）在计算过程中，下列几个问题必须注意11：（1）卫星之间的钟差是利用导航电文中给出的钟差改正数统一到UTC时间。这里，考虑的钟差是指卫星与接收机之间的钟差。（2）在计算中采用了接收机的概略坐标，第一次计算出的结果是不精确的。因此，必须反复迭代计算，直到满足规定的限差为止。（3）在一般导航型接收机中，都是采用这一数学模型计算位置的。现有的接收机都能同时跟踪四个以上卫星，但在计算中仍然利用四个卫星，不过是结果挑选的四个卫星。为此，按卫星的星座分布分成若干组，计算其PDOP，最后选择和利用一组其PDOP为最小的卫星作为计算数据，以得到最高

17、的定位精度。在测地型接收机和高质量的导航接收机中，都具有8个以上的通道，能同时跟踪7颗以上的卫星。为了提高定位精度，在计算位置过程中，利用了所有的卫星观测值。在这样情况下，出现了多余观测，观测值的个数超过了未知数的个数，使得式（21）的右端不等于零 （23）式中，为残差向量。根据最小二乘法的原理，最后得到接收机的位置解为 （24）其精度为 （25）式中，为伪距测量中误差，为权系数阵。这种多余观测的优点在于消除了卫星定位的系统误差。过去，我们经常发现在仅用4颗卫星的差分定位中，当中间更换卫星时，位置会出现较大的偏移，等过了数秒后又逐渐回到原位。定位精度越高，这一现象越明显，当应用4颗以上的卫星定

18、位时，这一现象就不存在了。4 程序设计4.1 VB介绍 Visual Basic,简称VB，是Microsoft公司推出的一种Windows应用程序开发工具。是当今世界上使用最广泛的编程语言之一，它也被公认为是编程效率比较高的一种编程方法。无论是开发功能强大、性能可靠的商务软件，还是编写能处理实际问题的实用小程序，VB都是最快速、最简便的方法。4.2 绝对定位程序的总体设计 系统描述与需求分析基本功能要求：实现GPS L1 C/A码伪距单点定位功能，精度符合GPS标准定位服务SPS可选择输入/输出文件名，输入不正确时，可识别并提示（即星历和观测值文件不同步）可显示GPS时和北京时间的转换程序运

19、行稳定输入：原始数据为RINEX的O文件，星历为RINEX的N文件或SP3输出：每个历元下的单点定位值，方式为屏幕，自定义格式文件运行环境：Windows环境下的控制台程序或窗口程序且为单独可执行的文件 软件总体设计在需求分析的基础上通过抽象和分解将系统分解为模块，确定系统功能是否实现，即把软件需求转换为软件模块表示的过程。该过程是一个分解的过程输出计算输入图1软件的总体设计 设计软件系统结构采用某种方法，将系统按功能划分成模块。（划分）确定模块的功能。（功能）确定模块之间的调用关系。（调用）确定模块之间传递的信息。（接口）评价模块结构的质量。（质量） 数据结构与数据库设计1)主要的数据结构：

20、单历元观测值星历参数（广播星历、精密星历）坐标时间2)数据库设计：概念设计逻辑设计物理设计 程序中的模块1)完成各功能的模块：数据读取模块（RINEX读取） 开始读取RINEX卫星数据存于星历结构中读取一历元RINEX观测数据根据观测历元时刻和观测值，计算卫星发射时刻求解卫星发射时刻的卫星坐标修正卫星传播时间传播时间收敛否观测值进行各项改正组法方程系数，求解法方程计算坐标及精度信息与前一次计算值比较，收敛否输出该历元结果文件结束否结束图2程序流程图卫星位置计算模块 各项改正模块 数据显示输出模块2)基础模块：法方程解算模块（矩阵运算）时间模块坐标计算及转换模块 程序流程单点定位程序的流程是一个

21、反复迭代的过程（卫星位置、接收机位置的迭代），程序流程图如图2所示。5 计算实例及误差分析利用该软件打开星历文件和观测文件，计算出观测文件中连续观测到的卫星坐标、观测文件中GPS卫星的坐标，最终计算出接收机的位置坐标。步骤截图如下：图3软件主界面截图图4观测文件中GPS卫星的坐标图5观测文件中连续观测到的卫星坐标图6坐标及坐标改正数6 结论如何利用GPS的广播星历来获得比较准确的卫星位置, 进而为实时、精确的导航和定位服务是一件相当复杂的工作。如何在工作中考虑各种误差对定位精度的影响, 如何利用不同的算法来提高定位的精度, 一直是人们非常关心的课题。通过一个实例的计算, 详尽介绍了利用GPS广

22、播星历计算卫星位置和利用观测数据计算地面点坐标的方法和原理,简单的对计算结果的误差进行了分析, 证实了该算法的可行性。由于本作者水平有限，文中不足之处还有很多，如何进一步提高GPS定位的精度，还需要进一步的研究。致 谢本论文是在河南工程学院土木系刘小强老师的悉心指导下完成的。刘老师作为一名优秀的、经验丰富的教师，具有丰富的专业知识和教学经验，在整个论文实验和论文写作过程中，对我进行了耐心的指导和帮助，提出严格要求，引导我不断开阔思路，为我答疑解惑，鼓励我大胆创新，使我在这一段宝贵的时光中，既增长了知识、开阔了视野、锻炼了心态，又培养了良好的实验习惯和科研精神。在此，我向我的指导老师表示最诚挚的谢意！在论文即将完成之际，我的心情久久无法平静，从开始选题到顺利论文完成，有不知多少多少可敬的师长、同学、朋友给了我无数的帮助。感谢河南工程学院土木系提供的实验器材和测绘工具，感谢土木系全体老师给予我丰富的专业知识和各个方面的关心和帮助，感谢小组长的认真负责，感谢合作组员的热心协助。同时也要感谢土木系08级测绘专业的全体同学，正是由于