CN114444267A - 一种飞行计划航段水平航迹预测方法 - Google Patents

一种飞行计划航段水平航迹预测方法 Download PDF

Info

Publication number
CN114444267A
CN114444267A CN202111633894.4A CN202111633894A CN114444267A CN 114444267 A CN114444267 A CN 114444267A CN 202111633894 A CN202111633894 A CN 202111633894A CN 114444267 A CN114444267 A CN 114444267A
Authority
CN
China
Prior art keywords
turning
calculating
point
latitude
longitude
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.)
Pending
Application number
CN202111633894.4A
Other languages
English (en)
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.)
Xian Aeronautics Computing Technique Research Institute of AVIC
Original Assignee
Xian Aeronautics Computing Technique Research Institute of AVIC
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 Xian Aeronautics Computing Technique Research Institute of AVIC filed Critical Xian Aeronautics Computing Technique Research Institute of AVIC
Priority to CN202111633894.4A priority Critical patent/CN114444267A/zh
Publication of CN114444267A publication Critical patent/CN114444267A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Navigation (AREA)

Abstract

本发明提供一种飞行计划航段水平航迹预测方法,属于飞行管理系统设计领域,该方法包括:(1)椭球面坐标系定义,将位置点经纬度转化为椭球面坐标;(2)计算转弯半径R,计算应飞航向角(北纬、东经输入为正;南纬、西经输入为负进行计算),转弯起始点S距离中间航路点距离D计算;(3)将航路点、航路转弯点转换到球面坐标系中,通过共面建立方程,通过与距离公式联立解算转弯点;(4)将所有的转弯点平滑连接,完成水平航迹预测与生成并输出。本发明将转弯点与航路转换到球面坐标系中,通过三点共面与距离信息联立方程,从而计算出转弯点经纬度,得到预测航路,该方法思路清晰,抓住问题核心,效率高。

Description

一种飞行计划航段水平航迹预测方法
技术领域
本公开涉及机载飞行管理技术领域,尤其涉及一种飞行计划航段水平航迹预测方法。
背景技术
通用飞机越来越普及,广泛用于农、林、牧、渔产业﹑旅游业等。直飞功能是一种飞行计划的执行方式,飞机在飞行过程中,由于危险天气、空管限制等因素,需要偏离当前激活航线,直接飞至某一目标航路点或者目标截断航段上。
现有飞行管理系统的航迹预测方法为:1)先通过计算转弯半径R和当前飞行航向;2)再基于转弯半径R和当前飞行航向,确定圆心O坐标;3)以O点为圆心、半径为转弯半径R,沿着航向方向作圆弧L;4)根据截断航线角φc确定所述圆弧L的切线斜率k,以切线斜率k作圆弧L的切线,切点即为转弯结束点P,切线与目标截断航段的交点为理论截断点Q。这种方法是通过计算转弯圆心点O经纬度,再计算转弯点经纬度,因此必须计算转弯中心点坐标,导致过程复杂,效率低。
发明内容
有鉴于此,本公开实施例提供一种飞行计划航段水平航迹预测方法,区别于传统的转弯点经纬度预测方法,不需要计算转弯中心点经纬度,而是将转弯点与航路转换到球面坐标系中,通过三点共面与距离信息联立方程,从而计算出转弯点经纬度,方法直接有效。当所有的航路转弯点计算出来后,平滑连接所有的水平航路对应转弯点,形成一条预测航路。
为了实现上述目的,本发明提供如下技术方案:
一种飞行计划航段水平航迹预测方法,包括以下步骤:
步骤1:对地球建立椭球面坐标系:以地心为原点,建立O-XYZ坐标系;
步骤2:在飞行航段上确定三个航路点P1、P2、P3,获取每个航路点纬度和经度坐标P1(N1,E1),P2(N2,E2),P3(N3,E3),确定起始转弯点S、终止转弯点E和转弯中心点O,计算转弯半径R,再计算转弯起始点S与P2航路点之间的距离D;
步骤3:将航路点P1、P2和起始转弯点S转换到球面O-XYZ坐标系中,得到航路点P1、P2和起始转弯点S的球坐标,通过三点共面建立共面方程,结合距离计算公式求解得到起始转弯点S的经纬度;
步骤4:将起始转弯点S的经纬度带入共面方程,计算得到终止转弯点E的经纬度,将所有的转弯点平滑连接,得到水平航迹预测并输出。
进一步地,步骤2中,转弯半径R的计算方法为:
转弯半径R=(Vgs)2/g×tanφ,
其中,φ为当前横滚角;Vgs为预测地速,根据当前飞机飞行高度在飞机性能数据库中查取。
进一步地,步骤2中,转弯起始点S与P2航路点之间的距离D计算方法如下:
(1)将航路点P1,P2,P3经纬度转化为球面O-XYZ坐标系表示;
(2)计算航路点P1到P2、P2到P3应飞磁航向角;
(3)计算航路点P1_P2_P3航路夹角;
(4)结合转弯半径R和航路夹角计算出距离D。
进一步地,步骤3中,求解起始转弯点S的经纬度时取航路点P1和P2之间的值。
本发明提出的一种飞行计划航段水平航迹预测方法,在通过确定转弯半径后,假设将转弯起始点、转弯终止点已知处理。通过建立椭球面坐标模型,并基于该模型进行航段和转弯起始点,终止点模型建立,通过转弯点与航路点共面和距离计算信息,反向计算出转弯起始点和转弯终止点的坐标。避免了传统方法通过计算转弯圆心点O,再计算转弯点情形,更加直接有效率。
其中,在计算转弯点经纬度的过程中,通过最新的基于模型的系统工具进行建模,验证由仿真结果可得,该方法能够有效的计算转弯点经纬度,通过已知的航路点经纬度信息,转弯半径最后将三个点切换到球面坐标中,通过球面坐标的视角通过共面成功计算出转弯点经纬度。
附图说明
为了更清楚地说明本公开实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本公开的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例中椭球面坐标系图;
图2为本发明实施例中地球表面任意一点P投影到椭球面坐标图;
图3为本发明实施例中水平航路点航段示意图;
图4为本发明实施例中航路点预测转弯球面坐标系示意图;
图5为本发明实施例中平航迹预测平滑结果示意图;
图6为本发明实施例中水平航迹预测方法流程图。
具体实施方式
下面结合附图对本公开实施例进行详细描述。
以下通过特定的具体实例说明本公开的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本公开的其他优点与功效。显然,所描述的实施例仅仅是本公开一部分实施例,而不是全部的实施例。本公开还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本公开的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。基于本公开中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本公开保护的范围。
要说明的是,下文描述在所附权利要求书的范围内的实施例的各种方面。应显而易见,本文中所描述的方面可体现于广泛多种形式中,且本文中所描述的任何特定结构及/或功能仅为说明性的。基于本公开,所属领域的技术人员应了解,本文中所描述的一个方面可与任何其它方面独立地实施,且可以各种方式组合这些方面中的两者或两者以上。举例来说,可使用本文中所阐述的任何数目个方面来实施设备及/或实践方法。另外,可使用除了本文中所阐述的方面中的一或多者之外的其它结构及/或功能性实施此设备及/或实践此方法。
还需要说明的是,以下实施例中所提供的图示仅以示意方式说明本公开的基本构想,图式中仅显示与本公开中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
另外,在以下描述中,提供具体细节是为了便于透彻理解实例。然而,所属领域的技术人员将理解,可在没有这些特定细节的情况下实践所述方面。
本公开实施例的一种飞行计划航段水平航迹预测方法,如图6所示,具体实施方法如下:
(1)先对地球建立椭球面坐标系:以地心为原点,建立O-XYZ坐标系。赤道为基本面,垂直于赤道基本面经过地心的一条轴称为Z轴,指向北极为正。本初子午面与赤道面重合的轴称为X轴,向东为正。Y轴与XZ平面垂直构成右手定则,如图1所示。图2为地球表面任意一点P投影到椭球面坐标图。
(2)确定水平航路点航段P1(N1,E1),P2(N2,E2),P3(N3,E3)。N代表纬度,E代表经度。P1、P2、P3共同构成部分航路。如图3所示。起始转弯点S(Ns,Es)、终止转弯点E(Ne,Ee),转弯中心点O,转弯半径为R,起始转弯点S到P2航路点距离为D,如图3所示。北纬、东经为正;南纬西经为负;以示区分。
(3)将航路点P1、S、P2三个点转换到球面坐标系统中,如图4所示,再以球面坐标系视角进行分析。设三个航路点坐标P1(N1,E1),P2(N2,E2),P3(N3,E3),转弯中心点O,飞机起始转弯点为S,终止转弯点E。则P1,P2,S三点的球坐标为:
起始点P1:θ1=90-N1,
Figure BDA0003441857760000061
Figure BDA0003441857760000062
(西经)
航路点P2:θ2=90-N2,
Figure BDA0003441857760000063
Figure BDA0003441857760000064
(西经)
转弯点S:θ=90-Ns,
Figure BDA0003441857760000065
Figure BDA0003441857760000066
(西经)
分别将P1、P2、S点投影到OX、OY、OZ轴上;则:
Figure BDA0003441857760000067
Figure BDA0003441857760000068
Figure BDA0003441857760000069
其中r,r1,r2为经纬度S,P1,P2转换到地心坐标系的矢量形式。
由于三矢量共面条件可得:
r·r1×r2=0 (2)
(4)滚转角计算:依据DO-283A标准,设α为航迹改变角,单位为度。φ为滚转角,则滚转角和航迹改变角关系:
φ=min(0.5*α,23) (3)
(5)转弯半径R计算:
转弯半径R:计算转弯半径R,转弯半径R=(Vgs)2/g×tanφ,其中,φ为当前横滚角;Vgs为预测地速,根据当前飞机飞行高度在现有飞机性能数据库中查取;
(6)距离D计算:
1)将航路点P1,P2,P3经纬度转化为球面O-XYZ坐标系表示;
2)计算P1到P2、P2到P3应飞磁航向角;
3)再计算P1_P2_P3航路夹角;
4)结合转弯半径和夹角信息计算出距离D。
(7)又因:
D=Rω (4)
其中,ω为P2_O_S弧度夹角,进一步公式(2)可得:
cosω=cos(D/R) (5)
又因:cosω=cos N2 cos Ns cos(Es-E2)+sin Ns sin N2 (6)
由公式(2)和(6)联立可解出S点的经纬度。在这里需要注意的是,求解的S的经度应取P1、P2之间的值,在确定了Es后,再带入公式(2)求解S点的纬度信息Ns,tan Ns>0,Ns为北纬;tan Ns<0,Ns为南纬。计算转弯终止点E点同理。
其余航路点情况类似,等计算出来所有航路点对应的转弯起始点、转弯终止点后,将这些点平滑连接如图5所示,共同构成一条水平航迹预测航路。
本发明不同于传统方法必须计算转弯中心点坐标,本发明通过将问题化解到球面坐标系中,通过共面特点,矢量叉乘和距离信息联立,反向解算出航迹转弯点坐标,方法直接有效。当所有的航路转弯点计算出来后,平滑连接所有的水平航路对应转弯点,形成一条预测航路。
以上所述,仅为本公开的具体实施方式,但本公开的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本公开揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本公开的保护范围之内。因此,本公开的保护范围应以权利要求的保护范围为准。

Claims (4)

1.一种飞行计划航段水平航迹预测方法,其特征在于,包括以下步骤:
步骤1:对地球建立椭球面坐标系:以地心为原点,建立O-XYZ坐标系;
步骤2:在飞行航段上确定三个航路点P1、P2、P3,获取每个航路点纬度和经度坐标P1(N1,E1),P2(N2,E2),P3(N3,E3),确定起始转弯点S、终止转弯点E和转弯中心点O,计算转弯半径R,再计算转弯起始点S与P2航路点之间的距离D;
步骤3:将航路点P1、P2和起始转弯点S转换到球面O-XYZ坐标系中,得到航路点P1、P2和起始转弯点S的球坐标,通过三点共面建立共面方程,结合距离计算公式求解得到起始转弯点S的经纬度;
步骤4:将起始转弯点S的经纬度带入共面方程,计算得到终止转弯点E的经纬度,将所有的转弯点平滑连接,得到水平航迹预测并输出。
2.根据权利要求1所述的飞行计划航段水平航迹预测方法,其特征在于,步骤2中,转弯半径R的计算方法为:
转弯半径R=(Vgs)2/g×tanφ,
其中,φ为当前横滚角;Vgs为预测地速,根据当前飞机飞行高度在飞机性能数据库中查取。
3.根据权利要求2所述的飞行计划航段水平航迹预测方法,其特征在于,步骤2中,转弯起始点S与P2航路点之间的距离D计算方法如下:
(1)将航路点P1,P2,P3经纬度转化为球面O-XYZ坐标系表示;
(2)计算航路点P1到P2、P2到P3应飞磁航向角;
(3)计算航路点P1_P2_P3航路夹角;
(4)结合转弯半径R和航路夹角计算出距离D。
4.根据权利要求1所述的飞行计划航段水平航迹预测方法,其特征在于,步骤3中,求解起始转弯点S的经纬度时取航路点P1和P2之间的值。
CN202111633894.4A 2021-12-29 2021-12-29 一种飞行计划航段水平航迹预测方法 Pending CN114444267A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111633894.4A CN114444267A (zh) 2021-12-29 2021-12-29 一种飞行计划航段水平航迹预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111633894.4A CN114444267A (zh) 2021-12-29 2021-12-29 一种飞行计划航段水平航迹预测方法

Publications (1)

Publication Number Publication Date
CN114444267A true CN114444267A (zh) 2022-05-06

Family

ID=81365681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111633894.4A Pending CN114444267A (zh) 2021-12-29 2021-12-29 一种飞行计划航段水平航迹预测方法

Country Status (1)

Country Link
CN (1) CN114444267A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114895705A (zh) * 2022-05-07 2022-08-12 中国民航科学技术研究院 一种飞机在地面滑行阶段的转弯半径计算方法及装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114895705A (zh) * 2022-05-07 2022-08-12 中国民航科学技术研究院 一种飞机在地面滑行阶段的转弯半径计算方法及装置

Similar Documents

Publication Publication Date Title
US11479359B2 (en) UAV operation route planning method, UAV pesticide spreading method and device
CN109828607B (zh) 一种面向不规则障碍物的无人机路径规划方法及系统
CN103557867B (zh) 一种基于稀疏a*搜索的三维多uav协同航迹规划方法
US6529821B2 (en) Route planner with area avoidance capability
CN104833360B (zh) 一种二维坐标到三维坐标的转换方法
Wang et al. Real-time dynamic Dubins-Helix method for 3-D trajectory smoothing
CN108204814B (zh) 无人机三维场景路径导航平台及其三维改进路径规划方法
CN112097770B (zh) 一种多无人机协同全覆盖的路径规划方法、装置、存储介质及终端
CN107544067A (zh) 一种基于高斯混合近似的高超声速再入飞行器跟踪方法
CN106546245B (zh) 基于ads-b数据的飞机轨迹推断和平滑方法
CN111561933B (zh) 双重改进a星最短航路规划方法
CN109615936A (zh) 机载飞行管理系统中的直飞航迹预测方法和直飞方法
CN108225333B (zh) 一种用于航线规划的最优路径生成方法
US10424208B2 (en) Navigational aid method depending on meteorological conditions
CN111897364B (zh) 一种无人机区域轨迹规划方法、装置及可读存储介质
CN105844969B (zh) 根据气象条件来改进飞行器的飞行轨迹的方法
US20150355882A1 (en) Coordinate transformation device, non-transitory computer readable medium storing coordinate transformation program, and coordinate transformation method
CN112148033B (zh) 无人机航线的确定方法、装置、设备及存储介质
CN114444267A (zh) 一种飞行计划航段水平航迹预测方法
CN115562290A (zh) 一种基于a星惩罚控制优化算法的机器人路径规划方法
CN113361174B (zh) 基于stp模型的大型无人机碰撞概率计算方法
Si et al. TOM-Odometry: A generalized localization framework based on topological map and odometry
Stephan et al. Fast generation of landing paths for fixed-wing aircraft with thrust failure
CN112037330A (zh) 基于AirSim的无人机运行场景仿真方法
CN111121779A (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