基于MATLAB的GPS信号的仿真设计

摘要
扩频通信是近几年来迅速发展起来的一种通信技术。

在早期研究这种技术的主要目的是为提高军事通信的和抗干扰性能，因此这种技术的开发和应用一直是处于状态。

扩频技术在军事应用上的最成功例可以以美国和俄国的全球定位系统（GPS和GLONASS）为代表；在民用上GPS和GLONASS也都得到了广泛的应用，这些系统的基础就是扩频技术。

全球定位系统（GPS）用于对全球的民用与军用飞机、舰船、人员、车辆等提供实时导航定位服务。

GPS系统采用典型的CDMA体制，这种扩频调制信号具有低截获概率特性。

该系统主要利用直接序列扩频调制技术，采用的伪码有C/A码和P(Y)码两种。

本文讲述的是直接序列扩频通信技术在全球定位系统（GPS）中的应用。

主要介绍扩频通信中的伪码仿真，简要论述M序列和伪随机噪声码（P码和C/A码）与其产生，并使用MATLAB7.0仿真M序列、P码和C/A码的编码过程和仿真结果，介绍直扩频技术伪码的相关知识，重点介绍P码。

关键字：全球定位系统；直接扩频通信；伪码仿真
Abstract
Spread spectrum communication is a communications technology developed rapidly in recent years. In early studies the main purpose of this technology is to improve the military communications confidential and anti-jamming performance, therefore the development and application of this technology is always in secret state. Spread spectrum technology in the most successful military application examples are the United States and Russia could the global positioning system (GPS and GLONASS) for representative; In civil GPS and GLONASS also have been widely used，which foundation of system is the spread spectrum technology.
Global positioning system (GPS) is used to provide real-time navigation and positioning services for global civil and military aircraft, ships, personnel, vehicles and so on. GPS system adopts the typical CDMA system, which kind of spread spectrummodulation signals have low intercept probability characteristic. This system mainly used the direct sequence spread spectrum modulation technology, using the PRN code including C/A code, P codes and Y codes.
This article tells the direct sequence spread spectrum communication technology applied in global positioning system (GPS) .The article mainly introduces the pn code spread spectrum communication simulation, briefly discussing M sequence and pseudo random noise code (P yards and C/A yards) and its produce and use MATLAB7.0 simulate M series, P yards and C/A yards of encoding process and the simulation results, introducing pn code straight spread-spectrum technology knowledge, especially P yards.
Key: GPS; DS-SS;Pn code simulation
目录引言4
1GPS理论与其特性5
1.1GPS系统概述5
1.2GPS信号构成7
1.2.1M序列10
1.2.2C/A码13
1.2.3P码15
1.3小结19
2MATLAB软件19
2.1MATLAB软件简介20
2.2MATLAB应用概述21
2.2.1 MATLAB功能介绍21
2.2.2 MATLAB使用方法23
3 GPS卫星导航信号算法与其MATLAB仿真29
3.1 C/A码仿真代码与其仿真结果29
3.2 P码的仿真代码与其仿真结果32
3.3 结果分析与其相关性分析36
结论37
致38
参考文献39
附录A 英文原文40
附录B 中文翻译41
附录C C/A码源代码41
附录D P码源代码54
引言
全球卫星定位系统，简称GPS系统，可在全球围，全天候为用户连续地提供高精度的位置、速度和时间信息。

文中目的就是搭建一个GPS仿真平台，使得各种信号生成的算法能在该仿真平台上得以仿真实现，以验证算法的性能，从而为信号模拟器的研制提供理论依据。

因此对GPS的卫星信号的仿真必将推动中国自主研制的卫星导航系统的发展。

目前，以GPS为代表的卫星导航应用产业已成为当今国际公认的八大无线产业之一。

随着技术的进步、应用需求的增加，GPS以全天候、高精度、自动化、高效率等显著特点与其所独具的定位导航、授时校频、精密测量等多方面的强大功能，已涉足众多的应用领域，使GPS成为继蜂窝移动通信和互联网之后的全球第三个 IT经济新增长点。

GPS系统中P码的捕获通常是先捕获到C /A码，然后利用C /A 码调制的导航电文中的转接字(HOW )所提供的P码信息对P码进行捕获。

然而，C /A 码的码长短、码速率低，易受敌方干扰和欺骗, 在强干扰和欺骗的战争环境下，很难通过C /A 码来捕获到P码。

因此，直接捕获P码一直倍受美国军方的关注。

产生P码并对其特性进行分析对进一步研究直接P码的捕获有着重要的意义。

本设计所针对的是GPS卫星信号的伪码仿真，主要是M序列的生成方法和P码的编码原理与其仿真。

本课题研究思路是：
⑴研究GPS信号的生成原理和直扩频通信技术基础，M序列和P码算法原理。

⑵在理论研究的基础上，实现MATLAB软件生成M序列和P码，并对其做简单的相关性分析。

文章最后对仿真过程中实现的功能与出现的为题做出总结，总结论文所做的工作和需要更仔细研究的方向。

1GPS理论与其特性
1.1GPS系统概述
全球定位系统（GPS）是美国政府于20世纪70年代开始研制，于1994年全面建成的高精度、高动态的星际导航定位系统，该系统全天候地向全球围具有GPS接收机用户提供精确、连续的三维位置、三维运动和时间需要。

GPS信号分为民用的标准定位服务（SPS，Standard Positioning Service）和军规的标准精确定位服务（PPS，Precise Positioning Service）。

GPS目前处于良好的运行状态，并满足20世纪60年代所提出的最佳定位系统标准。

这个系统向有适当接收设备的全球围的用户提供精确、连续的三维位置和速度信息。

GPS也向全球广播世界协调时（UTC）。

组成卫星星座的24颗卫星被安排在6个轨道平面上，即每个平面上4颗。

这样的卫星星座配置确定了卫星某一时刻在轨道中的位置，而为hi与地平线以上的卫星数会随着时间和地点的不同而不同，最少可见到颗，最多可见到11颗，因此保证了在开放的天空下，地球上和近地空间任一点在任何时刻均可以同时观测到至少4颗GPS卫星，为准确定位提供可能。

GPS系统采用CDMA（Code Division Multiple Access，码分多址）技术将在轨的24颗卫星分开，每颗卫星使用同样的调制方式在一样的载波频率上调制各自唯一的伪随机码和数据信息，然后使用星载卫星天线发射信号。

卫星轨道分布如图1.1所示：
图1.1 卫星轨道分布
1.2GPS信号构成
GPS卫星向广大用户发送的用于导航定位的信号，是一种调制波，但有别于常用
的无线电广播电台发送的调频调幅信号，它是利用伪随机噪声码传送导航电文的调相信号。

GPS卫星信号是目前常用的两种违心导航定位信号之一，它包含有三种信号分量，即载波（L1和L2）、测距码（C/A码和P码）和数据码（D码，亦称基带信号或导航电文）。

而这所有这些信号分量都是在同一个基本频率f
=10.23MHz的控制下产生的。

GPS 卫星信号示意图如图1.2所示：
图1.2GPS卫星信号频率构成
GPS卫星发送的GPS卫星信号采用L波段的两种载频作载波，分别被称作L1的主频率和L2的次频率。

这些载波频率由扩频码（每一颗卫星均有专门的伪随机序列）和导航电文所调制。

所有卫星均在这两个一样的载波频率上发射，但由于伪随机码调制不同，
因此无明显的相互干扰。

GPS使用L频段的两种载频为（其中f
是卫星信号发生器的基准频率）：
L1载波：f
L1=154×f
=1575.42 MHz，波长λ
1
=19.032 cm；
L2载波：f
L2=120×f
=1227.6MHz，波长λ
2
=24.42 cm。

选择L波段的好处是：
(1）减少拥挤，避免“撞车”。

目前L波段的频率占用率低于其他波段，与其他工作频率不易发生“撞车”现象，有利于全球性的导航定位测量。

(2）适应扩频，传送宽带信号。

GPS卫星采用扩频技术发送卫星导航电文，其频带高达20 MHz左右，在占用率较低的L波段上，易于传送扩频后的宽带信号。

在载波L1上调制有C/A码、P码的数据码，而在载波L2上，只有调制有P码的数据码。

在无线通信技术中，为了有效地传播信息，一般均将频率较低的信号加载到频率较高的载波上，而这时频率较低的信号成为调制信号。

GPS信号是一种调制波，它不仅采用L波段的载波，而且采用扩频技术传送卫星导
航电文。

所谓“扩频”，是将原来打算发送的几十比特速率的电文变换成发送几兆甚至几十兆比特速率的由电文和伪随机噪声码组成的组合码。

采用扩频技术时，若信号功率仅为噪声功率的1/10，那么信号将深深地淹没在噪声之中而不易被他人捕获，从而使得信号具有极强的性。

GPS信号的调制波，是卫星导航电文和伪随机噪声码（Pseudo Random Noise Code，简称PRN码，或称伪噪声码）的组合码。

卫星导航电文是一种不归零二进制码组成的编码脉冲串，称之为数据码，记作D(t)，其码率为50 b/s。

对于距离地面20 000 km之遥的GPS卫星，扩频技术能有效地将很低码率的导航电文发送给用户。

其方法是用很低码率的数据码作二级调制（扩频）。

第一级，用50 Hz的D码调制一个伪噪声码，例如调制一个被叫做P码的伪噪声码，它的码率高达10.23 MHz。

D码调制P码的结果，便形成了一个组合码——P(t)D(t)，使得D码信号的频带宽度从50 Hz扩展到10.23 MHz，也就是说，GPS卫星从原来要发送50 b/s的D码，转变为发送10 230 b/s的组合码
P(t)D(t)。

在D码调制伪噪声码以后，再用它们的组合码去调制L波段的载波，实现D码的第二级调制，而形成向广大用户发送的已调波。

如图4-1所示，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载波上的调制过程与L1载波大致一样，不同的是L2载波可以用C/A(t)D(t)码、P(t)D(t)码或者P(Y)码来调制。

最后，卫星向地面发射这两种已调波L1和L2。

图1.3 GPS卫星信号的产生
需要注意的是，GPS信号虽然有几种分量（C/A易捕获码、 P精确码和D导航数据码），但是它们均来源于一个公共的10.23MHz的基准频率（见图1-1）。

它们的频率不仅与基准频率有一定的比例关系，而且相互之间也存在一定的比例关系，详细如表4-1所示。

这既有利于GPS卫星发送信号，又便于广大用户接收和测量GPS信号。

从表4-1中可以看出，在D码的一个码元，将有20 460个C/A码码元，204 600个P码码元，31 508 400个L1周期和24 552 000个L2周期。

表1.1 GPS信号的频率关系
相关频率基频F 载频f
L1载频f
L2
基准频率F 10.23MHz 154F 120F
C/A码的码频f
g F/10 f
L1
/1540 f
L2
/1200
P码的码频f
p F f
L1
/154 f
L2
/120
D码的码频f
d F/204600 f
L2
/31508400 f
L2
/24552000
1.2.1 M 序列
码是一种表达信息的二进制数与其组合，是一组二进制的数码序列。

例如，对0，1，2，3取两位二进制数的不同组合表示为：00，01，10，11。

这些二进制数的组合形式称之为码。

其中每一位二进制数称为1个码元或1比特（bit ）；每个码均含有两个二进制数，即两个码元或两个比特。

比特是码的度量单位，也是信息量的度量单位。

如果将各种信息，例如声音、图像以与文字等，按某种预定的规则表示为二进制数的组合形式，则这一过程就称为编码，也就是信息的数字化。

在二进制的数字化信息传输中，每秒所传输的比特数称为数码率，用以表示数字化信息的传输速度，其单位为bit/s(简写为b/s)。

码可以看作是以0和1为幅度的时间函数，用u(t)表示。

因此，一组码序列u(t) ，对于某个时刻t 而言，码元是0或1完全是随机的，但其出现的概率均为1/2。

这种码元幅值是完全无规律的码序列，称为随机噪声码序列。

它是一种非周期序列，无法复制。

但是，随机噪声序列却有良好的自相关性，GPS 测距码就是利用了其自身良好的自相关性才获得成功的。

这里，自相关性是指两个结构一样的码序列的相关程度，它由自相关函数描述。

为了说明这一问题，可将随机噪声码序列u(t)平移k 个码元，获得具有一样结构的新的码序列u(t)。

比较这两个码序列，假定它们的对应码元中，码值（0或1）一样的码元个数为S u ，而码元相异的码元个数为D u ，那么两者之差S u -D u 与两者之和S u +D u （即码元总
数）的比值，即定义为随机噪声码序列的自相关函数，用符号R(t)表示：
(4-1)
在实际应用中，可通过自相关函数R(t)的取值判断两个随机噪声码序列的相关性。

根据自相关函数R(t)的取值，即可确定两个随机噪声码序列是否已经“相关”，或者说，两个码序列的对应码元是否已完全“对齐”。

假设GPS 卫星发射一个随机序列u(t)，而GPS 信号接收机在收到信号的同时复制出结构与u(t)完全一样的随机序列u(t) ，由于信号传播延迟的影响，被接收的随机序列u(t)与u(t)之间产生了平移，即对应码元已错开，因而R(t )≈0。

若通过一个时间延迟器来调整，使它们的码元相互完全对齐，即有R(t)=1，那么就可以从GPS 接收机的时间延迟器中，测出卫星信号到达用户接收机的准确传播时间，再乘以光速便可确定卫星至观测站的距离。

所以，随机噪声码序列良好的自相关特性为GPS 测距奠定了基础。

M 序列是线行反馈移位寄存器产生的周期最长的序列，是多级移位寄存器或其他延迟元件通过线性反馈产生的最长的码序列。

下面以一个四级反馈移位寄存器组成的m 序列为例来说明，如图1.4所示。

在时钟脉冲的驱动下，每个存储单元的容，都按次序由上一级单元转移到下一单元，而最后一个存储单元的容便为输出。

同时，其中某两个存储单元，例如单元3和单元4的容进行模二相加后，再反馈输入给第一个存储单元。

u u u u
D S D S t R +-=)(
图1.4 四级反馈移位寄存器示意图
所谓的模二相加，是二进制数的一种加法运算，常用符号表示，其运算规则如下：
当移位寄存器开始工作时，置“1”脉冲使各级存储单元处于全“1”状态，此后在时钟脉冲的驱动下，移位寄存器经历15种不同的状态，然后再返回到“1”状态，从而完成一个周期（见表1.2）。

在四级反馈移位寄存器经历了上述15种状态的同时，其最末级存储单元输出了一个具有15个码元，且周期为15t
u
的二进制数码序列，称为m序列。

t u 表示时钟脉冲的时间间隔，即码元的宽度。

0,1
1,1
1
0，
1
1=
⊕
=
⊕
=
⊕
=
⊕
表1.2 四级反馈移位寄存器状态序列
状态编号 各级状态 模二加反馈 末级输出的二进制
数 ④ ③ ② ① ③ ⊕ ④ 1 1 1 1 1 0 1 2 1 1 1 0 0 1 3 1 1 0 0 0 1 4 1 0 0 0 1 1 5 0 0 0 1 0 0 6 0 0 1 0 0 0 7 0 1 0 0 1 0 8 1 0 0 1 1 1 9 0 0 1 1 0 0 10 0 1 1 0 1 0 11 1 1 0 1 0 1 12 1 0 1 0 1 1 13 0 1 0 1 1 0 14 1 0 1 1 1 1 15 0 1 1 1 1 0
由此可见，四级反馈移位寄存器所产生的m 序列，其一个周期可能包含的最大码元个数恰好等于24-1个。

因此，一般来说，一个r 级移位寄存器所产生的m 序列，在一个周期其码元的最大个数
（1.2） 与此相对应，这时m 序列的最大周期为：
（1.3） 式中，N u 也称为码长。

由于移位寄存器不容许出现全“0”状态，因此2r -1码元中，“1”的个数总比“0”的个数多一个。

这样，当两个周期一样的m 序列其对应码元完全对齐时，自相关系数R(t)=1，而在其他情况则有 （1.4）
当r 足够大时，就有R(t )≈0。

所以，伪随机噪声码与随机噪声码一样，具有良好的自相关性，而且是一种结构确定、可以复制的周期性序列。

GPS 信号接收机就是利用这一特征使所接收的伪随机噪声码和机产生的伪随机噪声码达到对齐同步，进而捕获和识别来自不同GPS 卫星的伪随机噪声序列。

由于受GPS 卫星至用户GPS 接收机的路径信号传播延迟的影响，被接收的伪随机码和复制的伪随机码之间产生了平移；如果通过一个时间延迟器来对复制的伪随机码进行移动，使两者的相关函数值为1，则可以从时间延迟器中测出对齐码元所用的时间，从而可以较准确地确定由卫星到接收机的距离。

由此可知，伪随机序列的良好的自相关特性，对于利用GPS 卫星的测距码进行精密测距是非常重要的。

12-=r
u N (21)r
u u u u
T t N t =-=121
1)(--=-=r u
N t R
m 序列有下列特性: （1）均衡性：在一个周期中，“1”与“0”的数目基本相等，“1”比“0”的数目多一个。

它不允许存在全“0”状态。

（2）游程分布：在序列中，一样的码元连在一起称为一个游程。

一般来说，长度为1的游程占总数的1/2，长度为2的游程占总数的1/4，依此类推。

连“1”的游程和连“0”的游程各占一半。

（3）移位相加特性：一个m 序列m P 与其经过任意次延迟移位产生的另一个序列m r 模二相加，得到的m S 仍是m 序列，即
S
r P m m m =⊕ （1.5）
（4）伪噪声特性：如果对随机噪声取样，并将每次取样按次序排成序列，可以发现其功率谱为正态分布。

由此形成的随机码具有噪声码的特性。

m 序列在出现概率、游程分布和自相关函数等特性上与随机噪声码十分相似。

正因为如此，我们将m 序列称为伪随机码，或人工复制出来的噪声码。

1.2.2 C/A 码
C/A 码（Coarse Acquisition Code ）是Gold 码，用于粗测距和捕获GPS 卫星信号。

它是由两个10级反馈移位寄存器组合产生的，其序列长度为1023（基数码）。

因为C/A 码的基码速率是1.023MHz ，因此伪随机序列的重复周期是1023/1.023*106或1ms 。

图1.5描述了GPS C/A 码发生器的结构方案。

图1.5 C/A 码发生器
两个移位寄存器于每星期日子夜零时，在置“1”脉冲作用下处于全“1”状态，同时在频率为f 1=f 0/10=1.023MHz 时钟脉冲驱动下，两个移位寄存器分别产生码长为N=210-1=1023、周期为1ms 的两个m 序列G 1(t)和G 2(t)。

这时G 2(t)序列的输出不是在该移位寄存器的最后一个存储单元，而是选择其中两个存储单元进行二进制相加后输出，由此得到一个与G 2(t)平移等价的m 序列G 21（即与延时等价）。

再将其与G 1(t)进行模二相加，将可能产生1023种不同结构的C/A 码。

C/A 码不是简单的m 序列，而是由两个具有一样码长与数码率，但结构不同的m 序列相乘所得到的组合码，称为戈尔德（Gold ）序列。

（1.6） 采用不同的it 0值，可能产生1023个G 2(t)，再加上G 1(t)和G 2(t)本身，共可能产生1025种结构不同的C/A 码供选用。

这些C/A 码具有一样的码长N=210-1=1023bit ，一样的码元宽t u =1/f 1=0.98 μs（相当于293.1 m ）和一样的周期T u =Nt u =1 ms 。

从这些G(t)码中选择32个码以PRN1，PRN2，…，PRN32命名各种GPS 卫星。

由于C/A 码长很短，只有1023比特，易于捕获。

在GPS 定位中，为了捕获C/A 码，以测定卫星信号的传播延时，通常需要对C/A 码逐个进行搜索。

若以50个码元每秒的速度搜索，对于只有1023个码元的C/A 码，搜索时间只要20.5 s 。

通过C/A 码捕获卫星后，即可获得导航电文，通过导航电文提供的信息，便可以很容易地捕获GPS 的P 码。

所以，C/A 码除了作为粗测码外，还可作为GPS 卫星信号P 码的捕获码。

120/()()()C A t G t G t it =⋅+
C/A 码的码元宽度较大。

假设两个序列的码元对齐误差为码宽的1/10～1/100，则此时相应的测距误差为29.3～2.93 m 。

随着现代科学技术的发展，使得测距分辨力大大提高。

一般最简单的导航接收机的伪距测量分辨力可达0.1m 。

C/A 码的码长、码元宽度、周期和数码率分别为：码长N u =210-1=1023bit ；码元宽度t u ≈0.977 52μs，相应的长度为293.1m ；周期T u =N u ·t u =1 ms ；数码率为1.023 Mb/s 。

不同的GPS 卫星所使用的C/A 码的上述四项指标一样，但编码规则不同，这样既便于复制又便于区分。

C/A 码具有以下特点：
（1）由于C/A 码的码长较短，易于捕获，而通过捕获C/A 码所得到的信息，又可以方便地捕获P 码，因此，通常称C/A 码为捕获码。

在GPS 导航和定位中，为了捕获C/A 码以测定卫星信号传播的时间延迟，通常对C/A 码进行逐个搜索，而C/A 码总共只有1023个码元，若以50码元每秒的速度搜索，仅需20.5 s 便可完成。

（2）C/A 码的码元宽度较大。

若两个序列的码元相关误差为码元宽度的1/10～1/100, 则此时所对应的测距误差可达29.3～2.9 m 。

由于其精度较低，所以称C/A 码为粗捕获码。

1.2.3 P 码
图1.6给出了在GPS 中用于实现码分多址技术的产生直接序列PRN 码的高层方框图。

每个合成的PRN 码由前面两个另外的码发生器导出。

在每种情况下，第二个码发生器的输出在其与第一个的输出由异或电路合并之前要相对于第一个进行延时，延时的量是可变的。

卫星的PRN 码与延时的量是相关联的。

在P 码的情况下，延时的基码整数与PRN 码一样。

对于C/A 码来说，对每颗卫星都有特别的延时。

表1.3列出了这些延时。

C/A 码延时可以由一种简单而有效的技术来实现，这种技术不需要使用延时寄存器。

图1.6 GPS 码发生器
表1.3 C/A 码和P 码的码相位分配和码序列初始
段
P 码（Precise Code ）是卫星
的精测码，码
速率为10.23
MHz ，它是由两
组各有两个12
级反馈移位寄存器结合产生的，其基本原理与C/A 码相似，但其线路设计细节远比C/A 码复杂，且严格。

这四个移位寄存器称为X1A ，X1B ，X2A ，X2B 。

图1.7示出了这种移位
卫星PRN 号码 C/A 码抽头 选择 C/A 码延时 /基码 P 码延时 /基码 前10个C/A 基码（八进制） 前12个C/A 基码（八进制） 1 2⊙6 5 1 1440 4444 2 3⊙7 6 2 1620 4000 3 4⊙8 7 3 1710 4222 4 5⊙9 8 4 1744 4333 5 1⊙9 17 5 1133 4377 6 2⊙10 18 6 1455 4355 7 1⊙8 139 7 1131 4344 8 2⊙9 140 8 1454 4340
寄存器方案的详细方框图。

图中未包括用于设定或读出移位寄存器和计数器相位状态所必须的控制。

图1.7 P码发生器
12级反馈移位寄存器产生的m序列的码元总数为212-1=4095，采用截短法将两个12级m 序列截短为一周期中码元数互为素数的截短码。

所有4个反馈移位寄存器的自然周期都是按如下方式截短的：X1A和X2A在4092个基码之后复位，去掉了它们自然的4095个基码序列的最后3个基码，寄存器X1B和X2B再4093个基码后复位，去掉了它们自然的4095个基码的最后2个基码，.这就导致X1B序列的相位相对于X1A学列在每一个X1A 寄存器循环都滞后一个基码。

结果在X1A和X1B之间有相对相位移动。

在X2A和X2B之间也发生类似的相位移动。

在GPS星期的起始点所有移位寄存器均置于初始状态。

在每个X1A时元的末尾，X1A移位寄存器也复位到其初始状态。

在每个X1B时元的末尾，X1B 移位寄存器复位到其初始状态。

在每个X2A时元的末尾，X2A移位寄存器复位到其初始状态。

在每个X2B时元的末尾，X2B移位寄存器复位到其初始状态。

A和B移位寄存器
的输出（第12级）用异或电路合并起来，由X1A⊕X1B导出X1序列，由X2A⊕X2B导出X2序列。

X2序列被延迟i个基码（相应与SVi）以形成X2i。

SVi的P码是Pi=X1⊕X2i。

在X1A 产生了4092 个码片的时候, X1A 完成了一个周期, 并产生了一个脉冲SETX1AEPOCH 使X1A 复位, X1B 则需要产生4093 个码片时才产生脉冲SETX1BEPOCH. 当X1B输出了3749 个周期的时候, 它停止移位, 等待X1A 输出了3750个周期后对它发送重新启动的信号. 在X1A 完成了3750个周期后, 它又产生
了一个叫做X1EPOCH 的脉冲. 所以, 在每一个X1EPOCH 周期中, X1B 要停止并等待X1A 4092- 3749= 343个码片的时间. X2A 和X2B 的工作过程与X1A 和X1B 的工作过程类似. 同时X2A每经历3750个周期后都要比X1A 滞后37个时钟周期. X1B, X2A, X2B这3个线性反馈移位寄存器都要参照X1A。

表1.4 GPS码发生器多项式和初始状态
在将X1A和X1B通过模二相加或波形相乘，得到周期为4092×4093的长周期码。

再对乘积码截短，截出周期为1.5s、码元数N1=15.345×106的X1。

同样的方法，在另外一组中，两个12级移位寄存器产生X2，只是X2码比X1码周期略长一些。

总之，两个子码N
1
和N
2
均是由两个12级移位寄存器产生的截短码。

其码速率均为
10.23 Mb/s，N
2
的码长比N
1
的多37个码元。

两个子码的码长分别为：
（4.7）
（4.8）
因此P码的码元数为：
（4.9）
相应的周期为：
（4.10）
66
1
10.2310 1.515.34510
N=⨯⨯=⨯
6
2
15.3451037()
N b
=⨯+
)
(
10
35
.214
2
1
b
N
N
N⨯
=
⋅
=
267()38(
P
P
N
T
f
=≈≈
天星期)
在乘积PN
1(t)·PN
2
(t+n
i
τ)，n
i
可取0，1，2，…，36。

这样可得到37种P码。

在
实际应用中，P码采用7天的周期，即在PN
1(t)·PN
2
(t+n
i
τ)中截取一段周期为7天的
P码，并规定每星期六午夜零点使P码置“1”状态作为起始点。

在这37个P码中，32个供GPS卫星使用，5个供地面站使用。

这样，每颗卫星所使用的P码便具有不同的结构，但码长和周期一样。

P码的特征是：码长N
u =2.35×1014bit；码元宽度t
u
≈0.097 752μs，相应长度为29.3m；
周期T
u =N
u
·t
u
≈267天；数码率为10.23 Mb/s。

因为P码的码长约为6.19×1012bit，所以如果仍采用搜索C/A码的办法来捕获P码，即逐个码元依次进行搜索，当搜索的速度仍为50码元每秒时，那将是无法实现的（约需14×105天）。

因此，一般都是先捕获C/A码，然后根据导航电文中给出的有关信息，便可捕获P码。

另外，由于P码的码元宽度为C/A码的1/10，这时若取码元的相关精度仍为码元宽度的1/10～1/100，则由此引起的距离误差约为2.93～0.29 m，仅为C/A码的1/10。

所以P 码可用于较精密的导航和定位，称为精码。

目前美国政府对P码，不提供民用，因此一般GPS用户实际只能接收到C/A码。

1.3小结
全球定位导航系统随着技术的发展和用户对系统功能的需求改变，也在不断进行现代化改进，本章主要介绍GPS的总体结构，整体分析GPS信号的组成和调制所用频率，继而全面分析了M序列、C/A码、P码的结构，详细解释了信号发生器的工作流程，为仿真的实现奠定基础。

2MATLAB软件。