PN码产生器的理论研究及MATLAB仿真毕业设计论文(1)

1、本科毕业设计论文本科毕业设计论文题题 目目 PN 码产生器的理论研究及 MATLAB 仿真实现 毕业 任务书一、题目 PN 码产生器研究及其 MATLAB 仿真二、研究主要内容1、 了解 PN 码的根本定义：PN 码Pseudo-Noise Code ，是一具有与白噪声类似的自相关性质的 0 和 1 所构成的编码序列，最广为人知的二位元 P-N Code 是最大长度位移暂存器序列，简称 m-序列， 他具有长 2 的 N 次方 - 1 个位元，m 级暂存器来产生。同时 PN 码分长码与短码，在 CDMA 中的担当不同的角色。2、主要研究内容：了解 GPS 全球卫星导航定位系统及扩频通信的相关概念

2、，重点是研究 GPS 信号结构，并设计 GPS 信号的 m 序列和 C/A 码及 P 码的生成算法，对 GPS 的信号结构及特性进行深入了解并明确仿真要到达的目标，再通过对直接序列扩频通信等理论的研究及分析来确定 GPS 信号的生成算法9，在此过程中，要掌握 PN 码相关特性，最终实现 GPS 信号的产生，最后在理论证明成立的情况下利用 MATLAB 仿真验证他们的随机特性.具体以 m 序列发生器为例。 三、主要技术指标了解 GPS 卫星的相关概念及原理，分析 GPS 信号的特性及其产生机制，并掌握扩频通信中伪随机序列的概念和应用；掌握移位存放器产生m序列的根本知识，并熟悉伪随机序列的特性；掌

3、握 C/A 码发生器的原理，以及 C/A 码的自相关性及其互相关性，并进行性能仿真分析；了解 P 码在 GPS 系统中的重要地位，掌握俩个子码 X1，X2 的产生原理，以及通过 X1，X2 复合构成的 P 码的原理。4、进度和要求 第 15 周 查阅相关资料文献，完成开题报告。设计论文 第 610 周 进一步解析m序列、C/A 码、P 码发生器的原理，完成方案设计并用 MATLAB 进行仿真验证；第 1116 周 完成课程设计的全部内容，按要求编写课程设计毕业论文， 制作电子演示文稿。五、主要参考书及参考资料1 刘基余，GPS 卫星导航定位原理及方法，第二版，北京：科学出版社，2021.2 徐

4、明远，邵玉斌.MATLAB 仿真在通信与电子信息工程中的应用. 西安电子科技大学出版社，2005.3，2021.4 JAMES BAO-YEN TSOUI 著 陈军 潘顶峰等 译.GPS 软件接收机根底.电子工业出版社.2007 .5 明浩GPS 系统中的 P 码产生与特性分析J沈阳理工大学学报，2005.6 王惠南GPS 导航原理与应用M北京：科学出版社，20037 Spaans E.J.The Munich Satellite Navigation Summit 2007J. European Journalof Navigation , 2007 .8 杨俊，武奇生编著. GPS 根本原

5、理及其 MATLAB 仿真. 西安电子科技大学出版社，2006.9MATLAB/Simulink，2021.10 JAMES BAO-YEN TSOUI 著 陈军 潘顶峰等 译.GPS，2007.11 樊昌信通信原理M第 5 版北京：国防工业出版社，200112.13 J.B.Lozow. Analysis of Direct P(Y)-Code Acquisition, Navigation, Spring 199714 A1ison BrownModeling and Simulation of GPS Using Software Signal Generation and Digita

6、l Signal Reconstruction, Proceedings of the ION National Technical Meeting, January 200015 Charles CCounselmanMultipathRejecting GPS Antennas，Proceedings of the IEEE，V01.87，N O.1，January 1999.16 李天文 GPS 原理及应用. 科学出版社，2003.17 陈涛.GPS接收机基带信号处理的研究和设计D.上海：上海交通大学，2021.18MatLab/Simulink的GPS系统仿真J.系统仿真学报，2006

7、.19 谢刚 . GPS 原理与接收机设计. 电子工业出版社 .2021.20 赵鸿图 矛艳 通信原理 MATLAB 仿真教程. 人民邮电出版社.2021.摘 要 GPSGlobal Positioning System是基于卫星的定位系统，被认定为当前定位导航设备中最重要的开展方向，在军事和民用中有着广泛的应用，并且在日常生活方面是一个难以用数字预测的广阔领域。GPS 像移动 ， 机，计算机互联网对我们生活的影响一样，人们的日常生活也已离不开它。GPS 系统采用典型的 CDMA 体制，这种扩频调制信号具有低截获概率特性1。该系统主要利用直接序列扩频调制技术，采用的伪码有 C/A 码和 P(Y

8、)码两种。目前，伪随机码在扩频通信中的应用非常广泛，这一技术也在卫星导航系统中得到了利用。在 GPS 中，导航卫星向用户发送的导航定位信号，就是利用伪随机码传送导航电文的调相信号。本文首先介绍了 GPS 卫星系统的组成，简要论述 M 序列和伪随机噪声码P 码和 C/A 码及其产生，重点研究 C/A 码的产生原理并用MATLAB 软件对 C/A 码的生成进行了仿真。关键词：GPS，伪随机码，MATLAB，C/A 码ABSTRACTABSTRACT GPS (Global Positioning System) satellite-based positioning system, is cons

9、idered the most important development direction in the current positioning and navigation equipment, has been widely applied in military and civilian. Everyday life is a difficult figure for broad areas. GPS like mobile phones, fax machines, computers and the Internet in our lives, peoples daily liv

10、es will do without it. The GPS system uses a typical CDMA system, such a spread spectrum modulation signal has a low probability of intercept characteristics. The system is mainly the use of direct sequence spread spectrum modulation techniques, the use of pseudo-code C / A code and P (Y) code two.

11、At present, the pseudo-random code is widely used in spread spectrum communications using this technology is also satellite navigation system. In the GPS navigation satellite navigation and positioning signals sent to the user, it is using to send the navigation message of pseudo-random code phase m

12、odulation signal. This paper first describes the composition of the GPS satellite system, and briefly discusses the M sequences and pseudo-random noise code (P code and C / A code) and it focus on the principle of the C / A code and the MATLAB software, C / A code generate the simulation.KEYKEY WORD

13、SWORDS: GPS，pseudo random noise code，MATLAB，C/A code目目 录录第一章 绪论.1 引言.11.2 课题研究的意义和目的.21.3 本文的内容和组织结构.3第二章 GPS 系统原理概述 .4 GPS 卫星信号 .4 GPS 卫星的载波信号 .7 GPS 卫星的导航电文 .8第三章 GPS 信号构成 .93.1 M 序列 .9 m序列产生的原理 .9 m序列的性质 .11 m序列特征多项式确定 .13 GOLD 组合码 .14 C/A 码 .15 C/A 码发生器原理 .16 C/A 码发生器实现结构 .17 C/A 码的自相关系数 .18 C/A

14、 码的互相关系数 .19 P 码 .20 P 码产生原理 .22 P 码特性 .24第四章 GPS 卫星导航算法及其 MATLAB 软件仿真 .26 MATLAB 软件简介 .26 MATLAB 功能介绍 .27 M 序列仿真及其仿真结果 .27 C/A 码仿真及其仿真结果 .29 P 码仿真及其仿真结果 .30 本章小结.32全文总结.33致谢.34参考文献.35毕设总结.37附录.38第 1 章 绪论 引言 全球定位系统简称 GPS 系统，可在全球范围内，全天候为用户连续提供高精度的位置、速度和时间信息。本课题的目的就是搭建一个 GPS 仿真平台，使得各种信号生成的算法能在该仿真平台上得以

15、仿真实现，以验证算法的性能，从而为信号模拟器的研制提供理论依据。因此对 GPS 的卫星信号的仿真必将推动中国自主研制的卫星导航系统的开展。GPS 是美国政府于 20 世纪 70 年代开始研制，于 1994 年全面建成的高精度、高动态的星际导航定位系统，该系统全天候地向全球范围内具有 GPS 接收机用户提供精确、连续的三维位置、三维运动和时间需要。GPS 信号分为民用的标准定位效劳SPS，Standard Positioning Service和军规的标准精确定位效劳PPS，Precise Positioning Service 。目前，以 GPS 为代表的卫星导航应用产业已成为当今国际公认的八

16、大无线产业之一。随着技术的进步、应用需求的增加，GPS 以全天候、高精度、自动化、高效率等显著特点及其所独具的定位导航、授时校频、精密测量等多方面的强大功能，已涉足众多的应用领域，使 GPS 成为继蜂窝移动通信和互联网之后的全球第三个 IT 经济新增长点。 在卫星定位系统领域，美国的技术可谓最常见和最成熟的，该系统是美国历时16年、耗资130亿美元，由军方控制并资助建成的。目前，全球众多国家都在使用这个系统对地面的汽车、海上的船只和天空的飞行器以及卫星、导弹进行全天候的和实时的准确定位，定位精度成10米。如今，GPS 已经是全球性的系统，而且满足二十世纪六十年代提出的最正确定位标准。它能向有适

17、当接收设备的全球范围用户提供精确、连续的三维位置和速度信息。而对于用户来说，最重要的就是 GPS 接收机，接收机根据同时收到的48颗卫星的位置信息，应用差分定位原理，每隔13秒向用户播报一次其位置(经纬度)、速度、高度和时间信息，以供用户或用户的系统使用112。GPS由3个区段组成:卫星星座，地面控制/监视网络和用户接收设备。GPS卫星星座由24颗分布于平均高度为20210km的6个轨道面的卫星组成，卫星运行周期为n小时58分(12恒星时)。地面监控系统由1个主控站、3个注入站和5个监测站组成，用于跟踪观测GPS卫星，计算卫星星历，监测和控制卫星的健康状况，保持精确的GPS时间系统，向卫星注入

18、导航电文和控制命令。GPS用户设备的核心是GPS接收机，其主要功能是接收并处理GPS卫星发播的导航信号，计算出用户接收机的位置、速度和时间信息3。目前，卫星导航定位系统中大多采用伪码扩频技术。GPS 卫星发射的导航信号是将基带信号先经伪随机噪声码也称伪随机码或者伪码扩频，再对载频进行 BPSK 调制形成新的信号，采用这种信号格式，可使系统具有较高的抗干扰能力和保密性，其关键在于使用了自相关特性很强的伪随机码扩频技术，GPS系统中使用了两种伪随机码，一种是时钟速率为 10.23MHz 用于精密测距的精确军用码，简称 P 码；另一种是时钟速率为 1.023MHz 用于分址，搜捕卫星信号的粗捕获民用

19、码，简称 C/A 码。它在一个周期中有 1023 个码位，周期为1ms。GPS 的 P 码序列捕获要求首先捕获 C/A 码序列导航电文中的转换码，以辅助完成捕获 P 码。利用 C/A 码快速、粗略的进行测距，来引导 P 码进行精确测距。然而，C /A 码的码长短、码速率低，易受敌方干扰和欺骗， 在强干扰和欺骗的战争环境下，很难通过 C/A 码来捕获到 P 码，直接捕获 P 码一直倍受美国军方的关注，产生 P 码并对其特性进行分析对进一步研究 P 码的捕获有着重要的意义。1.2 课题研究的意义和目的GPS 系统采用的是直接序列扩频通信体制，在两个 L 波段上 ，分别用伪随机噪声码 P 码和 C

20、/A 码扩频调制了导航电文 D 码。L1 的中心频率为1575.42MHz，L2 的中心频率为 1227.6MHz。载波信号 L1 上调制有 P 码、C /A码和 D 码，而载波信号 L2 上只调制 P 码和 D 码6。 GPS 系统中 P 码的捕获通常是利用 C/A 码来完成的，用户首先捕获到 C/A码，然后利用 C/A 码调制的导航电文中的握手字(HOW-handover word)所提供的P 码信息对 P 码进行捕获。由于 P 码在战争中显得十分重要, 直接捕获 P 码一直倍受美国军方的关注。产生 P 码并对其特性进行分析对进一步研究 P 码的捕获有着重要的意义，而 C/A 码在民用中也

21、发挥了很重要的作用。 1.3 本文的内容和组织结构 本文研究的题目是 GPS 卫星系统中 PRN 码伪随机码发生器的理论原理，主要研究的的是 C/A 码和 P 码的产生原理。本文的组织结构如下： 第一章 绪论，介绍了GPS全球定位系统的概况，概述了课题的背景、研究意义和目的、国内外研究现状以及伪随机码在GPS卫星系统的重要程度，给出了本文的主要内容和组织结构。 第二章 GPS系统原理概论，主要介绍了GPS系统的组成和GPS信号的组成，并对GPS载波信号和导航电文做了简单介绍，阐述了伪随机码信号。 第三章 GPS信号的构成。主要介绍构成GPS信号的m序列、CA码、P码；通过了解其性质、原理及结构

22、，深化对GPS信号的进一步认识。 第四章 对MATLAB软件的运用做了简单的论述，并通过了解m序列的产生，研究C/A码及P码的产生原理及性质，在MATLAB中对其进行了仿真，这是本章重点。 第二章 GPS 系统原理概述 GPS 卫星信号 GPS 卫星向广阔用户发送的用于导航定位的信号，是一种调制波，但有别于常用的无线电播送电台发送的调频调幅信号，它是利用伪随机噪声码传送导航电文的调相信号。GPS 卫星信号是目前常用的两种违心导航定位信号之一，它包含有三种信号分量，即载波L1 和 L2、测距码C/A 码和 P 码和数据码D 码，亦称基带信号或导航电文。而这所有这些信号分量都是在同一个根本频率 f

23、0113。GPS 卫星信号示意图如图 2.1 所示： 基本频率f0=10.23MHzL11575.42MHzC/A码1.023MHzP码10.23MHzL21227.6MHzP码10.23MHz50BPS数据码（导航电文、或D码）1/10154120 图 GPS 卫星信号频率构成示意图GPS 卫星发送的 GPS 卫星信号采用 L 波段的两种载频作载波，分别被称作 L1 的主频率和 L2 的次频率。这些载波频率由扩频码每一颗卫星均有专门的伪随机序列和导航电文所调制。所有卫星均在这两个相同的载波频率上发射，但由于伪随机码调制不同，因此无明显的相互干扰16。GPS 使用 L 频段的两种载频为其中 f

24、0是卫星信号发生器的基准频率： L1 载波：fL1=154f0=1575.42 MHz，波长 1=19.032 cm； L2 载波：fL2=120f0 =1227.6MHz，波长 2=24.42 cm。选择 L 波段的好处是：(1减少拥挤，防止“撞车。目前 L 波段的频率占用率低于其他波段，与其他工作频率不易发生“撞车现象，有利于全球性的导航定位测量。(2适应扩频，传送宽带信号。GPS 卫星采用扩频技术发送卫星导航电文，其频带高达 20 MHz 左右，在占用率较低的 L 波段上，易于传送扩频后的宽带信号。在载波 L1 上调制有 C/A 码、P 码的数据码，而在载波 L2 上，只有调制有 P码的

25、数据码。在无线通信技术中，为了有效地传播信息，一般均将频率较低的信号加载到频率较高的载波上，而这时频率较低的信号成为调制信号。GPS 信号是一种调制波，它不仅采用 L 波段的载波，而且采用扩频技术传送卫星导航电文。所谓“扩频，是将原来打算发送的几十比特速率的电文变换成发送几兆甚至几十兆比特速率的由电文和伪随机噪声码组成的组合码。采用扩频技术时，假设信号功率仅为噪声功率的 1/10，那么信号将深深地淹没在噪声之中而不易被他人捕获，从而使得信号具有极强的保密性。GPS 信号的调制波，是卫星导航电文和伪随机噪声码PseudoRandom Noise Code，简称 PRN 码，或称伪噪声码的组合码。

26、卫星导航电文是一种不归零二进制码组成的编码脉冲串，称之为数据码，记作 D(t)，其码率为 50 b/s。对于距离地面 20000 km 之遥的 GPS 卫星，扩频技术能有效地将很低码率的导航电文发送给用户。其方法是用很低码率的数据码作二级调制扩频。第一级，用 50Hz 的 D 码调制一个伪噪声码，例如调制一个被叫做 P 码的伪噪声码，它的码率高达 10.23 MHz。D 码调制 P 码的结果，便形成了一个组合码P(t)D(t)，使得 D 码信号的频带宽度从 50Hz 扩展到 10.23 MHz，也就是说，GPS 卫星从原来要发送 50b/s 的 D 码，转变为发送 10230 b/s 的组合码

27、 P(t)D(t)。在 D 码调制伪噪声码以后，再用它们的组合码去调制 L 波段的载波，实现D 码的第二级调制，而形成向广阔用户发送的已调波。D 码的数据首先同伪噪声码 C/A 码和 P(Y)码模二相加后，形成组合码 C/A(t)D(t)和 P(t)D(t)，然后才调制 L1 载波。需要注意的是，组合码 C/A(t)D(t)和 P(t)D(t)是通过相移键控BPSK调制到 L1 载波上的。在 L1 载波上，C/A(t)D(t)调制和 P(t)D(t)调制在相位上是正交的。因此在这两个合并的 L1 载波频率上的 C/A(t)D(t)调制和 P(t)D(t)调制之间有 90的相移。L2 载波上的调

28、制过程与 L1 载波大致相同，不同的是 L2 载波可以用 C/A(t)D(t)码、P(t)D(t)码或者 P(Y)码来调制。最后，卫星向地面发射这两种已调波 L1 和 L2。 L1载波C/A码导航电文P码L2载波卫星时钟F=10.23MHz=90模2和加法器；调制器；信号合成器 图 GPS 信号的产生框图需要注意的是，GPS 信号虽然有几种分量C/A 易捕获码、 P 精确码和 D 导航数据码，但是它们均来源于一个公共的 10.23MHz 的基准频率。它们的频率不仅与基准频率有一定的比例关系，而且相互之间也存在一定的比例关系，详细如表 2-1 所示。这既有利于 GPS 卫星发送信号，又便于广阔用

29、户接收和测量GPS 信号。从表 2-1 中可以看出，在 D 码的一个码元内，将有 20460 个 C/A 码码元，204600 个 P 码码元，31508400 个 L1 周期和 24552000 个 L2 周期。相关频率基频 F载频 fL1载频 fL2基准频率 F154F120FC/A 码的码频 fgF/10fL1/1540fL2/1200P 码的码频 fpFfL1/154fL2/120D 码的码频 fdF/204600fL2/31508400fL2/24552000表 GPS 信号的频率关系 GPS 卫星的载波信号 GPS 卫星的测距码信号和导航电文信号都属于低频信号，其中 C/A 码和

30、P码的数码率分别为 1.023Mbit/s 与 10.23Mbit/s,而 D 码数据码的数码率仅为 50bit/s。GPS 卫星离地面约 2104km，其电能又非常紧张，因此很难将上述数码率很低的信号传输到地面。解决这一难题的方法，就是另外发射一种高频信号，并将低频的测距码信号和导航电文信号加载到这一高频信号上，构成一高频的调制波发射给地面。GPS 卫星采用 L 频带的两种不同频率的电磁波作为高频信号，分别称为 L1载波与 L2载波6。其中：L1载波的频率f1=1575.42MHz，波长1=19.03cm，其上调制 C/A 码、P 码以及导航电文；L2载波的频率f2=1227.6MHz，波长

31、 2=24.42cm，其上仅调制 P 码和导航电文。GPS 卫星发射信号的频率，都要受卫星上原子钟的基准频率的控制。GPS 微星原子钟基准频率的 1/10，而 L1载波的频率f1为基准频率154 倍，L2载波的频率f2那么取基准频率f0的 120 倍。 GPS 卫星的导航电文1 2 3 4 530s1 2 3 4 5 6 7 8 9 106s0.6s0.02s子帧4、5含25页1个帧含5个子帧1个子帧含10个字一个子帧含30bit 导航电文的组成 导航电文是由 GPS 卫星向用户播发的一组包含卫星在空间的位置、卫星状态、卫星时钟的修正参数、Z 计数等重要数据的二进制代码，它是用户使用 GPS进

32、行定位时必不可少的数据。导航电文的传输速率为 50bit/s，它以帧为单位向外发送，每个主帧 1500bit，每帧 30s1。如图 2-3 所示，每个主帧分为 5个子帧，每个子帧包含 300bit，长度为 6s。每个子帧又可分为 10 个字，每个字又是由 30bit 组成的。其中第四子帧、第五子帧各有 25 个页面，需要 750s才能全部播发完毕，也就是每 30s 翻转一页，12.5min 完整地播发一次，其内容仅仅在卫星注入新的导航数据后才得以更新。第一子帧、第二子帧、第三子帧每隔 30s 重复一次，其内容每小时更新一次，所以在卫星定位后可以根据已有的导航电文预测下一段时间导航电文6。 第三

33、章 GPS 信号构成3.1 M 序列伪随机码又称伪随机序列，是一种可以预先确定并重复产生和复制具有随机统计特性的二进制码序列。在现代工程实践中，伪随机信号在移动通信、导航、雷达和保密通信、通信系统性能的测量等领域中有着广泛的应用。例如，在连续波雷达中可用作测距信号，在遥控系统中可用作遥控信号，在多址通信中可用作地址信号，在数字通信中可用作群同步信号，还可用作噪声源以及在保密通信中的加密作用等。伪随机发生器在测距、通信等领域的应用日益受到人们重视。伪随机信号与随机信号的区别在于，随机信号是不可预测的，它在将来时刻的取值只能从统计意义上去描述。伪随机序列实质上不是随机的，而是收发双方都知道确实定性

34、周期信号。之所以称其为伪随机序列，是因为它表现出白噪声采样序列的统计特性，在不知其生成方法的侦听者看来像真的随机序列一样。m序列作为一种根本的序列，具有很强的系统性、规律性和相关性。m序列是伪随机序列中最重要的一种，是最长线性移位存放器序列，m序列易于实现，它具有优良的自相关特性，是在直接扩频通信系统中用于扩展要传递的信号。即m序列是由多级移位存放器构成，假设 N 为移位存放器的级数，n级移位存放器共有个状态，除去全 0 状态外还剩下种状态，因此它能2n21n产生的最大长度的码序列为位。产生m序列的线性反响移位存放器称作最21n长线性移位存放器。产生m序列的移位存放器的电路结构，其反响线连接不

35、是随意的，m序列的周期P也不能取任意值，而必须满足，式中，n是移21np 位存放器的级数。m序列码发生器是一种反响移位型结构的电路，它由n位移位存放器加异或反响网络组成，其序列长度，只有一个多余状态即全 0 状态，所以21nM 称为最大线性序列码发生器5。 m序列产生的原理移位存放器的后续状态可以用当前状态及特定矩阵来表示，这个矩阵是阶矩阵，称为 A A 矩阵，A A 矩阵的第r行对应移位存放器第r级反响输入状n n态。对于一个n级移位存放器序列产生器，其 A A 矩阵的第一个元素必定为1na1，否那么，该序列发生器就必然退化为级数小于n的移位存放器序列发生器，一个n级简单线性移位存放器序列发

36、生器，其 A A 矩阵有如下形式12311100000100000010nn nccccA在给定移位存放器的初始状态后，可由 A A 矩阵求出后续状态，即111111212221222212( )(1)(1)( )(1)(1)( )(1)( )(1)(1)nnnnnnnnnxjxjxjaaaxjxjaaaxjX jA X jAaaaxjxjxj因此，且1( )(1)nirrrxja xj()( )kX jkAX j当单位矩阵时，有，即移位存放器中的内容在第kAI()( )X jkX j个状态和第个状态是相同的，即序列发生器从第个状态开始，经过jjkj次状态转移后，又回到了第个状态，产生的序列长

37、度就为。因此，对于最kjk大长度线性移位存放器序列发生器，必然有 3.121nNAAI 对于阶矩阵 A A，假设x为其特征值，那么有，由于在二进制n n0AxI系统中，整理化简得11 0AxI 3.2121231310nnnnnxC xC xC xCx 定义特征方程和特征多项式分别为 3.30( )nn rrrF xC x 3.40( )nrrrf xC x其中。01,1nCC m序列的性质扩展频谱通信要求扩频序列应具有较好的随机特性，而序列既具有一定m的随机特性，又具有一定的周期性，因而它是一种伪随机序列，序列应具有m线性叠加特性、均衡特性、游程分布特性和自相关特性611。一个m序列与其经任

38、意次延迟移位后产生的另一个不同序列na模 2 相加，得到是仍是该m序列的延迟移位序列。例如 1110100 与向kna右移 3 位后的序列 1001110 逐位模 2 相加后的序列为 0111010 ，相当于原序列向右移 1 位后的序列，仍是m 序列。2.均衡特性在m序列的每个周期中， “1码元出现的数目为次， “0码元出21n12n现的数目为次，即“0的出现个数比“1的出现个数少一个。121n伪随机序列的平衡性是指序列中“1的数目只比“0 的数目多l。码的平衡性由码序列中的直流分量决定。平衡性好，那么载波抑制度大，从而有利扩频通信出现一些稳定的信号，从而易于被检测而导致保密能力的降低。分布特

39、性在序列当中，相同的码元连在一起，称为一个游程，一般来说，长度为 1的游程占总数的 1/2，长度为 2 的游程占总数的 1/4，依次类推，连“1的游程和连“0的游程各占一半。1、m序列的一个周期p=2n-1 中游程总数为 2n-1 。2、当 1游标长度kn-1,游程数目占总数 2-k； 当,连 1 和连“0游程各占一半。12kn 游程长度游程游程数41111130001211，00210，1，0，14 表 m序列的游程分布示意图 4.自相关特性 在一般情况下，假设有两个相同周期T的序列，那么两者之间的互相关系数定义为 3.5 121dx t x ttT 式中相对于的时间延迟； 2xt当时，就是

40、自相关系数，通俗地说，相关系数表示 12x txt t t序列和之间的相似程度。对于时元的离散采样的自相关系数为 1x t 2xtkt 3.6 11nkkitX tX tLP式中按前述二进制信号波形乘积法那么进行运算求得。对于m序列而言，它的自相关系数为： 3.7 01,0, 1, 21tiTP ittjjjT 等于除部位0和LP的整数倍以外的任何数对于一个周期为的m序列取值 1 或 0 ，其自相关函数如下列21np nana图。由图可见，当时，m序列的自相关函数出现峰值 1；当偏离 00( )R时，相关函数曲线很快下降；当，相关函数值为；当时，11p1 pp又出现峰值；如此周而复始。当周期

41、P 很大时，m序列的自相关函数与白噪声类似。这一特性很重要，相关检测就是利用这一特性，在有或无信号相关函数值的根底上识别信号，检测自相关函数值为 1 的码序列。从m 序列的自相关系数可以看出,m 序列是一个狭义伪随机码。当m序列的移位值为其周期的整数倍时，其自相关值取得最大值为 1，移位值取其他值时，其自相关值恒为。1 p R R1p-110P-1P+1pt图 m序列的自相关函数 m序列特征多项式确定n级线性移位存放器的如图 2.4 所示，图中表示反响线的两种可能连0iC 接方式，表示连线接通，第n-i 级输出参加反响中。表示连接线断1iC 0iC 开，第级输出未参加反响。ni c1c2c3C

42、n-1an-1an-2an-3a1a0C0=1Cn=1输出 图 n级线性移位存放器结构因此，一般形式的线性反响逻辑表达式为 3.8112201nnnnnin iiaC aC aC aC a将等式左面的移至右面，并将代入上式，那么上式可改写为 na001nnaC aC 3.9110niniC a定义一个与上式相对应的多项式 3.10 1niiiF xC x其中 x 的幂次表示元素的相应位置。该式称为线性反响移位存放器的特征多项式，特征多项式与输出序列的周期有密切关系。当 F(x)满足以下三个条件时，就一定能产生m序列：(1) F(x)是不可约的，即不能再分解多项式； (2) F(x)可整除,这里

43、；1px 21np (3)F(x)不能整除，这里。1qx qp可以确定m序列的特征多项式： 3.11 011nniniif xcc xc xc x GOLD 组合码在扩频通信系统中，大都采用线性或非线性移位存放器产生伪随机序列作为扩频序列，例如m序列、Gold 序列，因为它们不仅具有理想的随机特性，而且较容易产生。其中，最简单、最常用的伪随机序列是m序列。m序列具有锋利的自相关特性，有较小的互相关值，码元平衡。但序列数目不多，序列复杂度不大，这对于采用码分多址技术的 GPS 系统是不利的。1976 年 RGold 提出了新的一类序列：Gold 序列。它的自相关旁瓣值和互相关值与m序列的互相关值

44、一样，序列数目大大增加，序列复杂度也略有改善，是又一重要的伪随机序列，且在扩频通信等现代通信领域得到广泛应用811。 Gold 码是由两个周期和速率相同的码元构成不同的m序列组合而成。例如，现在具有相同周期LP=2n-1的两个m序列：和，那么它们构成的 Gold X Y序列为 3.12 GXY式中：为向左移动了j个码元的序列，而j= jY0,1,2,22n从式可以看出，由两个m序列和以及左移j个码元的，可以 X Y jY构成，共个 XY 1XY 2XY 21nXY21nGold 序列，连同原序列和，总共有个序列，成为 Gold 序列族。 X Y21n从上可见，Gold 码具有以下特点。 1)

45、Gold 码的速率和周期与构成它的m序列相同。 2) Gold 码不是m序列，它的互相关值可以用 Anderson 导出的下述公式进行估算 3.13 12112 12LPLP 式中：LP 为 Gold 码的长度周期。 3) Gold 码的结构简单，调整方便，具有大量的可用码型，适宜于码分多址的大量用户需求。N 级移位存放器能够产生独立的m序列数目为 3.1421nnJn式中：为尤拉函数，其数值等于数值中与为互21n0,1,2,22n21n素的正整数之个数。 C/A 码目前，卫星导航定位系统中大多采用伪码扩频技术。GPS 卫星发射的导航信号是将基带信号先经伪随机噪声码(也称伪随机码或伪码)扩频，

46、再对载频进行 BPSK 调制形成的信号，采用这种信号格式，可使系统具有较高的抗干扰能力和保密性，其关键在于使用了自相关特性很强的伪随机码扩频技术。GPS 系统中使用了两种伪随机码，一种是时钟速率为 10.23MHz 用于精密测距的精确军用码，简称 P 码。另一种是时钟速率为 1.023MHz 用于分址、搜捕卫星信号的粗捕获民用码，简称 C/A 码。它在一个周期中有 l023 个码位，周期为 1ms。GPS的 P 码序列捕获要求首先捕获 C/A 码序列导航电文中的转换码，以辅助完成捕获 P 码。利用 C/A 码快速、粗略的进行测距，来引导 P 码进行精确测距。所以，C/A 码在 GPS 的导航定

47、位中起着至关重要的作用816。GPS 系统中 P 码的捕获通常是利用 C/A 码来完成的,用户首先捕获到 C/A 码,然后利用 C/A 码调制的导航电文中的握手字(HOW 2hand overword)所提供的 P 码信息对 P 码进行捕获。由于 P 码在战争中显得十分重要,而且 C/A 码在民用中也发挥了很重要的作用,所以研究并实现 C/A 码,就有很重要的意义了。C/A 码是用于粗测和捕获的伪随机码。它是m序列的优选对组合码形成的 Gold 码(简称G 码)。G 码是长度相同而且互相关极大值最小的m序列逐位进行模二相加所构成的。改变产生它的两个m序列的相对位,就可以得到不同的码。对于长度为

48、N=2K-1 的m序列,每两个码就可以用这种方法产生 N 个不同的 G 码。在这 N 个码中,任何两个互相关最大值等于构成它们的两个m序列的互相关最大值。G 码的自相关旁瓣有起伏,但是,它的峰值不超过相关最大值。G 码最主要的优点在于广泛用于多址通信,这就是 GPS 中 C/A 码采用的主要原因。 C/A 码发生器原理GPS 中 C/A 码是 Gold 码，其序列长度为 1023 位。C/A 码是由两个 10 级移位存放器产生的m序列 G1 和 G2，经模二和后产生的。Gold 码的周期和速率与构成它的m序列是一致，但是，改变两个m序列之间的相位关系，可以组合成一种新的 Gold 码，以致第

49、j 颗卫星的 C/A 码为 (3.15) 1210jjpGtG t GtN式中：G1G2为构成 C/A 码的m序列； 为下述 P 码的码元宽度，且知=1/10.23MHz。ppNj为第j颗 GPS 卫星 C/A 码的两个m序列 G1G2之间的相位偏差系数，其值为正整数，以致不同的 GPS 卫星都有不同的 C/A 码，便于广阔用户做导航定位测量时，识别 GPS 卫星，捕获和跟踪到所需的 GPS 卫星的导航定位信号。10级线性反馈移位寄存器G1全置“1”电路时间基准源10分频器10级线性反馈移位寄存器G2平移等价序列形成电路C/A码G11.023MHZ时钟脉冲G2f10.23MHZ图 C/A 码的

50、生成框图z，它们的周期均为 1ms。在一个 C/A 码周期内具有 1023 个码元，换言之，C/A 码是由两个 1023bit 的m序列构成的，其发生器如下列图。图中两个 10 级线性反响移位存放器的特征多项式分别为 31011Gxxx ) 236891021Gxxxxxxx 依据上列两个特征多项式，便可生成两个确定的m序列，其一经过平移等价序列选择电路也称为相位选择器 ，致使不同的 GPS 卫星 C/A 发生器，能够产生不同的 C/A 发生器，亦即，变更为子序列 G2 的输出模二和。例如，对PRN01 卫星而言，采用第 2.6 级输出模二和，其输出假设令为，那么 PRN012,62G卫星的

51、C/A 码为 ()2,62011PRNC AGG对 PRN01 卫星而言，采用第 3、7 级输出的模二和，故知 PRN02 卫星的 C/A码为 (3.18)3,72021PRNC AGG由于 G1,G2 两个m序列均有 1023 个码元，G2 的平移等价序列便多达 1023个，加上不平移的 G2 和 G1 的组合，便可生产 1024 个周期为外围 1msec 和1023bit 的 C/A 码，实际上只选用其中的 37 个 C/A 码8。 C/A 码发生器实现结构C/A 码发生器是由两个 10 级的反响移位存放器产生的 G 码所构成的,其产生的码长都为 N=210-1=1023，在这两个反响移位

52、存放器分别对应两个m序列 G1和 G2，它们的特征多项式如式 3.16 所示。其中它们在统一时钟驱动下工作，当全局复位信号到来时,它们被赋于统一初相，各个存放器使用一个二输入异或门或六输入异或门通过线性反响方式连接8。通过在 G2 存放器对应位抽取抽头，并使各抽头通过异或产生 G2 序列的不同平移等价序列。其中,抽头的设定方式根据 GPS 卫星所规定的 C/A 码初相分配表来实现,比方 5 号 GPS 卫星规定相位选择为:19，那么对应 G2 存放器就设定一号、九号存放器的抽头相互异或,以得+到 5 号 GPS 卫星的 G2 序列的平移等价序列。得到的序列结果再和 G1 存放器输出结果相异或，

53、那么输出的就是 GPS 卫星的 C/A 码了。具体结构如图 3.4 所示。采用不同的抽头方式。可以得到的 1023 个 G2 序列,同时,再加上 G1、G2 本身,就可以得到 1025 种结构不同的 C/A 码，这些 C/A 码都具有相同的码长：1023bits，相同的码宽： 3.19usfT98. 01从这些码中可以选出 32 种，用以分配给 GPS 卫星。由于 C/A 码码长短，每秒可以搜索 1000 次,所以 C/A 码除了可以做粗测码,也可以作为 GPS 卫星的捕获码。C/A码发生器系统工作时钟RST复位信号1 2 3 4 5 6 7 8 9 10计数器Counter=1023码时钟R

54、ET复位信号1 2 3 4 5 6 7 8 9 10G1G2OUTCOUNTER RST图 3.4 C/A 码发生结构 C/A 码的自相关系数自相关函数P（t）Pt-RcRc0Gold码频谱0-20-30-40-50123456789边带数|S(f)|2/PsP=1023码元的C/A码频谱|S(f)| C/A 码的自相关系数和频谱C/A 码的自相关系数等于它及其时延序列乘积的积分平均值，亦即 jjGt Gt 3.20 1212G t Gtk G tGtk式中：横线表示对时间的平均值 k表示相位偏差系数且知10jpkN当时延时，C/A 码的自相关系数时，C/A 码的自相关系数0 3.21 121

55、201gGt Gtk Gt Gtk 当时延时，C/A 码的自相关系数等于构成它的两个m序列的互相关0系数，亦即： 3.22 12gGt Gtn C/A 码的互相关系数 C/A 码的互相关系数等于两个 C/A 码序列乘积的积分平均值，亦即 jjgGt Gt 3.23 21211tGGtk GtGt 当时延时，C/A 码的互相关系数为0 3.24 10gLP 当时延时，C/A 码的互相关系数等于它的自相关系数，即等于m序列0的互相关系数。综上可见，C/A 码具有一千余个可用码型，能够给相应数量的卫星分配确定的各自独立的 C/A 码。而且所有的 C/A 码均具有相同的 1ms 周期和 1023 个码

56、元，后者可以是 GPS 信号接收机仅以较短的时间搜索和捕获到 GPS 卫星发送的C/A 码，快速实现首次导航定位测量17。 P 码GPS 系统中 P 码的捕获通常是先捕获到 C/ A 码，然后利用 C/ A 码调制的导航电文中的转接字提供的 P 码信息对其进行捕获。然而 C/A 码的码长短、码速率低，易受敌方干扰和欺骗，在强干扰和欺骗环境下，很难通过 C/ A 码来捕获到 P 码；因此，直接捕获 P 码一直倍受美国军方的关注；产生 P 码并对其特性进行分析对进一步研究直接对 P 码捕获的接收机设计有着重要的意义5。P 码是由两个载波发送给 GPS 用户的另一个伪噪声码，是一个具有2.35E+1

57、4 个码元的特长序列，是美国军方严格控制使用的保密军用码。P 码的序列长度为 2.354695927651014,码速率为 10.23MHz ,序列周期 266. 41 天或 38.05 星期；其中，每颗 GPS 卫星使用序列的一个星期，故每个序列长度是 6.1871012,该码元长度远远大于 C/A 码 1023 的码元长度。P 码是由两个子码 X1，X2 构成的复合伪噪声码；每一颗 GPS 卫星采用各自不同的 P 码，其区别在于第二个子码 X2 存在一个相位偏差系数。其值为jn037 的正整数假设了一个不存在的零号卫星 。因此，第j颗 GPS 卫星的 P码其生成如图 3.6 所示为 3.2

58、5 120jjPtXt Xtn用于同步导航电文、子帧、C/A码和P码X1子码周期输出212级线性移位寄存器X1212级线性移位寄存器X2移位电路X1P码X2全置1电路时间基准源F=10. 23M HZ每星期六午夜零点置初态1. 5s周期脉冲时钟脉冲 P 码发生器的原理框图且规定，当子码 X1 的长度周期为 15345000bit 时，子码 X2 的最大长度周期那么为 15345037bit，因此，P 码的最大长度周期为: 3.2615345000 15345037235469592765000pLPbitP 码的码率为 10.23MHZ，故其码元宽度；P 码的最大61 10.23 10ps时间

59、周期为: 3.27266 9 45min55.5pppTPLPd hs由上可见，P 码是一个长达 2.351014个码元的伪噪声码，其时间周期到达266d 多。为了捕获和跟踪到如此之长的 P 码，对它采用了下述两项措施。1P 码的时间周期“截短。每个星期日的子夜零点作为截短周期的起点，每个星期六的午夜 24 时作为截短周期的重点，亦即 P 码的使用周期为 7d,而每颗 GPS 卫星又具有不同的 P 码。2采用分布捕获法。P 码的时间周期虽然被截短称为 7d，但 GPS 信号接收机仍不易捕获到 7d 周期的截短 P 码，而采用了二步捕获法：首先捕获和跟踪到一个仅有 1023 个码元的 C/A 码

60、，解译出它所传送的卫星导航电文，依据该电文转换码所提供的 P 码捕获信息，用户可以较快的捕获到截短周期为 7d 的 P 码，以使用 P 码进行导航定位测量19。 P 码产生原理P 码的两个子码 X1，X2 也是一种复合伪噪声码，它们均由 2 个 12 级线性移位存放器产生的m序列模二和相加而成5。对于子码 X1 而言,它的 2 个 12 级线性移位存放器的特性多项式分别为 68111211afXxxxx 3.281258910111211bfXxxxxxxxx 对于子码 X2 而言，它的 2 个 12 级线性移位存放器的特性多项式分别为 134578910111221afXxxxxxxxxxx