基于航空器性能的ADS-B航迹去噪方法

1.一种基于航空器性能的ADS-B航迹去噪方法，其特征在于，包括：建立航班时刻信息、航班ADS-B轨迹数据和民用航空器机型性能数据的数据库；获取指定航班的关键时刻信息和ADS-B航迹数据；对该航班的ADS-B航迹数据进行逻辑去噪；对该航班的ADS-B航迹数据进行数值去噪；所述航班时刻信息包括航空器机型信息、实际起飞时间ATOT和实际落地时间ALDT；所述民用航空器机型性能数据包括机型基本性能参数如最小操作重量W op,min 、参考重量W ref 、起飞构型失速速度V s,ref_TO 以及最大巡航压力高度HP max ；所述获取指定航班的关键时刻信息和ADS-B航迹数据的方法包括：根据航迹点监视时间进行升序排列；剔除具有相同监视时间的其他轨迹点，即在一个监视时间下只保留一个航迹点；对于场面航迹点数据，即字段OnGround＝1，进行删除；所述对该航班的ADS-B航迹数据进行逻辑去噪的方法包括：依据ADS-B航迹数据获取航班的ADS-B航迹点队列{A 0 ,A 1 ,A 2 ,A 3 ,…A n }，其中，A n 表示第n个ADS-B航迹点；航迹点A i ，判断其监视时间MTime i 是否在实际起飞时间ATOT到实际着陆时间ALDT范围内，即，MTime i ∈[ATOT,ALDT]，若不满足此条件，则航路点A i 判为噪点，直接从ADS-B航迹点队列中剔除，其中，A i 表示第i个ADS-B航迹点；根据机型性能定义最低飞行高度HP min 与最高巡航高度HP max ，对ADS-B航迹点A i ，判断其压力高度HP i 是否在最低飞行高度与最高巡航高度之间，即，HP i ∈[HP min ,HP max ]，若不满足此条件，则航路点A i 判为噪点，直接从ADS-B航迹点队列中剔除；对航迹点A i ，根据该点的压力高度HP i ，先计算出此高度上国际标准大气条件ISA下的大气温度T i ，最后计算出音速a i,ISA ，即：其中空气绝热系数k＝1.4，气体常数R＝287.05287m 2 /(K·s 2 )；判断地速GS i 是否在最小失速速度V s,min 与此高度上的音速a i,ISA 之间，即，GS i ∈[V s,min ,a i,ISA ],若不满足此条件，则航路点A i 判为噪点，直接从ADS-B航迹点队列中剔除。 2.如权利要求1所述的基于航空器性能的ADS-B航迹去噪方法，其特征在于，所述最小失速速度V s,min 的计算方法为：根据改航班机型的W op,min 、W ref 和V s,ref _ TO 计算得出该机型在起飞构型下的失速速度，即最小失速速度V s,min ： 3.如权利要求2所述的基于航空器性能的ADS-B航迹去噪方法，其特征在于，所述对该航班的ADS-B航迹数据进行数值去噪的方法包括：对航路点A i ，定义压力高度变化率R i,i+1 与其阈值函数F i,i+1 ，形式为：其中，RC max 和RD max 分别为最大上升率和最大下降率，且均为正数；计算出R 0,1 、R 0,2 和R 0,3 ，并求出F 0,1 +F 0,2 +F 0,3 ，如果结果值大于或等于2，则认为A 0 点为噪点，即A 0 高度值相对于A 1 和A 2 为噪点，此时将A 0 剔除以保证航迹点队列的第一个点为正常轨迹点，删除之后A 1 改为A 0 ，A 2 改为A 1 ，……，A n 改为A n-1 以此类推；计算出F i-1,i (i＝1,2,…,n)，如果F i-1,i ＝1，则轨迹点A i 判为噪点并将之剔除，依次类推，最后结束数值去噪过程。