CN113655807B - 一种计算空中匀变速运动目标飞行轨迹的方法 - Google Patents
一种计算空中匀变速运动目标飞行轨迹的方法 Download PDFInfo
- Publication number
- CN113655807B CN113655807B CN202110805000.9A CN202110805000A CN113655807B CN 113655807 B CN113655807 B CN 113655807B CN 202110805000 A CN202110805000 A CN 202110805000A CN 113655807 B CN113655807 B CN 113655807B
- Authority
- CN
- China
- Prior art keywords
- target
- point
- flight
- data
- calculating
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000008859 change Effects 0.000 title abstract description 3
- 230000008569 process Effects 0.000 claims abstract description 15
- 229930186657 Lat Natural products 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 230000003287 optical effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Navigation (AREA)
Abstract
本发明的一个实施例公开了一种计算空中匀变速运动目标飞行轨迹的方法,包括:获取目标的航段数据,根据航段数据,计算目标飞行过程中任意时刻的飞行数据。本申请将目标路径规划分解为匀变速运动的子线段,各子线段连接成连续的航迹,根据每段航迹的航迹数据,计算目标飞行过程中任意时刻的飞行数据,即可实现目标路径规划。无需多次坐标转换,计算过程简单,有广泛的应用前景。
Description
技术领域
本发明涉及飞行航迹领域。更具体地,涉及一种计算空中匀变速运动目标飞行轨迹的方法。
背景技术
在空情模拟的应用场景中,受地球曲率影响,实现空中目标的匀变速飞行往往涉及多次坐标转换,例如从地理坐标系转换到地心直角坐标系或其他直角坐标系,计算过程比较复杂。
发明内容
有鉴于此,本发明的第一个实施例提供一种计算空中匀变速运动目标飞行轨迹的方法,其特征在于,包括:
获取目标的航段数据,
根据航段数据,计算目标飞行过程中任意时刻的飞行数据。
在一个具体实施例中,将飞行轨迹划分为N个航段,每个航段的航段数据包括起始点数据和终止点数据,
在一个具体实施例中,所述起始点数据包括起始点经度、起始点纬度、起始点高度和起始点速度;
所述终止点数据包括:终止点经度、终止点纬度、终止点高度和终止点速度。
在一个具体实施例中,所述计算目标飞行过程中任意时刻的飞行数据包括:
计算目标平均速度:
Spd=(Spds+Spde)/2
其中,Spds为起始点速度,Spde为终止点速度;
计算起始点、终止点与地心夹角:
其中,Lats为起始点纬度,Late为终止点纬度,Lngs为起始点经度,Lnge为终止点经度;
根据起始点、终止点与地心夹角计算此航段目标飞行长度:
其中,Re为地球半径;
根据所述此航段目标飞行长度计算飞行时间:
T=S/Spd;
根据所述飞行时间计算目标飞行过程中任意时刻的飞行数据。
在一个具体实施例中,所述任意时刻的飞行数据包括:任一时刻目标所在位置的经度、纬度、高度和速度。
在一个具体实施例中,所述任意时刻t目标所在位置的高度为:
其中,Alts为起始点高度,Alte为终止点高度。
在一个具体实施例中,,所述任意时刻t目标所在位置的纬度为:
如果Lat>90,Lat=180-Lat,目标通过北极;
如果Lat<-90,Lat=-180-Lat,目标通过南极。
在一个具体实施例中,所述任意时刻t目标所在位置的经度Lng为:
如果目标在南北极点,此时设Lng=0;
如果绝对值大于1,此时设Lng=0;
如果Lnge>Lngs,
如果Lnge<Lngs,
其中,任意时刻t目标所在位置点、初始点与地心夹角θt为:
本发明的第二个实施例提供一种计算机设备,包括处理器及存储在存储有计算机程序的存储器,其特征在于,所述处理器执行所述程序时实现如第一个实施例中任一项所述的方法。
本发明的第三个实施例提供一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如第一个实施例中任一项所述的方法。
本发明的有益效果如下:
本申请将目标路径规划分解为匀变速运动的子线段,各子线段连接成连续的航迹,根据每段航迹的航迹数据,计算目标飞行过程中任意时刻的飞行数据,即可实现目标路径规划。无需转换到地心直角坐标系或其他直角坐标系,计算过程简单,有广泛的应用前景。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出根据本发明一个实施例的计算空中匀变速运动目标飞行轨迹的方法流程图
图2示出根据本发明一个实施例的航段组成示意图。
图3示出本发明的另一个实施例的计算机设备的结构示意图。
具体实施方式
为使本发明的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
如图1所示,一种计算空中匀变速运动目标飞行轨迹的方法,包括:
获取目标的航段数据,
根据航段数据,计算目标飞行过程中任意时刻的飞行数据。
将飞行轨迹划分为N个航段,前一航段的终止点经度、纬度、高度、速度和后一航段的起始点经度、纬度、高度、速度相同,以保证航迹的位置和速度的连续性。
在每一航段中,目标从起始点到终止点为匀变速运动,并且目标运动轨迹在起始点、终止点和地心三点所构成的平面内。构建由起始点、终止点、地心以及飞行轨迹的位置关系如图2所示,目标沿图中Ps-Pe曲线飞行。
每个航段的航段数据包括起始点数据和终止点数据,
所述起始点数据包括起始点经度、起始点纬度、起始点高度和起始点速度;
所述终止点数据包括:终止点经度、终止点纬度、终止点高度和终止点速度。
在一个具体实施例中,所述计算目标飞行过程中任意时刻的飞行数据包括:
计算目标平均速度:
Spd=(Spds+Spde)/2
其中,Spds为起始点速度,Spde为终止点速度;
计算起始点、终止点与地心夹角:
其中,Lats为起始点纬度,Late为终止点纬度,Lngs为起始点经度,Lnge为终止点经度;
根据起始点、终止点与地心夹角计算此航段目标飞行长度:
其中,Re为地球半径;
根据所述此航段目标飞行长度计算飞行时间:
T=S/Spd;
根据所述飞行时间计算目标飞行过程中任意时刻的飞行数据。
在一个具体实施例中,所述任意时刻的飞行数据包括:任一时刻目标所在位置的经度、纬度、高度和速度。
所述任意时刻t目标所在位置的高度为:
其中,Alts为起始点高度,Alte为终止点高度。
所述任意时刻t目标所在位置的纬度为:
如果Lat>90,Lat=180-Lat,目标通过北极;
如果Lat<-90,Lat=-180-Lat,目标通过南极。
所述任意时刻t目标所在位置的经度Lng为:
如果目标在南北极点,此时设Lng=0;
如果绝对值大于1,此时设Lng=0;
如果Lnge>Lngs,
如果Lnge<Lngs,
其中,任意时刻t目标所在位置点、初始点与地心夹角θt为:
本发明的另一个实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现,在实际应用中,所述计算机可读存储介质可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于电、磁、光、电磁、红外线、或半导体的***、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本实施例中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行***、装置或者器件使用或者与其结合使用。
如图3所示,本发明的另一个实施例提供的一种计算机设备的结构示意图。图3显示的计算机设备12仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。
如图3所示,计算机设备12以通用计算设备的形式表现。计算机设备12的组件可以包括但不限于:一个或者多个处理器或者处理单元16,***存储器28,连接不同***组件(包括***存储器28和处理单元16)的总线18。
总线18表示几类总线结构中的一种或多种,包括存储器总线或者存储器控制器,***总线,图形加速端口,处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说,这些体系结构包括但不限于工业标准体系结构(ISA)总线,微通道体系结构(MAC)总线,增强型ISA总线、视频电子标准协会(VESA)局域总线以及***组件互连(PCI)总线。
计算机设备12典型地包括多种计算机***可读介质。这些介质可以是任何能够被计算机设备12访问的可用介质,包括易失性和非易失性介质,可移动的和不可移动的介质。
***存储器28可以包括易失性存储器形式的计算机***可读介质,例如随机存取存储器(RAM)30和/或高速缓存存储器32。计算机设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机***存储介质。仅作为举例,存储***34可以用于读写不可移动的、非易失性磁介质(图3未显示,通常称为“硬盘驱动器”)。尽管图3中未示出,可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器,以及对可移动非易失性光盘(例如CD-ROM,DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下,每个驱动器可以通过一个或者多个数据介质接口与总线18相连。存储器28可以包括至少一个程序产品,该程序产品具有一组(例如至少一个)程序模块,这些程序模块被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块42的程序/实用工具40,可以存储在例如存储器28中,这样的程序模块42包括但不限于操作***、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块42通常执行本发明所描述的实施例中的功能和/或方法。
计算机设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24等)通信,还可与一个或者多个使得用户能与该计算机设备12交互的设备通信,和/或与使得该计算机设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡,调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口22进行。并且,计算机设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。如图3所示,网络适配器20通过总线18与计算机设备12的其它模块通信。应当明白,尽管图3中未示出,可以结合计算机设备12使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID***、磁带驱动器以及数据备份存储***等。
处理器单元16通过运行存储在***存储器28中的程序,从而执行各种功能应用以及数据处理,例如实现本发明实施例所提供的一种计算空中匀变速运动目标飞行轨迹的方法。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (8)
1.一种计算空中匀变速运动目标飞行轨迹的方法,其特征在于,包括:
获取目标的航段数据,
根据航段数据,计算目标飞行过程中任意时刻的飞行数据;
将飞行轨迹划分为N个航段,每个航段的航段数据包括起始点数据和终止点数据;
所述计算目标飞行过程中任意时刻的飞行数据包括:
计算目标平均速度:
Spd=(Spds+Spde)/2
其中,Spds为起始点速度,Spde为终止点速度;
计算起始点、终止点与地心夹角:
其中,Lats为起始点纬度,Late为终止点纬度,Lngs为起始点经度,Lnge为终止点经度;
根据起始点、终止点与地心夹角计算此航段目标飞行长度:
其中,Re为地球半径,Alts为起始点高度,Alte为终止点高度;
根据所述此航段目标飞行长度计算飞行时间:
T=S/Spd;
根据所述飞行时间计算目标飞行过程中任意时刻的飞行数据。
2.根据权利要求1所述的方法,其特征在于,所述起始点数据包括起始点经度、起始点纬度、起始点高度和起始点速度;
所述终止点数据包括:终止点经度、终止点纬度、终止点高度和终止点速度。
3.根据权利要求1所述的方法,其特征在于,所述任意时刻的飞行数据包括:任一时刻目标所在位置的经度、纬度、高度和速度。
4.根据权利要求1所述的方法,其特征在于,所述任意时刻t目标所在位置的高度为:
其中,Alts为起始点高度,Alte为终止点高度。
5.根据权利要求4所述的方法,其特征在于,所述任意时刻t目标所在位置的纬度为:
如果Lat>90,Lat=180-Lat,目标通过北极;
如果Lat<-90,Lat=-180-Lat,目标通过南极。
6.根据权利要求5所述的方法,其特征在于,所述任意时刻t目标所在位置的经度Lng为:
如果目标在南北极点,此时设Lng=0;
如果绝对值大于1,此时设Lng=0;
如果Lnge>Lngs,
如果Lnge<Lngs,
其中,任意时刻t目标所在位置点、初始点与地心夹角θt为:
7.一种计算机设备,包括处理器及存储在存储有计算机程序的存储器,其特征在于,所述处理器执行所述程序时实现如权利要求1-6中任一项所述的方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-6中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110805000.9A CN113655807B (zh) | 2021-07-16 | 2021-07-16 | 一种计算空中匀变速运动目标飞行轨迹的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110805000.9A CN113655807B (zh) | 2021-07-16 | 2021-07-16 | 一种计算空中匀变速运动目标飞行轨迹的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113655807A CN113655807A (zh) | 2021-11-16 |
CN113655807B true CN113655807B (zh) | 2024-02-02 |
Family
ID=78489436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110805000.9A Active CN113655807B (zh) | 2021-07-16 | 2021-07-16 | 一种计算空中匀变速运动目标飞行轨迹的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113655807B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000023967A1 (en) * | 1998-10-16 | 2000-04-27 | Universal Avionics Systems Corporation | Flight plan intent alert system and method |
CN104794213A (zh) * | 2015-04-27 | 2015-07-22 | 清华大学 | 基于飞行计划与导航数据库的航迹预测方法 |
CN111351488A (zh) * | 2020-03-03 | 2020-06-30 | 南京航空航天大学 | 飞行器智能轨迹重构再入制导方法 |
CN112231389A (zh) * | 2020-10-16 | 2021-01-15 | 中国民用航空华东地区空中交通管理局 | 基于航迹的可视化冲突模型构建方法、装置、电子设备及存储介质 |
WO2021071279A2 (ko) * | 2019-10-09 | 2021-04-15 | 주식회사 에스360브이알 | 지리상의 위치를 특정하는 방법과 이를 이용하는 데이터베이스 및 데이터베이스의 데이터베이스 |
CN112697146A (zh) * | 2020-11-19 | 2021-04-23 | 北京电子工程总体研究所 | 一种基于稳健回归的航迹预测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3025920B1 (fr) * | 2014-09-15 | 2016-11-04 | Thales Sa | Procede de calcul temps reel d'une trajectoire planifiee, notamment de plan de vol, combinant une mission, et systeme de gestion d'une telle trajectoire |
US11270596B2 (en) * | 2018-10-11 | 2022-03-08 | Xwing, Inc. | Autonomous path planning |
-
2021
- 2021-07-16 CN CN202110805000.9A patent/CN113655807B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000023967A1 (en) * | 1998-10-16 | 2000-04-27 | Universal Avionics Systems Corporation | Flight plan intent alert system and method |
CN104794213A (zh) * | 2015-04-27 | 2015-07-22 | 清华大学 | 基于飞行计划与导航数据库的航迹预测方法 |
WO2021071279A2 (ko) * | 2019-10-09 | 2021-04-15 | 주식회사 에스360브이알 | 지리상의 위치를 특정하는 방법과 이를 이용하는 데이터베이스 및 데이터베이스의 데이터베이스 |
CN111351488A (zh) * | 2020-03-03 | 2020-06-30 | 南京航空航天大学 | 飞行器智能轨迹重构再入制导方法 |
CN112231389A (zh) * | 2020-10-16 | 2021-01-15 | 中国民用航空华东地区空中交通管理局 | 基于航迹的可视化冲突模型构建方法、装置、电子设备及存储介质 |
CN112697146A (zh) * | 2020-11-19 | 2021-04-23 | 北京电子工程总体研究所 | 一种基于稳健回归的航迹预测方法 |
Non-Patent Citations (2)
Title |
---|
ADS-B历史飞行轨迹数据清洗方法;王兵;;交通运输工程学报(第04期);223-232 * |
考虑地球扁率的RNAV航路规划;方学东;;测绘科学(第06期);130+150-151 * |
Also Published As
Publication number | Publication date |
---|---|
CN113655807A (zh) | 2021-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220020175A1 (en) | Method for Training Object Detection Model, Object Detection Method and Related Apparatus | |
CN113345228B (zh) | 基于拟合轨迹的行驶数据的生成方法、装置、设备及介质 | |
CN110677815A (zh) | 停留点识别方法、装置、计算机设备及存储介质 | |
CN109218562B (zh) | 时钟同步方法、装置、设备、存储介质和车辆 | |
CN110647675B (zh) | 停留点识别及预测模型训练方法、装置及存储介质 | |
CN110244765B (zh) | 一种飞行器航线轨迹生成方法、装置、无人机及存储介质 | |
CN110017840B (zh) | 导航数据的处理方法和装置、电子设备、计算机可读介质 | |
CN113655807B (zh) | 一种计算空中匀变速运动目标飞行轨迹的方法 | |
CN107729944B (zh) | 一种低俗图片的识别方法、装置、服务器及存储介质 | |
CN110737690A (zh) | 用户标签挖掘方法、装置、计算机设备和存储介质 | |
CN112419739A (zh) | 车辆定位方法、装置和电子设备 | |
CN109348408B (zh) | 一种车位确定方法、装置、电子设备、车辆及存储介质 | |
CN108268821A (zh) | 一种室内外场景识别方法及装置 | |
CN114356475B (zh) | 一种显示处理方法、装置、设备及存储介质 | |
CN112507532B (zh) | 一种平面空间二对一追踪的区域分类及制导控制方法 | |
CN111862052B (zh) | 检测缝隙的方法、装置、设备及介质 | |
CN112445242B (zh) | 航线的跟踪方法、装置、设备及存储介质 | |
CN112507529B (zh) | 一种智能体燃料节省制导追踪方法 | |
CN113759331B (zh) | 一种雷达定位方法、装置、设备及存储介质 | |
CN110992494B (zh) | 一种角色模型的消散方法、装置、设备及存储介质 | |
CN102375325B (zh) | 一种真三维仿真角度描述及直接投影显示方法 | |
CN113640845B (zh) | 位置数据延迟改正优化的方法、装置、设备及存储介质 | |
CN109919121B (zh) | 一种人体模型的投影方法、装置、电子设备及存储介质 | |
CN111091848B (zh) | 一种用于预测头部姿态的方法和装置 | |
CN116051925B (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 |