gps软件接收机的仿真与实现毕业论文

毕业论文GPS软件接收机的仿真与实现学院地质工程与测绘学院专业测绘工程摘要随着GPS的升级和新的卫星导航系统的发展，相比较传统GPS接收机，GPS软件接收机具有的成本低、灵活性高等优点越来越突出。它使用软件方法和少量硬件即可实现信号接收处理，可以直接由运行在微处理器上的MATLAB程序完成信号处理，因此具有良好的灵活性、可移植性及可扩展性。因此，研究GPS软件接收机的仿真平台具有重要意义。本文重点对GPS软件接收机的捕获和跟踪部分进行了研究，并在MATLAB中进行了定位解算。本文在掌握GPS软件接收机原理的基础上，实现了对信号的仿真、捕获、跟踪及定位。捕获部分为了提高GPS软件接收机的定位速度和定位精度，选用了在MATLAB环境下执行时间短、性能高的并行码相位搜索捕获算法。跟踪部分将码跟踪环和载波跟踪环组合在一起，降低了跟踪环路的复杂度。载波跟踪环路则选用了对1800相位转换不敏感的COSTAS环，以保证载波跟踪环路对信号的正确跟踪。最后在MATLAB环境下，编写了捕获、跟踪和数据处理等程序，用软件方式实现了对用户的定位，并对定位结果进行了分析概括，验证了所有算法的可行性，讨论了不足之处，为后续软件接收机的相关研究工作奠定了良好的基础。关键词GPS，软件接收机，仿真，捕获，跟踪，同步目录摘要IABSTRACTII第一章绪论111GPS发展概况与组成112GPS的组成1121GPS空间卫星星座部分2122地面控制部分2123用户设备部分313GPS接收机的发展概况314软件接收机的特点及国内外研究现状4141软件接收机的结构特点4142软件接收机的发展现状615课题研究的意义616论文研究的主要内容7第二章GPS信号的产生和结构921GPS信号的产生922GPS信号结构10221载波信号10222C/A码和P码10223导航电文13224GPS卫星信号的调制1423本章小结14第三章GPS信号的捕获1531GPS信号捕获原理1532GPS软件接收机捕获算法15321串行搜索捕获算法15322并行频域搜索捕获算法16323并行码相位搜索捕获算法1733本章小结19第四章GPS信号的跟踪2041解调过程2042锁相环原理2143载波跟踪2244码跟踪2445本章小结24第五章GPS软件接收机的MATLAB实现2651并行码相位搜索捕获算法的MATLAB实现及捕获结果2652GPS信号跟踪的MATLAB实现及跟踪结果2953软件接收机的定位结果3454本章小结35总结与展望37致谢38参考文献39附录41第一章绪论本章简略介绍了GPSGLOBALPOSITIONINGSYSTEM的相关知识，介绍了软件定义无线电技术，阐述了软件接收机国内外的发展现状，最后一节介绍了本文的主要研究内容及研究思路。11GPS发展概况与组成GPS全球定位系统，其全称为定时和测距地导航卫星NANSTARNAVIGATIONSYSTEMTIMINGANDRANGING。是美国从上世纪70年代开始研制，历时21年，耗资200亿美元，于1994年全面建成，具有在海、陆、空全方位实时三维导航与定位能力的新一代卫星导航与定位系统18。它具有全球覆盖、连续全天候工作、高精度和能为高动态平台提供服务等特点，可以提供军用和民用两种服务，并与1995年初达到了全运行能力1。1999年1月，美国政府宣布了在新的GPS卫星信号中增加两个民用信号L2C和L5的提案，这些新增信号将会大大增加民用用户的精度，为全球的民用、商用和科研用户提供帮助，截止目前，新的GPSIII计划又提出了重新评估整个GPS体系结构，以满足直到2030年的民用和军用用户的需求。相对于GPS，目前存在的其他卫星导航系统还包括俄罗斯的全球导航卫星系统、欧盟的伽利略系统和我国的北斗导航系统，它们分别为不同的政府和部门提供服务2。12GPS的组成GPS是由美国国防部为军事目的建立的、陆海空三军联合研制的新一代空间卫星导航定位系统。该系统是以卫星为基础的无线电导航定位系统，具有全能性陆地海洋航空和航天、全球性、全天候连续性和实时性的导航定位和定时的功能，能为各类用户提供精密的三维坐标、速度和时间。GPS定位系统主要由三部分组成空间星座部分、地面监控部分、用户设备三部分。GPS联合计划办公室对这些组成部分的正式计划性术语分别是空间段、控制段和用户设备区段。空间星座部分包括在轨卫星，它们向用户设备提供测距信号和数据电文。地面控制部分更新卫星的时钟校正量、星历，以及其他用于确定用户位置、速度和时间的参数。此外还对卫星进行跟踪和维护，维持卫星的轨道布局，并监测卫星的健康状况和信号的完好性3。用户接收设备部分完成导航、授时等功能。三者各有独立的功能和作用，但又是有机配合、缺一不可的整体系统。图11显示了GPS定位系统的三个组成部分及相互关系。图11GPS系统组成121GPS空间卫星星座部分GPS空间卫星星座的基本配置由24颗分布在6个等间隔轨道平面上的卫星组成。每个轨道平面上有4颗卫星，它们按与地球赤道面成55度的相同方向运行。卫星的这种特定分布保证了在地球表面和近地空间任一点、任何时刻均可同时观测到至少4颗GPS卫星，情况好时最多可观测到N颗。通过测量这些卫星播发的信号到达用户接收机的时间，用户可以用4颗卫星确定出自己的位置。每颗GPS卫星都携带有高稳定度的两台伽原子钟和两台艳原子钟，为用户提供准确定时，并保持与GPS系统时同步。卫星上的S波段接收机用于接收地面注入站的数据和指令，L波段双频发射机用于向GPS用户发射包括卫星星历时钟校正参数、电离层传输延迟校正参数等信息的导航信号4。GPS系统采用码分多址体制将24颗卫星分开，每颗卫星以相同的频率和同样的调制方式使用各自唯一的伪随机码发射信号。伪随机码是一种周期性的可复制的、具有良好自相关特性的二进制伪随机序列。122地面控制部分地面监控部分由一个主控站，5个全球监测站和3个地面控制站组成。每个监测站配有GPS接收机，监测站的主要任务是对每颗卫星进行观测，并向主控站提供观测数据。主控站采集各个监测站传送来的数据，根据采集的数据计算每一颗卫星的星历、时钟校正量、状态参数、大气校正量等，并按一定格式编辑成导航电文传送到注入站。地面控制站也称作地面天线，地面控制站与卫星之间有通信链路。由主控站传来的卫星星历和时钟参数以S波段射频上行注入到各个卫星19。123用户设备部分用户设备部分的核心是GPS接收机，处理由卫星发射的L波段信号。GPS接收机可以捕获到用户当前位置所能观测到的待测卫星信号，并跟踪这些卫星的运行，对所接收到的GPS信号进行变换、放大和处理，以便测量出GPS信号从卫星到接收机天线的传播时间，提取出GPS卫星所发送的导航电文。对获得的导航电文中的广播星历、卫星时钟修正等参数进行处理而实时地计算出用户的三维位置，甚至三维速度和时间。13GPS接收机的发展概况GPS接收机，作为GPS系统中的地面用户设备，是与人们获取GPS信号并应用于自己生活的方式与媒介，所以接收机的发展则是GPS系统发展不可缺少的一部分。从第一台商用GPS接收机面世至今己有近30年。早期的GPS接收机基于模拟信号处理，采用微处理器做应用计算工作。由于使用了大量的模拟器件，早期的GPS接收机体积非常庞大。随着1989年第二阶段卫星投入使用，GPS接收机实现了用微处理器和一个附加的集成电路来完成原来所有由模拟信号处理来完成的工作。随着集成电路技术的成熟，现代GPS接收机采用特定用途的ASIC进行信号处理，并采用高速微处理器进行应用计算，实现了由模拟器件到数字固态器件的转变。这种利用集成电路制造的GPS接收机功能强大而且体积比较小，具有许多优点。上述GPS接收机均采用了大量模拟器件或集成电路来实现，通常称这种采用硬件来实现的接收机为传统GPS接收机5。从目前国内外GPS接收机应用的情况来看，GPS接收机己经嵌入到了我们日常生活的许多用品中，例如蜂窝电话、个人数字助理和汽车等，其应用范围从为计算机网络提供同步的参考时间源涉及到了无人驾驶车辆的导航。同时，对于商用和军用方面，GPS接收机在陆用、航空、海事和空间引导等各个领域的应用越来越广泛，为用户定位、导航和监测发挥着巨大的作用。伽利略计划是欧洲提出的民用全球卫星导航系统研制计划。伽利略计划的最终目标是建立一个独立的，性能优于GPS，与现有全球卫星导航系统具有互用性的民用全球卫星导航系统。该系统计划从2008年开始投入正式运行，2005年12月28日，伽利略计划首颗实验卫星“GIOVEA”顺利发射升空。中国与欧盟于2004年10月9日在北京签署中欧伽利略计划技术合作协议，成为参加该计划的第一个非欧盟成员国6。2003年5月25日，我国在西昌将第三颗“北斗一号”送入太空，与2000年发射的前两颗一起构成了我国完备的卫星导航定位系统，即北斗卫星导航系统（COMPASSNAVIGATIONSATELLITESYSTEM），简称CNSS7。目前改进型的“北斗二号”正在研制中，计划2007年初发射两颗北斗导航卫星，2008年左右满足中国及周边地区用户对卫星导航系统的需求，并进行系统组网和试验，逐步扩展为全球卫星导航系统。可为我国陆地、海洋、空中和空间的各类军事和民用提供多种业务保障，尤其对提高我国国防现代化有着重要的意义。无论是GPS还是计划中的伽利略系统和北斗导航系统，用户部分的核心都是导航接收机，可以预计，随着新一代全球卫星导航定位系统的建立，需要开发适用于新一代导航卫星的接收机。新一代导航接收机与现有的GPS接收机是不兼容的，用户向下一代导航接收机的过渡将面临巨大的升级费用。现代数字GPS接收机多是基于数字信号处理专用芯片和高速微处理器构成的，采用ASIC处理数字中频信号来实现跟踪及解调等工作，其速度较快，但算法固定且不容易改变硬件架构，当算法改变时，必须重新进行芯片设计，代价高昂。软件无线电跟GPS接收机的结合孕育出了GPS软件接收机，采用软件无线电思想的软件接收机，基于相对通用的硬件平台，通过加载不同的运用程序实现不同的功能，可以很好的解决上述问题11。14软件接收机的特点及国内外研究现状141软件接收机的结构特点接收机是卫星导航系统的重要组成部分，是卫星导航系统与用户之间的唯一接口。传统接收机大多都是基于数字信号处理的专用芯片和用于快速计算的微处理器构成的。其主要结构如图12所示。天线接收的射频信号经过射频前端进行放大、混频、滤波、下变频，形成中频信号，然后经模拟/数字转换器得到数字中频信号。此信号经过基带信号处理模块和导航解算模块之后输出定位结果，发送给用户。在硬件接收接中基带信号处理采用专用集成电路（APPLICATIONSPECIFICINTEGRATEDCIRCUIT,ASIC实现信号的捕获和跟踪。天线RF定位结果图12硬件接收机结构图硬件接收机对信号的处理通过ASIC实现，其运算速度快，但算法固定且不容易改变硬件架构。当出现新的发射信号，或者卫星发射的信号结构发生变化时，需要由芯片厂商重新设计并更换接收机内部的专用芯片才可实现，这样不仅代价高昂，更严重影响了接收机的灵活性。另外，硬件接收机对用户来说就像一个黑盒子，用户只能获得最终处理的数据结果，无法控制其内部信号处理过程，一旦出现问题，用户无法查找原因。为了解决硬件接收机的弊端，将软件定义无线电引入到卫星导航接收机的研究中，实现接收机的软件化设计。其全新的设计思想和体系结构，打破了以硬件为核心的传统接收机功能单一、可扩展性差的设计局限13。软件接收机的结构如图13所示天线数字RFIF图13软件接收机结构图与硬件接收机不同的是，软件接收机除了射频前端和A/D转换器用硬件实现外，其他功能全部由软件来实现，包括基带信号处理过程中对卫星信号的捕获与跟踪。当出现新的卫星信号时，只需要修改软件，增加对新信号的处理模块即可，不需要更换或升级硬件。若信号出现失锁成不能正常捕获时，可通过软件调整不同的捕获和跟踪参数等方法来解决，方便寻找和分析原因。利用软件接收机通过多个导航定位系统的数据组合可以加快定位速度，提高定位精度。另外，通过集成软件接收机与相关的测试分析工具构建一个完整的平台。在此平台上，可以方便的测试和优化接收机的结构、功能、性能，研究信号特性，探讨信号处理方法，设计和验证接收机各组成单元一些关键参数的选择等。总之，软件接收机与硬件接收机相比，成本低，易于实现，具有更高的灵活性、通用性和开放性等优点，成为卫星导航接收机新的发展方向。软件部分硬件部分IFRF前端基带信号处理导航解算用户接口A/D转换软件部分硬件部分IFRF前端A/D转换基带信号处理用户接口导航解算142软件接收机的发展现状近年来越来越多的学者开始关注GPS软件接收机，软件接收机研究成为导航领域的一个热点。受元器件发展缓慢的影响，早期的软件接收机的研究大多只实现了跟踪和导航解算软件化，而基带信号的相关处理仍需要采用硬件实现。随着IC技术的发展，软件接收机的研究越来越趋向于真正的软件化。射频前端输出数字中频信号后，信号的捕获、跟踪和导航数据的解调与解算全部通过软件实现。2001年，美国DATAFUSIONCORPORATION研发了一套基于MATLAB/C的完整的GPS单频软件接收机开发测试平台。2001年，美国斯坦福大学与瑞典吕勒奥科技大学共同合作开发了四通道的GPS实时软件接收机，之后又合作研究了直接采样的射频前端，并设计了GPS/GALILEO多频接收机117。2003年，美国CENTERFORREMOTESENSINGINC开发出一套开放式的GPS软件接收机开发和测试平台，通过图形用户界面可以对该平台各个软件模块进行连接和配置。2003年，美国CORNELL大学开发出一套12通道单频GPS实时软件接收机，在普通PC上实现了12通道的实时处理。2004年，丹麦ALBORG大学GPS中心研制出基于MATLAB的GPS软件接收机测试样机，并在2005年的美国导航学会上提出了单频GPS和GALILEO软件接收机的研制方案2。2006年，瑞典的NORDNAV公司开发山GPS软件接收机R30。在该软件接收机中RF前端输出的数字中频信号通过USB20接口输入至PC机进行处理，最多可并行处理24个通道的信号，每个信号通道的捕获门限、相关器数量及载波相位和码相位的跟踪环路参数都可独立设置，而且接收机还支持实时和后处理两种工作方式，在PC上可并行实时处理12通道的数据1。近几年我国对接收机的研究和开发也有了很大进步。北京航空航天大学的金天等人实现了基于PENTIUM/AMD的处理器的12通道GPS实时软件接收机。中国科学院光电研究院的廖炳瑜等人实现了基于FPGAXILINXXC2V3000和DSPTMS320C6416的星载双频GPS软件接收机2。此外，越来越多的高校和企业也开始了对软件接收机进行研究，这将促进我国导航技术的发展。15课题研究的意义GPS软件接收机较之传统的GPS接收机是一个巨大的进步，它具有以下许多优点。1实用性使用软件方法和少量硬件进行GPS信号接收处理。GPS软件接收机成本低，开发速度快，具有高度的可配置性。它只需要进行软件修改，就可以对接收机进行修改和升级。当出现新的GPS信号频率或使用新的伪随机码序列时，只需对接收软件进行更改就可以适应新的信号。2灵活性GPS软件接收机也可以在进行修改后用于接收处理新的GPS信号和其他卫星定位系统的信号。可以根据环境使用不同的门限值或者调整和修改参数，保证接收机在任何情况下都能快速准确地实现定位。随着美国政府新的GPS计划的提出和其他新的导航系统的出现，GPS软件接收机的灵活性越来越显著。3高性价比适应于各种卫星导航系统；软件模块可重复使用；可供教育、研究机构研究测试新的软件算法，评估接收机性能；还可以用来实现比较复杂的算法。4强适应性若新的卫星导航系统出现后，只要对软件稍做调整，便可以适应不用的导航系统，而不需要另外购买新的接收机。5模块化设计部分软件模块可以进行复用，从而可以降低开发的成本8。基于以上优点，GPS软件接收机的实现必将为其带来广阔的应用前景。较之传统的GPS接收机，软件接收机能快速适应现代化的GPS信号和新的卫星导航系统，它的出现和壮大将为GPS定位技术的进一步发展注入新的活力，是卫星导航系统现代化的重要保障。另外，我国正在建设北斗卫星导航系统，GPS软件接收机修改后便可以应用于北斗定位算法的验证和开发。因此，深入研究GPS软件接收机对我国的卫星导航事业也具有促进作用。16论文研究的主要内容本文的研究的主要内容为GPS软件接收机的仿真与实现。在课题内容中重点研究了GPS软件接受机的算法，并用编程语言进行了实现，以及GPS软件接收机的具体设计与实现。论文主要内容分为6个部分，其中各章内容分别为第1章绪论主要介绍了GPS的发展概况和结构组成、GPS接收机的发展概况，阐述了课题的研究意义，给出了本文的主要内容和组织结构。第2章GPS信号的产生和结构本章详细分析了GPS信号的产生过程和GPS信号包含的具体内容和结构。第3章GPS信号的捕获讨论了时域滑动相关算法、基于快速傅里叶变换（FFT）的并行频域搜索和并行码相位搜索捕获算法，对其做了性能对比分析。第4章GPS信号的跟踪介绍了锁相环原理。码跟踪采用非相干延迟锁定环（DLL），载波跟踪采用经典的COSTAS锁相环。第5章GPS软件接收机的MATLAB实现介绍了用编程语言实现GPS软件接收机的算法，并对各个部分进行了成果展示与分析。结论与展望本章对本文的工作做了总结，并指出下一步的研究方向。第二章GPS信号的产生和结构GPS卫星向广大用户发送的用于导航定位的信号，是一种调制波，但有别常用的无线电广播电台发送的调频调幅信号，它是利用伪随机码传送导航电文的调相信号。GPS卫星信号是目前常用的两种卫星导航定位信号之一，它包含有载波信号（L1和L2、测距码（C/A码和P码）和数据码（D码，也称导航电文）。21GPS信号的产生GPS卫星播发的信号在UHF（ULTRAHIGHFREQUENCY）频段的两个频率部分L1和L2上传播（UHF频段覆盖了从500MHZ到3GHZ的所有频率）。下图21描述了GPS信号的产生原理图21GPS信号的产生在图21的最左边，时钟信号F01023MHZ分别和154和120相乘产生LL和L2载波信号。同时时钟信号F0在左下角限幅器作用后输入到PY码生成器和C/A码生成器。在图21的最底部，数据生成器产生导航数据。码生成器和数据生成器通过PY码生成器产生的信号X1同步9。在码生成器右边，PY码生成器和C/A码生成器分别产生的PY码和C/A码分别通过模2加法器和导航数据组合。产生的PY码导航数据和C/A导航数据输入到LL频率支路的两个调制器，通过二进制相移键控BPSK方法调制到载波LL上。PY部分信号衰减3DB后，PY支路和C/A支路相加产生LL信号。另外一路PY码导航数据输入到L2频率支路的调制器，通过二进制相移键控方法调制到载波L2上，这部分信号衰减6DB后产生L2信号。22GPS信号结构221载波信号GPS卫星测距码信号和导航电文的传输速率都很低，C/A码、P码和导航电文的传输速率分别为L023MBIT/S、1023MBIT/S和50BIT/S。且GPS卫星离地面约20200KM，其电能非常紧张，因此很难将测距码和导航电文速率很低的信号传输到地面。解决这一难题的方法就是将测距码信号和导航电文调制到L频段的两种载波频率FL1和FL2上发射，其中FL1154F0157542MHZ，波长11903CM，其上调制C/A码、P码和导航电文；FL2120F0122760MHZ，波长22442CM，其上仅调制C/A码或者P码其中一种，或者都不调制。222C/A码和P码每个GPS卫星发送的测距码有两种精密测距码PY码和粗捕获码C/A码，两种均为伪随机码。伪随机码全称伪随机噪声码PRNPSEUDORANDOMNOISECODE，简称伪码，是一种可以预先确定又可重复产生和复制，具有类似于白噪声随机统计特性的二进制码序列。伪码的产生方式很多，GPS系统中采用的是由M序列即最长线性反馈移位寄存器序列产生的复合码。M序列是由多级反馈移位寄存器产生的，不同级数的反馈寄存器以及不同的反馈抽头都将产生不同的M序列。1C/A码的生成原理GPS信号中的C/A码是序列长度为1023位码片数的GOLD码。GOLD码是M序列的复合码，由两个码长相等、码时钟速率相同的M序列优选对模2和构成。改变产生它的两个M序列的相对相位，就能得到一个新的GOLD序列。GOLD码有较优良的自相关和互相关特性，构造简单，产生的序列数多，广泛的应用于码分多址系统10。本GPS软件接收机使用带有C/A码的L1频段的GPS信号进行定位。GPS信号中C/A码的码速率为1023MHZ，周期为1MS。其中每个基码码片的宽度为1MS/10239775NS。C/A码的产生原理如下图22所示图22C/A码的产生原理根据C/A码的产生原理可知，G1的输出为第10位线性移位寄存器的输出即MLS输出，G2的输出则依赖于相位选择器，相位选择器选择两个位的输出经过一个模2加法器产生的码序列为G2的输出。C/A码发生器将G1和G2的输出作为输入，在用一个模2加法器产生C/A码。每到1023个周期，移位寄存器的值会重置为“1”，开始产生下一周期的码。由图25可知，G1和G2的驱动时钟为1023MHZ。G1寄存器的生成多项式为211031XXG由式21可知，寄存器的第3位和第10位的状态反馈给输入。同样，G2寄存器的生成多项式为2210986322XXXGPS星座中总共有37颗卫星，因此要产生37个不同的PRN码序列。如何产生不同的C/A码主要取决于G2序列的相位选择器。表21列出了不同卫星的对应的不同相位选择器的选择。表21码相位分配卫星ID编号GPSPRN信号编号码相位选择码延迟码片C/A码前10个码片11265144022376162033487171044598174455191711336621018145577181391131882914014549931014116231010232511504111134152164212125625417501313672551764141478256177215158925717751616910258177617171446911561818254701467191936471163320204747217152121584731746222269474176323231350910632424465121706252557513174326266851417612727795151770282881051617742929168591127303027860145331313886116253232498621712表中共有5列，第1列是卫星的ID（IDENTITY）号，第2列是PRN信号号码，第3列是码相位选择，用于形成G2产生器的输出，第4列是以码元数表示的码延迟，这个延迟就是MLS和G2输出之差。一旦码相位定了，延迟数也就定了。最后一列为每颗卫星的C/A码的前10位的八进制表，用于检查产生的C/A码是否有误11。C/A码之所以用作GPS信号传输，是因为它有以下良好的自相关和互相关特性。1互相关特性对于卫星I和卫星K对应的C/A码CIK和CKK，其互相关函数如下任意M230102MLMRLIIK2自相关特性对于卫星K对应的C/A码CKK，其自相关函数如下240102MLMRKLK125123102LKLK由上述特性可以看出，所有的C/A码之间是几乎互不相关的。除了完全对齐的情况，所有的C/A码自身也是几乎不相关的。这种特性使得很容易找出两个完全对齐的相同的C/A码。2P码的生成原理P码的码速率为1023MHZ，每个基码码片约宽977NS，对应的空间距离为293M。由于其基码码片宽度为C/A码的1/10，所以其测距精度远远高于C/A码，故称P码为精密测距码。P码的产生原理与C/A码相似，由两个码长互素的子码X1和X2组成的模2和复码产生。两个子码X1和X2是由24级移位寄存器产生的截短码，X1和X2的码速率均为1023MHZ。不同的卫星分别对应于不同的P码。223导航电文导航电文是用户用来定位和导航的数据基础，主要包括卫星星历、时钟改正、电离层时延改正、工作状态信息已经C/A码转到捕获P码的信息，这些数据以二进制的形式发送给用户，因此导航电文又称为数据码，或称之为D码。图23导航电文的组成格式导航电文的格式如图23所示，它的基本单位是1500BITS的一个主帧，导航电文的传输速率为50BPS，因此，每30S才能够传送完一个主帧。每一个主帧包含5个子帧，前三个子帧各有10个字码，每个字码30BITS；第4、5子帧各有25个页面，共有37500BITS，长达125MIN。第4、5子帧与第1、2、3子帧不同，它需要750S才能传送完毕全部信息量。224GPS卫星信号的调制传统的GPS信号属于直接序列扩频信号，信号的调制过程只是周期性地改变载波的相位，而载波的频率没有变化。各模块都是在以频率F01023MHZ为基准时钟的基础上产生的，基准频率F0经过限幅器为C/A码和P码生成器提供了稳定的时钟，导航数据和C/A码是由P码发生器生成的X1历元控制同步和叠加而成。叠加后的C/A码和P码再通过调制器以BPSK相互正交的调制方式调制到载波L1上，即分别调制在载波L1的同相和正交通道上。P码部分的信号衰减3DB后，与C/A码部分信号相加就是L1频段上发射的信号。L2频段上的GPS信号的结构与L1频段上信号基本一致，但是它通过切换开关仅提供P码与导航数据异或、C/A码与导航数据异或和单独P码中的一种12。由此，GPS卫星发射信号的结构分别为L1载波上26COSSIN1011TTDPATTDCATSLLLL2载波上272CO2PB图24信号调制结果图23本章小结本章介绍了有载波信号、测距码、D码组成的GPS信号结构，并详细介绍了C/A码的生成原理及GPS卫星信号的调制，为后续捕获和跟踪技术奠定了基础。第三章GPS信号的捕获在传统GPS接收机中，捕获部分通常是用专用的集成电路ASIC来实现的。而在GPS软件接收机中，捕获部分是通过软件来实现的。本章首先介绍了几种基本的捕获算法。随后给出了一种能在软件接收机中有效实现捕获的算法，并在MATLAB中利用算法对一组实际采集的信号进行了捕获。31GPS信号捕获原理在软件接收机中，捕获是利用伪随机码良好的相关特性，检测本地C从码和卫星信号伪随机码的相关输出。对某一颗可见卫星，当本地C/A码的码相位、本地载波的频率和输入信号中的码相位以及载波频率相匹配时，有最大的相关值。本地C/A码的码相位与输入信号中的码相位在任何一边的偏移超过一个码元时，有最小的相关值。当检测到最大相关值时便达到了捕获的目的。通过对GPS信号的捕获，可以确定犯颗卫星中的可见卫星，并可计算出可见卫星的载波频率和OA码的码相位的粗略值。接收到的GPS信号通常是几颗可见卫星播发的GPS信号的组合，如下式3121TSTSTSN上式中表示卫星播发的GPS信号，是接收机接收到的所有颗可见TSNNN卫星播发的GPS信号的组合。当捕获卫星时，输入信号和本地所产生的不同KTC/A码序列、不同本地载波序列相乘。利用C/A码良好的相关特性，通过检测到的最大相关值，找到可见卫星和对应卫星的载波频率和C/A码的码相位的粗略值，达到捕获的目的。其中信号捕获输出的C/A码的码相位和载波频率是用于跟踪的初始条件13。32GPS软件接收机捕获算法对于GPS卫星信号有多种捕获方案，以下几节将分别介绍三种捕获算法串行搜索捕获算法，并行频域搜索捕获算法和并行码相位搜索捕获算法。321串行搜索捕获算法串行搜索捕获算法在码分多址CDMACODEDIVISIONMULTIPLEACCESS系统中是一种常用的方法。GPS就是一种CDMA系统。下图31给出了串行搜索捕获算法的基本原理框图。图31串行搜索捕获算法原理图由图31可以看出，串行搜索捕获算法基于本地产生的C/A码序列、本地振荡器产生的载波与输入信号之间的相乘运算。对于一个特定的卫星，C/A码产生器产生一个本地C/A码序列。这个本地C/A码序列的码相位变化范围是01022个码元。输入信号首先和本地产生的C/A码序列相乘，然后再和本地振荡器产生的载波相乘。其中和本地振荡器产生的载波相乘后产生支路，和本地载波偏移900相位后相乘产生I正交的支路。然后支路和支路分别经过积分运算和平方运算后相加作为输出。QIQ理论上，由于C/A码只在信号上调制，信号的能量应该在支路。但是由于不II知道接收到的GPS信号的相位，解调后的信号不一定等于卫星产生的信号。所以I为了保证能够检测到信号，对支路和支路分别做积分运算是必要的。最后的输出信号是输入信号和本地产生信号的相关值。如果输出大于预先给定的门限值，那么本地产生的载波频率和C/A码的相位就是正确的。然后将该载波频率值和C/A码的相位值传递给跟踪部分，作为跟踪部分的初始值。串行搜索捕获算法中有两次不同的扫描运算过程。一次是对中频IF10KHZ范围内所有频率以一定步长的频率搜索，本论文采用500HZ为步长；一次是对1023个不同的码相位搜索。这两次的搜索过程总共有41943次运算。如下式32419302152103频率码相位显然，41943次运算的运算量非常大。这种耗时的搜索是串行搜索捕获方法的主要弱点。但是它同时也具有算法简单，实现起来简单直接的优点14。322并行频域搜索捕获算法如果能从搜索过程中消除对其中一个参数的搜索，或者能使对两个参数的搜索同时进行，算法性能将会大大提高。这一小节介绍的并行频域搜索捕获算法，就是对一个参数的搜索。这种方法利用了傅立叶变换将信号从时域变换到频域。并行频域搜索捕获算法的原理图如下图32所示图32并行频域搜索捕获算法原理图图32中，输入信号首先和本地产生的C/A码序列相乘。对于一个特定的卫星，C/A码的相位介于01022个码元之间。然后对相乘后的信号进行傅里叶变换，使其变换到频域。进行傅里叶变换时，可以通过离散傅里叶变换DFT或快速傅里叶变换FFT实现。其中FFT比DFT的运算速度快，但是FFT要求输入序列的长度等于2的正整数次幂15。323并行码相位搜索捕获算法相比较以500HZ为步长的载波频率搜索总量41，码相位的搜索总量1023是比较大的。上节介绍的并行频域搜索捕获方法中，通过将信号在频域的搜索并行化方式，消除了对41个可能的载波频率进行必要的搜索过程。如果能使捕获过程在相位域并行化，那么相比较并行频域搜索捕获方法中的1023次搜索，这种情况下只需要搜索41次。下面介绍的并行码相位捕获算法就是用到了以上提到的并行化码相位搜索过程的办法。捕获的目的就是为了对输入信号和一个C/A码序列做相关运算。相比较串行搜索捕获方法中输入信号分别和1023个不同码相位的C/A码序列相乘，对输入信号和具有不同码相位偏移的C/A码序列做循环互相关更方便。下面是一种通过傅立叶变换实现循环相关的方法。两个长度同为N的有限长序列和的离散傅立叶变换和计算如NXYKXY下3610/2NNNKNJEKX37/KJYY两个长度同为N的有限长序列和的循环互相关计算如下XNZ381010NMNMYXNYNZ下面的分析中均省略了中的缩放因子，的N点离散傅立叶变换Z计算如下KZNKNJNMNEYXKZ/21010/210/2NNMKJNMKJE39YX式中是的复共轭。KX由以上分析可以看出，两个序列和的循环互相关的傅立叶变换XYNZ可以用这两个序列的频域表达式的复共扼和表示。然后再对KKXY取反傅立叶变换，就可以得到序列和的循环互相关的表达式KZN3101ZFZ下图34是并行码相位搜索捕获算法的原理图图34并行码相位搜索捕获算法原理图图34中，输入信号首先和不同相位的两路本地载波相乘，和本地载波相乘后的产生支路信号，和本地载波偏移900相位后的信号相乘产生支路信号。然后支IQI路和支路作为一路复输入信号进行傅立叶变换。本地C/A码经QNJINX过傅立叶变换后取共轭，然后其共轭和输入信号经过傅立叶变换得到的信号相X乘。通过反傅立叶变换将这部分相乘后的信号转换到时域，最后通过取绝对值得到输入信号和C/A码之间的相关值。绝对值的大小代表了输入信号和C/A码相关性的强弱，通过找出结果中的一个最大值，确定出输入信号和C/A码相关性最好时C/A码的各种状态信息。相比较前面的捕获方法，并行码相位搜索捕获方法将搜索量减少到了对41个不同载波频率的搜索。每次捕获时，C/A码的傅立叶变换只需要计算一次，41个不同的载波频率分别要进行傅立叶变换和反傅立叶变换。所以并行码相位搜索捕获算法的计算效率主要依赖于上述过程的实现。并行码相位搜索捕获算法参数估计的频率精度和串行搜索捕获算法相同。对于C/A码相位，由于并行码相位搜索捕获算法对每个抽样的码相位都给出了相关值，其精度大于以上两种方法。因为如果抽样频率是10MHZ，一个抽样后的1MS的C/A码序列有10000个点，所以对于一个码相位的计算结果有10000个值。相比较前面两种方法的1023个值，其码相位的精度是明显增大的。和其他捕获方法一样，并行码相位搜索捕获算法可以直接通过方框图34实现。相比较前面的两种捕获算法的方框图，并行码相位搜索捕获算法的原理图中几乎没有添加新的模块，这样许多元件就可以在这种方法的实现中复用16。33本章小结本章介绍了信号捕获的基本原理和常用的几种基本的捕获方法，并对捕获中的有关参数进行了分析。随后给出了在本文研究的软件接收机中有效实现捕获的算法，并在MATLAB中利用该算法对一组实际采集的信号成功进行了捕获。31节详细介绍了进行GPS信号捕获的基本原理。32节介绍了三种捕获算法，其中串行搜索捕获算法基于本地产生的C/A码序列、本地振荡器产生的载波与输入信号之间的相乘运算。串行搜索捕获算法中有两次不同的扫描运算过程。一次是对中频IF10KHZ范围内所有频率以一定步长的频率搜索，本论文采用500HZ为步长，一次是对1023个不同的码相位搜索。这种耗时的搜索是串行搜索捕获方法的主要弱点，但是它同时也具有算法简单，实现起来简单直接的优点。并行频域搜索捕获算法利用了傅立叶变换将信号从时域变换到频域，通过将信号在频域的搜索并行化，消除了对41个可能的载波频率的搜索过程，只对1023个不同的码相位进行了搜索，减少了搜索量，但计算效率没有明显提高。并行频域搜索捕获算法中载波频率的精度是1HZ，比串行搜索捕获算法500HZ的精度高。并行码相位捕获算法通过傅立叶变换实现循环相关，使得捕获过程在相位域并行化，将搜索量减少到了对41个不同载波频率的搜索，运算效率较高，适合软件接收机设计。这种捕获算法对每个抽样的码相位都给出了相关值，其精度明显大于串行搜索捕获算法和并行频域搜索捕获算法，其频率精度和串行搜索捕获算法相同。第四章GPS信号的跟踪完成信号捕获以后，得到的载波频率和C/A码相位只是对GPS信号的粗略估计。跟踪部分的主要目的就是得到更精确的载波频率和C/A码相位，保持对信号的跟踪和解调出特定卫星的导航数据。下图41是一个基本的解调原理图。图41解调原理图图41给出了从输入信号中解调出导航信息的基本原理。首先，输入信号和一个载波信号相乘，相乘后可以消除输入信号中的载波。然后将消除载波后的信号和C/A码相乘，相乘后输出导航信息。由以上过程可以看出，跟踪模块必须产生载波信号和C/A码信号来完成准确的跟踪和卫星信号的解调。下面是对GPS信号解调过程的详细描述。41解调过程设和分别是卫星播发的L1载波频率和L2载波频率，是卫星对应的1LF2KKCC/A码序列，是卫星对应的PY码序列，是导航数据序列，、，和KPKDP1L是C/A部分信号和PY码部分信号的功率。卫星播发的GPS信号可以表2PLTSK示为2SIN2COS111TFTPTFTDCTSLKLLKK41I22TFTDK在经过含有滤波部分和下变频部分的前端接收设备后，接收到的如下TSK42INCOS21TPTTPTSIFKLIFKCK其中输入信号中的载波经过前端接收设备后下变频成中频。IF该信号然后经过A/D转换被采样。由于在经过窄带滤波器时，P码失真变形。所以在经过A/D转换后的信号中以噪声表示。式（42）中的信号经过A/D转换后NE如下43COSNEDCSIFKK其中表示该信号在时域是离散的，大小为的正整数倍。NF1要从经过A/D转换后的信号中提取出导航数据，必须将该信号转换到基带。NDK图41中所示，首先通过输入信号和本地载波相乘去除输入信号中的载波。如果本地载波的频率和相位和输入信号中载波的频率和相位相同，它们相乘后如下COSCOSNNCNIFIFKIFK442121DCDKIK上式中第一项是导航信息和C/A码相乘，第二项中含有频率为2倍中频的载波信号。式（44）中相乘后得到的信号通过一个低通滤波器后第二项被滤除，得到信号4521NCK然后通过将去除载波后的信号和本地C/A码做相关运算去除信号中的码。如果NCK本地C/A码和输入信号的中的码序列和码相位完全相同，相关运算后的输出如下4610NNDNKNNKK其中是导航信息和信号中点的总数的乘积。NNDK以上只是对一颗卫星的一个信号的解调分析，这样减少了公式的复杂度，给出了解调过程的一种简单思路。在实际信号中，接收到的信号是由许多可见卫星播发的信号组成的。图41中，解调时需要产生本地载波和本地C/A码两个信号。要产生这两个完全准确的信号，就必须对其加一些反馈。通常称产生本地载波的反馈环路为载波跟踪环，产生本地C/A码的反馈环路为码跟踪环。42锁相环原理锁相环是一个相位负反馈误差闭环控制系统，主要用于调节本地晶体振荡器的输出频率来和输入信号的频率相匹配的环路，由鉴相器、环路滤波器和压控振荡器三部分组成。鉴相器是一个误差检测单元，比较输入相位与输出相位后得到误差信号。鉴相器特性多种多样，不同的锁相环的区别主要在于鉴相器的选取。环路滤波器实质上低通滤波器，它滤除误差信号的高频分量，输出一个信号去控制压控振荡器的频率。环路滤波器的阶数既决定了锁相环的性能又决定环路的阶数，若选用N阶滤波器，则其锁相环路为N1阶。压控振荡器是一个线性调频振荡器，其振荡频率在工作区域内随控制信号作线性变化。由环路滤波器输出的控制信号控制振荡器的输出频率，使之与输入频率同步。基本锁相环如图42所示，实际应用中的各种形式的环路都是由基本环路演变而来的。AB图42基本锁相环A时域结构BS域结构图42A是锁相环的时域结构，图，42B是经过拉普拉斯变换得到的S域结构。由于相位与频率之间存在积分的关系，为了便于分析，设输入相位为，压控振TI荡器的输出相位为，相位比较器测量两者的相位差，放大器为相位比较器的TF0K增益，低通滤波器滤除环路中的噪声。43载波跟踪载波跟踪环的主要作用是对码跟踪环的输出信号进行解调，得到导航电文数据，同时得到载波DOPPLER频移观测量，用于对接收机进行高精度测速。要成功的解调出导航数据，必须产生完全准确的本地载波信号。通常采用相位锁定环PLLPHASELOCKEDLOOP和频率锁定环FLLFREQUENCYLOCKEDLOOP来跟踪一个载波信号。PLL是通过测量本地载波相位和输入载波之间的相位差来保持对输入信号相位的跟踪。FLL是通过测量本地载波频率和输入载波之间的频率差来保持对输入信号频率的跟踪。本文采用PLL实现载波跟踪，其基本原理图如下图43所示图43PLL实现载波跟踪原理图图43中，前两个乘法器用于消除输入信号中的载波和C/A码，其中C/A码的消除需要用到码跟踪延迟锁相环中的输出。环路鉴相器用于找出跟踪时关于本地载波PI的相位差。载波环路鉴相器输出的相位差或者相位差的函数先通过环路滤波器，然后滤波后的信号作为反馈控制载波发生器来调整本地载波的频率大小。经过以上过程后，载波发生器会产生与输入信号中的载波相同的本地载波。图43是一个一般的相位锁定环，这种锁定环对1800的相位转换比较敏感。GPS信号中的导航数据位在经过一个20MS的周期后存在相位跳变，所以GPS接收机中需要的PLL应该对1800的相位转换不敏感。本文中研究的GPS软件接收机采用了COSTAS环。COSTAS环很重要的一个特性就是对1800的相位转换不敏感，因此它对导航数据位的跳变也会不敏感。GPS软件接收机中常用它作为载波跟踪环路。下图44是COSTAS环的基本原理图图44COSTAS环原理图图44中COSTAS环内包括两个乘法器。一个是本地载波和与C/A码相乘后的输入信号相乘形成支路，一个是经过900相移的本地载波和与C/A码相乘后的输入信号相I乘形成支路。COSTAS环的主要目的是将能量保持在支路，这点是通过对载波发生QI器加一些反馈来实现的。载波环路鉴别器用于找到反馈给载波相位振荡器的反馈因子。本文在研究GPS软件接收机时，选用本地载波和输入信号的相位差作为反馈项，该相位差计算如下436TANCOS21INDIQKK437TAKI由式（437）可以看出，当支路为零I支路取最大值时，本地载波和输入信号的相位差最小。这种反正切鉴别器是COSTAS环常用的鉴别器中最精确的，但也相对是最耗时的。下表41介绍了几种常用的鉴别器表41几种常用的鉴别器鉴别器特征鉴别器的输出与成比例KQISIGNDSIN鉴别器的输出与成比例K2鉴别器的输出相位差KI1TA当输入信号和本地载波的相位差为0或是时，表41中所示的这三种鉴相器018的输出均为零，所以COSTAS环也就对导航数据位跳变引起的1800的相位转换不敏感，也就是在导航数据位发生跳变的时候，载波跟踪环依然可以保持对信号的跟踪。44码跟踪伪码跟踪目的是进一步提高本地伪码与信号伪码的相关程度，更好地完成对扩频码的解扩，从而提高伪距精度，得到准确的导航数据。在捕获阶段，本地C/A码的码元误差被限制在半个采样间隔之内，在跟踪阶段，要将这个误差进一步地缩小并持续地保持在误差允许的范围以内。常用的接收扩频信号有两种方式1相干接收在本地产生一个相干载波，用相干载波调制本地伪码，然后和接收信号相乘，实现载波分离及信号解扩；2非相干接收用一个超前滞后延迟锁定环直接从GPS信号中提取码延迟误差，完成对码的跟踪。相干接收需要复现同步的本地载波和伪码，很难实现。因此，本文选用非相干超前滞后结构形式的延迟锁定环DLL作为码跟踪环。其结构如图45所示。图45码跟踪（DLL）原理图45本章小结利用捕获完成后得到的载波频率和C/A码相位的粗略值，跟踪部分可以得到更精确的载波频率和C/A码相位，保持对信号的跟踪和解调出特定卫星的导航数据。本章主要对GPS信号跟踪部分的载波跟踪部分和码跟踪部分做了深入研究。41节分析了对接收到的卫星信号进行解调的基本原理。输入信号和完全准确的本地载波和本地C/A码信号相乘后输出导航信息。其中产生本地载波的反馈环路为载波跟踪环，产生本地C/A码的反馈环路为码跟踪环。由于载波跟踪部分和码跟踪部分都有线性锁相环部分，在GPS信号跟踪部分，通常采用跟踪