本发明专利技术属于导航接收机设备领域，涉及一种GNSS信号载波跟踪和导航解算紧组合的滤波方法，具体为：步骤1，接收机将接收到GNSS信号变为数字中频信号r(t)；步骤2，各通道NCO分别生成两路正交信号；步骤3，相关器进行相干积累；步骤4，将相关器的输出结果输至跟踪通道中的鉴别器进行处理；步骤5，本地载波跟踪滤波器进行滤波处理；步骤6，提取观测量；步骤7，导航解算滤波器进行滤波；步骤8，根据接收机中导航滤波器的解算结果，预测下一时刻各个通道中的NCO控制字；步骤9，将本地载波跟踪滤波器的滤波结果与导航解算滤波器的预测值进行紧组合，获得最终通道NCO的控制参数,将紧组合的结果用于控制通道NCO，并同步反馈到本地载波跟踪滤波器中。
【技术实现步骤摘要】
本专利技术属于导航接收机设备领域，具体涉及一种GNSS信号载波跟踪和导航解算紧组合的滤波方法，其可运用在卫星导航系统中接收终端类设备的研制中。
技术介绍
随着卫星导航系统应用范围的不断扩大，导航接收机可能面临特殊应用环境，如丛林，复杂城区等，这就对导航信号的接收提出了更高的要求，目前对导航信号的接收跟踪方法主要有基于标量结构的信号跟踪方法和基于矢量结构的信号跟踪方法，标量信号跟踪方法对利用本地信号跟踪滤波器对每颗卫星进行单独跟踪，其跟踪性能灵敏度及连续性不如矢量跟踪方法，矢量跟踪方法将所有可以收到的卫星信号在一个跟踪环路中进行联合跟踪，可以看作是一种特殊的本地信号跟踪滤波器和导航接收机中的导航滤波器间的一种深组合方法，这种深组合方式仅仅保留了本地信号跟踪滤波器中的误差估计环节，并将其作为导航解算滤波器的输入，对导航信号进行联合跟踪处理，该方法可以有效提高跟踪门限，并且可以对较弱信号进行连续跟踪，该方法的稳健性较差，且系统复杂度较高，不利于接收机的实现。
技术实现思路
针对上述介绍的两种方法不能较好的解决导航接收机对弱信号跟踪性能与系统复杂性及稳健性间的矛盾，本专利技术提供了一种GNSS信号载波跟踪和导航解算紧组合的滤波方法，可以在提高导航接收机对弱信号跟踪性能的同时，保留接收机中的本地载波跟踪滤波器，有效降低系统的复杂性，提高接收机的稳健性。具体技术方案如下：一种GNSS信号载波跟踪和导航解算紧组合的滤波方法，包括以下步骤：步骤1，接收机通过天线接收到GNSS信号后，GNSS信号经过接收机中的射频前端、AD转换器后变为数字中频信号r(t)，t表示接收时刻；步骤2，接收机中为每颗卫星分配一个跟踪通道，N颗卫星对应N个跟踪通道，每个跟踪通道中均含有一个本地载波生成装置NCO(NCO:NumericallyControlledOscillator)，用于生成频率控制字为fNCO的两路正交信号，分别为同相信号sI(t)和正交信号sQ(t)，以其中任意一个通道为例，上述两路信号具体为：sI(t)＝cos(2πfNCO·t)sQ(t)＝-sin(2πfNCO·t)将进入该通道的数字中频信号分为相同的两路信号，一路信号与同相信号相乘混合为sI(t)·r(t)，另一路信号与正交信号相乘混合为sQ(t)·r(t)；步骤3，接收机跟踪通道中的相关器，用于将本地复制生成的信号和接收信号进行相干积累，相干积分时间为Tc，k表示跟踪环路中第k个跟踪历元，每个历元对应的时长为Tc，故输出信号的积分区间为(k-1)·Tc到k·Tc，得到相关值Ik和Qk，具体为：步骤4，接收机跟踪通道中的鉴频器装置，用于获得本地复制信号和接收信号间的误差估计参数，其鉴频器输出结果εk为：其中atan2表示四象限反正切函数；步骤5，接收机中的本地载波跟踪滤波器进行滤波，其步骤如下：本地载波跟踪滤波器的系统方程为Yk+1＝ΦY·Yk+wk其中Yk＝[fk,αk]T为本地载波跟踪滤波器的系统状态向量，fk,αk分别表示信号的多普勒和多普勒变化率，单位分别为Hz，Hz/s；T表示向量或矩阵的转置符号，wk＝[ωrf·wd；(ωrf/c)·wa]T为系统噪声，wd为由接收机中晶体振荡器引起的频率噪声，其噪声谱密度为qd；wa是系统频率变化率噪声，其功率谱密度为qa，ωrf是载波频率，c为光速；ΦY是本地载波跟踪滤波器的系统状态转移矩阵，具体为：QY是wk对应的过程噪声协方差矩阵，具体为：鉴频器的输出结果εk为本地载波跟踪滤波器的新息，其协方差矩阵为RY，具体为：其中C/N0为信号载噪比，表示信号强度的大小；π为圆周率。本地载波跟踪滤波器的量测矩阵为HY，具体为：HY＝[1-Tc/2]本地载波跟踪滤波器的具体实施步骤描述为：Step51：计算系统状态向量预测值Yk-及其预测协方差矩阵其中，PY,k-1表示Yk-1对应的协方差矩阵；Step52：计算本地载波跟踪滤波器增益矩阵KY：Step53：根据新息更新系统状态向量及其协方差矩阵：其中I表示单位矩阵；最终第k个历元NCO的控制参数为：fNCO,k＝Yk(1)其中Yk(1)表示系统状态向量的第一个元素；步骤6，提取观测量，即从N个通道中获取导航解算所需要的多普勒测量结果，生成观测矢量其中为卫星i在第k个历元的多普勒频率测量结果，N为卫星总数，i取值范围为1,2,…,N；的协方差矩阵为Rz，具体为：其中为第i颗卫星信号本地载波跟踪滤波器中的系统状态协方差矩阵，通过步骤5中的step53获取，其第一行第一列对应的元素为多普勒测量结果的精度，diag表示对角矩阵符号；步骤7，导航解算滤波器，导航解算滤波器通常采用一个ExtendKalmanFilter(EKF)实现，其系统方程为Xk+1＝Φ·Xk+ωk其中Xk＝[vx,ax,vy,ay,vz,az,δf]为第k个历元时刻接收机的运动状态矢量，vx,vy,vz为ECEF坐标系下的三维速度结果，ax,ay,az为ECEF(Earth-Centered,Earth-Fixed)坐标系下的三维加速度结果，δf为接收机上时钟的频率误差，状态转移矩阵为Φ，具体表示为：其中，Tb是导航解算滤波器的解算间隔；ωk是系统过程噪声，其协方差矩阵为Q，具体为：其中Qf＝Sf·TbSa为加速度噪声功率谱密度，Sf为时钟频率变化噪声功率谱密度。根据观测量步骤6中的观测量得到系统的测量方程为：其中H为测量矩阵，它由接收机到卫星的空间几何构型决定，为第k个历元时刻卫星的运动状态，上标s表示所有卫星的集合，其中为所有卫星在x方向上的速度值的集合，其元素个数为N个，同理可获得δfS为卫星上时钟的频率误差；导航解算滤波器测速解算的具体过程描述如下：Step71,计算接收机的运动状态矢量预测结果及其协方差矩阵其中Pk-1表示Xk-1的协方差矩阵。Step72，计算导航解算滤波器的增益矩阵Gk：Step73，更新接收机的运动状态矢量及其协方差矩阵：步骤8，通道NCO预测器，根据接收机中导航滤波器的解算结果，预测下一时刻各个通道中的NCO控制字，具体为：其中，Xk-1为第k-1个历元时刻的接收机的运动状态矢量。计算预测精度为：这样对于第i颗卫星而言，下一时刻NCO参数的预测结果及精度为：其中表示向量的第i个元素，表示矩阵的第i行第i列元素值；步骤9，本地载波跟踪滤波器的滤波结果和导航解算滤波器的预测值进行紧组合，即将步骤5所得结果Yk(1)和步骤8所得结果进行线性加权，从而获得最终通道NCO的控制参数fNCO,k，具体为其中λ为组合权重，它由两者的精度决定，具体为：PY,k(1,1)为对应的通道中本地载波跟踪滤波器的系统状态向量第1行第1列元素值，若是第i个通道，则为需要说明的是紧组合的结果fNCO,k再控制通道NCO的同时，还需要同步反馈到本地载波跟踪滤波器中，以保持本地载波跟踪滤波器的状态和通道NCO参数的一致性。上述步骤中描述了本地载波跟踪滤波器和导航解算滤波器更新频度一致的情况，实际上，两者的更新频度可能出现不一致，但通常情况为导航解算滤波器的更新间隔为本地载波跟踪滤波器更新间隔的M倍，M为整数，在这种情况下，本方法中本地载波跟踪滤波器和导航解算滤波器的紧组合将发生导航解算滤波器进行更新本文档来自技高网...
【技术保护点】
一种GNSS信号载波跟踪和导航解算紧组合的滤波方法，其特征在于，包括以下步骤：步骤1，接收机通过天线接收到GNSS信号后，GNSS信号经过接收机中的射频前端、AD转换器后变为数字中频信号r(t)，t表示接收时刻；步骤2，接收机中为每颗卫星分配一个跟踪通道，N颗卫星对应N个跟踪通道，每个跟踪通道中均含有一个本地载波生成装置NCO，用于生成频率控制字为fNCO的两路正交信号，分别为同相信号sI(t)和正交信号sQ(t)，每个通道中的两路信号具体为：sI(t)＝cos(2πfNCO·t)sQ(t)＝‑sin(2πfNCO·t)将数字中频信号分为相同的两路信号，一路信号与同相信号相乘混合为sI(t)·r(t)，另一路信号与正交信号相乘混合为sQ(t)·r(t)；步骤3，接收机跟踪通道中的相关器进行相干积累，设相干积分时间为Tc，k表示跟踪环路中第k个跟踪历元，每个历元对应的时长为Tc，故输出信号的积分区间为(k‑1)·Tc到k·Tc，得到相关值Ik和Qk，具体为：Ik=∫(k-1)·Tck·TcsI(t)·r(t)dt]]>Qk=∫(k-1)·Tck·TcsQ(t)·r(t)dt]]> 步骤4，将相关器的输出结果输至跟踪通道中的鉴别器进行处理，得到信号间的误差估计参数εk：ϵk=a tan2(Ik-1Qk-Qk-1Ik,Ik-1Ik+Qk-1Qk)Tc]]> 其中atan2表示四象限反正切函数； 步骤5，本地载波跟踪滤波器进行滤波处理，具体过程为： 本地载波跟踪滤波器的系统方程为， Yk+1＝ΦY·Yk+wk其中Yk＝[fk,αk]T为本地载波跟踪滤波器的系统状态向量，fk,αk分别表示信号的多普勒和多普勒变化率，单位分别为Hz，Hz/s；wk＝[ωrf·wd；(ωrf/c)·wa]T为系统噪声，wd为由接收机中晶体振荡器引起的频率噪声，其噪声谱密度为qd；wa是系统频率变化率噪声，其功率谱密度为qa，ωrf是载波频率，c为光速；ΦY是本地载波跟踪滤波器的系统状态转移矩阵，具体为：ΦY=1Tc01]]> wk是系统过程噪声，QY是wk对应的过程噪声协方差矩阵，具体为：QY=ωrf2qdTc000+(ωrfc)2qaTc3/3Tc2/2Tc2/2Tc]]> 鉴频器的输出结果εk为本地载波跟踪滤波器的新息，其协方差矩阵为RY：RY=(12πTc)2·1C/N0·Tc·(1+1C/N0·Tc)]]> 其中C/N0为信号载噪比；本地载波跟踪滤波器的量测矩阵为HY，具体为：HY＝[1 ‑Tc/2]本地载波跟踪滤波器的滤波步骤如下：Step51：计算系统状态向量预测值及其预测协方差矩阵Yk-=ΦY·Yk-1]]>PY,k-=ΦY·PY,k-1·ΦYT+QY]]>其中，PY,k‑1表示Yk‑1对应的协方差矩阵；Step 52：计算本地载波跟踪滤波器增益矩阵KY：KY=PY,k-·HYT(HY·PY,k-·HYT+RY)-1]]>Step53：根据新息更新系统状态向量及其协方差矩阵：Yk=Yk-+KY·ϵk]]>PY,k=(I-KY·HY)·PY,k-]]>其中I表示单位矩阵；第k个历元NCO的控制参数fNCO,k为：fNCO,k＝Yk(1)其中Yk(1)表示系统状态向量的第一个元素；步骤6，提取观测量，从N个通道中获取导航解算所需要的多普勒测量结果，生成观测矢量其中为卫星i在第k个历元的多普勒频率测量结果，N为卫星总数，其协方差矩阵为Rz：Rz=diag([PY,k1(1,1),PY,k2(1,1),...,PY,kN(1,1)])]]>其中为第i颗卫星信号本地载波跟踪滤波器中的系统状态协方差矩阵，通过步骤5中的step 54获取，其第一行第一列对应的元素为多普勒测量结果的精度；步骤7，导航解算滤波器进行滤波，具体过程如下：导航解算滤波器系统方程为，Xk+1＝Φ·Xk+ωk其中Xk＝[vx,ax,vy,ay,vz,az,δf]为第k个历元时刻接收机的运动状态矢量，vx,vy,vz为ECEF坐标系下的三维速度结果，ax,ay,az为ECEF坐标系下的三维加速度结果，δf为接收机上时钟的频率误差，状态转移矩阵为Φ，具体为：Φ=&Ph...
【技术特征摘要】
1.一种GNSS信号载波跟踪和导航解算紧组合的滤波方法，其特征在于，包括以下步骤：步骤1，接收机通过天线接收到GNSS信号后，GNSS信号经过接收机中的射频前端、AD转换器后变为数字中频信号r(t)，t表示接收时刻；步骤2，接收机中为每颗卫星分配一个跟踪通道，N颗卫星对应N个跟踪通道，每个跟踪通道中均含有一个本地载波生成装置NCO，用于生成频率控制字为fNCO的两路正交信号，分别为同相信号sI(t)和正交信号sQ(t)，每个通道中的两路信号具体为：sI(t)＝cos(2πfNCO·t)sQ(t)＝-sin(2πfNCO·t)将数字中频信号分为相同的两路信号，一路信号与同相信号相乘混合为sI(t)·r(t)，另一路信号与正交信号相乘混合为sQ(t)·r(t)；步骤3，接收机跟踪通道中的相关器进行相干积累，设相干积分时间为Tc，k表示跟踪环路中第k个跟踪历元，每个历元对应的时长为Tc，故输出信号的积分区间为(k-1)·Tc到k·Tc，得到相关值Ik和Qk，具体为：Ik=∫(k-1)·Tck·TcsI(t)·r(t)dt]]>Qk=∫(k-1)·Tck·TcsQ(t)·r(t)dt]]>步骤4，将相关器的输出结果输至跟踪通道中的鉴别器进行处理，得到信号间的误差估计参数εk：ϵk=atan2(Ik-1Qk-Qk-1Ik,Ik-1Ik+Qk-1Qk)Tc]]>其中atan2表示四象限反正切函数；步骤5，本地载波跟踪滤波器进行滤波处理，具体过程为：本地载波跟踪滤波器的系统方程为，Yk+1＝ΦY·Yk+wk其中Yk＝[fk,αk]T为本地载波跟踪滤波器的系统状态向量，fk,αk分别表示信号的多普勒和多普勒变化率，单位分别为Hz，Hz/s；wk＝[ωrf·wd；(ωrf/c)·wa]T为系统噪声，wd为由接收机中晶体振荡器引起的频率噪声，其噪声谱密度为qd；wa是系统频率变化率噪声，其功率谱密度为qa，ωrf是载波频率，c为光速；ΦY是本地载波跟踪滤波器的系统状态转移矩阵，具体为：ΦY=1Tc01]]>wk是系统过程噪声，QY是wk对应的过程噪声协方差矩阵，具体为：QY=ωrf2qdTc000+(ωrfc)2qaTc3/3Tc2/2Tc2/2Tc]]>鉴频器的输出结果εk为本地载波跟踪滤波器的新息，其协方差矩阵为RY：RY=(12πTc)2·1C/N0·Tc·(1+1C/N0·Tc)]]>其中C/N0为信号载噪比；本地载波跟踪滤波器的量测矩阵为HY，具体为：HY＝[1-Tc/2]本地载波跟踪滤波器的滤波步骤如下：Step51：计算系统状态向量预测值及其预测协方差矩阵Yk-=ΦY·Yk-1]]>PY,k-=ΦY·PY,k-1·ΦYT+QY]]>其中，PY,k-1表示Yk-1对应的协方差矩阵；Step52：计算本地载波跟踪滤波器增益矩阵KY：KY=PY,k-·HYT(HY·PY,k-·HYT+RY)-1]]>Step53：根据新息更新系统状态向量及其协方差矩阵：Yk=Yk-+KY·ϵk]]>PY,k=(I-KY·HY)·PY,k-]]>其中I表示单位矩阵；第k个历元...
【专利技术属性】
技术研发人员：唐小妹，林红磊，徐博，欧钢，吴健，李蓬蓬，彭竞，伍微，雍玲，
申请(专利权)人：中国人民解放军国防科学技术大学，
类型：发明
国别省市：湖南;43
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。