著录项信息
专利名称 | 一种矢量形式信息分配系数的组合导航系统自适应联邦滤波方法 |
申请号 | CN201110266447.X | 申请日期 | 2011-09-09 |
法律状态 | 授权 | 申报国家 | 暂无 |
公开/公告日 | 2012-02-15 | 公开/公告号 | CN102353378A |
优先权 | 暂无 | 优先权号 | 暂无 |
主分类号 | G01C21/24 | IPC分类号 | G;0;1;C;2;1;/;2;4查看分类表>
|
申请人 | 南京航空航天大学 | 申请人地址 | 江苏省南京市白下区御道街29号
变更
专利地址、主体等相关变化,请及时变更,防止失效 |
权利人 | 南京航空航天大学 | 当前权利人 | 南京航空航天大学 |
发明人 | 陈计辉;熊智;王融;张丹;刘建业;张旭;张承;刘伟霞;姚小松 |
代理机构 | 南京经纬专利商标代理有限公司 | 代理人 | 许方 |
摘要
本发明公开了一种矢量形式信息分配系数的自适应联邦滤波方法,该方法在分析GPS、天文(CNS)和SAR三种导航传感器的工作环境和工作特性基础上,建立地理系下惯性导航系统与各导航传感器的线性化量测方程,构建滤波子系统;然后分别对滤波子系统进行卡尔曼滤波,结合子系统估计误差协方差阵特征值和子系统可观测矩阵奇异值,设计完成了一种矢量信息分配系数计算方案,针对矢量分配系数会导致协方差阵不对称问题,设计了新的信息分配方程;最后采用联邦滤波方法对惯性导航的误差状态进行了最优估计。本发明导航精度高,能充分发挥多传感器组合导航系统对动态环境下惯性导航系统误差状态量的估计作用。
1.一种矢量形式信息分配系数的组合导航系统自适应联邦滤波方法,其特征在于包括以下步骤:
(1)通过建立惯性导航系统INS的误差状态量方程,得到对惯性导航系统误差状态量的数学描述,惯性导航系统误差状态量X定义为:
其中:φN,φE,φD分别表示惯性导航系统误差状态量中的北向平台误差角状态量、东向平台误差角状态量和地向平台误差角状态量;δvN,δvE,δvD分别表示惯性导航系统误差状态量中的北向速度误差状态量、东向速度误差状态量和地向速度误差状态量;
δL,δλ,δh分别表示惯性导航系统误差状态量中的纬度误差状态量、经度误差状态量和高度误差状态量;εbx,εby,εbz分别表示惯性导航系统误差状态量中的X轴、Y轴、Z轴方向陀螺常值漂移误差状态量;εrx,εry,εrz分别表示惯性导航系统误差状态量中的X轴、Y轴、Z轴方向陀螺一阶马尔可夫漂移误差状态量; 分别表示惯性导航系统误差状态量中的X轴、Y轴和Z轴方向加速度计零偏,上标T为转置;
(2)建立地理系下各子系统的量测方程,包括GPS/INS量测方程、CNS/INS量测方程和SAR/INS量测方程;
(3)将步骤(2)所述的各子系统的量测方程中子系统误差状态量进行KF滤波,将子系统KF滤波结果送入联邦滤波器;
(4)根据各导航子系统的协方差矩阵以及可观测性矩阵,求取一种矢量形式的联邦滤波信息分配系数,使每个子系统的每个状态变量得到不同的信息分配系数,根据该矢量形式的联邦滤波信息分配系数,建立系统的过程信息在各子滤波器之间的分配原则;
其中,求取一种矢量形式的联邦滤波信息分配系数的具体步骤如下:
基于子系统协方差矩阵的特征值求取矢量信息分配系数Αi,基于子系统可观测矩阵奇异值求取矢量信息分配系数γi,然后融合Αi、γi信息得到联邦滤波信息分配系数Βi,具体步骤为:
1)Αi的求解
将子系统协方差矩阵Pi按特征值分解为:
式中,Λi=diag{λi1,λi2,…,λin},λi1,λi2,…,λin为Pi的特征值,n为Pi阵的阶数,Li为矩阵Pi进行特征值分解后得到的矩阵,上标T为转置;
对第i个局部滤波器的误差状态量Xi的每一个状态分量xij,xij表示第i个局部滤波器误差状态估计中的第j个分量,独立进行信息分配系数计算,得到信息分配系数为:
i=1,2,…,N;j=1,2,…,n
式中,λij为状态变量xij对应的特征值,N为子系统个数;
Xi对应的信息分配系数为矩阵形式:
2)γi的求解
p×q
设某时间段动态系统的可观测性矩阵为Q,Q∈R ,对Q进行奇异值分解,得到T
Q=USV
式中:U=[u1 u2 … up],为p×p阶正交矩阵,称为Q的左奇异矩阵,u1、u2…up分别为左奇异列向量,V=[v1 v2 … vq],为q×q阶正交矩阵,称为Q的右奇异矩阵,v1、v2…vq分别为右奇异列向量; Λr×r=diag(σ1,σ2,…,σr),其中r、σi分别为矩阵Q的秩和奇异值,i=1,2,…,r;
对矩阵V进行分析,得到σi为与其对应的右奇异向量vi中取得最大绝对值的状态变量所对应的奇异值:
对每个子系统的可观测性矩阵分别进行奇异值分解,进一步对Xi的每一个分量xij独立进行信息分配系数计算,得到信息分配系数为:
i=1,2,…,N;j=1,2,…,n
式中,σij为状态变量xij对应的奇异值,N为子系统个数;
则Xi对应的信息分配系数为矩阵形式:
3)Βi的求解
即
同时,联邦滤波信息分配系数满足信息守恒原理:
所述的根据该矢量形式的联邦滤波信息分配系数,建立系统的过程信息在各子滤波器之间的分配原则,表达式如下:
式中, 表示k时刻i子滤波器的状态估计信息量, 表示k时刻联邦滤波
器的总的状态估计信息量, 表示k时刻i子滤波器的过程噪声信息量, 表示k时刻联邦滤波器的总的过程噪声信息量, 表示k时刻i子滤波器的状态估计值,表示k时刻联邦滤波器的状态估计值,Βi表示总的状态估计信息量 的信息分配系数,ΒiQ为Βi中后9个对角线元素,表示总的过程噪声信息量 的分配系数;
(5)联邦滤波器对步骤(3)中子系统送来的滤波结果进行数据融合,输出全局最优估计值,从而对惯性导航系统的导航误差进行修正;
其中所述融合的方法为:
式中,Ρg表示联邦滤波后的误差协方差阵,Ρi表示i子滤波器误差协方差阵, 表示联邦滤波器的状态估计值, 表示i子滤波器的状态估计值,上标“-1”表示求逆。
一种矢量形式信息分配系数的组合导航系统自适应联邦滤\n波方法 \n技术领域\n[0001] 发明涉及一种矢量形式信息分配系数的自适应联邦滤波方法,属于飞行器组合导航技术领域,可应用于高空长时间飞行的航空飞行器导航参数的确定,用于提高导航系统的导航精度。 \n背景技术\n[0002] 随着航空航天技术的发展,对导航与制导系统的性能要求越来越高,组合导航系统已成为当前重要的导航系统方案。在导航多传感器信息融合领域,分散化滤波技术日益受到人们重视,Carlson提出的分散化联邦滤波器得到越来越广泛的应用。信息分配是联邦滤波器研究和设计中的关键,直接影响联邦滤波器的精度和容错性等性能。 \n[0003] Carlson最先提出的信息分配原则是固定比例的,通常根据正常工作子系统的个数进行固定平均分配。考虑到在实际的高动态导航环境中,各子滤波器的性能和估计质量都是不断变化的,为进一步提高联邦滤波器的性能,国内外学者提出了多种动态的信息分配方法,主要包括基于估计协方差特征值法、基于协方差阵Fronenius范数方法、基于可观测矩阵条件数法等,这些动态信息分配方法从一定程度上实现了对系统状态的实时跟踪,滤波效果较固定比例的分配方法有所改善。 \n[0004] 经分析不难发现,以上研究的信息分配系数均为标量形式。即把状态变量作为一个整体考虑来进行信息反馈,其实质是把状态变量的每个分量的特性看作是相同的,具有完全相同的估计精度和收敛速度。但这种观点是不合理的,根据实际情况知,每个传感器功能特性和精度的不同使得每个局部滤波器建立的状态方程和量测方程是不同的,而状态方程的特性又决定了状态的每个分量具有不同的可观测度和收敛速度,在基于惯性导航的组合模式中,平台的方位误差角收敛速度明显低于其它的状态分量。再者,不同的传感器其输出特性也是不同的,如SAR图像匹配仅能输出经度、纬度信息,大气数据系统仅能提供高度及地向速度信息,从而导致不同的状态分量具有不同的估计精度。 \n[0005] 因此,现有标量形式信息分配系数的联邦滤波方法难以具体实时跟踪反映每个状态变量的变化特性,导航精度不高,不能充分满足动态长航时飞行器对导航精度的要求。 发明内容\n[0006] 本发明所要解决的技术问题是克服标量形式信息分配系数的联邦滤波不能充分具体反映 每个状态量变化特性的缺陷,提供一种矢量形式信息分配系数的自适应联邦滤波方法。 \n[0007] 本发明为解决上述技术问题采用以下技术方案: \n[0008] 本发明基于系统误差协方差阵特征值和系统可观测矩阵奇异值的矢量形式信息分配系数的联邦滤波方法,包括以下步骤: \n[0009] (1)通过建立惯性导航系统INS的误差状态量方程,得到对惯性导航系统误差状态量的数学描述,惯性导航系统误差状态量X定义为: \n[0010] \n[0011] 其中:φN,φE,φD分别表示惯性导航系统误差状态量中的北向平台误差角状态量、东向平台误差角状态量和地向平台误差角状态量;δvN,δvE,δvD分别表示惯性导航系统误差状态量中的北向速度误差状态量、东向速度误差状态量和地向速度误差状态量;\nδL,δλ,δh分别表示惯性导航系统误差状态量中的纬度误差状态量、经度误差状态量和高度误差状态量;εbx,εby,εbz分别表示惯性导航系统误差状态量中的X轴、Y轴、Z轴方向陀螺常值漂移误差状态量;εrx,εry,εrz分别表示惯性导航系统误差状态量中的X轴、Y轴、Z轴方向陀螺一阶马尔可夫漂移误差状态量; 分别表示惯性导航系统误差状态量中的X轴、Y轴和Z轴方向加速度计零偏,上标T为转置; \n[0012] (2)建立地理系下各子系统的量测方程,包括GPS/INS量测方程、CNS/INS量测方程和SAR/INS量测方程; \n[0013] (3)将步骤(2)所述的各子系统的量测方程中子系统误差状态量进行KF滤波,将子系统KF滤波结果送入联邦滤波器; \n[0014] (4)根据各导航子系统的协方差矩阵以及可观测性矩阵,求取一种矢量形式的联邦滤波信息分配系数,使每个子系统的每个状态变量得到不同的信息分配系数,根据该矢量形式的联邦滤波信息分配系数,建立系统的过程信息在各子滤波器之间的分配原则; [0015] (5)联邦滤波器对步骤(3)中子系统送来的滤波结果进行数据融合,输出全局最优估计值,从而对惯性导航系统的导航误差进行修正。 \n[0016] 本发明采用以上技术方案与现有技术相比,具有以下技术效果: \n[0017] 本发明克服了现有技术无法实时具体反映状态变量特性的不足,构建了一种适用于动态长时间飞行航空飞行器的组合导航系统联邦滤波方法,它具有以下优点: [0018] (1)根据每个滤波子系统每个状态变量的估计精度分别自适应进行分配系数调整; \n[0019] (2)根据每个滤波子系统每个状态变量的可观测性分别自适应进行分配系数调整; \n[0020] (3)根据矢量形式信息分配系数建立的过程信息在各子滤波器之间的分配原则,确保系 统的误差协方差阵为对称阵,从而有效改善联邦滤波器的稳定性和精度。 附图说明\n[0021] 图1为本发明的矢量信息分配系数联邦滤波的一种实施例子的流程图; \n[0022] 图2为仿真的一条飞行航迹; \n[0023] 图3为本发明的导航经度误差与传统滤波导航经度误差的仿真对比图; \n[0024] 图4为本发明的导航纬度误差与传统滤波导航纬度误差的仿真对比图; \n[0025] 图5为本发明的导航高度误差与传统滤波导航高度误差的仿真对比图。 \n具体实施方式\n[0026] 下面结合附图对本发明的技术方案做进一步的详细说明: \n[0027] 如图1所示,本发明的原理是:从地理系导航的角度入手,依据系统状态方程和各子系统的线性化量测方程,构成滤波子系统,从各滤波子系统中提取估计协方差矩阵Pi和可观测性矩阵Qi,综合两者信息求取矢量形式的信息分配系数Bi,从而进一步完成信息分配、最优融合等步骤,实现对组合导航误差状态量的最优估计。具体实施方法如下: [0028] 一、建立惯性导航系统的误差状态量方程 \n[0029] 选择导航坐标系为北东地地理水平坐标系(OnXnYnZn),采用线性卡尔曼滤波器进行组合,系统的状态方程为惯性导航系统的误差状态量方程,通过对惯性导航系统的性能及误差源的分析,可以获得惯导系统的误差状态量方程为: \n[0030] \n[0031] 式中,\n[0032] 其中φN,φE,φD为平台误差角;δvN,δvE,δvD为速度误差;δL,δλ,δh为纬度、经度和高度误差;εbx,εby,εbz,εrx,εry,εrz分别为陀螺常值漂移误差和一阶马尔可夫漂移误差; 为加速度计零偏,上标T为转置。 \n[0033] 二、建立地理系下各子系统的量测方程 \n[0034] ①GPS/INS量测方程 \n[0035] \n[0036] 式中,vnINS veINS vdINS表示为惯性导航系统的北、东、地方向速度,vnGPS,veGPS,vdGPS表示为GPS接收机的北、东、地方向速度,LINS,λINS,hINS表示为惯性导航系统的纬度、经度、高度测量值,LGPS,λGPS,hGPS表示为GPS接收机的纬度、经度、高度测量值,Rn,Re分别表示北向、东向地球曲率半径,L表示航迹纬度值,MnGPS、MeGPS、MdGPS表示为GPS接收机的北、东、地方向的测速误差,NnGPS,NeGPS,NdGPS为GPS接收机沿北、东、地方向的位置测量误差,均考虑为白噪声。 \n[0037] ②CNS/INS量测方程 \n[0038] \n[0039] 其中 θ、ψ表示为飞行器的真实俯\n仰角和航向角。式中γrINS,θpINS,ψhINS分别表示惯性导航系统的横滚、俯仰、航向角,γrCNS,θpCNS,ψhCNS分别表示天文星敏感器的横滚、俯仰、航向角,OrCNS、OpCNS、OhCNS表示为天文星敏感器对横滚角、俯仰角与航向角在地理系下的等效测量误差,均考虑为白噪声。 [0040] ③SAR/INS量测方程 \n[0041] \n[0042] 其中 式(4)中LINS,λINS分别表示惯性导航系统的纬度、\n经度测量值,LSAR,λSAR分别表示SAR景象匹配导航的纬度、经度测量值,NnSAR,NeSAR为SAR景象匹配导航沿北、东方向的位置测量误差。 \n[0043] 三、子系统KF(Kalman Filter)滤波,估计子系统的误差状态量 \n[0044] 状态方程和量测方程的离散化及卡尔曼滤波器 \n[0045] 当采用线性卡尔曼滤波器时,需要对上面连续形式的系统状态方程(1)和量测方程(2)、(3)、(4)进行离散化,从而获得离散形式的系统方程。其离散化形式如下: [0046] \n[0047] 式中 Td为迭代周期。 \n[0048] 从而可以获得系统的线性化卡尔曼滤波器方程如下: \n[0049] \n[0050] \n[0051] \n[0052] \n[0053] \n[0054] 上式中, 表示k-1时刻的状态对k时刻的状态的最优估计值,又称一步预测估值,Φk,k-1表示k-1时刻到k时刻的状态转移矩阵, 表示k-1时刻的系统状态估计值, 表示k时刻的系统状态估计值,Kk表示增益矩阵,Zk表示k时刻的观测矢量,Hk表示k时刻的观测系数矩阵,Pk|k-1表示最优预测估计误差协方差阵,Pk-1表示k-1时刻的系统误差协方差阵,Qk-1表示k-1时刻的噪声方差矩阵,Γk-1表示k时刻的噪声矢量对k+1时刻状态矢量影响的噪声系数矩阵,Rk表示k时刻的量测方差矩阵,Pk表示k时刻的系统误差协方差阵,I为单位矩阵。 \n[0055] 四、矢量形式信息分配系数求取及信息分配 \n[0056] 基于子系统协方差矩阵特征值的矢量系数Ai的求解 \n[0057] 子系统协方差矩阵Pi按特征值分解为: \n[0058] \n[0059] 式中,Λi=diag{λi1,λi2,Λ,λin},λi1,λi2,L,λin为Pi的特征值,n为Pi阵的阶数。 \n[0060] 对Xi的每一个分量xij(表示第i个局部滤波器状态估计中的第j个分量)独立进行信息分配系数计算,信息分配系数为: \n[0061] \n[0062] i=1,2,L,N;j=1,2,L,n \n[0063] 式中,λij为状态变量xij对应的特征值,N为子系统个数。 \n[0064] Xi对应的信息分配系数为矩阵形式: \n[0065] \n[0066] 基于子系统可观测矩阵奇异值的矢量系数Υi的求解。 \n[0067] 设某时间段动态系统的可观测性矩阵为Q(Q∈Rp×q),对Q进行奇异值分解,得 [0068] Q=USVT (8) \n[0069] 式中:U=[u1 u2 L up],V=[v1 v2 L vq]都是正交矩阵;\n[0070] Λr×r=diag(σ1,σ2,L,σr),其中r,σi(i=1,2,L,r)分别为矩阵Q的秩和奇异值。对矩阵V进行分析,得到σi为与其对应的右奇异向量V中取得最大绝对值的状态变量所对应的奇异值。若σi的值较大,则相应的系统状态变量具有较好的观测性,可以获得较高精度的估计;若σi的值较小,则相应的系统状态变量可能会出现奇异,落入不可观测区间。 \n[0071] 对每个子系统的可观测性矩阵分别进行奇异值分解,进一步对Xi的每一个分量xij(表示第i个局部滤波器状态估计中的第j个分量)独立进行信息分配系数计算,信息分配系数为: \n[0072] \n[0073] i=1,2,L,N;j=1,2,L,n \n[0074] 式中,σij为状态变量xij对应的奇异值,N为子系统个数。 \n[0075] 则Xi对应的信息分配系数为矩阵形式: \n[0076] \n[0077] 融合子系统协方差矩阵特征值和可观测矩阵奇异值的矢量系数Υi的求解 [0078] \n[0079] 即 \n[0080] \n[0081] 信息分配系数满足信息守恒原理: \n[0082] \n[0083] 主滤波器信息分配: \n[0084] 根据矢量形式信息分配系数将系统的过程信息在各子滤波器之间的进行分配,表达式如下: \n[0085] \n[0086] 五、联邦主滤波器信息融合 \n[0087] 将各子滤波器估计信息进行融合,得到全局最优估计。 \n[0088] 融合算法为: \n[0089] \n[0090] 图2为仿真的一条飞行航迹,包括滑跑、加速、上升、滚转、转弯等机动方式。 [0091] 经过仿真得到,对相同的状态变量,不同的子系统获得了不同的信息分配系数;对同一个子系统,不同的状态变量也得到了不同的信息分配系数;分配系数曲线均随时间动态变化。该种矢量信息分配方案能够实时具体反映每个状态量的变化特性。 \n[0092] 图3~图5的仿真结果表明,该方法能够根据导航子系统中各个状态变量的精度特性和可观测特性进行实时自适应调节,滤波精度高,对飞行环境变化有更强的适应性。
法律信息
- 2022-08-19
未缴年费专利权终止
IPC(主分类): G01C 21/24
专利号: ZL 201110266447.X
申请日: 2011.09.09
授权公告日: 2013.08.21
- 2013-08-21
- 2012-03-28
实质审查的生效
IPC(主分类): G01C 21/24
专利申请号: 201110266447.X
申请日: 2011.09.09
- 2012-02-15
引用专利(该专利引用了哪些专利)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 |
1
| |
2010-10-20
|
2010-06-03
| | |
2
| |
2010-10-13
|
2010-05-28
| | |
3
| | 暂无 |
1996-11-20
| | |
4
| |
2008-05-14
|
2007-12-12
| | |
被引用专利(该专利被哪些专利引用)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 | 该专利没有被任何外部专利所引用! |