一种多传感器融合的无人机复杂天候飞行控制方法

1.本发明涉及无人机飞行控制系统的技术领域。


背景技术：

2.当前世界涌起了发展无人机的热潮，无论是军用领域，还是民用领域，无人机的地位愈加重要，应用范围也愈加广泛。然而目前的无人机在复杂天候情况下会出现诸多不利于安全航行的问题，例如，降雨情况下无人机的各个舵面、机翼、尾翼的气动性能会受到影响，其飞行稳定性与操纵性会出现恶化；结冰情况下无人机外形尤其是机翼会发生改变，其甚至可能造成机翼升力的骤然减少，引发重大安全事故。
3.目前对于无人机在天候情况较为复杂时，往往因为各种因素的干扰而无法得到正确的位置信息。例如，在降雨、阵风等影响下无人机的可能会出现掉高现象，但因为天候影响，单一传感器反馈得到的无人机高度信息存在不准确，误差较大情况，从而便会影响到对应无人机的高度控制。因此研究复杂天候情况下无人机的飞行控制对于无人机的应用发展具有重要意义。


技术实现要素：

4.本发明的目的是针对目前无人机在复杂天候条件下无法得到正确的高度信息，使其在复杂天候条件下出现飞行控制困难的问题等缺陷，提供一种可在复杂天候情况下进行稳定的高度控制的无人机飞行控制方法。
5.本发明的技术方案如下：
6.一种多传感器融合的无人机复杂天候飞行控制方法，其包括：
7.s1根据地轴系与机体轴系转换关系建立无人机运动学方程；
8.s2建立无人机中与高度值相关的多个传感器的观测模型；
9.s3根据所述观测模型及所述无人机运动学方程，建立无人机高度状态方程和量测方程；
10.s4基于所述无人机高度状态方程和量测方程，通过卡尔曼滤波法得到多个传感器高度数据融合后的无人机飞行高度值；
11.s5构建无人机高度控制回路，根据所得数据融合后的飞行高度值调控无人机飞行高度。
12.根据本发明的一些优选实施方式，所述观测模型包括大气数据计算机观测模型、无线电高度表观测模型及差分gps观测模型。
13.根据本发明的一些优选实施方式，所述无人机运动学方程建立如下：
[0014][0015][0016]
[0017]
u＝v
0 cosαcosβ
[0018]
v＝v
0 sinβ
[0019]
w＝v
0 cosβsinα
[0020]
其中，地轴系为ox
e
y
e
z
e
，为地轴系下沿x轴、y轴、z轴无人机速度，u,v,w分别为无人机速度沿机体轴系x轴、y轴和z轴的分量；θ,φ,ψ分别为俯仰角，滚转角和偏航角；v0,α,β分别为无人机速度，无人机迎角和无人机侧滑角。
[0021]
根据本发明的一些优选实施方式，所述观测模型包括：
[0022]
如下的大气数据计算机观测模型：
[0023]
h1＝h+b
h
+v1，
[0024]
其中，h1为大气数据计算机测量的高度值，h为无人机真实高度，b
h
为大气数据计算机测高常值偏差，v1为大气数据计算机观测方差；
[0025]
如下的无线电高度表观测模型：
[0026]
h2＝h+v2，
[0027]
其中，h2为无线电高度表测量的高度值，v2为无线电高度表的观测方差；
[0028]
如下的差分gps的观测模型：
[0029]
h3＝h+v3，
[0030]
其中，h3为差分gps测量的高度值，v3为差分gps的观测方差。
[0031]
根据本发明的一些优选实施方式，所述状态方程如下：
[0032]
x＝[h b
h
]
t
[0033]
所述量测方程如下：
[0034]
z＝[h
1 h
2 h3]
t
，且：
[0035][0036][0037]
其中，w1、w2为系统噪声。
[0038]
根据本发明的一些优选实施方式，所述卡尔曼滤波法中，卡尔曼滤波器的时间更新方程为：
[0039][0040][0041]
p
(k+1|k)
＝p
(k|k)
+q
(k|k)
，
[0042]
和/或，
[0043]
量测更新方程为：
[0044]
[0045][0046]
p
(k+1)
＝(i
‑
k
(k+1)
c
(k+1)
)p
(k+1|k)
[0047]
其中：
[0048][0049][0050][0051]
其中，k
(k+1)
为卡尔曼增益；c
(k+1)
为观测输出矩阵；p
(k+1)
为估计均方差矩阵；q
(k+1)
为系统噪声的协方差矩阵；r
(k+1)
为量测噪声的协方差矩阵。
[0052]
根据本发明的一些优选实施方式，所述高度控制回路采用pid控制。
[0053]
根据本发明的一些优选实施方式，所述pid控制中，以期望飞行高度作为输入，以所述数据融合后的飞行高度值作为反馈。
[0054]
本发明具备以下有益效果：
[0055]
本发明所提供的一种多传感器融合的无人机复杂天候飞行控制方法，首先通过卡尔曼滤波方法融合多传感器测得高度数据，从而得到复杂天候情况下无人机真实飞行高度，再通过控制系统对无人机进行高度控制使得无人机能够在复杂天候环境下进行安全稳定飞行。该发明所提出的方法相比于现有方法，优势在于能够在复杂天候情况下通过多传感器对无人机进行高度测量，而不仅仅是对单一数据，通过数据融合便可以消除因为复杂天候带来的高度测量不准确问题，得到可靠的无人机飞行高度数据，然后再通过控制系统便可以控制无人机达到期望飞行高度，为无人机在复杂天候情况下安全飞行提供参考。
附图说明
[0056]
图1为一种具体的复杂天候下无人机飞行控制方法的流程示意图。
[0057]
图2为一种具体的无人机高度控制回路。
[0058]
图3为实施例中真实高度与大气数据计算机测量结果对比曲线。
[0059]
图4为实施例中真实高度与无线电高度表测量结果对比曲线。
[0060]
图5为实施例中真实高度与差分gps测量结果对比曲线。
[0061]
图6为实施例中真实高度与多传感器数据融合结果对比曲线。
具体实施方式
[0062]
以下结合实施例和附图对本发明进行详细描述，但需要理解的是，所述实施例和附图仅用于对本发明进行示例性的描述，而并不能对本发明的保护范围构成任何限制。所有包含在本发明的发明宗旨范围内的合理的变换和组合均落入本发明的保护范围。
[0063]
参照图1，以某型无人机在复杂天候下飞行情况为例，一种多传感器融合的无人机复杂天候飞行控制方法，包括以下步骤：
[0064]
第一步，根据地轴系与机体轴系转换关系建立无人机运动学方程，如下所示：
[0065]
[0066][0067][0068]
u＝v0cosαcosβ
[0069]
v＝v0sinβ
[0070]
w＝v0cosβsinα
[0071]
其中，地轴系为ox
e
y
e
z
e
、是固定于地球中心随地球一起转动的坐标系，原点o
e
即为地球中心；定义体轴系ox
b
y
b
z
b
、是固连于无人机飞行器上的坐标系，原点o
b
取无人机质心，ox
b
平行于飞行器的纵向中心线并指向运动方向，oy
b
垂直于ox
b
z
b
平面并指向飞行器右侧，oz
b
位于对称平面指向下方，构成右手系，z
e
为无人机纵向位移，u,v,w分别为无人机速度沿x轴、y轴和z轴的分量；θ,φ,ψ分别为俯仰角，滚转角和偏航角；v0,α,β分别为无人机速度，无人机迎角和无人机侧滑角。
[0072]
第二步，建立无人机中各高度传感器的观测模型，其中，所述观测模型包括：大气数据计算机观测模型、无线电高度表观测模型及差分gps观测模型；
[0073]
更具体的，各观测模型可包括：
[0074]
大气数据计算机的观测模型：
[0075]
h1＝h+b
h
+v1[0076]
其中，h1为大气数据计算机测量的高度值，h为无人机真实高度，b
h
为大气数据计算机测高常值偏差，v1为大气数据计算机观测方差。
[0077]
无线电高度表的观测模型：
[0078]
h2＝h+v2[0079]
其中，h2为无线电高度表测量的高度值，v2为无线电高度表的观测方差。
[0080]
差分gps的观测模型：
[0081]
h3＝h+v3[0082]
其中，h3为差分gps测量的高度值，v3为差分gps的观测方差。
[0083]
第三步，根据所述观测模型及所述无人机运动学方程，建立无人机高度状态方程和量测方程；
[0084]
其中，所述无人机高度状态方程的建立可具体如：
[0085]
将与无人机纵向位移z
e
的关系为h＝
‑
z
e
无人机真实高度h代入无人机纵向运动方程可得：
[0086][0087]
在三个传感器都正常工作情况下，系统的状态方程和量测方程如下所示：
[0088][0089]
z＝cx+v
[0090]
其中系统状态x＝[h b
h
]
t
，系统观测量z＝[h
1 h
2 h3]
t
，a为系统状态矩阵，c为量测矩阵，w为系统噪声矩阵，v为量测噪声矩阵，
[0091]
将其中各项代入，可以得到状态方程和量测方程具体形式如下所示：
[0092][0093][0094]
第四步，设计kalman滤波器，得到在多传感器数据融合下的无人机飞行高度值；
[0095]
更具体的，根据卡尔曼滤波法，可以得到时间更新方程为：
[0096][0097][0098]
p
(k+1|k)
＝p
(k|k)
+q
(k|k)
[0099]
量测更新方程为：
[0100][0101][0102]
p
(k+1)
＝(i
‑
k
(k+1)
c
(k+1)
)p
(k+1|k)
[0103]
其中：
[0104][0105][0106][0107]
其中，k
(k+1)
为卡尔曼增益；c
(k+1)
为传感器观测输出矩阵；p
(k+1)
为估计均方差矩阵；q
(k+1)
为系统噪声的协方差矩阵；r
(k+1)
为量测噪声的协方差矩阵。
[0108]
由此可以得到无人机飞行高度的无偏估计值该值与无人机真实高度h非常接近，可以认为相等。
[0109]
第五步，构建无人机高度控制回路，根据所得飞行高度值调控无人机飞行高度。
[0110]
其中，更具体的，参照图2，所述高度控制回路可采用pid控制方法，以期望飞行高度作为输入，所得多传感器数据融合的真实高度值作为反馈，调节无人机飞行高使无人机真实飞行高度达到输入的期望飞行高度，最终完成复杂天候条件下无人机的飞行高度控制。
[0111]
在一些具体实施例中，根据上述具体实施方式所得真实高度与大气数据计算机测量结果曲线如图3所示，真实高度与无线电高度表测量结果曲线如图4所示，真实高度与差分gps测量结果曲线如图5所示，真实高度与多传感器数据融合结果曲线如图6所示，可以看出，通过本发明的多传感器融合的方式所得的高度数据与真实高度最为一致。
[0112]
以上实施例仅用于展示本发明的部分实施方式，本发明的保护范围并不仅局限于上述实施例。凡属于本发明思路下的技术方案均属于本发明的保护范围。应该指出，对于本
技术领域的普通技术人员来说，在不脱离本发明原理的前提下的改进和润饰，这些改进和润饰也应视为本发明的保护范围。