[优秀毕业论文]GPS导航解算仿真系统

GPS导航解算仿真系统作者姓名： 专业班级： 指导老师： 摘 要全球卫星定位系统是具有全球性、全能性(陆地、海洋、航空与航天)、全天候优势的导航定位、定时和测速系统。它汇集了当代最先进的空间技术、通讯技术及微电子技术，以其定位精度高，全天候获取信息，用户接收设备轻巧、价格比较低廉等诸多优点而被世人注目，在很多领域获得了应用。GPS用户部分的核心是GPS接收机。其主要由基带信号处理和导航解算两部分组成。导航解算部分主要包括根据导航数据中的星历参数实时进行各可视卫星位置计算。此次设计是针对导航解算仿真系统的Matlab仿真实现，因要求不高，所以对卫星运动做了理想化处理，摄动力对卫星的影响忽略不计（所以为无摄运动），采用开普勒定律及最小二乘法计算其轨道参数，对其运动规律进行简略分析，并使用Matlab编程仿真实现了卫星的运功轨道平面、运动动态、可见卫星的分布及利用可见卫星计算出用户位置。关键字：GPS导航解算 无摄运动 伪距 matlab仿真Simulation On GPS Navigation Solution System Using MATLABAbstract：The global positioning system is a global, versatility (land, sea and air and space),all-weather advantage of navigation and positioning, timing and speed system. It collects the contemporary most advanced technology, the communication technology and the space of microelectronics technology, with its positioning accuracy, all-weather information, users receive lighter, the price is cheap equipment advantages of men, and in many fields.GPS is the core of the user part of the GPS receiver. The main by the baseband signal processing and navigation solution of two parts. Navigation solution, including some of the major navigation data based on parameters in the ephemeris satellite positions in real-time computation of the visual. The navigation solution designed for matlab Simulation Simulation System, for less demanding, so do the idealized satellite motion processing, the impact of the satellite perturbation is negligible (so as not perturbed motion), the use of Kepler Law and the least squares method to calculate the orbital parameters, its law of motion for a brief analysis and simulation using matlab programming the pump power to achieve the satellite orbital plane, movement dynamics, the distribution of visible satellites and the use of visible satellites to calculate the location of the user.Key words：GPS navigation solution；Non-disturbed motion；Pseudorange；Matlab simulation目录第1章 前言11.1 课题背景11.2 本课题研究的意义和方法21.3 GPS的发展简史及系统概述3第2章 GPS的坐标、时间系统和卫星的运动52.1 地球坐标系52.2 时间系统72.3 卫星的运动82.3.1开普勒定律82.3.2 无摄卫星运动的轨道参数102.3.3 真近点角的概念及其求解112.3.4 卫星瞬时位置的求解12第3章 GPS导航解算153.1 GPS导航系统的基本原理153.2 GPS定位中主要误差及消除算法173.3 GPS星历结构及解算过程19第4章 GPS导航解算系统的MATLAB仿真234.1 卫星可见性的估算234.2 GPS卫星运动的MATLAB仿真24结 论32致谢34参考文献35附录37 44成都理工大学毕业设计(论文)第1章 前言1.1 课题背景 GPS系统的前身为美军研制的一种子午仪卫星定位系统(Transit），1958年研制，64年正式投入使用。该系统用5到6颗卫星组成的星网工作，每天最多绕过地球13次，并且无法给出高度信息，在定位精度方面也不尽如人意。然而，子午仪系统使得研发部门对卫星定位取得了初步的经验，并验证了由卫星系统进行定位的可行性，为GPS系统的研制埋下了铺垫。由于卫星定位显示出在导航方面的巨大优越性及子午仪系统存在对潜艇和舰船导航方面的巨大缺陷。美国海陆空三军及民用部门都感到迫切需要一种新的卫星导航系统。1973年12月 ,美国国防部批准它的陆海空三军联合研制新的卫星导航系统: NAVSTAR/GPS。它是英文“Navigation Satellite Timing and Ranging/Global Positioning System” 的缩写词。其意为 “卫星测时测距导航/全球定位系统”,简称 GPS。这个系统向有适当接受设备的全球范围用户提供精确、 连续的三维位置和速度信息 ,并且还广播一种形式的世界协调时(UTC) 。通过遍布全球的（21+3）GPS导航卫星，向全球范围内的用户全天候提供高精度的导航、跟踪定位和授时服务。目前，GPS已在地形测量，交通管理，导航，野外勘探，空间宇宙学等诸多领域得到了广泛的应用。目前全球共有4大GPS系统，分别是：美国 GPS ，由美国国防部于20世纪70年代初开始设计、研制,于1993 年全部建成。1994 年,美国宣布在10年内向全世界免费提供 GPS使用权,但美国只向外国提供低精度的卫星信号。欧盟 “伽利略”，1999 年 欧洲提出计划 ,准备发射 30 颗卫星 ，组成 “伽利略” 卫星定位系统。俄罗斯“格洛纳斯”，尚未部署完毕。始于上世纪70年代需要至少 18 颗卫星才能确保覆盖俄罗斯全境;如要提供全球定位服务 ，则需要 24 颗卫星。中国“北斗”2003 年我国北斗一号建成并开通运行，不同于 GPS,“北斗”的指挥机和终端之间可以双向交流。四川大地震发生后 ,北京武警指挥中心和四川武警部队运用 “北斗” 进行了上百次交流。北斗二号系列卫星今年起将进入组网高峰期 ，预计在 2015 年形成由三十几颗卫星组成的覆盖全球的系统。1.2 本课题研究的意义和方法GPS系统是一个很庞大的系统，包含了天文，地理，计算机，电磁学，通信学，信息学等等。通过本文对GPS的学习研究，最重要的还是要学习其原理：卫星运动原理；卫星定位原理；卫星跟踪原理；卫星导航解算原理等等。通过基础原理的学习，一方面，可以使我们更进一步的理解卫星运动，定位的实现方法及导航的解算原理和方法；通过仿真，进一步了解简单定位的方法及其在仿真平台上的实现途径；另一方面，也可以培养我们自学的能力，训练仿真模拟的技巧和方法。至今，基本上完成了课题的要求，通过不断的注入既定参数，可以更加详细，直观的理解基本的导航解算原理和实现方法！1.3 GPS的发展简史及系统概述 图1-1 GPS系统结构图GPS导航定位以其定位精度高、观测时间短、测站间无需通视、可提供三维导航、操作简便、全天候作业、功能多、应用广泛等特点著称。组成卫星星座的24颗卫星被安排在6个轨道平面上，即每个平面上4颗。一个分布在全世界的地面控制/监视网监视着卫星的运行状态。用GPS信号可以进行海、空和陆地的导航、导弹的制导、大地测量和工程测量的精密定位、时间的传递和速度的测量等。对于测绘领域，GPS卫星定位技术已经用于建立高精度的全国性的大地测量控制网，测定全球性的地球动态参数；用于建立陆地海洋大地测量基准，进行高精度的海岛陆地联测以及海洋测绘；用于检测地球板块运动状态和地壳形变；用于工程测量，成为建立城市与工程控制网的主要手段。用于测定航空航天摄影瞬间相机位置，实现仅有少量的地面控制或无地面控制的航测快速成图，导致地理信息系统、全球环境遥感监测的技术革命4。目前，GPS、GLONASS、INMARSAT等系统都具备了导航定位功能，形成了多元化的空间资源环境。这一多元化的空间资源环境，促使国际民间形成了一个共同的策略，即一方面对现有系统充分利用，一方面积极筹建民间GNSS系统，待2011年左右，GNSS纯民间系统建成，全球 将形成GPS/GLONASS/GNSS三足鼎立之势，才能从根本上摆脱对单一系统的依赖，形成国际共有、国际共享的安全资源环境。世界才可以将卫星导航作为单一导航手段的最高应用境界。国际民间的这一策略，反过来又影响和迫使美国对其GPS使用政策作出更开放的调整。多元化的空间资源环境的确立，给GPS的发发展应用创造了一个前所未有的良好的国际环境。GPS技术广泛应用于武器导航：精确制导导弹、巡航导弹，车辆导航：车辆调度、监控系统，船舶导航：远洋导航、港口/内河引水，飞机导航：航线导航、进场着陆控制，星际导航：卫星轨道定位及个人导航：个人旅游及野外探险等诸多方面。 第2章 GPS的坐标、时间系统和卫星的运动2.1 地球坐标系GPS定位测量当中，要用到两种坐标系，即天球坐标系和地球坐标系。天球坐标系是指坐标原点和各坐标轴的指向在空间是保持不变的，可以很方便的描述卫星的运动和状态。而地球坐标系则是与地球体相关联的坐标系，用于描述地面测量站的位置。所以在导航解算中，我们将采用地球坐标系。要描述一个物体的位置必须要有相关联的坐标系,地球表面的GPS接收机的位置是相对于地球而言的.因此,要描述GPS接收机的位置，需要采用固联于地球上随同地球转动的坐标系、即地球坐标系作为参照系。地球坐标系有两种几何表达形式，即地球直角坐标系（平地球坐标系）和地球大地坐标系（真地球坐标系）。地球直角坐标系的定义是：原点O与地球质心重合，Z轴指向地球北极，X轴指向地球赤道面与格林威治子午圈的交点(即0经度方向)，Y轴在赤道平面里与XOZ构成右手坐标系(即指向东经90度方向)。地球大地坐标系的定义是：地球椭球的中心与地球质心重合，椭球的短轴与地球自转轴重合。地球表面任意一点的大地纬度为过该点之椭球法线与椭球赤道面的夹角 ，经度为该点所在之椭球子午面与格林威治大地子午面之间的夹角 ，该点的高度h为该点沿椭球法线至椭球面的距离。设地球表面任意一点P在地球直角坐标系内表达为P(x，y，z)，在地球大地坐标系内表达为P ( ， ，h)。则两者互换关系为：大地坐标系变为直角坐标系： （2-1）式中：n为椭球的卯酉圈曲率半径，e为椭球的第一偏心率。若椭球的长半径为a，短半径为b，则有 直角坐标系变为大地坐标系，可由下述方法求得由叠代法获得c为地心纬度， ep为椭圆率可设初始值=c进行叠代，直到|i=1-i| 小于某一门限为止。 （2-2）这两种坐标系在定位系统中经常交叉使用，必须熟悉两种坐标系之间的转换关系。图2-1 世界地心坐标系综合上面的，可以得到以下的结论：在GPS 定位系统所用的空间坐标系统当中，我们一般采用天球坐标去研究卫星的空间运动，而采用地球坐标去研究地面监控站点。2.2 时间系统 时间系统是卫星定位测量过程中的一个重要概念。现时的GPS测量的方法是通过接收和处理GPS卫星发射的无线电信号，以确定用户接收机和观测卫星间的距离，然后通过一定的数学方法以确定接收机所在的具体位置，为得到接收机和卫星的准确距离，必须获得无线电信号从卫星传输至接收机这一过程中的精确时间，因而利用卫星技术进行精密的定位和导航，必须要获得高精度的时间信息，这需要一个精确的时间系统。现行的卫星定位测量中与之紧密相关的时间系统有三种:世界时，原子时和动力学时。 世界时是以地球自转为基准的一种时间系统。根据不同的空间参考点，又可分为恒星时，太阳时，世界时三种。原子时以物质内部原子跃迁时所辐射和吸收的电磁波频率来定义的，具有很高的稳定性和复现性，是现时段最为理想的时间系统。在目前的导航定位系统中，均采用了原子时作为其高精度的时间基准。动力学时(DynamicTime，DT)是天体力学中用以描述天体运动的时间单位。当以太阳系质心建立起天休运动方程时，所采用的时间参数称为质心力学时(BaryeeniricDynamicTime，TDB);当以地球质心建立起天体运动方程时，所采用的时间参数称为地球力学时（Terrestrial Dynamic Time, TDT）。2.3 卫星的运动 卫星在空间绕地球运动的时候，除了受到地球重力场引力的作用外，还受到了太阳，月亮和其它的天体引力以及太阳光压，大气的阻力和地球潮汐力的影响。卫星的实际运动轨道非常的复杂，很难用非常精确的数学模型加以描述。在各种力作用对卫星影响的过程当中，以地球的引力场的作用最大，而其它力的影响则相对的小得多。通常把作用到卫星上的力按其影响的大小分成两部分：一类是中心力；一类是摄动力，也称为非中心力。假定地球为均匀球体的地球引力，称为在心力，它决定了卫星运动的基本规律和基本特征，由此决定地球的轨道，可以视为理想的轨道。非中心力包括地球非球形对称的作用力，日、月引力，大气阻力，光辐射压力以及地球的潮汐力等。摄动力的作用，使卫星偏离了既定的理想轨道。而在它影响下，卫星的运动称为卫星的受摄运动。而上述理想状态的卫星运动则称为无摄运动。卫星在地球的引力场当中所做的无摄运动，也称为开普勒运动，其规律可以由开普勒三大定律来描述。2.3.1开普勒定律 开普勒第一定律：卫星运动的轨道是个椭圆，而该椭圆的一个焦点和地球的质心重合。这一个定律表明了，在中心引力的作用下，卫星绕地球轨道运行的轨道面，是一个通过地球质心的精致平面。轨道椭圆一般称期为开普勒椭圆，其形状和大小都不变。在轨道上，卫星离地球质心远的一点叫做远地点，近的一点就做近地点。轨道图形可以表示为如图2-2： bs图2-2 卫星的椭圆运行轨道 Ms fs asM 远地点P 近地点P 卫星绕地球质心运动的轨道方程为：R= (2-3) 在该式当中，R是卫星的地心距离；as为开普勒椭圆的长半径；es为开普勒椭圆的偏心率；fs为真近点角，它描述了任意时刻，卫星在轨道上面相对于近地点的位置，是时间的函数，其定义见上图所示。 开普勒定义定律阐述了卫星运动轨道的基本形态及其与地心的关系。C开普勒第二定律：卫星的地心向径，即地球质心与卫星质心间的距离向量，在相同的时间内所扫过的空间面积是相等的。（如图2-3所示） 图2-3 卫星地心向径在相同的时间间隔内扫描的面积DPOO P与任何其它的运动物体一样，在轨道上面运动的卫星，也具有两种的能量：位能和动能。位能就是指仅仅受到地球重力场的影响，其大小和卫星的在轨高度有关。在近地点其位能最小，而在远地点其位能最大。卫星在任一个时刻t所具有的位能为 (G为万有引力常量，M为地球的质量,ms为卫星的质量)。动能则是由卫星的运动所引起的，其大小是卫星的运动速度的函数。如果取卫星的运动的速度为vs,则其动能为。根据能量守恒定律，卫星的势能与动能的总量是不变的，即 =常量 (2-4)因此，当卫星运行到近地点的时候，其动能最大；在远地点的时候，其动能最小，由此，开普勒第二定律所包涵的内容是：卫星在椭圆轨道上的运行速度是不断变化的，在近地点处的速度最大，而在远地点的速度最小。开普勒第三定律：卫星运动周期的平方与轨道椭圆长半径的立方之比为一个常数，而该常数等于地球引力常数和地球质量的乘积GM的倒数。其数学表达式为： (2-5)在这个式子当中，为卫星的运行周期.如果我们假设卫星的平均角速度为N,则有:N= (2-6)于是，开普勒第三定律（2-4）可以写成： (2-7)或者表示为常用的形式：N= (2-8)显然，当开普勒的长半径确定了过后，卫星运动的平均角速度就得到了确定，且保持不变。2.3.2 无摄卫星运动的轨道参数卫星的无摄运动，一般的可以由下面的6个参数（如图2-4）来描述：图2-4开普勒轨道参数 As-卫星轨道的长半径 Es -卫星轨道的偏心率 -升交点的赤径 i-卫星轨道面的倾角 s-近地点角距，即升交点与近地点的夹角 fs-卫星的真近点角，在轨道平面上为卫星与进地点的地心角距。当这6个参数一旦确定后，卫星在任意瞬时的相对于地球的空间位置及其速度，就被唯一的确定了！2.3.3 真近点角的概念及其求解在描述卫星无摄运动的6个参数当中，只有fs是关于时间的函数，其他的都是一般的参数。所以，计算卫星瞬时的位置的关键，计算出参数fs，并由此确定卫星的空间位置及其和时间的关系。为此，需要引进两个参数Es和Ms去计算真近点角。Es:偏近点角，如果定义过卫星质心做平行与椭圆短半轴的直线，M为该直线与近地点到椭圆中心连线的交点，则椭圆平面上近地点P到M的圆弧所对应的圆心角就是Es。Ms:平近点角。它是一个假设量，如果卫星在轨道运行的平速度为n，则平近点角定义为：Ms=n(t-t0) (2-9)t0为卫星过近地点的时刻，t为观察卫星的时刻。由上面的式子知道，卫星的平近点角仅仅为卫星平均速度的时间的函数，对于一个确定的卫星来说，这个参数可以认为是常数。3其中Ms与Es有关系如下：Ms=Es-essinEs (2-10)为了计算卫星的瞬时速度， 需要确定卫星运行的真近点角fs。由于有以下的关系成立：ascosEs=rcosfs+ases (2-11)于是将上式带入到 (2-3)中就得到：Cosfs= (2-12)或者得到以下常用的形式：Tan()= (2-13)2.3.4 卫星瞬时位置的求解 对于任意的观测时刻，根据卫星的平均运行速度n,根据2-9，2-10，2-11，便可以唯一确定真近点角fs。这样，卫星于任一观测历元t，相对于地球瞬间空间的位置便可以随之确定。若以直角坐标的原点 与地心M重合，轴指向近地点且垂直于轨道的平面，轴在轨道平面上垂直轴构成右手关系。于是，卫星任意时刻的轨道坐标可以表示成为： cosfs = sinfs (2-14) 0 而由上面的分析，可以得到： = cosE- = SinE (2-15) = 0而要把这个轨道坐标系坐标表示成为天球坐标的话，由于他们的坐标原点都是地球的质心，但是坐标轴的指向是不相同的，为了使他们的坐标轴相同，应该将坐标系（,）依次做下面的变化旋转：（1） 绕轴顺时针旋转角度s，使轴的指向由近地点变为升交点。（2） 绕轴顺时针旋转角度i，使轴与Z轴相同。（3） 绕轴顺时针旋转角度，使轴指向春分点。实现上述三步的旋转矩阵分别为R1，R2，R3，即： Coss -Sins 0R1= Sins Coss 0 (2-16) 0 0 1 1 0 0 R2= 0 Cosi -Sini (2-17) 0 Sini Cosi Cos -Sin 0 R3= Sin Cos 0 (2-18) 0 0 1于是得到了在天球坐标系下面的卫星位置坐标可以表示为： X CosE-es Y =R3R2R1 =R3R2R1as SinE (2-19) Z 0利用转换关系Rz(g)得到相应的地球坐标系的坐标： Cosg Sing 0而Rz(g) = -Sing Cosg 0 (2-20) 0 0 1 图2-5 描述卫星的运行及其轨道的参数当我们在一个地方同时可见的卫星如果多于四颗（GPS卫星的轨道设计和运动时间的安排使得用户在地球的任意位置（两极个别地点除外），都能够看到411颗的卫星）的时候，我们可以用最小二乘法去解算未知数，这样，充分的利用了已知的数据信息，使得结果的偏差最小化。第3章 GPS导航解算GPS用户部分的核心是GPS接收机。其主要由基带信号处理和导航解算两部分组成。其中基带信号处理部分主要包括对GPS卫星信号的二维搜索、捕获、跟踪、伪距计算、导航数据解码等工作。导航解算部分主要包括根据导航数据中的星历参数实时进行各可视卫星位置计算；根据导航数据中各误差参数进行星钟误差、相对论效应误差、地球自转影响、信号传输误差(主要包括电离层实时传输误差及对流层实时传输误差)等各种实时误差的计算，并将其从伪距中消除；根据上述结果进行接收机PVT（位置、速度、时间）的解算；对各精度因子(DOP)进行实时计算和监测以确定定位解的精度。图3-1 GPS接收机的基本构成框图本章中重点讨论GPS接收机的导航解算部分，基带信号处理部分可参看有关资料。本文讨论的假设前提是GPS接收机已经对GPS卫星信号进行了有效捕获和跟踪，对伪距进行了计算，并对导航数据进行了解码工作。3.1 GPS导航系统的基本原理GPS导航系统的基本原理是测量出已知位置的卫星到用户接收机之间的距离，然后综合多颗卫星的数据就可知道接收机的具体位置。16要达到这一目的，卫星的位置可以根据星载时钟所记录的时间在卫星星历中查出。而用户到卫星的距离则通过纪录卫星信号传播到用户所经历的时间，再将其乘以光速得到(由于大气层电离层的干扰，这一距离并不是用户与卫星之间的真实距离，而是伪距（PR）：当GPS卫星正常工作时，会不断地用1和0二进制码元组成的伪随机码（简称伪码）发射导航电文。GPS 系统使用的伪码一共有两种，分别是民用的C/A码和军用的P(Y)码。C/A码频率1.023MHz,重复周期一毫秒，码间距1微秒，相当于300m；P 码频率10.23MHz，重复周期266.4天，码间距0.1微秒，相当于30m。而Y码是在P码的基础上形成的，保密性能更佳。导航电文包括卫星星历、工作状况、时钟改正、电离层时延修正、大气折射修正等信息。它是从卫星信号中解调制出来，以50b/s调制在载频上发射的。导航电文每个主帧中包含5个子帧每帧长6s。前三帧各10个字码；每三十秒重复一次，每小时更新一次。后两帧共15000b。导航电文中的内容主要有遥测码、转换码、第1、2、3数据块，其中最重要的则为星历数据。当用户接受到导航电文时，提取出卫星时间并将其与自己的时钟做对比便可得知卫星与用户的距离，再利用导航电文中的卫星星历数据推算出卫星发射电文时所处位置，用户在WGS-84大地坐标系中的位置速度等信息便可得知。可见GPS导航系统卫星部分的作用就是不断地发射导航电文，如图3-2。然而，由于用户接受机使用的时钟与卫星星载时钟不可能总是同步，所以除了用户的三维坐标x、y、z外，还要引进一个t即卫星与接收机之间的时间差作为未知数，然后用4个方程将这4个未知数解出来。所以如果想知道接收机所处的位置，至少要能接收到4个卫星的信号。图3-2导航电文的组成GPS接收机可接收到可用于授时的准确至纳秒级的时间信息；用于预报未来几个月内卫星所处概略位置的预报星历；用于计算定位时所需卫星坐标的广播星历，精度为几米至几十米（各个卫星不同，随时变化）；以及GPS系统信息，如卫星状况等。GPS接收机对码的量测就可得到卫星到接收机的距离，由于含有接收机卫星钟的误差及大气传播误差，故称为伪距。对0A码测得的伪距称为UA码伪距，精度约为20米左右，对P码测得的伪距称为P码伪距，精度约为2米左右。GPS接收机对收到的卫星信号，进行解码或采用其它技术，将调制在载波上的信息去掉后，就可以恢复载波。严格而言，载波相位应被称为载波拍频相位，它是收到的受多普勒频移影响的卫星信号载波相位与接收机本机振荡产生信号相位之差。一般在接收机钟确定的历元时刻量测，保持对卫星信号的跟踪，就可记录下相位的变化值，但开始观测时的接收机和卫星振荡器的相位初值是不知道的，起始历元的相位整数也是不知道的，即整周模糊度，只能在数据处理中作为参数解算。相位观测值的精度高至毫米，但前提是解出整周模糊度，因此只有在相对定位、并有一段连续观测值时才能使用相位观测值，而要达到优于米级的定位 精度也只能采用相位观测值。按定位方式，GPS定位分为单点定位和相对定位（差分定位）。单点定位就是根据一台接收机的观测数据来确定接收机位置的方式，它只能采用伪距观测量，可用于车船等的概略导航定位。相对定位（差分定位）是根据两台以上接收机的观测数据来确定观测点之间的相对位置的方法，它既可采用伪距观测量也可采用相位观测量，大地测量或工程测量均应采用相位观测值进行相对定位。在GPS观测量中包含了卫星和接收机的钟差、大气传播延迟、多路径效应等误差，在定位计算时还要受到卫星广播星历误差的影响，在进行相对定位时大部分公共误差被抵消或削弱，因此定位精度将大大提高，双频接收机可以根据两个频率的观测量抵消大气中电离层误差的主要部分，在精度要求高，接收机间距离较远时（大气有明显差别），应选用双频接收机。 3.2 GPS定位中主要误差及消除算法 GPS定位中的主要误差有：星钟误差，相对论误差，地球自转误差，电离层和对流层误差。211）星钟误差星钟误差是由于星上时钟和GPS标准时之间的误差形成的，GPS测量以精密测时为依据，星钟误差时间上可达1ms，造成的距离偏差可达到300Km，必须加以消除。一般用二项式表示星钟误差。 （3-1）GPS星历中通过发送二项式的系数来达到修正的目的。经此修正以后，星钟和GPS标准时之间的误差可以控制在20ns之内。2）相对论误差由相对论理论，在地面上具有频率 的时钟安装在以速度 运行的卫星上以后，时钟频率将会发生变化，改变量为： （3-2）即卫星上时钟比地面上要慢，要修正此误差，可采用系数改进的方法。GPS星历中广播了此系数用以消除相对论误差，可以将相对论误差控制在70ns以内。3）地球自转误差GPS定位采用的是与地球固连的协议地球坐标系，随地球一起绕z轴自转。卫星相对于协议地球系的位置(坐标值)，是相对历元而言的。若发射信号的某一瞬间，卫星处于协议坐标系中的某个位置，当地面接收机接收到卫星信号时，由于地球的自转，卫星已不在发射瞬时的位置坐标值)处了。也就是说，为求解接收机接收卫星信号时刻在协议坐标系中的位置，必须以该时刻的坐标系作为求解的参考坐标系。而求解卫星位置时所使用的时刻为卫星发射信号的时刻。这样，必须把该时刻求解的卫星位置转化到参考坐标系中的位置。设地球自转角速度为 we，发射信号瞬时到接收信号瞬时的信号传播延时为t ，则在此时间过程中升交点经度调整为则三维坐标调整为 （3-3）地球自转引起的定位误差在米级，精密定位时必须考虑加以消除。4）电离层和对流层误差电离层是指地球上空距地面高度在50-1000km 之间的大气层。电离层中的气体分子由于受到太阳等天体各种射线辐射，产生强烈的电离，形成大量的自由电子和正离子。电离层误差主要有电离层折射误差和电离层延迟误差组成。其引起的误差垂直方向可以达到50米左右，水平方向可以达到150米左右。目前，还无法用一个严格的数学模型来描述电子密度的大小和变化规律，因此，消除电离层误差采用电离层改正模型或双频观测加以修正。对流层是指从地面向上约40km范围内的大气底层，占整个大气质量的99%。其大气密度比电离层更大，大气状态也更复杂。对流层与地面接触，从地面得到辐射热能，温度随高度的上升而降低。对流层折射包括两部分：一是由于电磁波的传播速度或光速在大气中变慢造成路径延迟，这占主要部分；二是由于GPS卫星信号通过对流层时，也使传播的路径发生弯曲，从而使测量距离产生偏差。在垂直方向可达到2.5米，水平方向可达到20米。对流层误差同样通过经验模型来进行修正。GPS星历中通过给定电离层对流层模型以及模型参数来消除电离层和对流层误差。实验资料表明，利用模型对电离层误差改进有效性达到75%，对流层误差改进有效性为95%。3.3 GPS星历结构及解算过程 要得到接收机的位置，在接收机时钟和GPS标准时严格同步的情况下，则待求解位置是3个未知变量，需要3个独立方程来求解。但是实际情况中，很难做到接收机时钟和GPS标准时严格同步，这样，我们把接收机时间和GPS标准时间偏差也作为一个未知变量，这样，求解就需要4个独立方程，也就是需要有4颗观测卫星。图3-3 GPS定位示意图（未考虑时间偏差）假设接收机位置为（xu，yu，zu） ，接收机时间偏差为 tu，则由于时间偏差引起的距离偏差为为得到的伪距观测值。我们可以得到联立方程 （3-4）将上式线性化，即在真实位置（xu，yu，zu）进行泰勒级数展开，忽略高次项，得到 （3-5）其中，式（3-5）即为实际计算的叠代公式，叠代终止条件是真实位置 （xu，yu，zu）的变化量小于某一个阈值，最终得到可以作为调整接收机时间偏差的依据，计算一般采用矩阵方式求解。要求解该方程，我们还需要预先知道4颗卫星的位置 （xj，yj，zj），而卫星位置可以从该卫星的星历中获得。GPS卫星星历给出了本星的星历，根据星历可以算出卫星的实时位置，并且星历中给出了消除卫星星钟误差、相对论误差、地球自转误差、电离层和对流层误差的参数，根据这些参数计算出的卫星位置，可以基本上消除上述误差。求解卫星位置的基本步骤为：计算卫星运行平均角速度计算归化时间；计算观测时刻的平近点角；计算偏近点角；计算卫星矢径；计算卫星真近点角；计算升交点角距；计算摄动改正项；计算经过摄动改正的升交距角、卫星矢径、轨道倾角；计算观测时刻的升交点经度；计算卫星在地心坐标系中的位置。特别值得指出的是，在计算卫星真近点角Vk时，应采用公式 （3-6）其中，e为偏心率， Ek为卫星偏近点角。有部分参考书籍计算卫星真近点角的公式有误，会导致卫星真近点角 的象限模糊问题，从而无法得到卫星正确位置。 进行上述计算后，再根据星历中广播的各误差参数进一步消除各项误差。这样，我们就得到一个完整的利用GPS星历进行导航定位解算的过程。第4章 GPS导航解算系统的matlab仿真4.1 卫星可见性的估算当初，卫星星座的设计要求在全球范围内任何时候，任何位置都必须保证至少四颗以上的卫星导航信号。换句话说，并不是所有的卫星都能够被一个用户所看见。一颗卫星信号能否被接收与下列因素是有关系的7：1 地球是否影响了该GPS卫星信号的传播。图4-1说明了地球对GPS信号的影响。如果卫星处于图中的阴影部分，则对图中的飞机是不能够接收到该卫星的信号的!图4-1 可见卫星的测量原理 E 飞机 h Rn 2 GPS接收机是否位于该GPS卫星发射天线的范围内。这种情况主要针对的是航天器上面的GPS接收机。GPS卫星信号的发射张角大约为21.3度，大于卫星到水平面的张角13.9度，这样就保证了一些飞行高度较高的航天器在高空可以更多的接收到GPS卫星的信号。但是对于那些超出发射角的飞行器就收不到信号了。3利用卫星地心用户之间的张角，可以估算出来可卫星信号能否被收到。如果这个角度小于90度的话，就可以收到，反之不能够收到。特别的，当这个角度刚刚为90度的时候，我们一般认为是收不到的。4.2 GPS卫星运动的MATLAB仿真首先在进行MATLAB仿真时要创建函数，可以避免代码的复制，提高代码重用性，以减少整体代码的复杂性，增强可读性和可管理性；此外也可以提高移植性，更易于调试和隔离错误，每一段程序最优化了，程序整体的性能也就提高了。下面的5个程序是对本章节相关内容仿真实现的子程序，通过主程序的调用完成相关仿真。主程序见附录。 %程序一：利用不同颜色的曲线绘制三维图像函数function plot3c(x,y,z,color)%封装了Matlab的Plot3函数switch (color)%r代表红色，g代表绿色，b代表蓝色，c代表青绿色，%m代表品红色，y代表黄色，k代表黑色，w代表白色%plot3语句为Matlab中绘制三维曲线的库函数 case 0 plot3(x,y,z,w-); case 1 plot3(x,y,z,r-); case 2 plot3(x,y,z,g-); case 3 plot3(x,y,z,c-); case 4 plot3(x,y,z,m-); case 5 plot3(x,y,z,y-); case 6 plot3(x,y,z,b-); case 7 plot3(x,y,z,k-); otherwiseend%程序二 绘制立方体，x,y,z代表立方体中心，而Lx,Ly,Lz为立方体的长宽高function boxplot3(x,y,z,Lx,Ly,Lz,color)x0=x-Lx/2;y0=y-Ly/2;z0=z-Lz/2;x=x0 x0 x0 x0 x0+Lx x0+Lx x0+Lx x0+Lx;y=y0 y0 y0+Ly y0+Ly y0 y0 y0+Ly y0+Ly;z=z0 z0+Lz z0+Lz z0 z0 z0+Lz z0+Lz z0;index=zeros(6,5);%定义index向量index(1,:)=1 2 3 4 1;index(2,:)=5 6 7 8 5;index(3,:)=1 2 6 5 1;index(4,:)=4 3 7 8 4;index(5,:)=2 6 7 3 2;index(6,:)=1 5 8 4 1;for k=1:6 plot3c(x(index(k,:),y(index(k,:),z(index(k,:),color)hold on %后续图形曲线叠加到前面的图形窗口中显示end%程序三 绘制两个立方体（表示卫星）function drawsatellite(moveX,moveY,moveZ,color)Length=500;Width=500;Height=500;boxplot3(moveX,moveY,moveZ,Length,Width,Height,color);boxplot3(moveX,moveY,moveZ,Length,Width*5,Height/10,color);%程序四 绘制卫星轨迹function drawsatelliteorbita=26560;% 卫星轨道的长半轴.e=0.02;%e是椭圆的偏心率.E=0:0.1:2*pi;x=a*(cos(E)-e);y=a*sqrt(1-e2)*sin(E);z=0*E;DtoR=2*pi/360;A1=32.8 92.8 152.8 212.6 272.8 332.8;%卫星星座数据.参看图4-2for k=1:6 A=A1(k)*DtoR;%升交点的经度 B=55*DtoR;%轨道的倾角 C=pi/100;%近地点幅角 %总共有6个卫星轨道平面 R1=cos(A) -sin(A) 0;%分别依据公式（2-16.2-17.2-18) sin(A) cos(A) 0; 0 0 1; R2=1 0 0; 0 cos(B) -sin(B); 0 sin(B) cos(B); R3=cos(C) -sin(C) 0; sin(C) cos(C) 0; 0 0 1; L1=length(E); R312=R1*R2*R3; Ans=R312*x;y;z;%将卫星轨道坐标系转换到地球赤道坐标 x1=Ans(1,:); y1=Ans(2,:); z1=Ans(3,:); plot3c(x1,y1,z1,k);%调用程序一 hold on; axis equal;%纵、横坐标轴采用等长刻度 axis off;%取消坐标轴end 图4-2 GPS星座的平面投影%程序五绘制地球function drawearth()%time 是参数%利用这个参数,可以绘制一个看起来是旋转的地球r=6400;time=0;j1=0:pi/10:2*pi;w1=-pi/2:pi/10:pi/2;L1=length(w1);L2=length(j1);for n=1:L1 z=ones(L2,1); z=z*r*sin(w1(n); temp=r*cos(w1(n); x=temp*sin(j1); y=temp*cos(j1); plot3(x,y,z);% 绘制三维图形显示 hold on; grid;end%figure(3);unit=ones(1,1);z0=ones(1,1);x0=ones(1,1);y0=ones(1,1);for n=1:L2 %n=7; for m=1:L1 temp=w1(m); temp2=j1(n)+time*pi/12; z=r*sin(temp); x=r*cos(temp)*sin(temp2); y=r*cos(temp)*cos(temp2); z1=unit*z; x1=unit*x; y1=unit*y; z0=z0 z1; x0=x0 x1; y0=y0 y1; end z0(:,1)=; x0(:,1)=; y0(:,1)=; plot3(x0,y0,z0); axis equal; axis off; hold on; boxplot3(0,0,0,100,100,100,7); tempx=6400; tempy=0; tempz=0; cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,1); tempx=0; tempy=6400; tempz=0; cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,2); tempx=0; tempy=0; tempz=6400; cube=100; boxplot3(tempx,tempy,tempz,cube,cube,cube,4); end 仿真主程序一：模拟运动的卫星，如图4-3程序见附录：图4-3 模拟运动的卫星因为程序仿真结果为动态仿真结果，获得影像动画播放12次影像动画，每秒显示2帧.由此可知卫星的天空瞬时位置是随着时间的变动而发生变化的。对于我们在地球上面的一个用户来说，一天当中的不同时刻看到的卫星是不相同的。这个仿真程序的功能实际上就是仿真了在一个轨道上面的卫星在不同的时候