中频数据模拟软件开发报告.docx

《中频数据模拟软件开发报告.docx》由会员分享，可在线阅读，更多相关《中频数据模拟软件开发报告.docx（23页珍藏版）》请在课桌文档上搜索。

1、中频数据模拟软件开发报告作者，李子月2012-7-10第1章绪论1、1软件编写背景、意义卫星导航和惯性导航是现在运用最为广泛的两种导航模式，两种导航模式具有性能互补特性，GPS/1NS组合导航是种比较志向的导航模式。松组合、紧组合中两种系统相互协助不够亲密，并且理论上存在确定的误差，系统性能不够稳定，1999年被提出的GPS/INS深组合导航是现在国际上的探讨热点。GPS/INS深组合导航理论上不同于松、紧组合模式，在数据同步、数据处理、误差分析、状态估计方面都存在难点，硬件实现特别困难。现在国外从事深组合探讨的除了些先进探讨所以与军工单位搭建了硬件测试平台，大部分都是处于半实物仿真阶段：国内

2、探讨由于技术上受到种种限制，科研主要以仿真为主，其中，国防科技高校、哈工大、北航等导航基础较强的高校有的试验室搭建了简易的半实物仿真平分。由于深组合导航算法困难，须要用到I、Q相关数据和惯导数据进行数据融合，并且数据同步难以实现，所以对数据要求也很高，真实数据有时难以满意软件须要，特殊是在理论探讨的初期，仿真数据的优势远大于真实采集数据。除此之外仿真中频数据还有其他方面的优势：（1）可以仿真现实中难以获得的运动轨迹数据。以高动态导航为探讨背景的科研1：作者很难获得高速、高动态载体的卫星数据，软件仿真中频数据可以仿真随意飞行轨迹参数卜的数据。（2）可以随意设定信噪比，以与可以捕获到的卫星数目，这

3、就便利了弱信号跟踪、抗干扰等技术的探讨。（3）仿真的中频数据由干全部参数都是人为设定，科研人员就可以清晰的知道试验的理论值，便利了误差分析。不但在GPS/INS深组合探讨中仿真中频数据具有很大的优势，在许多领域仿真中频数据都是个很好的选择。国外很早就有了GPS中频数据仿真的ma1.1.ab工具箱，但是价格特别昂贵，国内的许多试验室都开发了自己的仿真软件。1、2软件的功能和系统概况本软件是全部采纳matIab语言编写，基本功能就是仿真出可供软件接收机识别并处理、解算出位置、速度等信息的数据。基于木软件有些参数和数据可以随意选取，具体如下：（1）仿真中频数据必需采纳R1.NEX格式的导航电文文件，

4、其中包含r固定时间段之内的数据解算卫星信息，这种数据文件时可以随意选取的。(/)表示数据码，是载波1.1.的角频率，叫是载波1.1.信号的初始相位。去除掉P(Y)码，本软件中的信号表达式为：SJ“(，)=AGa)。(小in3/+J+,2.2)式中巧表示信号的噪声。从公式22中可以看出，仿真的中频信号包含理论上包含伪码、导航电文和载波信号，这三种信号调制在起再加上噪声就构成了GPS信号。伪码用来捕获和计算伪距，导航电文中包含/定位所须要的卫星轨道参数以与信号中的时间参数，载波中的多普勒频移可以计兑速度，载波相位可以用来精确定位。三者在调制过程中的关系如下：19cm我波1.1.：=1575.42M

5、HZC/A码：1.023MCPS1.J、3码片1S40周我泡ZZ数据码：50bps、每比特20周期C/A码/，图2.2载波、C/A码、导航电文调制2.1.2C/A码信号C/A码也就是PRN码，是种由0、1数字组成的伪随机码。GPS官方文档ICDGPS-200供应TC码的产生原理，如图2.3所示:图2.3C/A码生成原理C/A内码片速率为1.023MHz,每个C/A科长度为1023码片,对应的C/A码周期为Ims,码片宽度大约为伪码具有良好的自相关特性，自相关函数如下：%(r)=yCBga+r)dt,rG(-T,T)(2.3)由公式2.3可得，只有当同一个C/A码并且码相位误差很小的时候勺才能得

6、到个较大的峰值。GPS接收机每个通道中的C/A码与输入信号做相关运算，利用C/A码良好的自相关性就可以检测到输入信号中所包含哪几颗卫星。C码不但可以用来捕获卫星，还可以用来测量载体与卫星之间的距离，也就是伪距。通过捕获和跟踪可以精确的知道输入信号的码相位，然后计算出信号放射时间，接收时fE可以从本地时钟获得，做差就可以得到伪距。每一个码片对应的距离大约为300米，通过C/A码求得伪距误差为几十米左右，假如想获得更高精度的伪距，可以借助于载波测距，精度可达厘米级。导航电文导航电文数据率为50bps,对应的每个数据位长度为20ms,调制在GPS信号中导航电文是0、1组成的二进制数，这些二进制数代表

7、了计算卫星轨道信息所须要的参数。导航电文每一帧包含5个子帧，每一子帧又包含10个字，每一个字包含30数据位，每一数据位占据20ms。导航电文结构图如下：GPS导航电文结构从图中可以看出，每一个字0.6s,每一子帧6s,30s传送完一帧。前三子帧包含当前时刻全部的卫星星历信息，历史信息保存的4、5子帧，25帧才组成完整的历书，所以完整的卫星星历须要25*6=12.5min才能发送完毕。由于前三子帧已包含计算所须要的星历数据，18s的卫星数据就可以进行定位。2.2GPS信号误差卫星信号由卫星放射到接收机接收，这一过程须要信号传输2万多千米，这一过程会给信号带来许多误差。首先，卫星放射的卫星信号中本

8、身就包含H星星历误差以与卫星时钟带来的误差：其次H星信号传输须要穿越大气层，这一过程带来的误差可以归结为电离层延时和对流层延时；另外，信号传输到地面被接收机接收还会带来多路径效应、电磁干扰以与接收机噪声等误差。总之，接收机接收到的信号包含了许多噪声，假如要精确的定位计算就必须要对这些噪声模型有一个清晰地相识，下面将分别对这些误差进行介绍。时钟误差由于卫星上的原子钟时间偏移和频率漂移，造成卫星时钟与标准的GPS系统时间相比，t时刻北星时钟所包含的误差可以表示如下：&=%式fJ(2.4)式中a、。匚为导航电文中供应的误差参数，J为导航电文供应的时间参数,TimeofC1.OCk的缩写，表示时钟时刻

9、。除此之外，卫星时钟误差还包含相对论校正量和群延迟校正，总的时钟误差为：St=At+Mr-Tgd(2.(5)心由导航电文供应，相对论校正量M可以表示为：/,=Fe(7sinE1(2.(6)其中，e,为卫星的轨道偏心率，见表示轨道长度半径，耳为偏近点角，这三个参数都可以从导航电文获得。常熟产的值为：F=二=-4.442807633X10“】s(2.(7)和C为常数，可以查表获得。卫星星历误差通过地面监控部分计算出的卫星轨道参数用来描述里星的运行轨道，这些参数理论上是正确的，但是卫星在运行过程中受到的作用力特别困难，无法精确建模。通过星历参数计算出的卫星模型与真实的卫星运行轨道之间确定存在着误差，

10、但是这个误差没有精确的模型，由于地球转动特别缓慢，以与信号传输时间特别短，星历误差对定位造成的膨响很小。电离层和对流层误差电离层是高度约为70、100km出的大气层，电离层中充溢了电离子和电子，当电磁波穿过时波的传输速度和方向会发生变更，电离层对电磁波的影响与其频率有关，当GPS信号穿过电离层时电离层会对伪码和载波造成相反的影响。GPS信号穿过电离层伪码的传输速率会变慢，而载波的速率却会变快，二者变更的量相等，方向相反。多频接收机中电离层误差可以通过不同频率波之间的延时大小不同精确估计出延时，单频接收机中可以依据误差模型进行计算。对流层位于离地面约为IOkm的高度，对流层中包含了大气层中99%

11、的质量，其中氧气、氮气、水蒸气是造成对流层延时的主要缘由，与电离层相比，对流层对信号造成的误差要小得多，由于对流层和气象有关，实际中不便利实时获得精确的气象资料，一般的矫正都会采纳简化的对流层模型，而这种模型有多种形式的近似。这里列举种：T=247Sine+0.0121(2.(8)2. 2.4其他误差卫星信号传输到地面，经过四周的建筑物反射到达天线，这种经过反射的信号会对正确的信号造成影响，对给定位结果带来误差。另外，天线、放大器以与各部分的电子器件热噪声、信号量化误差、定位算法、软件接收机中的数值精度误差等都会给定位带来误差。这些误差都比较小，在实际结算过程中可以归结为一项进行估计。第三章G

12、PS中频数据仿真2.1 中频信号仿真模型GPS卫星放射的1.1.载波信号频率为1575.42MHz的高频信号，难以对其进行仿真，接收机的射频前端接收到射频信号之后仃先对其进行下变频，使信号载波频率变换到中频，此时的卫星数据成为中频数据。本软件仿真的对象即为GPS中须数据。考虑全部的误差因素，仿真中领信号的表达式为：S=尤必QHr三1.c,tp-&tm,-+33v-&/-)cos(.+wrfXr-7；,-以“,_+iv-tup)+a1.,+n1.(3.1)其中i表示可视卫星编号，N表示可见星的总数I，户表示信号的功率，。,表示导航电文，C,表示伪码，表示接收机时间，表示传输时间延时，血出S、4心

13、应分别代表电离层、对流层、卫星时钟漂移、多径效应所带来的误差，/表示中频信号角频率，?表示多普勒频移角频率，,表示信号噪声。本仿真软件的应用背景是高动态下的导航与定位，重点不在于电离层、对流层、多径效应等误差因素的影响，为了是问题更有利于分析算法的有效性，也为了软件仿真更有效率，本软件只考虑了卫星时钟误差带来的影响。简化之后的中频信号表达式如下：s”=回一乙+/e,-乙+&)cos3*+4)(/-7；+怎，)+仰+M1(3.2)在仿真过程中设置电离层、对流层、名径效应的影响为0,噪声设置为白噪声。在这个表达式中高动态体现在心上，表示的是多普勒频移用频率，载体与卫星的相对运动速度越快，?数值越大

14、，载体的机动性越强,例变更越快。当Q数值很大并且变更很快时，载波的频率就会变更的很快，造成跟踪环路无法锁定卫星信号，导致卫星信号失锁，无法导航定位。对应于3.2式中的中频信号模型，本软件的仿真框架图如下：图3.1中频数据仿真框架图对应于上式3.2中描述的中频信号模型，下面对中频信号的仿真过程做具体介绍。C/A码的生成依据2.1中介绍方法生成即可；导航电文可以通过网络卜载IGS精密星历数据文件获得，首先，读取文件中的参数，然后依据确定的规则把这些十进制的参数乘以一些系数之后转化为二进制数据，最终把这些二进制数据依据导航电文的格式编排，最终一部中还须要用到奇偶校验位的计算.；载波仿真中只要确定了中

15、频信号的频率和多普勒频移就可以干脆生成，中频频率是软件中设定的参数，这步中软件计算的重点就集中于多普勒频移的计算。2.2 导航电文仿真仿真导航电文须要先清晰其中包含的卫星参数，以与这些参数在导航定位中的作用。导航电文中的参数（软件接收机基础75页）前三子帧中包含了定位解算所须要的导航电文，本仿真试验中仿真了18秒数据，包含r前三子帧卫星轨道参数。第一子帧中的参数：WW（星期数）、用户位宜精度、卫星健康状态、/ODC卫（星时钟的数据龄期）、心，（群延迟估计）、卫星时钟修正因子。其次子帧中的参数：IODE（卫星星历数据的数据龄期）、Cn（卫星轨道半径的正弦调和修正值）、n（计算值的平均移动误差、M

16、O（参考时间的平近点角）、C“（纬度辎角的余弦调和修正项的幅度）、e,（椭圆离心率）、C,（纬度相角的正弦调和修正项的幅度）、R（卫星轨道长半轴平方根）、1.（时间星历的参考时间）。第三子帧中的参数：Ck（倾斜角的余弦调和修正项的幅度）、Q（卫星轨道的升交点经度）、C,（轨道倾角的正弦调和修正值）、（，（参考时刻的11星轨道倾角）、G（轨道半径的余弦调和修正值）、&（轨道近地点角距）、；（轨道倾角变更率）。以上GPS导航数据是导航解算必需的参数，仿真试验中这些参数可以通过网络下载RINEX格式的导航电文文件获得，下面对导航电文文件做简要介绍。RI三导航电文文件结构RINEX（ReceiverI

17、ndependentExchange）,共有三种格式观测文件、导航文件和气象文件（meteoro1.ogica1.）本软件所须要的是导琉文件，导航文件分为前一部分的文件头和后面的数据记录部分。文件头中的参数具体如F：标签第61-80列，从零起先计数描述格式RINEXVERSION/TYPE版本(2.10)-文件类型(N表示导航文件)9.2,1IX,A1.,19XPGM/RUNBY/DATE- 创建当前文件的程序名称- 创建当前文件的机构名称- 创建文件的时间A20,A20,A20*COMMENT注择语句A60*I0NA1.PHA历元的电离层参数A0-A32X,4D12.4IONBETA历元的电

18、离层参数B0-B32X,4D12.4DE1.TA-UTC:AO,A1.1T,W用于计算UTC的参数AO,A1.:多项式系数T:UTC参考忖间W:UTC参考周3X,2D19.12,2191.EAPSECONDS跳秒16ENDOFHEADER文件头部分的最终一行60X导航文件数据记录部分具体内容如下:观测记录描述格式PRN/EPOCH/SVC1.K-卫星编号12,-星历:Toc-时钟时间IX,12.2,年（2字节，依据须要可补IX,12,零）IX,12,nIX,12,日IX,12,时F5.1,分3D19.12秒-SV钟偏（SeCondS）SV钟漂（sec/SeC）-SV钟漂率（secsec2）BR

19、OADCASTORBIT-I-IODE数据分发历书3X,4DI9.12Crs（米）-De1.tan（弧度/秒）-MO（弧度）BROADCASTORBIT-2Cuc3X,4D19.12（弧度）-e第一偏心率-Cus（弧度）SqrI（八）（sqrt（m）BROADCASTORBIT-3-Toe历元时间（GPS周内的秒）Cic（弧度）-OMEGA（弧度）CIS（弧度）3X,4D19.12BROADCASTORBIT-4i（弧度）-Crc（米）-omega（弧度）-OMEGADOT（d1.&度/秒）3X,4D19.12BROADCASTORBIT-5-IDOT3X,4D19.12（弧度/秒）- 1.2

20、载波上的码- GPS周号（和TOE一起运用）- 1.2P码数据标记BROADCASTORBIT-6-SV精度（米）-SV健康度-TGD（秒）-IODC数据分发时钟3X,4D19.12BROADCASTORBIT-7-信息传递时间（GPS周秒，从HOW字的Z计数计算得到一间隔（小时）-备用-备用3X,4D19.12导航电文生成通过RINEX格式的导航电文文件可以获得导航电文中的参数，本软件的Reae1.eph.m函数从导航电文文件WUhn2890.IOn中读取星历数据（这个导航电文文件可以随意选取），读取的数据是十进制数，要转化为导航电文必需把这些数据转化为确定长度的二进制数。PerfeCtde

21、C2bin.m函数的作用是把导航电文中的数据转化为二进制位，区分于dec2bin.m,木函数的可以转化包括负数、小数在内的实数。GPS_naVdaIa_generaIe.m函数的作用是生成每一颗卫星的导航电文,本函数首先调用Rea(1.eph.m,读取导航电文文件中的参数到da1.a变量中。data中包含的导航电文中的参数经过PerfeCtdeC2bin.m函数转变为二进制位，这些二进制位依据导航电文格式编排导航电文中每一个字的前24位，这样生成的导航电文每帧为10*24的格式。每个字的后6位为奇偶校验位，奇偶校验位是利用前24计算生成，本软件中计算奇偶校验位的函数为Paritym0生成的6位

22、奇偶校验位附加在24位导航电文后面组成完整的导航电文。生成导航电文是生成中频数据中重要的步，之后还要计算中频数据中的其他参数，码相位、载波相位、多普勒频移、信号传输时间等这些参数在卜一节做具体介绍。3. 3中频数据中叁数计算卫星参数计算上文中已经计算出r导航电文中的计算卫星轨道所须要的参数，木软件中$21_口。5支2氏111函数的作用是计算指定卫星的卫星速度、位置。软件中首先调用Rea(1.eph.m函数，读取的卫星参数保存在变jnavda1.a中，CaICUIate_sata1.ite.m函数调用navdata,计算得出卫星的位置、速度，Ca1.CU1.aIe_Saudi1.em输出的变量经

23、过整理，变为所须要的位置、速度格式作为函数sajposition.m的输出量。卫星位置、速度的计算过程具体如下。卫星的轨道计算须要建立在对卫星轨道理论分析的基础之上，须要用到开普勒定律、开普勒方程式、坐标系之间的转换、真近点角和平近点角的含义。除了理论分析之外，计算过程还须要知道所需的16项星历数据参数。具体计算过程如下所示：(1)计算归化时间。因为卫星的星历数据都是相对于参考时刻而言的，所以须要将观测时刻做归一化处理1.F(3.(3)归化之后的时间(限制在-302400到302100之间。(2)计算卫星的平均角速度n。将导航电文中的扰动修正项带入%书(3.(4)式中为地心引力常数，为椭圆半长

24、轴，由星历数据供应，在利用星历数据中供应的修正项函得到平均角速度=%+r?(3.(5)(3)计算信号放射时刻的平近点角虬M=Mii+叫(3.(6)(4)计算卫星在信号放射时刻的偏近点角线E=/+etsinE1.(3.(7)式中e,为卫星椭圆轨道的偏心率，这是个超越方程，须要用迭代法求解，一般来说只需迭代十次即可。(5)计算信号放射时刻的真近点角匕。(6)计算卫星的地心向径r=a(-etcosE1)(3.(8)(7)计算信号放射时刻的升交点角距次。利用H星星历中供应的卫星轨道的近地点角距。，代入如下公式：件=VA+&(3.(9)(8)计算信号放射时刻摄动修正项为4,沌和西。升交点角距修正项：6氏

25、=Cm8s2+Cmsin2.(3.10)卫星地心向径修正项：西=Crrcos2+Cnsin2伪(3.11)星轨道倾角修正项：区=Ckcos2+C11sin21.(3.12)Cv、G、C、,、C-CI都来自于卫星星历数据。得到的这三个修正项用来校正内，和(：,=次+&Z(3.13)=r+r(3.14)(=+i(+(3.15)(9)计算卫星信号放射时刻在椭圆轨道平面直角坐标系中的位置。(3.16)yj=rtsinzt(3.17)(10)计算信号放射时刻的升交点赤经,计匏公式如卜.：C,=。0+(C-Qm-Cav(3.18)式中5、C为星历数据供应的卫星参数，Q,等于7.292115146105ra

26、ds,是地球自转角速度。(11)计算卫星在WGS-84地心地固直角坐标系中的位置，公式如下:.vt=;CoSQA-yiCOSiSini(3.19)-=)；sin-yjCOSi*cosi(3.20)ZA=Ksin(3.21)计算出的是在WGS-84地心地固坐标系中的坐标，假如想要得到最为广泛应用的大地坐标系，即经纬高，还须要转换到IVGS-84大地坐标系。计算出卫星的坐标之后对位置函数求导数即可得到计算卫星的速度。到这里已经计算出卫星的速度、位置，假如想要得到卫星信号仿真所须要的多普勒频移、码相位、载波相位还必需知道载体的轨迹，卜面介绍载体轨迹的仿真过程。3.3.2 栽体轨迹仿真载体的运动轨迹要

27、参照导航电文数据的采集地点，以保证在生成的运动轨迹环境卜的可视卫星与仿真的卫星吻合，假如采纳的导航电文文件包含32颗卫星，运载体的轨迹可随意生成，然后依据载体轨迹计算特定时间段之内的可视卫星，然后对可视卫星进行仿真。木仿真试验采纳的导航电文文件是接收机定点采集的，文件中只包含了可视卫星的星历数据，所以载体轨迹也应当参照星历采集地点的经纬度，然后依据轨迹生成算法进行计算。为了证明仿真的中频数据正确性本试验中采纳了三个轨迹，定点，直线，圆周运动。3.3.3 信号传输时间计第上文中已经计兑出了卫星的速度、位置和载体的速度位置，接卜.来要计算载波相位和伪码相位，相位的计算必须要利用载体与卫星之间的距离

28、信息，也就是伪距的计算。伪距是接收时刻载体的位置与信号放射时刻卫星的位置之间的距离，计算伪距可以通过计算信号传播时间乘以光速获得，具体过程如下。相对于R星地心地固坐标系时刻在转动，计算的每一时刻的丑星位置都是相对当前的地心地固坐标系。以卫星信号放射时刻地心地固坐标系为参考，接收机位置已知，接收时刻卫星的位置已知，放射时刻卫星的位置。给7；赋初值为67ms,Tt=T1.-Tp(3.22)其中r,表示信号接收时间，7；表示信号放射时间，7；表示信号传输时间。依据Z和卫星星历可以计算出卫星在信号放射时刻地心地固坐标系下的坐标%,计算伪距必须要在同一个坐标系下，考虑地球白转，对地心地固坐标系进行转换：

29、COs(纬*Qsin(u*7j,)0%=-SinM*7；)cos(yk*T;,)OOOI(3.23)坐标系转换之后可以计算得到新的：?一一(3.24)进行迭代计算，得到，其中3.23和3.24中的参数定义如下。%,：相对于frame为接收时刻地心地固坐标系卜的接收时刻接收机位置。%.：相对J-frame为接收时刻地心地固坐标系下的放射时刻卫星位置。%：相对于frame为放射时刻地心地固坐标系卜的放射时刻卫星位置。3.3.4 多普勒频移以与其它参数的计算利用以上的计算结果，伪距等于当前时刻载体位置和信号放射时刻卫星位置之间的距离，多普勒频移等于接收时刻载体的速度和卫星速度在伪距方向的重量之和，计

30、算过程如下：“-(v,-vw),%)+1.)rrr(3.25)r=.W-Sf+()”)+(Zifr-Z“/(3.26)=-Ac(3.27)其中参数定义如下：vrf：卫星信号放射时刻速度与载体接收时刻速度在伪距方向上速度之差。vw:载体在X方向上的重量。vhi：载体在)，方向上的重量。Va,:载体在2方向上的重量。r:伪距。:多普勒频移。：载波频率。C：光速。木软件中PR_dop_Ca1.cu1.a1.e.m函数的功能就是计算每一个时间点的伪距、载波相位、伪内相位、多普勒频移。利用上文中叙述的算法，通过反熨迭代计算出信号传播时间，据此进一步解算出要求解的变量。3.4卫星信号生成以上介绍的内容已经

31、把生成卫星信号所须要的参数计算完毕，接卜来的内容就是如何利用这些参数生成卫星信号。生成卫星信号须要导航电文、伪码、以与我波，这些信号的生成方法上文中已作了介绍。在本软件的仿真过程中首先调用导航电文生成程序，把生成的导航电文保存在变量data,然后再把每一个导航为扩展为20个相同的数字，这样每一个数字对应1ms,也就是一个C/A码周期。其次步是仿真伪码，仿真伪内的过程中须要调用上文中计算的伪碍相位，生成的伪码和导航电文进行调制，调制过程中须要参照传输时间把对应J-C码的导航电文向前移动相应的量，然后二者相乘。第三步就是生成载波信号，载波信号是基于程序中设定的中频频率，在此基础上加上多普勒频移，组

32、成载波的频率，调用上文中计算的载波相位，就可以生成完整的载波信号，此处的载波信号是以来样频率为基准的一些点。在本软件中，程序生成的主程序是GPS_成_Signa1.SimUIation2.m,执行过程中首先谢用PRdopCa1.CU1.ate函数生成全部北星的伪距、多普勒频移、教波相位、伪码相位，把这些变量作为点用函数generate_perfect_GPS_signa1.的输入量。generate_perfect_GPS_signa1函数把仿真的20秒数据分为，10段执行，分40次调用genF_Signa1.3、gen_IF_Signa1_6gen_IF_Signa1._16gen_IF_Signa1._19函数,最终数据经过采样函数SamP1.e采样之后作为最终的中频数据保存。