本发明涉及gnss接收机,具体涉及一种联合gnss接收机多通道信号的最优鉴别和跟踪方法。
背景技术:
1、随着卫星全球导航系统(gnss)在各领域的深入应用,导航接收机使用性能越来越受到恶劣环境导致的挑战。导航定位模块作为导航接收机的核心部件,其主要包括捕获、跟踪和数据解算三部分,为增强导航定位模块复杂环境下的适应能力,开展了一系列的研究与改进。目前的研究的工作可分为两类:一是对导航定位模块的关键部件进行改造和升级,例如用卡尔曼滤波器代替传统dll(延迟锁定环)、pll(锁相环)/fll(锁频环)环路,以便提高性能;二是对导航定位模块的结构进行调整,比较典型的作法是将通道独立工作的标量跟踪接收机升级为通道融合的矢量跟踪接收机。
2、针对第一类改进,除了优化滤波器和nco(数字控制振荡器)外,也经常通过扩展线性牵入范围、增加相干积分和非相干积分时间类提高鉴别器性能。
3、针对第二类改进,矢量跟踪产生性能提升的实质是将信号处理和位置解算进行一体化设计,通过多通道融合的信息优势转化为跟踪阶段的处理优势。因此,矢量跟踪机制正常工作的前提是要得到先验的位置信息,才能进行后续的滤波融合操作。但对于无法持续跟踪的环境,如失锁重捕,此时没有建立完整的闭合环路,即不能发挥矢量跟踪的作用,需要重新回到标量模式中。
技术实现思路
1、本发明的目的是解决对导航定位模块的关键部件进行改造和升级的提升力度有限,需要加长积分时间改善性能,以及矢量跟踪只能在持续跟踪发挥作用的不足之处,而提供一种联合gnss接收机多通道信号的最优鉴别和跟踪方法。无论是否存在有效位置信息,都可通过通道信号融合的方式获得辅助增益,改善接收机整体处理信号的能力。
2、为了解决上述现有技术所存在的不足之处,本发明提供了如下技术解决方案:
3、一种联合gnss接收机多通道信号的最优鉴别和跟踪方法,其特殊之处在于,包括如下步骤:
4、步骤1、初始化滤波器的卡尔曼滤波的状态矩阵sn、测量矩阵δω′、状态估计误差协方差矩阵、卡尔曼滤波的系数矩阵f,并设置初始的预测状态估计值sm+1;根据状态矩阵sn、卡尔曼滤波的系数矩阵f,计算状态传输方程,其中初始的权矩阵定义为单位矩阵;
5、步骤2、使用卡尔曼滤波的系数矩阵f和当前时刻状态估计误差协方差矩阵,计算预测状态估计误差协方差矩阵;
6、步骤3、将卫星信号通过多个通道输入鉴别器,通过鉴别器获取接收信号向量,将接收信号向量与本地复现信号向量进行比较,得到参数残差向量;对参数残差向量进行抗差估计后,作为滤波器的输入值;
7、所述接收信号向量为信号幅值ai、码相位延迟τi、载波相位载波频率fd,i;所述参数残差向量为码相位差、载波相位差以及信号幅值差;
8、步骤4、将步骤3的参数残差向量作为预测状态实际值,与步骤1设置的初始的预测状态估计值sm+1比较,计算观测残差,并计算观测残差协方差矩阵;
9、步骤5、依据鉴别器输出的参数残差向量的载噪比,定义一个权因子函数,以在信号层面获得权因子;将权因子作为对角元素构成一个对角矩阵,然后与初始的权矩阵构成等价权矩阵
10、步骤6、通过等价权矩阵步骤4的观测残差协方差矩阵求解增益矩阵kk;
11、步骤7、通过步骤6的增益矩阵kk、步骤4的观测残差更新下一时刻状态估计值sm+1、当前时刻状态估计误差协方差矩阵;
12、返回步骤2,重复进行鉴别和跟踪。
13、进一步地,步骤1中,所述状态矩阵sn如下:
14、
15、其中,xuyuzu表示用户位置,表示用户速度,tu、表示钟偏、钟漂,表示信号幅值,表示信号载波相位;pu表示xuyuzutu,表示
16、所述测量矩阵δω′如下:
17、
18、δau表示信号幅度残差,表示信号载波相位残差,δpu表示经过伪码相位残差δτu处理折算后的用户位置以及钟偏残差。
19、进一步地,步骤1中,所述状态传输方程如下:
20、
21、
22、其中,m表示当前时刻,m+1表示下一时刻,sm表示当前时刻状态估计值,sm+1表示下一时刻的状态向量,ωm表示当前时刻参数向量,in×n表示单位矩阵,q表示噪声矩阵,fd表示信号频率,t表示积分时间,c表示光速,hecef表示地心地固坐标系的余弦矩阵。
23、进一步地,步骤3中,所述参数残差向量δω通过下式求解:,
24、
25、其中,x表示接收到的接收信号向量,表示本地复现信号向量,其是基于鉴别器产生的参考信号和接收信号向量进行匹配和同步后得到的向量,n表示噪声向量,g表示最小二乘系数矩阵,δx为接收到的信号向量与本地复现信号向量之间的残差,用于调整本地产生的参考信号的频率。
26、进一步地,所述步骤5具体如下:
27、步骤5.1、采取离散的方法求鉴别器输出的参数残差向量均值,进一步利用曲线拟合求得鉴别s曲线,从而分析鉴别器的线性工作区间;
28、步骤5.2、根据步骤5.1得到的线性工作区间,定义权因子函数如下:
29、
30、其中,表示第i个通道的权因子,cnri表示第i个通道的载噪比;
31、步骤5.3、将权因子作为对角元素构成一个对角矩阵,然后与初始的权矩阵构成等价权矩阵
32、进一步地,所述步骤5.1具体为:
33、根据鉴别器的输出函数f(δθ+φ)求解鉴别器的输出期望μ:
34、
35、其中,δθ为真实相位差,即载波相位估计值与载波相位真实值之差,φ表示热噪声引起的相位误差,p(φ)表示φ的概率密度函数;
36、由于鉴别器的输出包括两部分,即out=δθ+φ,则鉴别器输出的跟踪参数残差矩阵的期望可以表达为:
37、
38、利用鉴别器的输出、真实相位差δθ获得鉴别s曲线;
39、根据鉴别s曲线,获得偏差在x%内的线性工作区间lrx%为:
40、
41、kd表示鉴别增益,μ′表示载波相位真实值;
42、x取值为1或5或10。
43、进一步地,步骤6中,所述增益矩阵kk如下:
44、
45、其中,表示卡尔曼滤波的系数矩阵的转置,r表示观测残差协方差矩阵。
46、与现有技术相比,本发明的有益效果是:
47、(1)本发明一种联合gnss接收机多通道信号的最优鉴别和跟踪方法,在鉴别器鉴别阶段创造性的提出基于最小二乘估计的多通道鉴别方法,该多通道鉴别方法可以通过最小二乘一次性从获得信号幅值ai、码相位延迟τi、载波相位载波频率fd,i,进一步得到参数残差向量,然后将经过抗差处理之后的参数残差向量作为卡尔曼滤波的输入,从而得到滤波之后的参数;本发明克服了传统矢量接收机通过位置反馈提高处理增益、无法在失锁环境中持续跟踪的不足。
48、(2)本发明一种联合gnss接收机多通道信号的最优鉴别和跟踪方法,结合后续的抗差矢量跟踪算法,建立鉴别和滤波协调处理的理论联系,构建了以最终参量(用户位置、用户速度)、中间过程参量(信号幅值、信号载波相位)作为状态量的混合卡尔曼滤波框架,在鉴别阶段只输出幅值、码相位以及载波相位,其中码相位在卡尔曼滤波输入之前需要转化为位置和钟偏,并在此过程中利用构造的等价权矩阵求解增益矩阵,从而保证最终滤波输入的结果是经过抗差估计处理的,获得更加精确的原始观测数据。
49、(3)本发明一种联合gnss接收机多通道信号的最优鉴别和跟踪方法中,抗差环路在信号级别进行抗差处理,权重的选择主要依据载噪比的强度;本发明根据鉴别器的线性工作区间确定信号级的权值,以实现在短相干积分时间条件下从鉴别和滤波环节提高信号处理的增益,使得gnss接收机在连续和突变的恶劣环境下均能够取得更好的定位效果。