基于非合作导航卫星信号的无码多普勒定轨方法技术

本发明专利技术公开了基于非合作导航卫星信号的无码多普勒定轨方法，通过接收机的射频前端将导航卫星信号转换为数字中频信号，对中频信号进行消码处理后得到连续载波信号，并估计载波信号的初始频率；将该频率值作为锁相环路跟踪的初值，恢复载波并测量提取出各个导航星信号的多普勒频移；然后采用基于多普勒频移特征的匹配算法寻找出某频移对应的导航星；建立高精度测量模型，结合轨道动力学模型，采用卡尔曼滤波器估计出高精度轨道参数。本发明专利技术能在明码信号不可用的情况下，利用卫星导航系统发送的加密导航信号，实现高精度定轨。

【技术实现步骤摘要】
基于非合作导航卫星信号的无码多普勒定轨方法
本专利技术属于卫星导航
，特别涉及了基于非合作导航卫星信号的无码多普勒定轨方法。
技术介绍
卫星轨道精密测量是卫星应用的基础。目前，国内外采用的定轨方式主要是依赖地基测控网，结合星上应答机获得测距、测速信息后精密校正轨道并上注，但是这种定轨方式不但在恶劣环境时生存保障较为困难，而且定轨精度约为百米到十米量级并随着星上轨道预报后精度逐渐下降，有时难以满足精度要求。法国提出并研制的DORIS定轨系统，是一套星上自主实时定轨软件，其解算周期为10s，能够快速、准确地确定卫星轨道。中国首颗海洋动力环境卫星海洋二号(HY-2)主要用于观测全球海洋动力环境参数。卫星激光测距系统观测精度可达到毫米级，可以满足精度要求，但是全球观测台站少、分布不均，以及受观测区域和气候条件的限制，部分台站还无法实现常规白天测距，使观测数据较为有限，而且SLR设备较为复杂庞大，价格不菲。目前，GPS、GLONASS等国外的卫星导航系统是当前天基跟踪测量的主流手段，是地基定轨系统的重要补充，其中以GPS为主。但是GPS卫星导航系统是国外军方控制，在恶劣环境下，明码信号面临关闭的风险，这样就无法通过普通的接收机得到卫星伪距、伪距率，也就无法实现定位和定轨。我国的北斗系统目前是一个区域导航系统，尚不能实现全轨道覆盖。
技术实现思路
为了解决上述
技术介绍
提出的技术问题，本专利技术旨在提供基于非合作导航卫星信号的无码多普勒定轨方法，在明码信号不可用的情况下，利用卫星导航系统发送的加密导航信号，实现高精度定轨。为了实现上述技术目的，本专利技术的技术方案为：基于非合作导航卫星信号的无码多普勒定轨方法，包括以下步骤：(1)对接收到的导航卫星信号经数字下变频后，采样得到数字中频信号，消除数字中频信号的保密调制码，根据消码后的信号计算得到载波频率的估计值；(2)将步骤(1)得到的载波频率的估计值作为载波跟踪环跟踪的初始频率，产生载波跟踪环的本振信号，将步骤(1)中消码后的信号作为载波跟踪环的输入信号，恢复出载波信号，并测量提取出各个导航卫星信号的多普勒频移；(3)计算各个导航卫星信号的先验值，结合步骤(2)得到的多普勒频移的实测值，建立含有接收机频差的多普勒频移实测值和先验值的识别匹配算法，完成导航星的识别与匹配；(4)利用统计定轨算法实现轨道测量，利用轨道动力学方程建立状态方程，利用载波相位测量值建立相位差测量方程，再利用卡尔曼滤波器估计用户卫星的轨道参数。其中，步骤(1)中所述根据消码后的信号得到载波频率的估计值的具体步骤；(A)选取合适的搜索步长，产生若干候选频率，对每一个候选频率产生一个载波相位调整信号；(B)将步骤(A)产生的载波相位调整信号分别与消码后的信号相乘，对每一个相乘后的信号按照一定单位长度进行分块，累加分块数据，从而得到信号增强的新数据；(C)对步骤(B)得到的新数据进行FFT运算，检查频谱结果，寻找超过门限的尖峰，对于每一个尖峰，记录其对应的频率值，即为载波频率的估计值。其中，步骤(3)的具体步骤：(a)根据步骤(2)恢复出的载波信号，由用户卫星的初始轨道和周期上注的导航卫星星历计算出各导航卫星信号的多普勒频移先验值；(b)将步骤(2)得到的多普勒频移实测值和步骤(a)得到的多普勒频移先验值分别按照高低顺序排列成两组数据样本，引入待求解的整体平移量，该整体平移量用于将多普勒频移先验值作整体平移；(c)对各导航卫星信号的多普勒频移实测值与平移后的多普勒频移先验值的频差绝对值求和得到匹配指标；(d)求解在匹配指标取最小值时的整体平移量，将多普勒频移先验值作整体平移，实现多普勒频移实测值与多普勒频移先验值的匹配，实现导航星识别。其中，在步骤(d)中，采用三分法求解在匹配指标取最小值时的整体平移量，其具体步骤：(Ⅰ)假设某导航卫星信号的多普勒频移先验值作平移的频率区间是选取该区间的1/3、2/3处的坐标把整个区间分成三个小区间[m1,m2]，(Ⅱ)分别计算在m1、m2处匹配指标e(m1)、e(m2)，若e(m1)＞e(m2)，则舍弃区间令否则舍弃令(Ⅲ)重复步骤(Ⅰ)-(Ⅱ)，直到小于某一预设阈值时，即为整体平移量。其中，在步骤(1)中，采用平方法消除数字中频信号的保密调制码。其中，在步骤(1)中，采用交叉相关法消除数字中频信号的保密调制码。采用上述技术方案带来的有益效果：本专利技术可以提高利用卫星导航系统实现卫星定轨的可靠性和安全性，可用作卫星定轨接收机的备份工作模式。同时本专利技术可望达到较高定轨精度。附图说明图1是本专利技术的流程图；图2是射频信号下变频原理框图；图3是平方法消码示意图；图4是交叉相关法消码示意图；图5是锁相环结构示意图；图6是本专利技术三分法的3种情况示意图。具体实施方式以下将结合附图，对本专利技术的技术方案进行详细说明。如图1所示本专利技术的流程图，用户卫星射频前端将导航卫星信号转换为数字中频信号，对得到的中频信号进行消码处理得到连续载波信号；基于相位补偿、数据块累加及频谱分析来精确估计消码后的载波的频率；将频率估计值作为载波跟踪环路的频率初始值，载波环路对载波信号动态跟踪，测量出各个导航星信号的多普勒频移；然后采用基于多普勒频移特征的匹配算法寻找出某频移对应的导航星；将导航卫星当作一个频率已知的无线电信标，建立高精度测量模型，结合轨道动力学模型，采用扩展卡尔曼滤波算法估计出高精度轨道参数。其具体过程如下。1、中频信号获取及其消码处理在GPS接收机中，通常采用混频方式将射频信号从高频下变频得到频率较低的中频模拟信号，继而通过模数转换器转换为数字中频信号。射频信号下变频原理如图2所示。在图2中，射频信号r1(t)和本地载波信号r2(t)经过混频器混频，得到下变频的中频信号。便于信号处理，射频信号采用正交下变频得到正交的中频信号。混频后的信号分为两路，分别采样得到正交的卫星中频信号r(t)。导航卫星信号经正交采样后可表示为：r(t)＝ACDej2πft+n(t)(1)式中，A为信号幅度，C为未知的伪随机码，D为电文数据位，ej2πft是用复指数形式描述的正弦信号，其中f为中频载波频率，t是时间，π是圆周率，j表示复数的虚部，n(t)为白噪声。对中频信号r(t)进行平方或者交叉相关法进行消码处理，平方法和交叉相关法消码原理分别如图3、4所示。以平方法为例，对中频信号r(t)平方可得：r2(t)＝A2C2D2ej2π2ft+n2(t)+2ACDej2πft×n(t)(2)＝A2ej2π2ft+n′(t)处于±1状态的调制码，经平方后均为+1，而+1不改变信号相位和载波频率。平方后的载波信号变为连续的二倍频载波，其他两项频谱是噪声，n′(t)＝n2(t)+2ACDej2πft×n(t)。2、无码载波频率估计技术传统的捕获环节利用伪随机码良好的相关特性，完成对载波频率和码相位的粗略估计，实现跟踪通道初始化以及导航卫星识别。但是，消码处理后的信号削弱了信号强度，所以对该信号进行FFT变换，几乎检测不到峰值；而且在轨高动态条件下的多普勒频移搜索范围为几十kHz，若采用逐步搜索法，会耗费大量时间。所以载波频率初值的获取是需要解决的关键问题。由于信号叠加可增强信号强度，以平方消码后的信号为例，以T为单位分成L本文档来自技高网...

【技术保护点】
基于非合作导航卫星信号的无码多普勒定轨方法，其特征在于，包括以下步骤：(1)对接收到的导航卫星信号经数字下变频后，采样得到数字中频信号，消除数字中频信号的保密调制码，根据消码后的信号计算得到载波频率的估计值；(2)将步骤(1)得到的载波频率的估计值作为载波跟踪环跟踪的初始频率，产生载波跟踪环的本振信号，将步骤(1)中消码后的信号作为载波跟踪环的输入信号，恢复出载波信号，并测量提取出各个导航卫星信号的多普勒频移；(3)计算各个导航卫星信号的先验值，结合步骤(2)得到的多普勒频移的实测值，建立含有接收机频差的多普勒频移实测值和先验值的识别匹配算法，完成导航星的识别与匹配；(4)利用统计定轨算法实现轨道测量，利用轨道动力学方程建立状态方程，利用载波相位测量值建立相位差测量方程，再利用卡尔曼滤波器估计用户卫星的轨道参数。

【技术特征摘要】
1.基于非合作导航卫星信号的无码多普勒定轨方法，其特征在于，包括以下步骤：(1)对接收到的导航卫星信号经数字下变频后，采样得到数字中频信号，消除数字中频信号的保密调制码，根据消码后的信号计算得到载波频率的估计值；(2)将步骤(1)得到的载波频率的估计值作为载波跟踪环跟踪的初始频率，产生载波跟踪环的本振信号，将步骤(1)中消码后的信号作为载波跟踪环的输入信号，恢复出载波信号，并测量提取出各个导航卫星信号的多普勒频移；(3)计算各个导航卫星信号的先验值，结合步骤(2)得到的多普勒频移的实测值，建立含有接收机频差的多普勒频移实测值和先验值的识别匹配算法，完成导航星的识别与匹配；具体步骤如下：(a)根据步骤(2)恢复出的载波信号，由用户卫星的初始轨道和周期上注的导航卫星星历计算出各导航卫星信号的多普勒频移先验值；(b)将步骤(2)得到的多普勒频移实测值和步骤(a)得到的多普勒频移先验值分别按照高低顺序排列成两组数据样本，引入待求解的整体平移量，该整体平移量用于将多普勒频移先验值作整体平移；(c)对各导航卫星信号的多普勒频移实测值与平移后的多普勒频移先验值的频差绝对值求和得到匹配指标；(d)求解在匹配指标取最小值时的整体平移量，将多普勒频移先验值作整体平移，实现多普勒频移实测值与多普勒频移先验值的匹配，实现导航星识别；(4)利用统计定轨算法实现轨道测量，利用轨道动力学方程建立状态方程，利用载波相位测量值建立相位差测量方程，再利用卡尔曼滤波器估计用户卫星...

【专利技术属性】
技术研发人员：郁丰，卢欢，何真，刘莹，吴云华，华冰，康国华，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32