一种无人机信息融合定位方法与流程

[0001]
本发明涉及无人机定位技术领域，具体为一种无人机信息融合定位方法。


背景技术：

[0002]
无人机作为未来战场上重要的侦察工具和攻击武器受到越来越多的关注。无人机在执行任务时，其机载光电设备对无人机的定位精度有较高要求。目前对于无人机编队的精准定位已有一些解决方案；
[0003]
专利cn2020100538363提出一种基于激光雷达与定位向量匹配多无人机协同定位方法，但激光雷达作为有源器件受天候和大气的影响，若在大雨、下雪、浓雾等非晴朗气候条件下，红外线波的侦测能力会大幅衰減，此时存在检测概率低、虚警率高等问题。而基于捷联式惯性导航系统(sins)和全球定位系统(gps)的组合导航系统可对无人机观测数据进行初步定位，再结合无源测向交叉定位技术进行信息融合，最后解算出精准无人机位置信息，这种多源信息融合的定位方法能够使无人机编队在恶劣天气下工作，定位盲区更小，更加适合于对无人机进行定位。


技术实现要素：

[0004]
本发明的目的就在于为了解决现有无人机定位在非晴朗气候条件下，红外线波的侦测能力会大幅衰減，存在检测概率低、虚警率高的问题，而提出一种无人机信息融合定位方法，本发明可以使多个无人机在松组合导航系统和多机无源交叉测向定位下进行信息融合，从而在复杂气候环境中获得更好的定位效果。
[0005]
本发明的目的可以通过以下技术方案实现：一种无人机信息融合定位方法，该方法包括以下步骤：
[0006]
步骤一：通过测向交叉定位法和uav视觉传感器系统测得各无人机的相对位置信息；无人机为无人机；
[0007]
步骤二：通过sins/gps组合导航系统解算出无人机在地面坐标系中的初步位置坐标；具体为：在其中一时刻，通过sins/gps松组合导航系统实时地得到目标无人机的初步位置信息及姿态角信息(θ2γ2ψ)；
[0008]
步骤三：在多机交叉定位中解算出无人机位置坐标，结合sins/gps组合导航系统输出的无人机位置坐标分别作差计算，把得到的差值作为量测值，再经过kalman滤波，估算出sins/gps组合导航系统的误差并对其进行校正输出经过校正后目标无人机的位置坐标；其中，多机交叉定位的具体表现为：当有多架无人机时，各无人机的定位除组合导航自定位外，通过其余无人机实现，即将其余无人机作为测向基站，分别对某架目标无人机进行测向得到一组俯仰角和方位角；
[0009]
当多架无人机数量等于二时，则通过双机三维定位算法进行分析，具体分析过程为：以a(x1,y1,z1)、b(x2,y2,z2)两无人机为测向基站，分别对无人机目标c(x,y,z)进行测向，分别得到一组俯仰角和方位角，即为(α1,β1)和(α2,β2)，两者相互组合得到四组测量子
集(α1,β1,α2)、(α1,β1,β2)、(α2,β2,α1)和(α2,β2,β1)，其中α
i
(i＝1,2)为方位角，设定x轴为正北方向，β
i
(i＝1,2)为俯仰角，由方位角α1、α2和俯仰角β1、β2的三角关系可得：
[0010][0011]
将上述方程表示为矩阵形式：
[0012]
hx＝z
[0013]
其中
[0014][0015]
x＝(x,y,z)
t
[0016][0017]
由此可得目标无人机坐标为：
[0018]
x＝h-1
z
[0019]
选取上述方程(1)中3组进行组合解出目标无人机的位置。
[0020]
优选的，所述相对位置信息包括相对方位角和俯仰角；初步位置坐标包括经度l、纬度b和高度h。
[0021]
优选的，当多架无人机数量大于二时时，对于每个无人机均可得到n组方位角和俯仰角的观测值，通过两两交叉定位进行定位解算，选择其中一个无人机为测向主站，坐标为a(x0,y0,z0)；其它无人机作为辅站，坐标分别为a(x
i
,y
i
,z
i
)(i＝1,2,...,n-1)，然后通过主站与辅站两两组合对目标进行定位，由此得到n-1个目标估计位置，分别为c(x
i
,y
i
,z
i
)(i＝1,2,...,n-1)，通过聚类法该组数据进行处理，从而得到目标位置坐标的估计值，具体步骤为：
[0022]
s1：设每个样本自成一类，即建立n-1类：标号(0)表示聚类开始运算前的状态，样本即一个目标位置坐标；
[0023]
s2：计算各类间的距离，得到一个n-1维的方阵d
(0)
，初始时即为各样本间的距离；
[0024]
s3：设前一步聚类运算中已求得距离矩阵d
(n)
，n为逐次聚类的次数，则求d
(n)
中的最小元素，设其为和两类间的距离，则将与合并为一类由此建立新的分类：
[0025]
s4：计算合并后新类别间的距离，即计算与其他未发生合并的间的距离；通过重心法作为距离计算准则，得d
(n+1)
；重复计算及合并，直到距离矩阵中的最
小元素超过定位精度上限或所有元素都属于同一类则输出包含元素最多的类的重心作为最终输出结果，即目标位置坐标的估计值。
[0026]
优选的，所述kalman滤波的具体过程为：设随机线性离散系统的状态方程和量测方程分别是：
[0027][0028]
式中：t为迭代周期，φ
k,k-1
是k-1到k时刻的一步转移矩阵；γ
k-1
是系统噪声矩阵；h
k
是量测矩阵；v
k
是量测噪声序列；w
k-1
为k-1时刻的系统噪声矩阵；
[0029]
设定过程噪声和观测噪声有如下统计特性：
[0030][0031]
式中q
k
为系统噪声序列的方差矩阵；r
k
为量测噪声方差矩阵；
[0032]
线性离散系统的基本卡尔曼滤波(kf)方程；状态一步预测方程：
[0033][0034]
状态估计方程：
[0035][0036]
最优滤波增益方程：
[0037][0038]
最后得到预测均方误差方程：
[0039][0040]
上述是离散kf的基本方程，给出初始值和p0，再根据k时刻的量测值z
k
，递推出k时刻的状态估计值
[0041]
与现有技术相比，本发明的有益效果是：基于测向交叉定位法利用uav视觉传感器系统测得各无人机的相对位置信息，利用sins/gps组合导航系统解算出无人机在地面坐标系中的初步位置坐标；在多机交叉定位中解算出无人机位置坐标，结合sins/gps组合导航系统输出的无人机位置坐标分别作差计算，把得到的差值作为量测值，再经过kalman滤波，估算出sins/gps组合导航系统的误差并对其进行校正输出经过校正后目标无人机的位置坐标，可以使多个无人机在松组合导航系统和多机无源交叉测向定位下进行信息融合，从而在复杂气候环境中获得更好的定位效果。
附图说明
[0042]
为了便于本领域技术人员理解，下面结合附图对本发明作进一步的说明。
[0043]
图1为本发明的无人机定位示意图；
[0044]
图2为本发明的无人机定位流程图；
[0045]
图3为本发明的双机交叉定位示意图。
具体实施方式
[0046]
下面将结合实施例对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
[0047]
请参阅图1-3所示，一种无人机信息融合定位方法，无人机定位算法就是通过视觉传感器系统测得各无人机之间的相对方位角α、俯仰角β，结合sins/gps组合导航系统给出的目标无人机位置(经度λ、纬度及高度h)和姿态角(俯仰角θ、横滚角γ及航向角ψ)，利用齐次坐标转换法解算得到无人机在大地坐标系下的经度l、纬度b和高度h的过程；
[0048]
本方法包括三部分，第一部分是基于测向交叉定位法利用uav视觉传感器系统测得各无人机的相对位置信息(相对方位角、俯仰角)；第二部分是利用sins/gps组合导航系统解算出无人机在地面坐标系中的初步位置坐标(经度l、纬度b和高度h)；第三部分是在多机交叉定位中解算出无人机位置坐标，结合sins/gps组合导航系统输出的无人机位置坐标分别作差计算，把得到的差值作为量测值，再经过kalman滤波，估算出sins/gps组合导航系统的误差并对其进行校正输出经过校正后目标无人机的位置坐标；该方法的示意图和流程图如图1、图2所示；
[0049]
实施例1
[0050]
由于在本定位系统中主要关注的最终输出结果为无人机的位置坐标因而可对其余误差来源做简化处理，即略过无人机姿态角(θ,γ,ψ)的测量误差，同时对减振器振动角(δθ
ba
,δγ
ba
,δψ
ba
)、测向法测得无人机的相对位置信息(α,β)以及其余无人机测得的自身位置坐标、目标无人机相对方位角和俯仰角作如下近似假设：
[0051]
考虑到航迹位置信息可得到(α,β,r)的准确值，可假设无人机对目标的定位过程中，参数x(x＝α,β,r)的参数误差δx都服从于均值为准确值，方差是σ
δx
的正态分布；类似地，可假设参数x(x＝δθ
ba
,δγ
ba
,δψ
ba
)的参数误差δx都服从于均值为0，方差是σ
δx
的正态分布，基站方位角误差、俯仰角误差，测向基站的自身定位误差都是零均值的，并且标准差分别为σ
φ
、σ
ε
、σ
s
，并且可以认为测向基站自身定位的误差与基站所测的方位角误差、俯仰角误差之间都是相互独立的；具体流程如下：
[0052]
(1)在某一时刻，通过sins/gps松组合导航系统实时地得到目标无人机的初步位置信息及姿态角信息(θ,γ,ψ)；松组合方式就是分别获得sins与gps的量测信息，再分别进行做差运算，把得到的差值作为组合导航系统的量测值，再经过kalman滤波，估计出sins的误差，然后对sins进行校正；
[0053]
(2)通过多机交叉定位测得的目标无人机位置坐标对上述初步位置坐标进行校正，得到校正后的目标无人机位置坐标；多机交叉定位具体实现为当有多架无人机时，各无人机的定位除组合导航自定位外，还可通过其余无人机实现，即将其余无人机作为测向基
站，分别对某架目标无人机进行测向得到一组俯仰角和方位角；下面就双机三维定位算法进行分析；
[0054]
如图3所示，a(x1,y1,z1)、b(x2,y2,z2)两无人机为测向基站，分别对无人机目标c(x,y,z)进行测向，分别得到一组俯仰角和方位角为(α1,β1)和(α2,β2)它们相互组合可以得到四组测量子集(α1,β1,α2)、(α1,β1,β2)、(α2,β2,α1)和(α2,β2,β1)，其中α
i
(i＝1,2)为方位角(假设x轴为正北方向)，β
i
(i＝1,2)为俯仰角；由方位角α1、α2和俯仰角β1、β2的三角关系可得：
[0055][0056]
可将上述方程表示为矩阵形式：
[0057]
hx＝z
[0058]
其中
[0059][0060]
x＝(x,y,z)
t
[0061][0062]
由此可得目标无人机坐标为：
[0063]
x＝h-1
z
[0064]
易知，由上述4个方程任选3组就可以进行组合解出目标无人机的位置；当推广到n(n＞2)个无人机的情形时，对于每个目标均可得到n组方位角和俯仰角的观测值，在定位解算时我们采用两两交叉定位；选择其中一个无人机为测向主站，坐标为a(x0,y0,z0)；其它无人机作为辅站，坐标分别为a(x
i
,y
i
,z
i
)(i＝1,2,...,n-1)，然后通过主站与辅站两两组合对目标进行定位，由此可以得到n-1个目标估计位置，分别为c(x
i
,y
i
,z
i
)(i＝1,2,...,n-1)；采用聚类法对这组数据进行处理，从而得到目标位置坐标的估计值，具体流程为：设每个样本(即一个目标位置坐标)自成一类，即建立n-1类：(标号(0)表示聚类开始运算前的状态)，计算各类间的距离(初始时即为各样本间的距离)，得到一个n-1维的方阵d
(0)
；然后，假设前一步聚类运算中已求得距离矩阵d
(n)
(n为逐次聚类的次数)，则求d
(n)
中的最小元素，设其为和两类间的距离，则将与合并为一类由此建立新的分类：最后，计算合并后新类别间的距离，即计算与其他未发生合并的间的距离，此处我们采用重心法作为距离计算准则，得d
(n+1)
；重复计算及合并，直到距离矩阵d
(n)
中的最小元素超过定位精度上限或所有元素都属于同一类则
输出包含元素最多的类的重心作为最终输出结果，即目标位置坐标的估计值；
[0065]
(3)将sins/gps松组合定位和无源测向多机交叉定位结果进行作差计算，经过kalman滤波，估算出松组合导航系统定位差值并校正；具体kalman滤波在sins/gps组合导航系统和无源交叉定位过程为：
[0066]
设随机线性离散系统的状态方程和量测方程分别是：
[0067][0068]
式中：t为迭代周期；φ
k,k-1
是k-1到k时刻的一步转移矩阵；γ
k-1
是系统噪声矩阵；h
k
是量测矩阵；v
k
是量测噪声序列；w
k-1
为k-1时刻的系统噪声矩阵；假设系统的过程噪声和观测噪声有如下统计特性：
[0069][0070]
式中q
k
为系统噪声序列的方差矩阵；r
k
为量测噪声方差矩阵；下面给出线性离散系统的基本卡尔曼滤波(kf)方程；kf的基本思想是采用信号与噪声的状态空间模型，通过前一时刻的估计值和现在时刻的观测值来更新状态变量，进而计算出现在时刻的估计值；状态一步预测方程：
[0071][0072]
状态估计方程：
[0073][0074]
最优滤波增益方程：
[0075][0076]
最后得到预测均方误差方程：
[0077][0078]
上述是离散kf的基本方程，只要给出初始值和p0，再根据k时刻的量测值z
k
，就可以递推出k时刻的状态估计值
[0079]
本发明在使用时，基于测向交叉定位法利用uav视觉传感器系统测得各无人机的相对位置信息，利用sins/gps组合导航系统解算出无人机在地面坐标系中的初步位置坐标；在多机交叉定位中解算出无人机位置坐标，结合sins/gps组合导航系统输出的无人机位置坐标分别作差计算，把得到的差值作为量测值，再经过kalman滤波，估算出sins/gps组合导航系统的误差并对其进行校正输出经过校正后目标无人机的位置坐标，可以使多个无人机在松组合导航系统和多机无源交叉测向定位下进行信息融合，从而在复杂气候环境中获得更好的定位效果。
[0080]
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽
叙述所有的细节，也不限制该发明仅为的具体实施方式。显然，根据本说明书的内容，可作很多的修改和变化。本说明书选取并具体描述这些实施例，是为了更好地解释本发明的原理和实际应用，从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。