CN112199816B - 航班ads-b数据模拟方法 - Google Patents

航班ads-b数据模拟方法 Download PDF

Info

Publication number
CN112199816B
CN112199816B CN202010893875.4A CN202010893875A CN112199816B CN 112199816 B CN112199816 B CN 112199816B CN 202010893875 A CN202010893875 A CN 202010893875A CN 112199816 B CN112199816 B CN 112199816B
Authority
CN
China
Prior art keywords
flight
data
ads
aircraft
subfield
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
CN202010893875.4A
Other languages
English (en)
Other versions
CN112199816A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010893875.4A priority Critical patent/CN112199816B/zh
Publication of CN112199816A publication Critical patent/CN112199816A/zh
Application granted granted Critical
Publication of CN112199816B publication Critical patent/CN112199816B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40—Business processes related to the transportation industry
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Business, Economics & Management (AREA)
  • Primary Health Care (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Human Resources & Organizations (AREA)
  • General Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种航班ADS‑B数据模拟方法，包括：S1，采集全球真实航班数据，并根据航班数据，计算飞机的飞行轨迹、飞行速度和飞行高度；S2，更新仿真时间，计算当前仿真时间飞机和卫星的位置，飞机的航向，计算当前仿真时间卫星覆盖范围，统计卫星覆盖范围内的飞机；S3，根据卫星覆盖范围内的飞机信息生成ADS‑B报文数据；S4，将卫星覆盖范围内的报文数据发送至该卫星；S5，重复步骤S2‑S4直至仿真结束。本发明的航班ADS‑B数据模拟方法，可以实现全球航班ADS‑B数据模拟生成并发送，成功模拟星间组网和数据实时传输过程，通过对ADS‑B数据进行接收及分析，可以在无实物条件下完成星间网络研究，指导星间网络设计，保证有效的航班监控。

Description

航班ADS-B数据模拟方法
技术领域
本发明涉及航空安全与监视技术领域，尤其涉及一种航班ADS-B数据模拟方法。
背景技术
为了空管安全，广播式自动相关监视(ADS-B)是近年来广泛使用的航班监控技术，通过实时向外广播自身的身份ID、位置、速度、高度、航向等信息，飞机之间就可以了解到相互的信息。与雷达相比ADS-B技术具有位置准确、信息更加丰富的特点。通过构建低轨卫星星座，并搭载星基ADS-B载荷，可实现全球无缝实时覆盖，通过星间链路将数据回传，可实现全球任意一点ADS-B数据近实时的数据采集、下传。
广域航空安全与监视技术项目要求，ADS-B数据自飞机发出，经星间链路传输到落地的时间不超过10秒，才能保证有效的航班监控。星间网络的设计对于上述指标的实现至关重要。由于天基***研制周期长、风险大、成本高，因此搭建实际的卫星***开展设计研究和实验是不现实的，只能通过高可信的半实物仿真进行。开展星间网络研究，可信的数据输入是基本条件，因此必须能够模拟全球的民航班机及其ADS-B报文发送。
发明内容
为解决上述现有技术中存在的技术问题，本发明提供了一种航班ADS-B数据模拟方法。具体技术方案如下：
一种航班ADS-B数据模拟方法，所述方法包括：
S1，采集全球真实航班数据，根据航班数据，计算飞机的飞行轨迹、飞行速度和飞行高度；
S2，更新仿真时间，计算当前仿真时间飞机和卫星的位置，飞机的航向，计算当前仿真时间卫星覆盖范围，统计卫星覆盖范围内的飞机信息；
S3，根据卫星覆盖范围内的飞机信息生成ADS-B报文数据；
S4，将卫星覆盖范围内的报文数据发送至该卫星；
S5，重复步骤S2-S4直至仿真结束。
在一种可能的设计中，真实航班数据包括动态航班数据和静态航班数据，动态航班数据包括：航班号、起点、起点位置、终点、终点位置、起飞时刻和降落时刻，静态航班数据包括以下信息：ICAO地址、纬度值、经度值和初始航向值；根据动态航班数据获得以下飞机信息：航班号、时刻、飞行状态、飞机位置、航向和航速；根据静态航班数据获得以下飞机信息：飞机位置和初始航向。
在一种可能的设计中，飞行状态根据起飞时刻、降落时刻和当前时刻确定；
飞机位置根据起点位置、终点位置和当前时刻确定；
航向通过当前位置和终点位置确定；
航速通过起点位置、终点位置、起飞时刻和降落时刻确定；
飞行高度在民航合法飞行高度范围内随机确定。
在一种可能的设计中，确定卫星覆盖飞机范围，包括：
通过卫星轨道文件，确定星下点经纬度(Lo，Lα)、幅角α、卫星高度h、地球半径R，经纬度及幅角单位为度，高度及半径单位为米；
计算中心角
给定某飞机位置P点坐标(Po,Pα,h)，计算P点与星下点形成的球心角距离，并转化为角度制：
判断θ与β大小关系，当0≤θ≤β，认为P点在该卫星覆盖域内。
在一种可能的设计中，ADS-B报文长度为112bit，通用格式如下：
其中，DF＝17格式用于S模式应答机发射的ADS-B消息，此格式下，CA字段表征了S模式应答机的能力；AA字段为应答机24个比特ICAO地址；ME字段包含ADS-B消息的主体；PI字段为奇偶校验位；
DF＝18格式用于非S模式应答机发射的ADS-B消息或TIS-B消息，此格式下，CF字段表明ME字段是含有ADS-B消息或者TIS-B消息，当CF＝0或1时，此格式传输的是ADS-B消息。
在一种可能的设计中，CA字段编码含义如下表所示：
在一种可能的设计中，ME字段前5个比特包含TYPE Code子字段，对于TYPE Code子字段某个值，紧接着的3个比特包含相应的SUBTYPE Code子字段，TYPE Code子字段与SUBTYPECode子字段共同决定ME字段剩余部分传送的ADS-B或TIS-B消息内容，如下表所示：
在一种可能的设计中，空中位置信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为9～18；
监视状态子字段编码含义如下：
单天线子字段用于表明ADS-B发射子***是否运行时采用了单个天线，单天线子字段编码值为0，表示双天线通道运行，单天线子字段编码值为1，表示仅单天线通道运行；
高度子字段表示ADS-B发射子***的高度，高度子字段对高度数据进行编码，如下表所示：
MSG比特位 41 42 43 44 45 46 47 48 49 50 51 52
ME比特位 9 10 11 12 13 14 15 16 17 18 19 20
码元比特位 C1 A1 C2 A2 C4 A4 B1 “Q” B2 D2 B4 D4
比特位16被指定为Q比特，Q＝0时，增量为100英尺，高度组合顺序为D2D4A1A2A4B1B2B4 C1C2C4，D2～B4为格雷码编码，最小递增单位500英尺，C组码为5周期循环码，递增单位为100英尺，以-1200英尺为起点；Q＝1时，增量为25英尺，组合顺序为C1A1 C2A2 C4A4B1B2 D2 B4D4，C1～D4为格雷码编码，递增单位为25英尺；
时间子字段表示位置数据中适用时间是否与准确的0.2秒UTC时间点同步；
CPR格式子字段用于指出紧凑式位置报告(CPR)格式类型(“偶”或“奇”)对经度与纬度数据进行编码。
在一种可能的设计中，地表位置信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为5～8；
运行子字段编码的含义如下，并取运行编码＝0：
编码(十进制) 含义 量化
0 无可用运行信息
1 飞机停止(地面速度<0.2315km/h(0.125哩/h))
2-8 0.2315km/h(0.125哩/h)<地面速度<1.852km/h(1哩/h) 0.2315km/h(0.125哩/h))步长
9-12 1.852km/h(1哩/h)<地面速度<3.704km/h(2哩/h) 0.463km/h(0.25哩/h))步长
13-38 3.704km/h(2哩/h)<地面速度<27.78km/h(15哩/h) 0.926km/h(0.50哩/h))步长
39-93 27.78km/h(15哩/h)<地面速度<129.64km/h(70哩/h) 1.852km/h(1哩/h))步长
94-108 129.64km/h(70哩/h)<地面速度<185.2km/h(100哩/h) 3.704km/h(2哩/h))步长
109-123 185.2km/h(100哩/h)<地面速度<324.1km/h(175哩/h) 9.26km/h(5哩/h))步长
124 324.1km/h(175哩/h)<地面速度 0.2315km/h(0.125哩/h))步长
125 为飞机减速预留
126 为飞机加速预留
127 为飞机转向预留
航向/地面航迹状态子字段的编码含义如下：
编码 含义
0 航向/地面航迹数据无效
1 航向/地面航迹数据有效
航向/地面航迹子字段用于报告航向，或者自北开始顺时针方向(航向符号标志位＝0)ADS-B发射子***的运动，航向/地面航迹子字段的编码含义如下：
在一种可能的设计中，空中速度信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为19；
子类型子字段编码值为1；
意图不变标志子字段、IFR能力标志子字段、导航不确定度子字段编码值均为0；
东/西向标志位子字段表示速度在东/西向的速度矢量，编码值为0时表示东向，编码值为1时表示西向；
东/西速度子字段中编码值的含义如下：
北/南向标志位子字段表示速度在北/南向的速度矢量，编码值为0时表示北向，编码值为1时表示南向；
北/南速度子字段中编码值的含义如下：
本发明技术方案的主要优点如下：
本发明的航班ADS-B数据模拟方法，可以实现全球航班ADS-B数据模拟生成并发送，成功模拟星间组网和数据实时传输过程，通过对ADS-B数据进行接收及分析，可以在无实物条件下完成星间网络研究，指导星间网络设计，保证有效的航班监控。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明一实施例提供的航班ADS-B数据模拟方法的流程图；
图2为本发明一实施例提供的航向示意图；
图3为本发明一实施例提供的卫星覆盖范围示意图；
图4为本发明一实施例提供的动态航班数据发送过程流程图；
图5为本发明一实施例提供的静态航班数据发送流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明的一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
以下结合附图，详细说明本发明实施例提供的技术方案。
本发明实施例提供了一种航班ADS-B数据模拟方法，该方法包括：
S1，采集全球真实航班数据，并根据航班数据，计算飞机的飞行轨迹、飞行速度和飞行高度。其中，航班数据可以包括起飞机场、目的机场、起飞时间和降落时间等。
S2，更新仿真时间，获取当前仿真时间飞机和卫星的位置，飞机的航向，计算当前仿真时间卫星覆盖，统计卫星覆盖范围内的飞机。
S3，根据卫星覆盖范围内的飞机信息生成ADS-B报文数据。其中，飞机信息可以包括飞机的飞行轨迹、飞行速度和飞行高度以及位置和航向等，其可由航班数据经过计算获得。
S4，将卫星覆盖范围内的报文数据发送至该卫星。
S5，重复步骤S2-S4直至仿真结束。
本发明实施例提供的航班ADS-B数据模拟方法，可以实现全球航班ADS-B数据模拟生成并发送，成功模拟星间组网和数据实时传输过程，通过对ADS-B数据进行接收及分析，可以在无实物条件下完成星间网络研究，指导星间网络设计，保证有效的航班监控。
具体地，以下对本发明实施例提供的航班ADS-B数据模拟方法的各步骤进行详细阐述：
步骤S1中，全球真实航班数据分为动态航班数据和静态航班数据，动态航班数据为在预设统计时限内发生飞行的飞机的数据，静态航班数据为在预设统计时限内一直处于静止状态的飞机的数据，其数据统计与处理方法如下：
(a)动态航班数据
当前通过数据采集，已收集到全球超过1000个主要机场的超过170万条航班信息。每条动态航班数据包含以下内容：
各信息字段含义如下：
航班号No：航班的标识信息，与原始数据库包含信息一致，生成全球航班飞行数据库前根据航班号去重，与ADS-B中ICAO地址一一对应。
时刻t：仿真历元，本发明实施例中时间更新步长设置为1s。
飞行状态flag：表示飞行是否处于飞行状态，根据该航班起飞、降落时刻与当前时刻判断；当TSrc<t<TDst，flag＝1,表示当前时刻飞机处于飞行状态，否则flag＝0，表示飞机处于未飞行状态，假设飞机停泊于起点，且不发出ADS-B信号。
飞机位置(纬度lat,经度lon,高度h)，
对于处在飞行状态的飞机，起点、终点坐标分别为(Latsrc，Lonsrc)和(Latdst，Londst)，其在时刻t的位置为(Lat(t),Lon(t))，计算方法如下：
先将起点、终点的经纬度转换成弧度制，以下计算都采用弧度角度计算，两点经纬度差分别为，ΔLon＝Londst-Lonsrc，ΔLat＝Latdst-Latsrc。
a＝sin2(Δlat/2)+cos(Latsrc)cos(Latdst)sin2(ΔLon/2)
两点间地心角为：
笛卡尔单位向量：
v1＝[cos(Lonsrc)cos(Latsrc)sin(Lonsrc)cos(Latsrc)sin(Latsrc)]
v2＝[cos(Londst)cos(Latdst)sin(Londst)cos(Latdst)sin(Latdst)]
求单位正切向量：
n1＝v1×v2(其中×为向量叉乘运算)
t1＝n1×v1(其中×为向量叉乘运算)
将t1单位化,可得
则t时刻的飞机所在经纬度Lon(t)和Lat(t)为：
Lon(t)＝arctan(p2(t)/p1(t))
注意判断反正切函数角度象限：
Lat(t)＝arcsin(p3(t))
再将Lon(t)和Lat(t)换算成角度制输出。
由于高度信息对仿真影响不大，因此假设所有飞机处于6000～1200km范围内飞行，赋予随机值。
航向Az：
航向方位角Az取值范围为0～360°，由起点终点连线确定方向，并假设保持恒定。
假设t时刻飞机从起点A(Latsrc,Lonsrc)朝向终点B(Latdst,Londst)飞行，则通过下式可得到A和B点夹角c的正弦值。
cos(c)＝cos(π/2-Latdst)*cos(π/2-Latsrc)+sin(π/2-Latdst)*sin(π/2-Latsrc)*cos(Londst-Lonsrc)
因此方位角正弦值为：
利用反三角函数求解Az时，会出现象限问题。由于方位角是以正北为0度起点，由东向南向西顺时针旋转360度。因此，若A点固定于原点，则将Az转化为0～360°范围时，有如下判断：
航速v：
根据飞机起点和终点的经纬度坐标，以(Latsrc，Lonsrc)和(Latdst，Londst)标记为某航班飞行的起点和终点的经纬度。Re为地球平均半径，取6371km，计算过程中均用弧度表示。根据以下公式可以计算出起点与终点的距离d：
d＝REarccos(cos(Latsrc)cos(Latdst)cos(Londst-Lonsrc)+sin(Latsrc)sin(Latdst))
忽略飞机起飞、降落过程，假设该飞机作匀速运动，则其飞行速度大小为v：
各方向分量为：
(b)静态航班数据
静态数据来源于天拓三号ADS-B数据库中的地表面数据信息，目前地表面信息数据为25998条。每条静态航班数据包含内容如下:
由于该数据为地表面信息数据，则设定该飞机的飞行高度、航速均为0。
根据静态航班数据，可以直接获得处于静止状态的飞机的位置和初始航向。
步骤S2中，确定卫星覆盖飞机范围的过程如下：
通过卫星轨道文件，已知星下点经纬度(Lo，Lα)、幅角α、卫星高度h、地球半径R，经纬度及幅角单位为度，高度及半径单位为米。
中心角
如果波束与地面有交点，则上式求解出的β值一定在[0,pi/2]内，则不用判断象限。
给定某飞机位置P点坐标(Po,Pα,h)，计算P点与星下点形成的球心角距离，并转化为角度制：
上式cos的取值范围为[0,pi]。
判断θ与β大小关系，当0≤θ≤β，认为P点在该卫星覆盖域内。
步骤S3，ADS-B报文生成
ADS-B报文长度为112bit，通用格式如下：
其中，DF＝17格式用于S模式应答机发射的ADS-B消息，此格式下，CA字段表征了S模式应答机的能力；AA字段为应答机24个比特ICAO地址；ME字段包含ADS-B消息的主体；PI字段为奇偶校验位。
DF＝18格式用于非S模式应答机发射的ADS-B消息或TIS-B消息，此格式下，CF字段表明ME字段是含有ADS-B消息或者TIS-B消息，当CF＝0或1时，此格式传输的是ADS-B消息。
本发明实施例中假设只含有S模式应答机，因此设DF＝17。
CA代码：
CA字段编码含义如下表所示：
本发明实施例中由于模拟的数据源是民航飞机，因此CA部分只选取4种场景CA＝4、5、6、7。
上表中A通信表示地对空通信，B通信表示空对地通信。
ICAO地址：
ICAO地址为24位二进制码，是分配给S模式应答机的编码，是飞机在ICAO注册的唯一识别代码。为简化设计方案，本发明实施例中ICAO号与后续航班号一一对应。
ICAO号全为0或者1代表发射器通告发射失败。
ADS-B消息主体ME字段：
ME字段前5个比特包含TYPE Code子字段。对于TYPE Code子字段某个值，紧接着的3个比特包含相应的SUBTYPE Code子字段。TYPE Code子字段与SUBTYPECode子字段(若存在该子字段)共同决定了ME字段剩余部分传送的ADS-B或TIS-B消息内容。
ADS-B消息类型的确定
其余值同样有各自类型，但出于简化仿真的角度而删去，本发明实施例中只发送地表位置消息、空中位置消息及空中速度消息。
(a)空中位置信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为9～18。本发明实施例中采用统一值为9。
监视状态子字段编码含义如下：
单天线子字段用于表明ADS-B发射子***是否运行时采用了单个天线，单天线子字段编码值为0，表示双天线通道运行，单天线子字段编码值为1，表示仅单天线通道运行。
高度子字段表示ADS-B发射子***的高度，高度子字段对高度数据进行编码，如下表所示：
MSG比特位 41 42 43 44 45 46 47 48 49 50 51 52
ME比特位 9 10 11 12 13 14 15 16 17 18 19 20
码元比特位 C1 A1 C2 A2 C4 A4 B1 “Q” B2 D2 B4 D4
比特位16被指定为Q比特，Q＝0时，增量为100英尺，高度组合顺序为D2D4A1A2A4B1B2B4 C1C2C4，D2～B4为格雷码编码，最小递增单位500英尺，C组码为5周期循环码，递增单位为100英尺，以-1200英尺为起点；Q＝1时，增量为25英尺，组合顺序为C1A1 C2A2 C4A4B1B2 D2 B4D4，C1～D4为格雷码编码，递增单位为25英尺。
时间子字段表示位置数据中适用时间是否与准确的0.2秒UTC时间点同步。
CPR格式子字段用于指出紧凑式位置报告(CPR)格式类型(“偶”或“奇”)对经度与纬度数据进行编码，编码规则参照《ADS-B_1090MHz_ES标准.pdf》。
编码纬度子字段中的空中纬度位置数据应按照《Minimum OperationalPerformance Standards(MOPS)for 1090MHz Automatic Dependent Surveillance-Broadcast(ADS-B)》进行编码。
编码经度子字段中的空中经度位置数据应按照《Minimum OperationalPerformance Standards(MOPS)for 1090MHz Automatic Dependent Surveillance-Broadcast(ADS-B)》进行编码。
(b)地表位置信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为5～8，本发明实施例中采用统一值为5。
运行子字段编码的含义如下：
编码(十进制) 含义 量化
0 无可用运行信息
1 飞机停止(地面速度<0.2315km/h(0.125哩/h))
2-8 0.2315km/h(0.125哩/h)<地面速度<1.852km/h(1哩/h) 0.2315km/h(0.125哩/h))步长
9-12 1.852km/h(1哩/h)<地面速度<3.704km/h(2哩/h) 0.463km/h(0.25哩/h))步长
13-38 3.704km/h(2哩/h)<地面速度<27.78km/h(15哩/h) 0.926km/h(0.50哩/h))步长
39-93 27.78km/h(15哩/h)<地面速度<129.64km/h(70哩/h) 1.852km/h(1哩/h))步长
94-108 129.64km/h(70哩/h)<地面速度<185.2km/h(100哩/h) 3.704km/h(2哩/h))步长
109-123 185.2km/h(100哩/h)<地面速度<324.1km/h(175哩/h) 9.26km/h(5哩/h))步长
124 324.1km/h(175哩/h)<地面速度 0.2315km/h(0.125哩/h))步长
125 为飞机减速预留
126 为飞机加速预留
127 为飞机转向预留
本发明实施例中运行子字段的编码值取0，无情景信息。
航向/地面航迹状态子字段的编码含义如下：
编码 含义
0 航向/地面航迹数据无效
1 航向/地面航迹数据有效
本发明实施例中取航向/地面航迹状态子字段的编码值为1。
航向/地面航迹子字段用于报告航向，或者自北开始顺时针方向(航向符号标志位＝0)ADS-B发射子***的运动，航向/地面航迹子字段的编码含义如下：
其他子字段与空中位置信息定义一致，经纬度解码时消息编码Nb＝19，参见《ADS-B_1090MHz_ES标准.pdf》。
(c)空中速度信息ME字段含义如下：
其中，类型TYPE Code子字段编码对应值为19；
子类型子字段编码值为1。
意图不变标志子字段、IFR能力标志子字段、导航不确定度子字段编码值均为0。
东/西向标志位子字段表示速度在东/西向的速度矢量，编码值为0时表示东向，编码值为1时表示西向。
东/西速度子字段中编码值的含义如下：
北/南向标志位子字段表示速度在北/南向的速度矢量，编码值为0时表示北向，编码值为1时表示南向。
北/南速度子字段中编码值的含义如下：
其余子字段的取值为0。
基于上述，本发明实施例提供的航班ADS-B数据模拟方法可以分为动态航班数据发送和静态航班数据发送，分别参见附图4和5。
利用本发明实施例提供的航班ADS-B数据模拟方法进行星间网络设计的研究时，可以假定卫星网络模型，并建立相应地卫星轨道文件，确定每个卫星的星下点位置和覆盖范围等，覆盖范围内飞机发出的ADS-B报文可由ADS-B数据接收软件进行接收并分析，确定卫星网络模型的可行性。
需要说明的是，在本文中，诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。此外，本文中“前”、“后”、“左”、“右”、“上”、“下”均以附图中表示的放置状态为参照。
最后应说明的是：以上实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种航班ADS-B数据模拟方法，其特征在于，所述方法包括：
S1，采集全球真实航班数据，根据航班数据，计算飞机的飞行轨迹、飞行速度和飞行高度；
S2，更新仿真时间，计算当前仿真时间飞机和卫星的位置，飞机的航向，计算当前仿真时间卫星覆盖范围，统计卫星覆盖范围内的飞机；
S3，根据卫星覆盖范围内的飞机信息生成ADS-B报文数据；
S4，将卫星覆盖范围内的报文数据发送至该卫星；
S5，重复步骤S2-S4直至仿真结束；
其中，真实航班数据包括动态航班数据和静态航班数据，动态航班数据为在预设统计时限内发生飞行的飞机的数据，静态航班数据为在预设统计时限内一直处于静止状态的飞机的数据，动态航班数据包括：航班号、起点、起点位置、终点、终点位置、起飞时刻和降落时刻，静态航班数据包括以下信息：ICAO地址、纬度值、经度值和初始航向值；
根据动态航班数据获得以下飞机信息：航班号、时刻、飞行状态、飞机位置、航向和航速；
根据静态航班数据获得以下飞机信息：飞机位置和初始航向；
飞行状态根据起飞时刻、降落时刻和当前时刻确定；
飞机位置根据起点位置、终点位置和当前时刻确定；
航向通过当前位置和终点位置确定；
航速通过起点位置、终点位置、起飞时刻和降落时刻确定；
飞行高度在民航合法飞行高度范围内随机确定；
对于处在飞行状态的飞机，起点、终点坐标分别为(Latsrc，Lonsrc)和(Latdst，Londst)，其在时刻t的位置为(Lat(t),Lon(t))，纬度lat，经度lon，飞机的起飞时刻为Tsrc，降落时刻为Tdst，计算方法如下：
先将起点、终点的经纬度转换成弧度制，以下计算都采用弧度角度计算，两点经纬度差分别为，ΔLon＝Londst-Lonsrc，ΔLat＝Latdst-Latsrc；
a＝sin2(Δlat/2)+cos(Latsrc)cos(Latdst)sin2(ΔLon/2)
两点间地心角为：
笛卡尔单位向量：
v1＝[cos(Lonsrc)cos(Latsrc) sin(Lonsrc)cos(Latsrc) sin(Latsrc)]
v2＝[cos(Londst)cos(Latdst) sin(Londst)cos(Latdst) sin(Latdst)]
求单位正切向量：
n1＝v1×v2，其中×为向量叉乘运算
t1＝n1×v1，其中×为向量叉乘运算
将t1单位化,可得
则t时刻的飞机所在经纬度Lon(t)和Lat(t)为：
Lon(t)＝arctan(p2(t)/p1(t))
注意判断反正切函数角度象限：
Lat(t)＝arcsin(p3(t))
再将Lon(t)和Lat(t)换算成角度制输出。
2.根据权利要求1所述的航班ADS-B数据模拟方法，其特征在于，确定卫星覆盖飞机范围，包括：
通过卫星轨道文件，确定星下点经纬度(Lo，Lα)、幅角α、卫星高度h、地球半径R，经纬度及幅角单位为度，高度及半径单位为米；
计算中心角
给定某飞机位置P点坐标(Po,Pα,h)，计算P点与星下点形成的球心角距离，并转化为角度制：
判断θ与β大小关系，当0≤θ≤β，认为P点在该卫星覆盖域内。
3.根据权利要求2所述的航班ADS-B数据模拟方法，其特征在于，ADS-B报文长度为112bit，通用格式如下：
其中，DF＝17格式用于S模式应答机发射的ADS-B消息，此格式下，CA字段表征了S模式应答机的能力；AA字段为应答机24个比特ICAO地址；ME字段包含ADS-B消息的主体；PI字段为奇偶校验位；
DF＝18格式用于非S模式应答机发射的ADS-B消息或TIS-B消息，此格式下，CF字段表明ME字段是含有ADS-B消息或者TIS-B消息，当CF＝0或1时，此格式传输的是ADS-B消息。
4.根据权利要求3所述的航班ADS-B数据模拟方法，其特征在于，CA字段编码含义如下表所示：
5.根据权利要求3所述的航班ADS-B数据模拟方法，其特征在于，ME字段前5个比特包含TYPE Code子字段，对于TYPE Code子字段某个值，紧接着的3个比特包含相应的SUBTYPECode子字段，TYPE Code子字段与SUBTYPECode子字段共同决定ME字段剩余部分传送的ADS-B或TIS-B消息内容，如下表所示：
6.根据权利要求5所述的航班ADS-B数据模拟方法，其特征在于，空中位置信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为9～18；
监视状态子字段编码含义如下：
单天线子字段用于表明ADS-B发射子***是否运行时采用了单个天线，单天线子字段编码值为0，表示双天线通道运行，单天线子字段编码值为1，表示仅单天线通道运行；
高度子字段表示ADS-B发射子***的高度，高度子字段对高度数据进行编码，如下表：
MSG比特位 41 42 43 44 45 46 47 48 49 50 51 52 ME比特位 9 10 11 12 13 14 15 16 17 18 19 20 码元比特位 C1 A1 C2 A2 C4 A4 B1 “Q” B2 D2 B4 D4
比特位16被指定为Q比特，Q＝0时，增量为100英尺，高度组合顺序为D2D4 A1A2A4B1B2B4 C1C2C4，D2～B4为格雷码编码，最小递增单位500英尺，C组码为5周期循环码，递增单位为100英尺，以-1200英尺为起点；Q＝1时，增量为25英尺，组合顺序为C1A1 C2A2 C4A4B1B2 D2 B4D4，C1～D4为格雷码编码，递增单位为25英尺；
时间子字段表示位置数据中适用时间是否与准确的0.2秒UTC时间点同步；
CPR格式子字段用于指出紧凑式位置报告格式类型对经度与纬度数据进行编码，格式类型包括“偶”或“奇”。
7.根据权利要求5所述的航班ADS-B数据模拟方法，其特征在于，地表位置信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为5～8；
运行子字段编码的含义如下，并取运行编码＝0：
十进制编码 含义 量化 0 无可用运行信息 1 飞机停止，地面速度<0.2315km/h(0.125哩/h) 2-8 0.2315km/h(0.125哩/h)<地面速度<1.852km/h(1哩/h) 0.2315km/h(0.125哩/h)步长 9-12 1.852km/h(1哩/h)<地面速度<3.704km/h(2哩/h) 0.463km/h(0.25哩/h)步长 13-38 3.704km/h(2哩/h)<地面速度<27.78km/h(15哩/h) 0.926km/h(0.50哩/h)步长 39-93 27.78km/h(15哩/h)<地面速度<129.64km/h(70哩/h) 1.852km/h(1哩/h)步长 94-108 129.64km/h(70哩/h)<地面速度<185.2km/h(100哩/h) 3.704km/h(2哩/h)步长 109-123 185.2km/h(100哩/h)<地面速度<324.1km/h(175哩/h) 9.26km/h(5哩/h)步长 124 324.1km/h(175哩/h)<地面速度 0.2315km/h(0.125哩/h)步长 125 为飞机减速预留 126 为飞机加速预留 127 为飞机转向预留
航向/地面航迹状态子字段的编码含义如下：
编码 含义 0 航向/地面航迹数据无效 1 航向/地面航迹数据有效
航向/地面航迹子字段用于报告航向，或者自北开始顺时针方向ADS-B发射子***的运动，航向/地面航迹子字段的编码含义如下：
8.根据权利要求5所述的航班ADS-B数据模拟方法，其特征在于，空中速度信息ME字段的含义如下：
其中，类型TYPE Code子字段编码对应值为19；
子类型子字段编码值为1；
意图不变标志子字段、IFR能力标志子字段、导航不确定度子字段编码值均为0；
东/西向标志位子字段表示速度在东/西向的速度矢量，编码值为0时表示东向，编码值为1时表示西向；
东/西速度子字段中编码值的含义如下：
北/南向标志位子字段表示速度在北/南向的速度矢量，编码值为0时表示北向，编码值为1时表示南向；
北/南速度子字段中编码值的含义如下：
CN202010893875.4A 2020-08-31 2020-08-31 航班ads-b数据模拟方法 Active CN112199816B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010893875.4A CN112199816B (zh) 2020-08-31 2020-08-31 航班ads-b数据模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010893875.4A CN112199816B (zh) 2020-08-31 2020-08-31 航班ads-b数据模拟方法

Publications (2)

Publication Number Publication Date
CN112199816A CN112199816A (zh) 2021-01-08
CN112199816B true CN112199816B (zh) 2023-08-15

Family

ID=74005172

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010893875.4A Active CN112199816B (zh) 2020-08-31 2020-08-31 航班ads-b数据模拟方法

Country Status (1)

Country Link
CN (1) CN112199816B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113890609B (zh) * 2021-09-28 2022-08-05 南京先进激光技术研究院 基于ads-b接收机报文数据的激光信道安全预警***和方法
CN114143001B (zh) * 2021-12-01 2024-04-26 中国电子科技集团公司第二十八研究所 一种基于公钥密码算法的ads-b防欺骗方法
CN115209358B (zh) * 2022-07-15 2024-04-09 深圳市道通智能航空技术股份有限公司 自动相关监视广播信息的发送方法、装置和存储介质
CN115826464B (zh) * 2022-11-29 2024-03-22 航科院中宇（北京）新技术发展有限公司 一种远机位节点采集***的采集方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202549080U (zh) * 2012-03-16 2012-11-21 中国民用航空总局第二研究所 一种雷达数据、飞行计划数据与ads-b数据融合***
CN104777464A (zh) * 2015-04-20 2015-07-15 中国民航大学 利用ads-b数据实现多雷达数据半实物仿真的方法
CN108021144A (zh) * 2017-12-29 2018-05-11 中国地质大学（武汉） 一种无人机航迹规划及动态威胁规避仿真设备
CN108615411A (zh) * 2016-12-13 2018-10-02 北京航管科技有限公司 一种实现航班信息处理的方法及装置
CN110224743A (zh) * 2019-06-24 2019-09-10 中国电子科技集团公司第五十四研究所 一种基于动态场景的星载ads-b测试信号的形成方法
CN111160770A (zh) * 2019-12-30 2020-05-15 南京航空航天大学 一种进场航班动态协同排序方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7107146B2 (en) * 2003-12-10 2006-09-12 Honeywell International Inc. Methods and systems for generating a terrain elevation map in a cartesian format
US10659145B2 (en) * 2017-01-11 2020-05-19 Aireon Llc Simulating reception of transmissions
US10650687B2 (en) * 2018-05-18 2020-05-12 Aireon Llc Decoding position information in space-based systems
US11024182B2 (en) * 2018-12-17 2021-06-01 The Boeing Company System and method to forecast flight delay based on real-time data

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202549080U (zh) * 2012-03-16 2012-11-21 中国民用航空总局第二研究所 一种雷达数据、飞行计划数据与ads-b数据融合***
CN104777464A (zh) * 2015-04-20 2015-07-15 中国民航大学 利用ads-b数据实现多雷达数据半实物仿真的方法
CN108615411A (zh) * 2016-12-13 2018-10-02 北京航管科技有限公司 一种实现航班信息处理的方法及装置
CN108021144A (zh) * 2017-12-29 2018-05-11 中国地质大学（武汉） 一种无人机航迹规划及动态威胁规避仿真设备
CN110224743A (zh) * 2019-06-24 2019-09-10 中国电子科技集团公司第五十四研究所 一种基于动态场景的星载ads-b测试信号的形成方法
CN111160770A (zh) * 2019-12-30 2020-05-15 南京航空航天大学 一种进场航班动态协同排序方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于BDS的ADS-B报文生成***;倪育德;张心一;刘瑞华;马圆晨;刘鹏;;电讯技术(11);1271-1277 *

Also Published As

Publication number Publication date
CN112199816A (zh) 2021-01-08

Similar Documents

Publication Publication Date Title
CN112199816B (zh) 航班ads-b数据模拟方法
Sun The 1090 megahertz riddle: a guide to decoding mode S and ADS-B signals
US6182005B1 (en) Airport guidance and safety system incorporating navigation and control using GNSS compatible methods
US6314363B1 (en) Computer human method and system for the control and management of an airport
Magazu III Exploiting the automatic dependent surveillance-broadcast system via false target injection
US20070021908A1 (en) System and method for data mapping and map discrepancy reporting
US10659145B2 (en) Simulating reception of transmissions
MXPA03003479A (es) Enlace de datos integrados en un receptor de vigilancia.
CN1592919A (zh) 确定活动飞行器之间的碰撞路径的方法和相关***及计算机软件程序产品
CN100449329C (zh) 一种符合广播式自动相关监视报文处理要求的数据格式转换方法
CN104777464B (zh) 利用ads‑b数据实现多雷达数据半实物仿真的方法
CN113553983B (zh) 一种星载ads-b与遥感图像结合的异常目标监测方法
CN109523863A (zh) 一种基于ads-b的雷达训练模拟***
Mott Measurement of airport operations using a low-cost transponder data system
CN104898142B (zh) 一种飞行器或场地车辆运行状态的监视设备及其方法
CN111768653B (zh) 一种ads-b测试数据模拟的方法
Sun The 1090 Megahertz Riddle
Baker et al. Mode S system design and architecture
Varga et al. ADS-B based real-time air traffic monitoring system
WO2023030641A1 (en) Aircraft collision avoidance method and device
Ostroumov et al. Software-Defined Transmitter to Support Automatic Dependent Surveillance-Broadcast
US8836571B2 (en) Method for transmission of a geographic coordinate
Pereira Space Based ADS-B receiver for ISTsat-1
CN116778759B (zh) 一种飞机实时监控***及方法
Ostroumov Spatial Indexing of Airplane Trajectory with Open Location Code and Hexagonal Hierarchical Index

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Yang Lei

Inventor after: Zhang Fang

Inventor after: Zhao Yong

Inventor after: Song Xin

Inventor after: Guo Jianming

Inventor after: Liu Xianfeng

Inventor after: Chen Quan

Inventor after: Li Xuan

Inventor before: Yang Lei

Inventor before: Zhang Fang

Inventor before: Zhao Yong

Inventor before: Song Xin

Inventor before: Guo Jianming

Inventor before: Liu Xianfeng

Inventor before: Chen Quan

Inventor before: Li Xuan

GR01 Patent grant
GR01 Patent grant