CN112666978B - 一种无人机自适应着陆导航的方法及装置 - Google Patents

一种无人机自适应着陆导航的方法及装置 Download PDF

Info

Publication number
CN112666978B
CN112666978B CN202011561067.4A CN202011561067A CN112666978B CN 112666978 B CN112666978 B CN 112666978B CN 202011561067 A CN202011561067 A CN 202011561067A CN 112666978 B CN112666978 B CN 112666978B
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
fitting
fitting parameter
landing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011561067.4A
Other languages
English (en)
Other versions
CN112666978A (zh
Inventor
张奇智
雷瑛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Northern Sky Long Hawk Uav Technology Co ltd
Original Assignee
Beijing Northern Sky Long Hawk Uav Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Northern Sky Long Hawk Uav Technology Co ltd filed Critical Beijing Northern Sky Long Hawk Uav Technology Co ltd
Priority to CN202011561067.4A priority Critical patent/CN112666978B/zh
Publication of CN112666978A publication Critical patent/CN112666978A/zh
Application granted granted Critical
Publication of CN112666978B publication Critical patent/CN112666978B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提供了一种无人机的自适应着陆导航的方法及装置,其中,该方法包括:确定下滑阶段的给定高度拟合曲线Hg(x)公式;确定Hg(x)的一阶微分求导方程;获取无人机进入下滑窗口时允许的飞行高度Ha和下滑窗口距离理想着陆点的水平待飞距离x0;获取着陆接地时允许的飞行高度Hb和着陆接地时距离理想着陆点的水平距离x1;将Ha、x0、Hb和x1代入上述公式建立航线拟合参数组的方程组;求解航线拟合参数组的方程组得到航线拟合参数;将航线拟合参数代入Hg(x)得到符合机场和空管以及无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。本发明中无人机可适应不同机场要求,在保证接地升降速度的同时保证着陆的位置精度。

Description

一种无人机自适应着陆导航的方法及装置
技术领域
本发明涉及一种飞行控制领域,尤其涉及一种无人机自适应着陆导航的方法及装置。
背景技术
受到地理环境(如机场外高楼)、空管(如有其它无人机在同一空域有飞行任务,不允许在特定高度着陆)等限制,无人机在着陆前的空域会有不同的要求,通常横侧向的条件比较宽泛,不存在规划难题,主要是受飞行高度的影响,无人机只能限制在某一高度下飞行。对于飞行导航与控制来讲,在不通过操作员人工干预的情况下很难实现可以应对不同着陆高度要求的导航算法对飞行进行引导并完成定点着陆。
目前普遍的无人机在着陆下滑的导航设计都采用固定下滑角度高度控制的控制方式,这种方法逻辑较为简单,主要通过高度控制保证着陆位置精度。也有少数采取下滑+拉平的控制方式,在下滑阶段的控制方式与上述方式相同,在下滑阶段后加入拉平阶段,牺牲位置的控制精度对升降速率进行调整。如图1所示,具体方式为:K为空管允许的下滑初始点,B点为退出下滑阶段进入拉平阶段的转段点,H为进入拉平阶段时无人机的飞行高度,G为理想着陆点,也是地面坐标系的原点。地面坐标系中,X轴为水平位置坐标轴,向右为正(无人机由K窗口向G点着陆,待飞距离=(0-无人机水平坐标位置),该值随无人机下滑而逐渐减小),Y轴为高度坐标轴,向上为正。无人机在KB段为下滑阶段,纵向控制为高度控制,KB的连线即为给定下滑线,λ为下滑阶段的下滑角,即KB延长线与X轴的夹角。无人机在BG段为拉平阶段,纵向控制为升降速度控制,防止接地时的升降速率过大导致起落架损伤。此段BG的连线不是给定轨迹曲线,可能在实际飞行中受到风干扰影响向发生改变。
根据场地需求以某一下滑角度直线下滑的算法,可能会因下滑角度陡峭导致接地时下降速率过大(设计给定下滑线时考虑无风状态,下降速率=真空速*下滑角度的正弦值,下降速率和下滑角度成正比),但对于带有起落架采取轮式起降的无人机来说,接地时升降速率过大对于起落架的损伤较大,会减少使用寿命。少数无人机在着陆前可通过图1中拉平的方式增大俯仰角保证着陆姿态并同时减小着陆升降速率,但是拉平阶段的纵向控制需要放弃高度控制转为升降速率控制,无人机不按给定航线飞行,在受到常值风干扰的情况下无法保证着陆位置精度,最终导致着陆点的位置与理想接地点G相差较大。
因此,如何在满足不同的机场或空域的限定条件的情况下自适应地生成着陆航线,如何在保证升降速度的同时保证着陆的位置精准等问题成为本领域技术人员亟待解决的问题。
发明内容
本发明旨在解决上述问题之一。
本发明的主要目的在于提供一种无人机自适应着陆导航的方法。
本发明的另一目的在于提供一种无人机自适应着陆导航的装置。
为达到上述目的,本发明的技术方案具体是这样实现的:
本发明一方面提供了一种无人机的自适应着陆导航的方法,包括:
步骤1、以无人机当前位置的待飞距离为自变量确定下滑阶段的给定高度拟合曲线Hg(x):
Hg(x)=a1+a2x+a3x2+a4x3公式(1),
其中,x表示所述无人机的当前位置距离理想着陆点的水平待飞距离,Hg(x)表示所述无人机的当前位置距离理想着陆点的水平待飞距离为x时的飞行高度,a1为第一航线拟合参数,a2为第二航线拟合参数,a3为第三航线拟合参数,a4为第四航线拟合参数;
步骤2、确定所述下滑阶段的给定高度拟合曲线Hg(x)的一阶微分求导方程:
Figure BDA0002859408380000021
其中,
Figure BDA0002859408380000022
表示所述无人机在任意飞行高度处的曲线斜率;
步骤3、根据机场和空管的要求获取所述无人机进入下滑窗口时允许的飞行高度Ha的值和所述下滑窗口距离理想着陆点的水平待飞距离x0的值;根据所述无人机的机型获取着陆接地时允许的飞行高度Hb的值和着陆接地时距离理想着陆点的水平距离x1的值;
步骤4、将所述Ha的值、所述x0的值、所述Hb的值和所述x1的值依次代入所述公式(1)和所述公式(2)中,建立航线拟合参数组的方程组:
Figure BDA0002859408380000023
步骤5、求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4
步骤6、将求解得到的所述第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值代入所述公式(1)得到符合所述机场和空管以及所述无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。
可选的,所述着陆接地时距离理想着陆点的水平距离x1的值为0;
所述无人机着陆接地时允许的飞行高度Hb为所述无人机着陆时重心相对地面的高度。
可选的,当Hg=Ha,x=x0时,
Figure BDA0002859408380000031
的值为0;当Hg=Hb,x=x1时,
Figure BDA0002859408380000032
其中,Hyd表示所述无人机着陆时期望的升降速率,Vx表示所述无人机的水平飞行速度。
可选的,所述无人机的水平飞行速度Vx通过以下公式计算得到:
Figure BDA0002859408380000033
其中,Vg表示所述无人机的给定真空速,所述Hyd的值小于所述无人机的落架能承受的最大升降速率。
可选的,所述求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4,包括:
设矩阵A为:
Figure BDA0002859408380000034
(’代表转置),
矩阵B为:
Figure BDA0002859408380000035
设C矩阵为求解矩阵,C=inv(B)*A,得出的C矩阵四个值分别为所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4的值
本发明另一方面提供了一种无人机的自适应着陆导航的装置,包括:
给定高度拟合曲线构建模块,用于以无人机当前位置的待飞距离为自变量确定下滑阶段的给定高度拟合曲线Hg(x):
Hg(x)=a1+a2x+a3x2+a4x3公式(1),
其中,x表示所述无人机的当前位置距离理想着陆点的水平待飞距离,Hg(x)表示所述无人机的当前位置距离理想着陆点的水平待飞距离为x时的飞行高度,a1为第一航线拟合参数,a2为第二航线拟合参数,a3为第三航线拟合参数,a4为第四航线拟合参数;
所述给定高度拟合曲线构建模块,还用于确定所述下滑阶段的给定高度拟合曲线Hg(x)的一阶微分求导方程:
Figure BDA0002859408380000041
其中,
Figure BDA0002859408380000042
表示所述无人机在任意飞行高度处的曲线斜率;
获取模块,用于根据机场和空管的要求获取所述无人机进入下滑窗口时允许的飞行高度Ha的值和所述下滑窗口距离理想着陆点的水平待飞距离x0的值;根据所述无人机的机型获取着陆接地时允许的飞行高度Hb的值和着陆接地时距离理想着陆点的水平距离x1的值;
求解模块,用于将所述Ha的值、所述x0的值、所述Hb的值和所述x1的值依次代入所述公式(1)和所述公式(2)中,建立航线拟合参数组的方程组:
Figure BDA0002859408380000043
所述求解模块,还用于求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4;将求解得到的所述第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值代入所述公式(1)得到符合所述机场和空管以及所述无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。
可选的,所述着陆接地时距离理想着陆点的水平距离x1的值为0;
所述无人机着陆接地时允许的飞行高度Hb为所述无人机着陆时重心相对地面的高度。
可选的,当Hg=Ha,x=x0时,
Figure BDA0002859408380000044
的值为0;当Hg=Hb,x=x1时,
Figure BDA0002859408380000045
其中,Hyd表示所述无人机着陆时期望的升降速率,Vx表示所述无人机的水平飞行速度。
可选的,所述无人机的水平飞行速度Vx通过以下公式计算得到:
Figure BDA0002859408380000046
其中,Vg表示所述无人机的给定真空速,所述Hyd的值小于所述无人机的落架能承受的最大升降速率。
可选的,所述求解模块通过以下方式求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4
设矩阵A为:
Figure BDA0002859408380000047
(’代表转置),
矩阵B为:
Figure BDA0002859408380000051
设C矩阵为求解矩阵,C=inv(B)*A,得出的C矩阵四个值分别为所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4的值。
由上述本发明提供的技术方案可以看出,本发明提供了一种无人机自适应着陆导航的方法及装置,可以根据不同的机场或空管的适航要求以及各种无人机机型的适航要求自适应地拟合出最优的给定高度拟合曲线,保证了无人机可在不同机场适应不同的下滑窗口要求,并在保证接地升降速度的同时保证着陆的位置精度。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为现有技术提供的无人机着陆轨迹设计的示意图;
图2为本发明实施例1提供的无人机自适应着陆导航的方法的流程图;
图3为本发明实施例1提供的着陆导航航线示意图;
图4为本发明实施例1提供的一个应用示例中求解航线拟合参数后的给定高度拟合曲线和实际飞行高度的示意图;
图5为本发明实施例1提供的一个应用示例中求解航线拟合参数后的给定高度拟合曲线中的升降速度变化曲线示意图;
图6为本发明实施例1提供的无人机自适应着陆导航的装置的结构示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或数量或位置。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
本实施例中,面对不同的空管和场地限制的要求以及各种机型的无人机着陆时的位置精度要求和升降速率,在设计着陆轨迹时,主要针对以下几个方面考虑:(1)进入下滑航线前的飞行高度要求;(2)进入下滑航线时距离着陆点的距离要求(3)着陆接地时无人机的升降速率要求;(4)着陆接地时的着陆位置精度要求(定点着陆,着陆时无人机距离理想着陆点前后不能超过50m)。其中,第(1)点与第(2)点的要求通常由机场当地的空管提出;第(3)点中,无人机的升降速率对于不同的无人机机型和起落架强度来说有不同的要求,通常根据起落架设备的承受能力或其他飞行控制需求而定;第(4)点中,着陆时的位置精度通常要求越高越好,尽可能实现定点着陆。
因此,本实施例通过以下思路来确定下滑阶段的给定高度拟合曲线,即给定高度拟合曲线:设计一套参数可自适应机场场地和空管要求的给定高度拟合曲线Hg(x),在保证可以按照机场空管要求高度和待飞距离进入下滑窗口的同时,接地时的着陆精度和下降速率也符合设计需求。
由于无人机在进入下滑阶段时通常通过调整俯仰角使无人机跟随给定高度的方式保证无人机跟随下滑线飞行,因此需要设计给定高度的导航算法。
因此,本发明实施例提供了一种无人机自适应着陆导航的方法。如图2所示,该无人机自适应着陆导航的方法,包括以下步骤(S101-S106):
S101、以无人机当前位置的待飞距离为自变量确定下滑阶段的给定高度拟合曲线Hg(x):
Hg(x)=a1+a2x+a3x2+a4x3公式(1)
其中,x表示无人机的当前位置距离理想着陆点的水平待飞距离,Hg(x)表示无人机的当前位置距离理想着陆点的水平待飞距离为x时的飞行高度,a1为第一航线拟合参数,a2为第二航线拟合参数,a3为第三航线拟合参数,a4为第四航线拟合参数;理想着陆点即为图1中的G点。
本实施例中,无人机在进入下滑阶段着陆轨迹的给定高度值可以通过公式(1)表示,通过公式(1)中x和Hg一一对应的关系保证任意时刻无人机的位置精度。
S102、确定下滑阶段的给定高度拟合曲线Hg(x)的一阶微分求导方程:
Figure BDA0002859408380000071
公式(2)中对x进行求导的意义是:针对公式(1)中给定高度拟合曲线的给定高度Hg
Figure BDA0002859408380000072
的值表示在任意高度给定值“Hg”处的曲线斜率,即该点所对应下滑角度的正切值,反应了给定高度拟合曲线中无人机实时的下滑角度λ。
S103、根据机场和空管的要求获取无人机进入下滑窗口时允许的飞行高度Ha的值和下滑窗口距离理想着陆点的水平待飞距离x0的值;根据无人机的机型获取着陆接地时允许的飞行高度Hb的值和着陆接地时距离理想着陆点的水平距离x1的值;
本实施例中,根据机场和空管的要求,不同的机场和空管允许的下滑窗口的Ha的值和x0的值不同,需根据实际飞入的机场要求实时获取。此外,Hb的值和x1的值对于不同的无人机机型和起落架强度来说有不同的要求,通常根据起落架设备的承受能力或其他飞行控制需求而定。其中,作为本实施例中一种可选的实施方式,无人机着陆接地时允许的飞行高度Hb可以为无人机着陆时重心相对地面的高度,一般来说,无人机着陆时重心相对地面的高度略高于起落架离地高度。该数值与飞行器设计相关,不同机型的重心对地高度不同,通常在1m以上。
S104、将Ha的值、x0的值、Hb的值和x1的值依次代入公式(1)和公式(2)中,建立航线拟合参数组的方程组:
Figure BDA0002859408380000073
作为本实施例中的一种可选实施方式,在无人机进入下滑窗口时,Hg=Ha,x=x0,其中,Ha为进入下滑窗口时允许的飞行高度,x0为下滑窗口距离理想着陆点的水平待飞距离,由于机场和空管的要求不同,相应的,Ha和x0的取值不同。本实施例中,通过选取符合机场和空管要求的Ha和x0的取值,可以保证下滑阶段的给定高度拟合曲线Hg(x)可以按照机场和空管要求的高度和待飞距离进入下滑窗口。
由于无人机进入下滑窗口的点与无人机退出定高阶段(进入下滑窗口之前为定高平飞阶段)的点为同一点,可视为进入下滑阶段的初始角度等于定高平飞的下滑角度(0°),即下滑曲线的斜率为0,因此,上述航线拟合参数组的方程组中,
Figure BDA0002859408380000081
的值为0。
作为本实施例中的一种可选实施方式,无人机在着陆接地时,Hg=Hb,x=x1,其中,Hb为无人机着陆接地时的飞行高度,x1为无人机着陆接地时距离理想着陆点的水平距离的值。本实施例中,理想着陆点为G点(参见图1和图3),G点为本算法设计中地面坐标系原点,无人机在着陆接地时,当前位置在水平方向上与原点重合,因此,x1的值为0。由于无人机机型和起落架强度来说有不同的要求,相应的,Hb的取值通常根据起落架设备的承受能力或其他飞行控制需求而定,也会不同。本实施例中,通过选取符合无人机机型和起落架强度要求的Hb的取值,可以保证无人机在按照下滑阶段的给定高度拟合曲线Hg(x)进行着陆接地时的着陆精度和下降速率也符合设计需求。
参见图3,针对G点,
Figure BDA0002859408380000082
的值为着陆时期望的下滑线斜率,对应的下滑角λ也反应了无人机在该点飞行速度的方向,可通过给定真空速Vg与升降速度Hyd计算所得。作为本实施例中的一种可选方式,具体计算方法如下:
Figure BDA0002859408380000083
的值为下滑角λ的正切值,从数值上等于无人机在该点升降速率与水平方向速度的比值
Figure BDA0002859408380000084
(即给定真空速Vg在地面坐标系中纵向分量与水平分量的比值),其中,Hyd为无人机着陆时期望的升降速率,该数值需根据起落架强度考虑取值,通常要略低于起落架所能承受的最大升降速率(如起落架允许以最大1.5m/s的升降速率着陆,则该值可设置为1~1.2m/s)。Vx为无人机的水平飞行速度,作为本实施例中一种可选的所述方式,该值可以通过将给定真空速Vg在地面坐标系正交分解后勾股定理计算所得:
Figure BDA0002859408380000085
S105、求解航线拟合参数组的方程组,得到第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4
本实施例中,在步骤S104中将符合机场和空管要求的Ha和x0、符合无人机机型以及起落架设备的承受能力的Hb和x1以及
Figure BDA0002859408380000091
代入航线拟合参数组的方程组后,就可以求解到航线拟合参数(a1、a2、a3和a4),而且,求解到的航线拟合参数(a1、a2、a3和a4)可自适应地生成符合机场场地和空管要求以及满足着陆精度和下降速率要求的给定高度拟合曲线Hg(x),在保证可以按照要求高度和待飞距离进入下滑窗口的同时,接地时的着陆精度和下降速率也符合设计需求。
下面,作为本实施例一种可选的实施方式,本实施例提供一种具体的求解航线拟合参数组的方程组,得到第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的实现方式,但本发明并不限于此,本发明对于求解方程组的方式不做限制。
根据该航线拟合参数组的方程组的构成和特点,写成矩阵计算的形式便于计算。
设矩阵A为:
Figure BDA0002859408380000092
(’代表转置)。
矩阵B为:
Figure BDA0002859408380000093
设C矩阵为求解矩阵,C=inv(B)*A,得出的C矩阵四个值分别为a1、a2、a3和a4的值。
以某长航时无人机为例,无人机的重心距离高度Hb为1.02m,其给定真空速Vg为30.55m/s,起落架允许最大升降速率为1.2m/s,故设计无人机着陆时期望的升降速度Hyd为-0.6m/s。针对某机场的适航空管要求,下滑窗口距离理想着陆点水平距离x0为4800m,高度差Ha约为310m。
根据上述情况,Ha=310m,
Figure BDA0002859408380000094
Hb=1m,
Figure BDA0002859408380000095
由于待飞距离x的值相对较大,将公式中待飞距离的量采取0.01的比例尺处理。
因此:矩阵A为:[310 0 1-0.0117]’
矩阵B为:
Figure BDA0002859408380000096
通过计算求得:a1=1,a2=-0.0117,a3=0.4028,a4=-0.0056。
S106、将求解得到的第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值代入公式(1)得到符合机场和空管以及无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。
具体地,将上述具体方式中求解得到的a1=1,a2=-0.0117,a3=0.4028,a4=-0.0056代入公式(1)得到无人机下滑阶段的给定高度拟合曲线如下:
Hg(x)=1-0.0117x+0.4028x2-0.0056x3
参见图4和图5,图4为求解航线拟合参数后的给定高度拟合曲线和实际飞行高度的示意图。图5为该给定高度拟合曲线中的升降速度变化曲线示意图。
图4中横坐标为飞行距离,纵坐标为飞行高度。从实际效果来看,在使用本实施例提供的无人机自适应着陆导航的方法对无人机进行引导后,可以保证精准的着陆位置的精确性,当无人机的飞行高度为1.023m时(即无人机重心对地高度为1.023m时),飞行距离为4800m,无人机在理想着陆点着陆,实现了精准着陆控制。图5中展现了升降速度变化趋势,对比图1可知,下滑曲线前半段下滑角逐渐增大,因此,升降速度较大,后半段临近接地,下滑角逐渐减小,因此升降速度逐渐减小,最终接地时升降速度为-0.5843m/s,满足设计要求,随后无人机转入滑跑段,升降速度为0m/s。
通过本发明实施例提供的无人机自适应着陆导航的方法,可以根据不同的机场或空管的适航要求以及各种无人机机型的适航要求自适应地拟合出最优的着陆导航航线,保证了无人机可在不同机场适应不同的下滑窗口要求,并在保证接地升降速度的同时保证着陆的位置精度。
实施例2
本实施例提供了一种无人机的自适应着陆导航的装置。该装置可以安装于无人机内,可以实施实施例1中提供的无人机的自适应着陆导航的方法。本实施例仅对该无人机的自适应着陆导航的装置的结构做简要描述,其他未尽事宜请参见实施例1中的描述。
如图6所示,该无人机的自适应着陆导航的装置包括:给定高度拟合曲线构建模块、获取模块和求解模块,其中:
给定高度拟合曲线构建模块,用于以无人机当前位置的待飞距离为自变量确定下滑阶段的给定高度拟合曲线Hg(x):
Hg(x)=a1+a2x+a3x2+a4x3公式(1),
其中,x表示无人机的当前位置距离理想着陆点的水平待飞距离,Hg(x)表示无人机的当前位置距离理想着陆点的水平待飞距离为x时的飞行高度,a1为第一航线拟合参数,a2为第二航线拟合参数,a3为第三航线拟合参数,a4为第四航线拟合参数;
给定高度拟合曲线构建模块,还用于确定下滑阶段的给定高度拟合曲线Hg(x)的一阶微分求导方程:
Figure BDA0002859408380000111
其中,
Figure BDA0002859408380000112
表示无人机在任意飞行高度处的曲线斜率;
获取模块,用于根据机场和空管的要求获取无人机进入下滑窗口时允许的飞行高度Ha的值和下滑窗口距离理想着陆点的水平待飞距离x0的值;根据无人机的机型获取着陆接地时允许的飞行高度Hb的值和着陆接地时距离理想着陆点的水平距离x1的值;
求解模块,用于将Ha的值、x0的值、Hb的值和x1的值依次代入公式(1)和公式(2)中,建立航线拟合参数组的方程组:
Figure BDA0002859408380000113
求解模块,还用于求解航线拟合参数组的方程组,得到第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4;将求解得到的第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值代入公式(1)得到符合机场和空管以及无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。
作为本实施例中的一种可选实施方式,着陆接地时距离理想着陆点的水平距离x1的值为0;
无人机着陆接地时允许的飞行高度Hb为无人机着陆时重心相对地面的高度。
作为本实施例中的一种可选实施方式,当Hg=Ha,x=x0时,
Figure BDA0002859408380000114
的值为0;当Hg=Hb,x=x1时,
Figure BDA0002859408380000115
其中,Hyd表示无人机着陆时期望的升降速率,Vx表示无人机的水平飞行速度。
作为本实施例中的一种可选实施方式,无人机的水平飞行速度Vx通过以下公式计算得到:
Figure BDA0002859408380000116
其中,Vg表示无人机的给定真空速,Hyd的值为1~1.2m/s。
作为本实施例中的一种可选实施方式,求解模块通过以下方式求解航线拟合参数组的方程组,得到第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4
设矩阵A为:
Figure BDA0002859408380000117
(’代表转置),
矩阵B为:
Figure BDA0002859408380000121
设C矩阵为求解矩阵,C=inv(B)*A,得出的C矩阵四个值分别为第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值。
通过本发明实施例提供的无人机自适应着陆导航的装置,可以根据不同的机场或空管的适航要求以及各种无人机机型的适航要求自适应地拟合出最优的给定高度拟合曲线,保证了无人机可在不同机场适应不同的下滑窗口要求,并在保证接地升降速度的同时保证着陆的位置精度。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (10)

1.一种无人机的自适应着陆导航的方法,其特征在于,包括:
步骤1、以无人机当前位置的待飞距离为自变量确定下滑阶段的给定高度拟合曲线Hg(x):
Hg(x)=a1+a2x+a3x2+a4x3公式(1),
其中,x表示所述无人机的当前位置距离理想着陆点的水平待飞距离,Hg(x)表示所述无人机的当前位置距离理想着陆点的水平待飞距离为x时的飞行高度,a1为第一航线拟合参数,a2为第二航线拟合参数,a3为第三航线拟合参数,a4为第四航线拟合参数;
步骤2、确定所述下滑阶段的给定高度拟合曲线Hg(x)的一阶微分求导方程:
Figure FDA0002859408370000011
其中,
Figure FDA0002859408370000012
表示所述无人机在任意飞行高度处的曲线斜率;
步骤3、根据机场和空管的要求获取所述无人机进入下滑窗口时允许的飞行高度Ha的值和所述下滑窗口距离理想着陆点的水平待飞距离x0的值;根据所述无人机的机型获取着陆接地时允许的飞行高度Hb的值和着陆接地时距离理想着陆点的水平距离x1的值;
步骤4、将所述Ha的值、所述x0的值、所述Hb的值和所述x1的值依次代入所述公式(1)和所述公式(2)中,建立航线拟合参数组的方程组:
Figure FDA0002859408370000013
步骤5、求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4
步骤6、将求解得到的所述第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值代入所述公式(1)得到符合所述机场和空管以及所述无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。
2.如权利要求1所述的方法,其特征在于:
所述着陆接地时距离理想着陆点的水平距离x1的值为0;
所述无人机着陆接地时允许的飞行高度Hb为所述无人机着陆时重心相对地面的高度。
3.如权利要求1或2所述的方法,其特征在于:
当Hg=Ha,x=x0时,
Figure FDA0002859408370000021
的值为0;
当Hg=Hb,x=x1时,
Figure FDA0002859408370000022
其中,Hyd表示所述无人机着陆时期望的升降速率,Vx表示所述无人机的水平飞行速度。
4.如权利要求3所述的方法,其特征在于:所述无人机的水平飞行速度Vx通过以下公式计算得到:
Figure FDA0002859408370000023
其中,Vg表示所述无人机的给定真空速,所述Hyd的值小于所述无人机的落架能承受的最大升降速率。
5.如权利要求4所述的方法,其特征在于:
所述求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4,包括:
设矩阵A为:
Figure FDA0002859408370000024
’代表转置,
矩阵B为:
Figure FDA0002859408370000025
设C矩阵为求解矩阵,C=inv(B)*A,得出的C矩阵四个值分别为所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4的值。
6.一种无人机的自适应着陆导航的装置,其特征在于,包括:
给定高度拟合曲线构建模块,用于以无人机当前位置的待飞距离为自变量确定下滑阶段的给定高度拟合曲线Hg(x):
Hg(x)=a1+a2x+a3x2+a4x3公式(1),
其中,x表示所述无人机的当前位置距离理想着陆点的水平待飞距离,Hg(x)表示所述无人机的当前位置距离理想着陆点的水平待飞距离为x时的飞行高度,a1为第一航线拟合参数,a2为第二航线拟合参数,a3为第三航线拟合参数,a4为第四航线拟合参数;
所述给定高度拟合曲线构建模块,还用于确定所述下滑阶段的给定高度拟合曲线Hg(x)的一阶微分求导方程:
Figure FDA0002859408370000026
其中,
Figure FDA0002859408370000031
表示所述无人机在任意飞行高度处的曲线斜率;
获取模块,用于根据机场和空管的要求获取所述无人机进入下滑窗口时允许的飞行高度Ha的值和所述下滑窗口距离理想着陆点的水平待飞距离x0的值;根据所述无人机的机型获取着陆接地时允许的飞行高度Hb的值和着陆接地时距离理想着陆点的水平距离x1的值;
求解模块,用于将所述Ha的值、所述x0的值、所述Hb的值和所述x1的值依次代入所述公式(1)和所述公式(2)中,建立航线拟合参数组的方程组:
Figure FDA0002859408370000032
所述求解模块,还用于求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4;将求解得到的所述第一航线拟合参数a1、第二航线拟合参数a2、第三航线拟合参数a3和第四航线拟合参数a4的值代入所述公式(1)得到符合所述机场和空管以及所述无人机的机型的要求的无人机下滑阶段的给定高度拟合曲线。
7.如权利要求6所述的装置,其特征在于:
所述着陆接地时距离理想着陆点的水平距离x1的值为0;
所述无人机着陆接地时允许的飞行高度Hb为所述无人机着陆时重心相对地面的高度。
8.如权利要求6或7所述的装置,其特征在于:
当Hg=Ha,x=x0时,
Figure FDA0002859408370000033
的值为0;
当Hg=Hb,x=x1时,
Figure FDA0002859408370000034
其中,Hyd表示所述无人机着陆时期望的升降速率,Vx表示所述无人机的水平飞行速度。
9.如权利要求8所述的装置,其特征在于:所述无人机的水平飞行速度Vx通过以下公式计算得到:
Figure FDA0002859408370000035
其中,Vg表示所述无人机的给定真空速,所述Hyd的值小于所述无人机的落架能承受的最大升降速率。
10.如权利要求9所述的装置,其特征在于:
所述求解模块通过以下方式求解所述航线拟合参数组的方程组,得到所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4
设矩阵A为:
Figure FDA0002859408370000041
’代表转置,
矩阵B为:
Figure FDA0002859408370000042
设C矩阵为求解矩阵,C=inv(B)*A,得出的C矩阵四个值分别为所述第一航线拟合参数a1、所述第二航线拟合参数a2、所述第三航线拟合参数a3和所述第四航线拟合参数a4的值。
CN202011561067.4A 2020-12-25 2020-12-25 一种无人机自适应着陆导航的方法及装置 Active CN112666978B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011561067.4A CN112666978B (zh) 2020-12-25 2020-12-25 一种无人机自适应着陆导航的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011561067.4A CN112666978B (zh) 2020-12-25 2020-12-25 一种无人机自适应着陆导航的方法及装置

Publications (2)

Publication Number Publication Date
CN112666978A CN112666978A (zh) 2021-04-16
CN112666978B true CN112666978B (zh) 2021-11-26

Family

ID=75408998

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011561067.4A Active CN112666978B (zh) 2020-12-25 2020-12-25 一种无人机自适应着陆导航的方法及装置

Country Status (1)

Country Link
CN (1) CN112666978B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987772A (zh) * 2021-05-11 2021-06-18 北京三快在线科技有限公司 无人机近地面飞行控制方法、装置、介质、设备和无人机
CN113320716B (zh) * 2021-06-16 2022-08-09 上海交通大学 具备自主调节足端触地点的外星逆境着陆器

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101763116A (zh) * 2008-12-24 2010-06-30 中国科学院自动化研究所 一种基于侧向导引的抗侧风着陆航迹跟踪控制方法
CN102390543A (zh) * 2011-08-23 2012-03-28 北京航空航天大学 一种无人机的纵向着陆轨迹的设计方法
CN103116359A (zh) * 2013-01-28 2013-05-22 北京航空航天大学 一种无人机着陆下滑初始段的引导方法
CN103197682A (zh) * 2013-03-22 2013-07-10 北京航空航天大学 一种无人机进场着陆下滑通道设计与下滑轨迹调整方法
CN105159308A (zh) * 2015-08-28 2015-12-16 北京控制工程研究所 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101763116A (zh) * 2008-12-24 2010-06-30 中国科学院自动化研究所 一种基于侧向导引的抗侧风着陆航迹跟踪控制方法
CN102390543A (zh) * 2011-08-23 2012-03-28 北京航空航天大学 一种无人机的纵向着陆轨迹的设计方法
CN103116359A (zh) * 2013-01-28 2013-05-22 北京航空航天大学 一种无人机着陆下滑初始段的引导方法
CN103197682A (zh) * 2013-03-22 2013-07-10 北京航空航天大学 一种无人机进场着陆下滑通道设计与下滑轨迹调整方法
CN105159308A (zh) * 2015-08-28 2015-12-16 北京控制工程研究所 一种可重复使用飞行器着陆段制导与控制律一体化耦合设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于日盲区紫外成像的无人机着陆引导技术研究;张磊 等;《中国激光》;20160731;第43卷(第7期);第0704001-1至0704001-10页 *

Also Published As

Publication number Publication date
CN112666978A (zh) 2021-04-16

Similar Documents

Publication Publication Date Title
CN112666978B (zh) 一种无人机自适应着陆导航的方法及装置
CN110908396B (zh) 可重复使用运载器全阶段再入返回制导方法
CN110727287B (zh) 用于确定爬升剖面的方法和系统以及预测爬升剖面的方法
US8849480B2 (en) Aircraft gross weight and center of gravity validator
CN106525375A (zh) 一种用于无人机抗风能力检测的自适应系统
US8660722B2 (en) Method for optimizing aircraft landing on a runway
CN109814593B (zh) 一种可自主寻热的低空太阳能无人机飞行控制方法和系统
US20100222944A1 (en) Automatic pilot pitch angle compensation
CN114065398B (zh) 一种大展弦比柔性飞行器飞行性能计算方法
CN111522356B (zh) 一种倾转旋翼无人机强鲁棒全包线一体化控制方法
CN106940568A (zh) 飞行控制系统及方法
CN113985914A (zh) 一种无人设备降落控制方法、装置、设备及存储介质
CN116954239B (zh) 一种基于改进粒子群算法的无人机航迹规划方法及系统
CN116339387B (zh) 一种狭小空间复杂紊流影响下无人机安全距离保持方法
CN115981376B (zh) 一种基于气流模型的无人机长滞空飞行规划方法和装置
CN107688689B (zh) 一种基于分层加权的飞行程序噪声评估方法
CN113589838B (zh) 一种基于圆柱位置离散化的三维轨迹调度方法
CN115510910A (zh) 一种优化飞行模拟器动感模拟算法中的滤波器参数的方法
CN111580533B (zh) 基于空气动力学的无人机信息收集方法及装置
CN109240338B (zh) 基于超声波传感器阵列的无人机自主降落方法
CN113961011A (zh) 一种航迹倾角保持控制方法
CN109211225A (zh) 获取大椭圆轨道空间物体剩余轨道寿命的方法、系统及设备
CN117894227A (zh) 一种无人机模拟卫星轨迹规划方法
CN116661471A (zh) 飞行器的拉平模式控制方法、系统、存储介质及终端设备
CN113110575A (zh) 无人机导地线巡视航迹规划方法、系统及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant