GPS导航系统毕业论文.doc

第1 章 绪论1.1 导航的基本概念导航是引导运载体到达预定目的地的过程。导航分两类：(1)自主式导航，用飞行器或船舶上的设备导航，有惯性导航、多普勒导航和天文导航等；(2)非自主式导航，用于飞行器、船舶、汽车等交通设备与有关的地面或空中设备相配合导航，有无线电导航、卫星导航。在军事上，导航还要配合完成武器投射、侦察、巡逻、反潜和援救等任务。高效、高精度的导航系统更是我国这种发展中国家赶超发达国家的战略性资源和倍能器。在军用方面，随着新时期军事战略方针的转变及高新技术武器装备的发展，导航定位定向系统已经成为我军现代化建设中一项不可缺少的重要军事技术装备，其重要性表现在：它是信息战必不可少的基础设备，是建立战场统一坐标的前提，是快速、准确火力部署的保障，同时又是实现武器精确打击能力的必要条件。所以，导航定位定向系统对迅速提高我军的综合作战能力，加快数字化部队建设至关重要；在民用方面，国外的导航定位定向系统己在大地测量、定向钻并、隧道掘进、地面车辆导航、飞机进场着陆、航天航空遥感、机载重力测量、公路监测、地下油气管道监测、矿井监测、激光断面监测等方面得到广泛地的应用，并取得了巨大的经济效益。在日常生活中我们经常接触到的导航是车载导航，车载导航属于非自主式导航，车载导航是利用车载gps（全球定位系统）配合电子地图来进行的，汽车gps导航系统由两部分组成： 一部分由安装在汽车上的gps接收机和显示设备组成；另一部分由计算机控制中心组成，两部分通过定位卫星进行联系。1.2 惯性导航(ins)概述通常说的惯性技术，是惯性器件、惯性测量、惯性导航、惯性制导和惯性稳定等技术的统称。惯性技术既是一门学科，也是一门工程技术，在陆、海、空、天各个领域有着广泛应用。惯性器件（陀螺仪和加速度计）、惯性仪表、惯性导航系统都是以牛顿力学定律为基础的。惯性导航系统通过加速度计实时测量载体运动的加速度，经积分运算得到载体的实时速度和位置信息。惯性技术是对载体进行导航的关键技术之一，惯性技术是利用惯性原理或其它有关原理，自主测量和控制运载体运动过程的技术，惯性测量和惯性敏感器技。现代惯性技术在各国政府雄厚资金的支持下，己经从最初的军事应用渗透到民用领域。惯性技术在国防装备技术中占有非常重要的地位。对于惯性制导的中远程导弹，一般说来命中精度 70%取决于制导系统的精度。对于导弹核潜艇，由于潜航时间长，其位置和速度是变化的，而这些数据是发射导弹的初始参数，直接影响导弹的命中精度，因而需要提供高精度位置、 速度和垂直对准信号。目前适用于潜艇的唯一导航设备就是惯性导航系统。1.3 卫星导航概述卫星导航是采用导航卫星对地面、海洋、空中和空间用户进行导航定位的技术。卫星导航系统通过测定的载体到导航卫星的距离、距离差等参数，并结合获取的或计算机得到的导航卫星瞬时位置来确定载体位置。目前主要有美国的gps、俄罗斯的glonass、欧洲的galileo和我国的北斗导航卫星系统。这里主要介绍gps导航定位。gps系统主要由空间部分、地面部分及用户设备组成，gps的空间部分由24颗卫星组成，包括21颗工作星和三颗备份卫星，位于距地面20200km高度的近圆轨道上，24颗卫星均匀分布在6个轨道面上，轨道倾角为55度，这样可以确保地球上的任何一个地方都可以同时观测到611颗卫星。布置在轨道上的备份星，可随时进入工作状态。gps卫星导航系统具有全球覆盖、全天候工作、定位精度高和用途广的优点。地面控制部分由一个主控站、三个注入站和若干监测站组成，主控站通过注入站向所有gps卫星注入新的导航信息，确保卫星的导航数据和时钟信息的精确性。用户设备部分包括接收机硬件和机内软件以及gps数据的后处理软件包，其主要功能是能够捕获到按一定卫星截止角所选择的的待测卫星，并跟踪这些卫星的运行。当接收机捕获到跟踪的卫星信号后，就可测量出接收天线至卫星的伪距和伪距率，解调出卫星轨道参数等数据。1.4 组合导航概述组合导航是现代导航理论和技术发展的必然结果。每种单一的导航系统都存在着各自的独特性能和局限性，几种不同的单一系统组合在一起，利用多种信息源互相补充，构成一种具有更多余度和更高导航精度的多功能系统即为组合导航系统。根据不同的应用要求与目的，可以构成不同的组合导航系统，由于惯性导航系统的自主性，目前多以惯性导航系统作为主导航系统构成组合导航系统。1.5 组合导航系统平台开发的现状及发展趋势鉴于各种单一的各种导航系统有这样那样的缺点，上世纪70年代，现代组合导航系统在航海、航空与航天等领域随着现代高科技的发展应运而生。随着电子计算机技术特别是微机技术的迅猛发展和现代控制系统理论的进步，组合导航技术开始迅猛发展起来。过去单独使用的各种导航设备，通过微型电子计算机有机的组合到一起，发展各自的特点、扬长避短，组合导航成为目前导航技术发展的重要方向之一。从上个世纪80年代开始，英、美、法等国的军方和民用部门开始了对gps/ins组合导航平台的研究。80年代后期欧美各国在ms visual studio c+6.0仿真软件平台上，实现了ins/gps组合导航的解算，同时在windows2000 操作界面上实现了曲线和数据的同步显示。近年来我国各大高校采用visual c+6.0编程，在windows操作系统下实现了人机界面和直观的数据显示。该平台将组合导航的理论方案与导航设备相结合，可以进行室内试验系统的调试和研究，还可以提供直观的界面显示、数据图表等。1.6 课题研究的意义惯性导航系统( ins) 通过加速度计实时测量载体运动的加速度，经积分运算得到载体的实时速度和位置信息，是一种完全自主的导航系统,具有不依赖外界信息、隐蔽性好、抗辐射性强、全天候等优点,是机载设备中能提供多种导航参数的重要导航设备。惯性导航系统一般由惯性传感器模块、导航解算模块、电路系统、电源模块、滤波模块和外壳等部分组成。但它的定位误差随时间而积累, 长时间工作后会产生大的误差, 使得惯性导航系统不宜作远距离导航。全球定位系统(gps)继惯性导航以后导航技术的又一大进展，其由空间卫星、地面监控设施以及用户接收机三部分组成。具有较高的导航精度,但是该系统不能提供如载体姿态等导航参数,且在载体上使用时,由于载体的机动运动,常使接收机不易捕获和跟踪卫星的载波信号,甚至对已跟踪的信号失锁。因此,为了克服惯性导航与全球定位系统的缺点,多根据ins和gps的导航功能互补的特点,以适当的方法将两者组合来提高系统的整体导航精度及导航性能以及空中对准和再对准的能力gps 接收机在惯导位置和速度信息的辅助下,也将改善捕获、跟踪和再捕获的能力,并在卫星分布条件差或可见星少的情况下导航精度不致下降过大。由于优点显著,gps/ ins 组合系统被一致认为是载体最理想的组合导航系统。1.7 课题研究内容及方案平台开发内容的程序包括对微惯性传感器/组件的实时采集、曲线绘制、显示和存储。组合导航平台可以实现mems imu和gps等姿态、位置、速度等解算结果的数据显示，曲线显示，数据存储等功能。研究方案如下：1.捷联惯导系统的工作原理；2.推导捷联惯导系统姿态、位置、速度更新计算，采用matlab/vb编程测试；3. 构建ins/gps导航系统界面；4. 完成对组合导航系统的实时采集、曲线的绘制、显示和存储程序编制；5. 针对不同导航任务要求，编制系统测试界面与后台运行程序。第 2 章 惯性导航系统及gps基本原理2.1 惯性导航系统的组成及分类惯性导航系统一般由惯性传感器模块、导航解算模块、电路系统、电源模块、滤波模块和外壳等部分组成。其中，惯性传感器模块是惯性导航系统的核心，一般由三个陀螺仪和加速度计通过正交装配构成。根据构建导航坐标系方法和途径的不同，可将惯性导航系统分为两种类型：采用物理平台模拟导航坐标系的系统称为平台式惯性导航系统；采用数学算法确定导航坐标系的系统称为捷联式惯性导航系统。平台式惯性导航系统是用物理平台直接模拟导航坐标系，导航计算比较简单。陀螺稳定平台能够隔离载体的角运动，给惯性器件提供较好的工作环境，系统的精度较高，但是物理平台本身结构复杂、体积大、制造成本高。捷联式惯性导航系统由于没有实体物理平台，结构简单、体积小、维护方便，但惯性器件直接安装在运载体上，工作环境恶劣，对惯性器件的环境适应性要求很高。同时，由于加速度计输出的加速度分量是沿载体坐标系轴向的，需经计算转换到某种导航坐标系中去，计算量大。惯性导航完全不依赖外面的声、光、电、磁等传播信号，可以实时、高精度地输出所需要的全部导航参数信息，自主地进行定位、导航，不受地域的限制，不受自然环境和人为干扰的影响，隐蔽性好，不论外太空、空间、地面、地下、水面及水下都能全天候的可靠工作，这是其他导航技术，如天文导航、无线电导航与定位、卫星导航等无法实现的。这些独特的优点使其成为国防、航天、航空、船舶与海洋、陆地交通等领域十分重要、不可替代的导航手段。2.2 捷联式惯性导航系统的工作原理“捷联”（strapdown）这一术语的英文原意是“捆绑”因此所谓的捷联式惯性导航系统就是将惯性测量装置的敏感器（陀螺仪和加速度计）直接固联在载体上。陀螺仪用于测量载体坐标系的3个轴的角速度信息，并将信息传输给导航计算机，经过误差补偿计算后进行姿态矩阵计算。加速度计组件用于测量载体坐标系的3个轴的加速度信息，并将信息同样传输到导航计算机，经过误差补偿计算后进行比力坐标系的变换。加速度计测量的是载体坐标系(b系)相对于惯性空间的加速度在载体坐标系中的投影，该测试量也称为比力。而对于捷联惯导系统，导航计算机要在导航坐标系中完成，因此，首先要将机体系中的测试量转换导航坐标系中的物理量，即实现由机体坐标系到导航坐标系的坐标转换。这一转换由姿态矩阵完成，而是利用陀螺仪的输出即载体相对惯性空间转动的角速率在载体坐标系下的投影计算得到。姿态矩阵是随时间的变化而不断变化的。另外，从姿态矩阵中可以单值的确定飞行器的姿态角。捷联式惯导系统中需要实时地求取姿态矩阵，以便提取飞行器姿态角（俯仰角、滚动角、航偏角）以及变换比力。捷联惯导系统原理简图如图2.1所示图2.1 捷联惯导系统原理简图2.3 常用坐标系惯性导航中所采用的坐标系可分为惯性坐标系与非惯性坐标系两类。惯性导航区别于其它类型的导航方案(如无线电导航、天文导航等)的根本不同之处就在于其导航原理是建立在牛顿力学定律一一又可称为惯性定争一的基础上的，“惯性导航”也因此而得名。然而牛顿力学定律是在惯性空间内成立的，这就首先有必要引入惯性坐标系，作为讨论惯性导航基本原理的坐标基准。我们知道，对飞行器进行导航的主要目的就是要实时地确定其导航参数，如飞行器的姿态、位置、速度等。飞行器的导航参数就是通过各个坐标系之间的关系来确定的，这些坐标系是区别于惯性坐标系、并根据导航的需要而选取的。我们称之为非惯性坐标系，如地球坐标系、地理坐标系、导航坐标系、平台坐标系以及机体坐标系等等。在惯性导航中常用的坐标系有以下几种：(1)地心惯性坐标系(下标为i) oexiyizi惯性坐标系是符合牛顿力学定律的坐标系，即是绝对静止或只做匀速直线运动的坐标系。地心惯性坐标系的原点oe。选在地球的中心，zi轴选在沿地轴指向北极的方向上，而xi、 yi轴则在地球的赤道平面内，并指向空间的两颗恒星，xi yi zi构成右手坐标系。三个坐标轴指向惯性空间固定不动，此坐标系是惯性仪表测量参考基准。(2)地球坐标系(下标为“)oxeyeze地球坐标系是固连在地球上的坐标系，它相对惯性坐标系以地球自转角速率e旋转，e =15.041070/小时。地球坐标系的原点在地球中心oe，oeze轴与oezi轴重合，oexeye在赤道平面内，xe轴指向格林威治经线，ye轴指向东经90度方向。(3)地理坐标系(下标为t) oxiyz地理坐标系是在载体上用来表示载体所在位置的东向、北向和垂线方向的坐标系。地理坐标系的原点o选在载体重心处，xi指向东，yi指向北，zi沿垂线方向指向天(东北天)。(4) 导航坐标系(下标为n) oxnynzn导航坐标系是在导航时根据导航系统一z作的需要而选取的作为导航基准的坐标系。 当把导航坐标系选得与地理坐标系重合时，可将这种导航坐标系成为指北方位系统;为了适应在极区附近导航的需要往往将导航坐标系的zn轴仍选的与zi轴重合，而使xn与xi及yn与yi之间相差一个自由方位角或游动方位角a，这种导航坐标系可称为自由方位系统或游动自由方位系统。(5)机体坐标系(下标为b) oxbybzb机体坐标下是固连在机体上的坐标系。机体坐标系的坐标原点o位于飞行器的重心处，xb沿机体横轴指向右，yb沿机体纵轴指向前，zb垂直于oxb yb，并沿飞行器的竖轴指向上oxbybz。构成右手坐标系。机体坐标系相对地理坐标系的方位为飞机的姿态和航向。 (6)平台坐标系(下标为p) -oxpypzp 平台坐标系是用惯导系统来复现导航坐标系时所获得的坐标系。平台坐标系的坐标原点。位于飞行器的重心处。当惯导系统不存在误差时，平台坐标系与导航坐标系相重合;当惯导系统出现误差时，平台坐标系就要相对导航坐标系出现误差角，就是平台的姿态角、。2.4 参数说明(1)位置 ：当地纬度，：当地经度，h：当地高度。(2)姿态角 :为载体的航向角。载体纵轴在水平面上的投影与地理子午线n之间的夹角即为载体航向角。航向角的数值是以地理北向为起点逆时针方向计算的。 :为载体的横滚角(也称为倾斜角)。载体纵向对称面与纵向铅垂平面之间的夹角即为横滚角。横滚角从铅垂平面算起，右倾为正，左倾为负。 :为载体的俯仰角。载体纵轴和纵向水平轴之间的夹角即为俯仰角，向上为正，向下为负。(3)比力 比力是指单位质量受到的位移加速度和重力加速度代数和，即单位质量上所受到外力作用的代数和。 :导航坐标系下的比力:机体坐标系下的比力(4)角速度 角速度用带有上下标的符号表示，如:，其下标含义为b系(机体坐标系)相对于i系(惯性坐标系)的转动角速度，上标含义为此角速度在b系(机体坐标系)中的投影。其它角速度符号含义与此相似。(5)坐标系变换矩阵坐标系变换矩阵也用带有上下标的符号表示，如:c，其含义为n系(导航坐标系)到b系(机体坐标系)的变换矩阵。其它坐标系变换矩阵符号的含义与此相似。(6)地球半径 把地球看作一个椭球体时:地球长轴半径re=6378.393km，椭圆度e=1/298.257。把地球看作一个球体时:地球半径re=6317km。(7)地球自转角速度 =15.0411/h=7.29212 x10-5rad/s(8)重力加速度 在地球表面附近，如果忽略向心加速度的影响，重力加速度的大小与p点到地球球心间距离的平方成反比，可以近似用下式计算得到: g(h)=g0 x ( 1一2h/re)式中:g0为赤道表面上的重力加速度值，g0 =9.8m/s2; h为p点离地面的高度。2.5 gps导航原理 gps导航系统的基本原理是测量出已知位置的卫星到用户接收机之间的距离，然后综合多颗卫星的数据就可知道接收机的具体位置。要达到这一目的，卫星的位置可以根据星载时钟所记录的时间在卫星星历中查出。而用户到卫星的距离则通过记录卫星信号传播到用户所经历的时间，再将其乘以光速得到（由于大气层电离层的干扰，这一距离并不是用户与卫星之间的真实距离，而是伪距（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导航系统卫星部分的作用就是不断地发射导航电文。然而，由于用户接受机使用的时钟与卫星星载时钟不可能总是同步，所以除了用户的三维坐标x、y、z外，还要引进一个t即卫星与接收机之间的时间差作为未知数，然后用4个方程将这4个未知数解出来。所以如果想知道接收机所处的位置，至少要能接收到4个卫星的信号。gps接收机可接收到可用于授时的准确至纳秒级的时间信息；用于预报未来几个月内卫星所处概略位置的预报星历；用于计算定位时所需卫星坐标的广播星历，精度为几米至几十米（各个卫星不同，随时变化）；以及gps系统信息，如卫星状况等。gps接收机对码的量测就可得到卫星到接收机的距离，由于含有接收机卫星钟的误差及大气传播误差，故称为伪距。对0a码测得的伪距称为ua码伪距，精度约为20米左右，对p码测得的伪距称为p码伪距，精度约为2米左右。第 3 章 捷联惯性导航算法及组合导航原理3.1 捷联惯性导航数学编排 捷联惯导系统的数学模型编排方程主要包括导航位置方程和姿态方程。对于捷联式惯导系统来说，需要导航计算机来计算导航平台。对于使用地理坐标系的导航系统而言，首先要用捷联式陀螺仪测量的角速度和导航计算机计算的角速度来计算飞机的姿态矩阵，然后从姿态矩阵的元素中提取飞机的姿态和航向信息，并把捷联式加速度计测量的比力信息用姿态矩阵变换到导航坐标系的比力，这样就可以在导航坐标系中求解飞机的速度和位置信息。捷联惯导系统由于没有机械平台，不能直接测出姿态角和航向角，所以要得到姿态角和航向角，必须在导航计算机中建立“数学平台”，即通过机体坐标系相对地理坐标系的方向余弦阵来计算机体的姿态角和航向角。在对捷联惯导系统工作机理进行全面分析的基础上，可以画出系统的原理方框图，如图3-1所示。图3-1捷联惯性导航系统原理方块图3.2 捷联惯性导航数学模型算法公式 据参考文献1，图3.1中各个模块的计算公式如下: 1.四元数q的即时修正我们选择四元数法作为捷联矩阵即时修正的算法。四元数是指由一个实数单位1和三个虚数单位i, j,k组成并具有下列实元的数。 （3.1）机体坐标系相对平台坐标系的转动可以用四元数q来表示，即 （3.2）q的即时修正可通过解下面的四元数微分方程来实现: （3.3）在求解公式(4.2)时，需要用到四元数的初始值，否则上式无法求解。确定四元数初始值的方法，可以根据初始对准中确定的姿态矩阵初始值中的元素，并利用四元数和姿态矩阵各对应元素相等的关系来确定。2. 捷联矩阵t的计算对于平台惯导系统，由于平台的存在，比力分量可以由加速度计直接测得，但是对于捷联惯导系统，加速度计是沿机体坐标系安装的，它只能测量沿机体坐标系的比力分量，因此引入捷联矩阵来实现机体坐标系到平台坐标系的坐标转换，由于根据捷联矩阵的元素可以确定飞机的姿态角，所以又可以叫做飞机姿态矩阵;由于捷联矩阵起到了平台的作用，所有又可以叫做“数学平台”。 由式(3.3)得出的四元数后，根据式(3.4)即可计算出捷联矩阵。 （3.4）3. 四元数q的最佳归一化由于计算机的算法误差会导致捷联矩阵称为非正交矩阵，对捷联矩阵进行正交化处理可以消除引起非正交的算法误差源的影响。实现四元数的归一化也就是完成了捷联矩阵t的正交化。以欧几里德范数最小为指标的四元数最佳归一化可由下式获得： (3.5)4.比力的坐标转换加速度计测量的比力通过矩阵t可转换为，即： (3.6)5.速度(即)的即时修正惯性导航的基本方程为： 以上可以写成如下矩阵形式： (3.7)6. 位置矩阵的即时修正 位置矩阵c是由地球坐标系转换到平台坐标系的方向余弦矩阵，它是纬度、精度与游动方位角的函数。由于位置矩阵的改变是由平台坐标系相对地球坐标系运动的角速率(又称为位置速率)可所引起的。所以位置矩阵c=可以通过求解下列的矩阵微分方程而获得： （3.8）将位置矩阵c与姿态矩阵(即捷联矩阵)t的即时修正相比较可以看出，由于位置矩阵的变化比姿态矩阵慢得多，所以它要求的即时修正的频率也要慢得多，而且不存在明显的非正交化误差，无需进行正交化处理。采用方向余弦罚可以直接求出位置矩阵c，因此我们采用方向余弦法进行位置矩阵的更新。3.3 捷联惯性导航速度与位置计算毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作 者 签 名： 日 期： 指导教师签名： 日期： 使用授权说明本人完全了解 大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名： 日 期： 学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名： 日期： 年 月 日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权 大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期： 年 月 日导师签名： 日期： 年 月 日致 谢时间飞逝，大学的学习生活很快就要过去，在这四年的学习生活中，收获了很多，而这些成绩的取得是和一直关心帮助我的人分不开的。首先非常感谢学校开设这个课题，为本人日后从事计算机方面