基于MATLAB的GPS信号的仿真.docx
- 文档编号:26262232
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:73
- 大小:1.18MB
基于MATLAB的GPS信号的仿真.docx
《基于MATLAB的GPS信号的仿真.docx》由会员分享,可在线阅读,更多相关《基于MATLAB的GPS信号的仿真.docx(73页珍藏版)》请在冰豆网上搜索。
基于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
Spreadspectrumcommunicationisacommunicationstechnologydevelopedrapidlyinrecentyears.Inearlystudiesthemainpurposeofthistechnologyistoimprovethemilitarycommunicationsconfidentialandanti-jammingperformance,thereforethedevelopmentandapplicationofthistechnologyisalwaysinsecretstate.SpreadspectrumtechnologyinthemostsuccessfulmilitaryapplicationexamplesaretheUnitedStatesandRussiacouldtheglobalpositioningsystem(GPSandGLONASS)forrepresentative;IncivilGPSandGLONASSalsohavebeenwidelyused,whichfoundationofsystemisthespreadspectrumtechnology.
Globalpositioningsystem(GPS)isusedtoprovidereal-timenavigationandpositioningservicesforglobalcivilandmilitaryaircraft,ships,personnel,vehiclesandsoon.GPSsystemadoptsthetypicalCDMAsystem,whichkindofspreadspectrummodulationsignalshavelowinterceptprobabilitycharacteristic.Thissystemmainlyusedthedirectsequencespreadspectrummodulationtechnology,usingthePRNcodeincludingC/Acode,PcodesandYcodes.
Thisarticletellsthedirectsequencespreadspectrumcommunicationtechnologyappliedinglobalpositioningsystem(GPS).Thearticlemainlyintroducesthepncodespreadspectrumcommunicationsimulation,brieflydiscussingMsequenceandpseudorandomnoisecode(PyardsandC/Ayards)anditsproduceanduseMATLAB7.0simulateMseries,PyardsandC/Ayardsofencodingprocessandthesimulationresults,introducingpncodestraightspread-spectrumtechnologyknowledge,especiallyPyards.
Key:
GPS;DS-SS;Pncodesimulation
引言
全球卫星定位系统,简称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,StandardPositioningService)和军规的标准精确定位服务(PPS,PrecisePositioningService)。
GPS目前处于良好的运行状态,并满足20世纪60年代所提出的最佳定位系统标准。
这个系统向有适当接收设备的全球范围内的用户提供精确、连续的三维位置和速度信息。
GPS也向全球广播世界协调时(UTC)。
组成卫星星座的24颗卫星被安排在6个轨道平面上,即每个平面上4颗。
这样的卫星星座配置确定了卫星某一时刻在轨道中的位置,而为hi与地平线以上的卫星数会随着时间和地点的不同而不同,最少可见到颗,最多可见到11颗,因此保证了在开放的天空下,地球上和近地空间任一点在任何时刻均可以同时观测到至少4颗GPS卫星,为准确定位提供可能。
GPS系统采用CDMA(CodeDivisionMultipleAccess,码分多址)技术将在轨的24颗卫星分开,每颗卫星使用同样的调制方式在相同的载波频率上调制各自唯一的伪随机码和数据信息,然后使用星载卫星天线发射信号。
卫星轨道分布如图1.1所示:
图1.1卫星轨道分布
1.2GPS信号构成
GPS卫星向广大用户发送的用于导航定位的信号,是一种调制波,但有别于常用的无线电广播电台发送的调频调幅信号,它是利用伪随机噪声码传送导航电文的调相信号。
GPS卫星信号是目前常用的两种违心导航定位信号之一,它包含有三种信号分量,即载波(L1和L2)、测距码(C/A码和P码)和数据码(D码,亦称基带信号或导航电文)。
而这所有这些信号分量都是在同一个基本频率f0=10.23MHz的控制下产生的。
GPS卫星信号示意图如图1.2所示:
图1.2GPS卫星信号频率构成
GPS卫星发送的GPS卫星信号采用L波段的两种载频作载波,分别被称作L1的主频率和L2的次频率。
这些载波频率由扩频码(每一颗卫星均有专门的伪随机序列)和导航电文所调制。
所有卫星均在这两个相同的载波频率上发射,但由于伪随机码调制不同,因此无明显的相互干扰。
GPS使用L频段的两种载频为(其中f0是卫星信号发生器的基准频率):
L1载波:
fL1=154×f0=1575.42MHz,波长λ1=19.032cm;
L2载波:
fL2=120×f0=1227.6MHz,波长λ2=24.42cm。
选择L波段的好处是:
(1)减少拥挤,避免“撞车”。
目前L波段的频率占用率低于其他波段,与其他工作频率不易发生“撞车”现象,有利于全球性的导航定位测量。
(2)适应扩频,传送宽带信号。
GPS卫星采用扩频技术发送卫星导航电文,其频带高达20MHz左右,在占用率较低的L波段上,易于传送扩频后的宽带信号。
在载波L1上调制有C/A码、P码的数据码,而在载波L2上,只有调制有P码的数据码。
在无线通信技术中,为了有效地传播信息,一般均将频率较低的信号加载到频率较高的载波上,而这时频率较低的信号成为调制信号。
GPS信号是一种调制波,它不仅采用L波段的载波,而且采用扩频技术传送卫星导航电文。
所谓“扩频”,是将原来打算发送的几十比特速率的电文变换成发送几兆甚至几十兆比特速率的由电文和伪随机噪声码组成的组合码。
采用扩频技术时,若信号功率仅为噪声功率的1/10,那么信号将深深地淹没在噪声之中而不易被他人捕获,从而使得信号具有极强的保密性。
GPS信号的调制波,是卫星导航电文和伪随机噪声码(PseudoRandomNoiseCode,简称PRN码,或称伪噪声码)的组合码。
卫星导航电文是一种不归零二进制码组成的编码脉冲串,称之为数据码,记作D(t),其码率为50b/s。
对于距离地面20000km之遥的GPS卫星,扩频技术能有效地将很低码率的导航电文发送给用户。
其方法是用很低码率的数据码作二级调制(扩频)。
第一级,用50Hz的D码调制一个伪噪声码,例如调制一个被叫做P码的伪噪声码,它的码率高达10.23MHz。
D码调制P码的结果,便形成了一个组合码——P(t)D(t),使得D码信号的频带宽度从50Hz扩展到10.23MHz,也就是说,GPS卫星从原来要发送50b/s的D码,转变为发送10230b/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.3GPS卫星信号的产生
需要注意的是,GPS信号虽然有几种分量(C/A易捕获码、P精确码和D导航数据码),但是它们均来源于一个公共的10.23MHz的基准频率(见图1-1)。
它们的频率不仅与基准频率有一定的比例关系,而且相互之间也存在一定的比例关系,详细如表4-1所示。
这既有利于GPS卫星发送信号,又便于广大用户接收和测量GPS信号。
从表4-1中可以看出,在D码的一个码元内,将有20460个C/A码码元,204600个P码码元,31508400个L1周期和24552000个L2周期。
表1.1GPS信号的频率关系
相关频率
基频F
载频fL1
载频fL2
基准频率F
10.23MHz
154F
120F
C/A码的码频fg
F/10
fL1/1540
fL2/1200
P码的码频fp
F
fL1/154
fL2/120
D码的码频fd
F/204600
fL2/31508400
fL2/24552000
1.2.1M序列
码是一种表达信息的二进制数及其组合,是一组二进制的数码序列。
例如,对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)相同的码元个数为Su,而码元相异的码元个数为Du,那么两者之差Su-Du与两者之和Su+Du(即码元总数)的比值,即定义为随机噪声码序列的自相关函数,用符号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的内容进行模二相加后,再反馈输入给第一个存储单元。
图1.4四级反馈移位寄存器示意图
所谓的模二相加,是二进制数的一种加法运算,常用符号表示,其运算规则如下:
当移位寄存器开始工作时,置“1”脉冲使各级存储单元处于全“1”状态,此后在时钟脉冲的驱动下,移位寄存器经历15种不同的状态,然后再返回到“1”状态,从而完成一个周期(见表1.2)。
在四级反馈移位寄存器经历了上述15种状态的同时,其最末级存储单元输出了一个具有15个码元,且周期为15tu的二进制数码序列,称为m序列。
tu表示时钟脉冲的时间间隔,即码元的宽度。
表2-4四级反馈移位寄存器状态序列
状态编号
各级状态
模二加反馈
末级输出的二进制数
④
③
②
①
③
④
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)
式中,Nu也称为码长。
由于移位寄存器不容许出现全“0”状态,因此2r-1码元中,“1”的个数总比“0”的个数多一个。
这样,当两个周期相同的m序列其对应码元完全对齐时,自相关系数R(t)=1,而在其他情况则有
(1.4)
当r足够大时,就有R(t)≈0。
所以,伪随机噪声码与随机噪声码一样,具有良好的自相关性,而且是一种结构确定、可以复制的周期性序列。
GPS信号接收机就是利用这一特征使所接收的伪随机噪声码和机内产生的伪随机噪声码达到对齐同步,进而捕获和识别来自不同GPS卫星的伪随机噪声序列。
由于受GPS卫星至用户GPS接收机的路径信号传播延迟的影响,被接收的伪随机码和复制的伪随机码之间产生了平移;如果通过一个时间延迟器来对复制的伪随机码进行移动,使两者的相关函数值为1,则可以从时间延迟器中测出对齐码元所用的时间,从而可以较准确地确定由卫星到接收机的距离。
由此可知,伪随机序列的良好的自相关特性,对于利用GPS卫星的测距码进行精密测距是非常重要的。
m序列有下列特性:
(1)均衡性:
在一个周期中,“1”与“0”的数目基本相等,“1”比“0”的数目多一个。
它不允许存在全“0”状态。
(2)游程分布:
在序列中,相同的码元连在一起称为一个游程。
一般来说,长度为1的游程占总数的1/2,长度为2的游程占总数的1/4,依此类推。
连“1”的游程和连“0”的游程各占一半。
(3)移位相加特性:
一个m序列mP与其经过任意次延迟移位产生的另一个序列mr模二相加,得到的mS仍是m序列,即
(1.5)
(4)伪噪声特性:
如果对随机噪声取样,并将每次取样按次序排成序列,可以发现其功率谱为正态分布。
由此形成的随机码具有噪声码的特性。
m序列在出现概率、游程分布和自相关函数等特性上与随机噪声码十分相似。
正因为如此,我们将m序列称为伪随机码,或人工复制出来的噪声码。
1.2.2C/A码
C/A码(CoarseAcquisitionCode)是Gold码,用于粗测距和捕获GPS卫星信号。
它是由两个10级反馈移位寄存器组合产生的,其序列长度为1023(基数码)。
因为C/A码的基码速率是1.023MHz,因此伪随机序列的重复周期是1023/1.023*106或1ms。
图1.5描述了GPSC/A码发生器的结构方案。
图1.5C/A码发生器
两个移位寄存器于每星期日子夜零时,在置“1”脉冲作用下处于全“1”状态,同时在频率为f1=f0/10=1.023MHz时钟脉冲驱动下,两个移位寄存器分别产生码长为N=210-1=1023、周期为1ms的两个m序列G1(t)和G2(t)。
这时G2(t)序列的输出不是在该移位寄存器的最后一个存储单元,而是选择其中两个存储单元进行二进制相加后输出,由此得到一个与G2(t)平移等价的m序列G21(即与延时等价)。
再将其与G1(t)进行模二相加,将可能产生1023种不同结构的C/A码。
C/A码不是简单的m序列,而是由两个具有相同码长及数码率,但结构不同的m序列相乘所得到的组合码,称为戈尔德(Gold)序列。
(1.6)
采用不同的it0值,可能产生1023个G2(t),再加上G1(t)和G2(t)本身,共可能产生1025种结构不同的C/A码供选用。
这些C/A码具有相同的码长N=210-1=1023bit,相同的码元宽tu=1/f1=0.98μs(相当于293.1m)和相同的周期Tu=Ntu=1ms。
从这些G(t)码中选择32个码以PRN1,PRN2,…,PRN32命名各种GPS卫星。
由于C/A码长很短,只有1023比特,易于捕获。
在GPS定位中,为了捕获C/A码,以测定卫星信号的传播延时,通常需要对C/A码逐个进行搜索。
若以50个码元每秒的速度搜索,对于只有1023个码元的C/A码,搜索时间只要20.5s。
通过C/A码捕获卫星后,即可获得导航电文,通过导航电文提供的信息,便可以很容易地捕获GPS的P码。
所以,C/A码除了作为粗测码外,还可作为GPS卫星信号P码的捕获码。
C/A码的码元宽度较大。
假设两个序列的码元对齐误差为码宽的1/10~1/100,则此时相应的测距误差为29.3~2.93m。
随着现代科学技术的发展,使得测距分辨力大大提高。
一般最简单的导航接收机的伪距测量分辨力可达0.1m。
C/A码的码长、码元宽度、周期和数码率分别为:
码长Nu=210-1=1023bit;码元宽度tu≈0.97752μs,相应的长度为293.1m;周期Tu=Nu·tu=1ms;数码率为1.023Mb/s。
不同的GPS卫星所使用的C/A码的上述四项指标相同,但编码规则不同,这样既便于复制又便于区分。
C/A码具有以下特点:
(1)由于C/A码的码长较短,易于捕获,而通过捕获C/A码所得到的信息,又可以方便地捕获P码,因此,通常称C/A码为捕获码。
在GPS导航和定位中,为了捕获C/A码以测定卫星信号传播的时间延迟,通常对C/A码进行逐个搜索,而C/A码总共只有1023个码元,若以50码元每秒的速度搜索,仅需20.5s便可完成。
(2)C/A码的码元宽度较大。
若两个序列的码元相关误差为码元宽度的1/10~1/100,则此时所对应的测距误差可达29.3~2.9m。
由于其精度较低,所以称C/A码为粗捕获码。
1.2.3P码
图1.6给出了在GPS中用于实现码分多址技术的产生直接序列PRN码的高层方框图。
每个合成的PRN码由前面两个另外的码发生器导出。
在每种情况下,第二个码发生器的输出在其与第一个的输出由异或电路合并之前要相对于第一个进行延时
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 MATLAB GPS 信号 仿真
![提示](https://static.bdocx.com/images/bang_tan.gif)