首页 / 专利分类库 / 信号装置 / 一种基于意图高度的航空器质量估计方法

一种基于意图高度的航空器质量估计方法

申请号 CN202311665331.2 申请日 2023-12-06 公开(公告)号 CN117894210A 公开(公告)日 2024-04-16
申请人 中国电子科技集团公司第二十八研究所; 发明人 蒋伟煜; 刘洲洋; 肖英超; 马龙彪; 王晓雨;
摘要 本 发明 提出了一种基于意图高度的航空器 质量 估计方法。插值计算各航迹点 风 速以及风对航空器质量估计的影响;比较气压高度、意图高度的高度差与 阈值 的大小,过滤影响质量估计的异常航迹点;设置质量序列,计算不同质量下的质量误差,误差最小的质量作为航空器估计质量;设置高度层序列,计算不同高度层下的航空器估计质量,并进行数据校验,剔除异常质量,计算正常质量的均值得到稳定的航空器估计质量。本发明充分考虑风对航空器质量估计的影响,根据意图高度与高度缓冲时间对航迹点过滤,剔除影响质量结果的异常航迹点,通过实施多高度层航空器质量估计,获得更加稳定的航空器估计质量。本发明在航班四维轨迹预测领域有广阔的应用前景。
权利要求

1.一种基于意图高度的航空器质量估计方法,其特征在于,包括以下步骤:
步骤1,解析目标航空器的S模式雷达数据,得到该航空器在各航迹点的数据集TP={tpi,i=1,2,...,N},其中,tpi为第i个航迹点的数据;
解析气象数据,得到气象数据集合WD={windj,j=1,2,...,M},其中,windj表示第j个向量;
步骤2,插值计算航迹点数据tpi对应的UV坐标系中U方向和V方向的风速分量wind_ui和wind_vi;
步骤3,设置高度阈值ηh,计算航迹点数据中的气压高度Hp与意图高度MCP的高度差值Δh,比较高度差值Δh与高度阈值ηh的大小,添加缓冲高度标签flag;
步骤4,设置缓冲时间Δt,根据缓冲时间更新航迹点的数据集合TR;
步骤5,分解真空速TAS到U方向和V方向,计算风对目标航空器的影响;
步骤6,设置质量序列,计算不同质量对应的航空器的功率power、Q值以及单位质量下功率power和Q值的均方误差error,比较均方误差error数值大小获得目标航空器的估计质量massest;
步骤7,设置高度层序列,筛选各个高度层中的的航迹点数据,计算不同高度层对应的航空器估计质量massfl;
步骤8,获取目标航空器的最终估计质量:对所有高度层的航空器估计质量massfl进行数值校验,去除异常质量,得到有效质量序列,计算有效质量的平均质量作为目标航空器的最终估计质量。
2.根据权利要求1所述的基于意图高度的航空器质量估计方法,其特征在于,步骤1中所述的第i个航迹点的数据tpi,为一个向量,具体如下:
tpi=[ti,hpi,mcpi,tasi,θi]
其中,ti、hpi、mcpi、tasi和θi分别代表第i个航迹点对应的时间、气压高度、意图高度、真空速和航向,N为航迹点总数。
3.根据权利要求2所述的基于意图高度的航空器质量估计方法,其特征在于,步骤1中所述的第j个风向量windj,具体如下:
风向量windj由UV坐标系中的方向U和方向V两个方向的分量组成,记作:
windj=(wind_uj,wind_vj)
其中,每个方向的风向量都是一个基于时间与空间的四维矢量,分别记作:
wind_uj=(tj,lj,lonj,latj)
其中,tj、lj、lonj和latj分别表示第j个风向量对应的时间戳、高度层、经度和纬度。
4.根据权利要求3所述的基于意图高度的航空器质量估计方法,其特征在于,步骤2中所述的插值计算航迹点数据tpi对应的UV坐标系中U方向和V方向的风速分量wind_ui和wind_vi,具体方法如下:
步骤2.1,获取与航迹点数据tpi中的时间ti最近的风向量的时间戳tj对应的时间索引,记作indext;
步骤2.2,获取与航迹点数据tpi中的高度hpi最近的风向量的高度层lj对应的索引,其中,最近的上高度层索引记作indexh_top,最近的下高度层索引记作indexh_bottom;
步骤2.3,根据时间索引indext与高度层索引indexh_top和indexh_bottom,得到当前高度层下的风向量集合WDl,如下:
WDl={windt,l,lon,lat|t=indext,l∈[indexh_top,indexh_bottom]}
遍历上述集合中的向量,计算高度层风速windlon,lat,计算方法为:
其中,x和y分别表示当前高度层下的风向量windj的经度和纬度,Dx,y表示风向量windj与第i个航迹点的距离,windlon,lat为当前高度层下第i个航迹点处的风速,根据上述方法分别计算上下高度层对应的风速wind_toplon,lat和wind_bottomlon,lat;
步骤2.3,第i个航迹点的风向量为上下高度风速向量wind_toplon,iat与wind_bottomlon,lat的均值,记作windi=(wind_ui,wind_vi),其中,wind_ui和wind_vi表示其在U方向和V方向上的分量。
5.根据权利要求4所述的基于意图高度的航空器质量估计方法,其特征在于,步骤3中所述的添加缓冲高度标签flag的具体方法为:
比较高度差值Δh与高度阈值ηh的大小,如果Δh<ηh,则判定为异常航迹高度数据,添加缓冲高度标签flag为0,否则判定为正常航迹高度数据,添加缓冲高度标签flag为1。
6.根据权利要求5所述的基于意图高度的航空器质量估计方法,其特征在于,步骤4中殴的根据缓冲时间Δt更新航迹点的数据集合TR,具体方法为:
步骤4.1,设置前向缓冲时间Δtforward,找到航迹数据集中缓冲高度标签为0的航迹点数据tpi,移该除航迹点向前Δtforward个数据,得到新的航迹点的数据集TP′={tpi,i=1,
2,...,N′},且N′≤N;
步骤4.2,设置后向缓冲时间Δtbackward,找到新的航迹数据集TP′中缓冲高度标签为0的航迹点tpi,移除该航迹点向后Δtbackward个数据,得到新的航迹点的数据集TP″={tpi,i=
1,2,...,N″},且N″≤N′≤N。
7.根据权利要求6所述的基于意图高度的航空器质量估计方法,其特征在于,步骤5中所述的计算风对目标航空器的影响,具体过程为:
步骤5.1,在新的航迹点数据集TP″中,分解航迹点数据tpi中的真空速tasi到U方向和V方向,方法如下:
其中,真空速tasi在U方向上的分量记作tas_ui,真空速tasi在V方向上的分量记作tas_vi;
步骤5.2,计算风对目标航空器的影响effecti,方法如下:
effecti=(wind_ui+1‑wind_ui)tas_ui+(wind_vi+1‑wind_vi)tas_vi
其中,wind_ui+1和wind_ui为当前航空器在ti+1和ti时刻的U方向上的风速分量,wind_vi+1和wind_vi为当前航空器在ti+1和ti时刻的V方向上的风速分量。
8.根据权利要求7所述的基于意图高度的航空器质量估计方法,其特征在于,步骤6中所述的获得目标航空器的估计质量massest,具体过程如下:
步骤6.1,根据BADA性能模型获取目标航空器的最小质量massmin与最大质量massmax,在质量区间[massmin×ηmin,massmax×ηmax]内线性插值K个质量数值,其中,ηmin和ηmax分别是质量上限阈值和质量下限阈值,得到质量序列MA={massk,k=1,2,...K},其中,massk表示该序列中的第k个质量;
步骤6.2,遍历质量序列,分别计算质量序列中每个质量对应的功率power与Q值,根据BADA性能模型,计算航迹点数据tpi对应的功率power与Q值,记作poweri(massk)与Qi(massk);
步骤6.3,计算第k个质量massk时全部航迹点的均方误差errork,方法如下:
其中,tpi∈TR″,得到massk对应的均方误差序列ERROR={errork,k=1,2,...,K},其中,errork数值最小对应的质量就是目标航空器的估计质量massest。
9.根据权利要求8所述的基于意图高度的航空器质量估计方法,其特征在于,步骤7中所述的计算不同高度层对应的航空器估计质量massfl,具体方法如下:
步骤7.1,设置航空器的最低飞行高度层fllow与最低飞行高度层flhigh,在高度层区间[fllow,flhigh]内线性插值Z个高度层数值,得到高度层序列FL={flz,z=1,2,...Z},其中,flz表示第z个高度层;
步骤7.2,遍历该高度层序列,比较航迹点数据tpi中的高度hpi与高度层flz,找到最近的高度对应的航迹点索引indexfl;
步骤7.3,计算航迹点索引indexfl前k个质量massk对应的全部航迹点的均方误差errork,方法如下:
得到对应高度层的航空器估计质量massfl,以及所有高度层下的航空器估计质量集合M={massfl,fl∈FL}。
10.根据权利要求9所述的基于意图高度的航空器质量估计方法,其特征在于,步骤8中所述的获取目标航空器的最终估计质量,具体过程如下:
步骤8.1,对航空器估计质量集合M={massfl,fl∈FL}中的每个数值进行合理性校验,包括:去除空数值,删除超出[massmin,massmax]质量区间的数值,得到航空器的有效估计质量集合M′={massfl,fl∈FL′};
步骤8.2,计算有效估计质量集合M′中全部质量的平均质量mass=average{massfl|fl∈FL′},该平均质量mass即为目标航空器的最终估计质量。

说明书全文

一种基于意图高度的航空器质量估计方法

技术领域

[0001] 本发明涉及一种航空器质量估计方法,特别是一种基于意图高度的航空器质量估计方法。

背景技术

[0002] 飞机轨迹预测一直是航空运输中许多机载和地面应用的关键问题。随着基于轨迹运行概念的推广,当前的空中交通管理与控制系统正在向基于轨迹运行的方式转变,高精度地预测航空器轨迹已经成为空中交通管理与控制的关键问题。
[0003] 目前,大多数方法都是基于质点模型与航空动学来进行轨迹预测。这类方法将航空器模拟为一个有质量的点,运用顿第二定律进行航空器的受力分析,通过一组微分方程将飞机的推力、阻力、重力与惯性加速度联系起来,在知道航空器初始状态(质量、当前推力设置、位置、速度、倾等)、大气条件(力、温度)和飞机意图(推力曲线、速度曲线、航线)的情况下,通过积分运算来获取航空器未来的速度和高度。其中,质量参数是影响轨迹预测准确性的关键参数。由于许多航空公司将航空器质量视为重要的竞争参数,不允许机载端下传质量信息,地面系统无法获取准确的航空器起飞质量,因此目前的通用做法是采用BADA提供的参考质量或者通过最小二乘法利用前序航迹点计算估计质量作为航空器质量。BADA提供的参考质量是一个固定的数值,由航空器机型决定,不能动态反映不同航空器之间的质量偏差,而通过最小二乘方法计算得到的估计质量,受前序航迹点影响很大,结果波动明显,无法提供稳定的航空器质量。

发明内容

[0004] 发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种基于意图高度的航空器质量估计方法。
[0005] 为了解决上述技术问题,本发明公开了一种基于意图高度的航空器质量估计方法,包括以下步骤:
[0006] 步骤1,解析目标航空器的S模式雷达数据,得到该航空器在各航迹点的数据集TP={tpi,i=1,2,...,N},其中,tpi为第i个航迹点的数据;
[0007] 解析气象数据,得到气象数据集合WD={windj,j=1,2,...,M},其中,windj表示第j个风向量;
[0008] 步骤2,插值计算航迹点数据tpi对应的UV坐标系中U方向和V方向的风速分量wind_ui和wind_vi;
[0009] 步骤3,设置高度阈值ηh,计算航迹点数据中的气压高度Hp与意图高度MCP的高度差值Δh,比较高度差值Δh与高度阈值ηh的大小,添加缓冲高度标签flag;
[0010] 步骤4,设置缓冲时间Δt,根据缓冲时间更新航迹点的数据集合TR;
[0011] 步骤5,分解真空速TAS到U方向和V方向,计算风对目标航空器的影响;
[0012] 步骤6,设置质量序列,计算不同质量对应的航空器的功率power、Q值以及单位质量下功率power和Q值的均方误差error,比较均方误差error数值大小获得目标航空器的估计质量massest;
[0013] 步骤7,设置高度层序列,筛选各个高度层中的的航迹点数据,计算不同高度层对应的航空器估计质量massfl;
[0014] 步骤8,获取目标航空器的最终估计质量:对所有高度层的航空器估计质量massfl进行数值校验,去除异常质量,得到有效质量序列,计算有效质量的平均质量作为目标航空器的最终估计质量。
[0015] 进一步的,步骤1中所述的第i个航迹点的数据tpi,为一个向量,具体如下:
[0016] tpi=[ti,hpi,mcpi,tasi,θi]
[0017] 其中,ti、hpi、mcpi、tasi和θi分别代表第i个航迹点对应的时间、气压高度、意图高度、真空速和航向角,N为航迹点总数。
[0018] 进一步的,步骤1中所述的第j个风向量windj,具体如下:
[0019] 风向量windj由UV坐标系中的方向U和方向V两个方向的分量组成,记作:
[0020] windj=(wind_uj,wind_vj)
[0021] 其中,每个方向的风向量都是一个基于时间与空间的四维矢量,分别记作:
[0022] wind_uj=(tj,lj,lonj,latj)
[0023]
[0024] 其中,tj、lj、lonj和latj分别表示第j个风向量对应的时间戳、高度层、经度和纬度。
[0025] 进一步的,步骤2中所述的插值计算航迹点数据tpi对应的UV坐标系中U方向和V方向的风速分量wind_ui和wind_vi,具体方法如下:
[0026] 步骤2.1,获取与航迹点数据tpi中的时间ti最近的风向量的时间戳tj对应的时间索引,记作indext;
[0027] 步骤2.2,获取与航迹点数据tpi中的高度hpi最近的风向量的高度层lj对应的索引,其中,最近的上高度层索引记作indexh_top,最近的下高度层索引记作indexh_bottom;
[0028] 步骤2.3,根据时间索引indext与高度层索引indexh_top和indexh_bottom,得到当前高度层下的风向量集合WDl,如下:
[0029] WDl={windt,l,lon,lat|t=indext,l∈[indexh_top,indexh_bottom]}
[0030] 遍历上述集合中的向量,计算高度层风速windlon,lat,计算方法为:
[0031]
[0032] 其中,x和y分别表示当前高度层下的风向量windj的经度和纬度,Dx,y表示风向量windj与第i个航迹点的距离,windlon,lat为当前高度层下第i个航迹点处的风速,根据上述方法分别计算上下高度层对应的风速wind_toplon,lat和wind_bottomlon,lat;
[0033] 步骤2.3,第i个航迹点的风向量为上下高度风速向量wind_toplon,lat与wind_bottomlon,lat的均值,记作windi=(wind_ui,wind_vi),其中,wind_ui和wind_vi表示其在U方向和V方向上的分量。
[0034] 进一步的,步骤3中所述的添加缓冲高度标签flag的具体方法为:
[0035] 比较高度差值Δg与高度阈值ηh的大小,如果Δh<ηh,则判定为异常航迹高度数据,添加缓冲高度标签flag为0,否则判定为正常航迹高度数据,添加缓冲高度标签flag为1。
[0036] 进一步的,步骤4中殴的根据缓冲时间Δt更新航迹点的数据集合TR,具体方法为:
[0037] 步骤4.1,设置前向缓冲时间Δtforward,找到航迹数据集中缓冲高度标签为0的航迹点数据tpi,移该除航迹点向前Δtforward个数据,得到新的航迹点的数据集TP′={tpi,i=1,2,...,N′},且N′≤N;
[0038] 步骤4.2,设置后向缓冲时间Δtbackward,找到新的航迹数据集TP′中缓冲高度标签为0的航迹点tpi,移除该航迹点向后Δtbackward个数据,得到新的航迹点的数据集TP″={tpi,i=1,2,...,N″},且N″≤N′≤N。
[0039] 进一步的,步骤5中所述的计算风对目标航空器的影响,具体过程为:
[0040] 步骤5.1,在新的航迹点数据集TP″中,分解航迹点数据tpi中的真空速tasi到U方向和V方向,方法如下:
[0041]
[0042] 其中,真空速tasi在U方向上的分量记作tas_ui,真空速tasi在V方向上的分量记作tas_vi;
[0043] 步骤5.2,计算风对目标航空器的影响effecti,方法如下:
[0044] effecti=(wind_ui+1‑wind_ui)tas_ui+(wind_vi+1‑wind_vi)tas_vi
[0045] 其中,wind_ui+1和wind_ui为当前航空器在ti+1和ti时刻的U方向上的风速分量,wind_vi+1和wind_vi为当前航空器在ti+1和ti时刻的V方向上的风速分量。
[0046] 进一步的,步骤6中所述的获得目标航空器的估计质量massest,具体过程如下:
[0047] 步骤6.1,根据BADA性能模型获取目标航空器的最小质量massmin与最大质量massmax,在质量区间[massmin×ηmin,massmax×ηmax]内线性插值K个质量数值,其中,ηmin和ηmax分别是质量上限阈值和质量下限阈值,得到质量序列MA={massk,k=1,2,...K},其中,massk表示该序列中的第k个质量;
[0048] 步骤6.2,遍历质量序列,分别计算质量序列中每个质量对应的功率power与Q值,根据BADA性能模型,计算航迹点数据tpi对应的功率power与Q值,记作poweri(massk)与Qi(massk);
[0049] 步骤6.3,计算第k个质量massk时全部航迹点的均方误差errork,方法如下:
[0050]
[0051] 其中,tpi∈TR″,得到massk对应的均方误差序列ERROR={errork,k=1,2,...,K},其中,errork数值最小对应的质量就是目标航空器的估计质量massest。
[0052] 进一步的,步骤7中所述的计算不同高度层对应的航空器估计质量massfl,具体方法如下:
[0053] 步骤7.1,设置航空器的最低飞行高度层fllow与最低飞行高度层flhigh,在高度层区间[fllow,flhigh]内线性插值Z个高度层数值,得到高度层序列FL={flz,z=1,2,...Z},其中,flz表示第z个高度层;
[0054] 步骤7.2,遍历该高度层序列,比较航迹点数据tpi中的高度hpi与高度层flz,找到最近的高度对应的航迹点索引indexfl;
[0055] 步骤7.3,计算航迹点索引indexfl前k个质量massk对应的全部航迹点的均方误差errork,方法如下:
[0056]
[0057] 得到对应高度层的航空器估计质量massfl,以及所有高度层下的航空器估计质量集合M={massfl,fl∈FL}。
[0058] 进一步的,步骤8中所述的获取目标航空器的最终估计质量,具体过程如下:
[0059] 步骤8.1,对航空器估计质量集合M={massfl,fl∈FL}中的每个数值进行合理性校验,包括:去除空数值,删除超出[massmin,massmax]质量区间的数值,得到航空器的有效估计质量集合M′={massfl,fl∈FL′};
[0060] 步骤8.2,计算有效估计质量集合M′中全部质量的平均质量mass=average{massfl|fl∈FL′},该平均质量mass即为目标航空器的最终估计质量。
[0061] 有益效果:
[0062] 本发明充分考虑风速对航空器质量估计结果的影响,根据意图高度与高度缓冲时间对航迹点过滤,可以剔除影响质量结果的异常航迹点,通过实施多高度层的航空器质量估计,获得更加稳定的航空器估计质量。附图说明
[0063] 下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
[0064] 图1为本发明一种基于意图高度的航空器质量估计方法的流程图
[0065] 图2为本发明实施实例的航空器估计质量的误差曲线示意图。
[0066] 图3为本发明实施实例的多高度层航空器估计质量的误差曲线示意图。

具体实施方式

[0067] 结合图1,本发明的一种基于意图高度的航空器质量估计方法,包括以下步骤:
[0068] 步骤1:解析航空器S模式雷达数据,得到各航迹点的数据集TP={tpi,i=1,2,...,N},记第i个航迹点为tpi,向量tpi=[ti,hpi,mcpi,tasi,θi],ti、hpi、mcpi、tasi、θi分别代表当前航迹点对应的时间、高度、意图高度、真空速和航向角,N为航迹点总数;解析气象数据,得到气象数据集合WD={windj,j=1,2,...,M},其中,windj表示第j个风向量,M为气象数据总数,风向量windj由正东方向U和正北方向V两个方向的分量组成,记作windj=(wind_uj,wind_vj),每个方向的风向量都是一个基于时间与空间的四维矢量,分别记作wind_uj=(tj,lj,lonj,latj),wind_vj=(tj,lj,lonj,latj),tj、lj、lonj、latj分别表示第j个风向量对应的时间戳、高度层、经度与纬度;
[0069] 步骤2:插值计算航迹点tpi对应的U、V方向风速分量wind_ui、wind_vi,插值方法为:
[0070] 步骤2.1,获取与当前航迹点tpi的时间ti最近的风向量的时间戳tj对应的索引,记作indext;
[0071] 步骤2.2,获取与当前航迹点tpi的高度hpi最近的风向量的高度层lj对应的索引,距离最近的上高度层索引记作indexh_top,距离最近的下高度层索引记作indexh_bottom;
[0072] 步骤2.3,根据时间索引indext与高度层索引indexh_top、indexh_bottom,得到当前高度层下的风向量集合WDl={windt,l,lon,lat|t=indext,l∈[indexh_top,indexh_bottom]},遍历集合中各向量计算高度层风速,计算方法为:
[0073]
[0074] 其中,x、y分别表示当前高度层下的风向量windj的经度、纬度,Dx,y表示风向量windj与航迹点tpi的距离,windlon,lat为当前高度层下航迹点tpi处的风速,分别计算上下高度层对应的风速wind_toplon,lat、wind_bottomlon,lat;
[0075] 步骤2.3,当前航迹点tpi的风向量为上下高度风速向量wind_toplon,lat与wind_bottomlon,lat的均值,记作windi=(wind_ui,wind_vi)。
[0076] 步骤3:设置合适的高度阈值ηh,计算气压高度Hp与意图高度MCP的差值Δh,比较高度差值Δh与高度阈值ηh的大小,添加缓冲高度标签flag,添加缓冲高度标签的方法为:
[0077] 步骤3.1,比较高度差值Δh与高度阈值ηh的大小,如果Δh<ηh,为异常航迹高度数据,添加缓冲高度标签0;
[0078] 步骤3.2,比较高度差值Δh与高度阈值ηh的大小,如果Δh≥ηh,为正常航迹高度数据,添加缓冲高度标签1。
[0079] 步骤4:设置缓冲时间Δt,根据缓冲时间更新航迹点的数据集合TR,更新航迹点的数据集合TR的方法为:
[0080] 步骤4.1,设置前向缓冲时间Δtforward,找到航迹数据集中缓冲高度标签为0的航迹点tpi,移除航迹点tpi向前Δtforward个数据,得到新的航迹点的数据集TP′={tpi,i=1,2,...,N′},且有N′≤N;
[0081] 步骤4.2,设置后向缓冲时间Δtbackward,找到航迹数据集中缓冲高度标签为0的航迹点tpi,移除航迹点tpi向后Δtbackward个数据,得到新的航迹点的数据集TP″={tpi,i=1,2,...,N″},且有N″≤N′≤N。
[0082] 步骤5:分解真空速TAS到U、V方向,计算风对航空器的影响,计算风对航空器的影响的过程为:
[0083] 步骤5.1,分解航迹点tpi的真空速TAS到U、V方向,计算方法如下:
[0084]
[0085] 其中,真空速tasi在U方向上的分量记作tas_ui,真空速tasi在V方向上的分量记作tas_vi,tasi航迹点tpi的真空速,θi航迹点tpi的航向角;
[0086] 步骤5.2,计算风对航空器的影响,方法如下:
[0087] effecti=(wind_ui+1‑wind_ui)tas_ui+(wind_vi+1‑wind_vi)tas_vi
[0088] 其中,wind_ui+1、wind_ui为航空器在ti+1、ti时刻的U方向上风速,wind_vi+1、wind_vi为航空器在ti+1、ti时刻的V方向上风速。
[0089] 步骤6:设置质量序列,计算不同质量对应的航空器的功率power、Q值以及单位质量下power和Q值的均方误差error,比较error数值大小获得航空器估计质量mass,计算航空器质量的具体过程如下:
[0090] 步骤6.1,根据飞机性能模型获取航空器的最小质量massmin与最大质量massmax,在质量区间[massmin×ηmin,massmax×ηmax]内线性插值K个质量数值,ηmin与ηmax分别是质量上限阈值和质量下限阈值,得到质量序列MA={massk,k=1,2,...K};
[0091] 其中,在某些实施例中,上述飞机性能模型可以采用BADA模型或其他模型实现,BADA(Base of aircraft data)模型是由欧控组织(EUROCONTROL)发布的一套飞机性能模型及相关数据库
[0092] 步骤6.2,遍历质量序列,分别计算每个massk对应的功率power与Q值,根据BADA性能模型,计算航迹点tpi对应的功率power与Q值,记作poweri(massk)与Qi(massk);
[0093] 步骤6.3,计算massk下全部航迹点的均方误差,方法如下:
[0094]
[0095] 其中,tpi∈TR″,得到massk对应的均方误差序列ERROR={errork,k=1,2,...,K},errork数值最小对应的质量就是航空器的估计质量massest。
[0096] 步骤7:设置高度层序列,筛选各个高度层的航迹点数据,计算不同高度层对应的航空器估计质量massfl,计算不同高度层下航空器估计质量massfl的方法如下:
[0097] 步骤7.1,设置航空器的最低飞行高度层fllow与最低飞行高度层flhigh,在高度层区间[fllow,flhigh]内线性插值Z个高度层数值,得到高度层序列Fl={flz,z=1,2,...Z};
[0098] 步骤7.2,遍历高度层序列,比较航迹点tpi的高度hpi与高度层flz,找到最近的高度对应的航迹点索引indexfl;
[0099] 步骤7.3,计算索引indexfl前massk对应的全部航迹点的均方误差,方法如下:
[0100]
[0101] 得到对应高度层的航空器估计质量massfl,以及所有高度层下的航空器估计质量集合M={massfl,fl∈FL}。
[0102] 步骤8:对所有高度层的航空器估计质量massfl进行数值校验,去除异常质量,得到有效质量序列,计算有效质量的平均质量作为航空器的最终估计质量,具体过程如下:
[0103] 步骤8.1,对航空器估计质量集合M={massfl,fl∈FL}中的每个数值进行合理性校验,去除空数值,以及超出[massmin,massmax]质量区间的数值,得到航空器的有效估计质量集合M′={massfl,fl∈Fl′};
[0104] 步骤8.2,计算质量集合M′中全部质量的平均质量mass=average{massfl|fl∈FL′},那么平均质量mass即为航空器的最终估计质量。
[0105] 下面结合图2至图3,通过仿真实验的实施例及其效果评价来进一步说明本发明。
[0106] 在本实施例中,如图2和图3中的曲线分别是航空器估计质量的误差曲线以及多高度层航空器估计质量的误差曲线。实验目标就是通过基于意图高度的航空器质量估计方法计算稳定的航空器估计质量。图2对应的质量的误差曲线中的三角形标志表示当选取三角形位置对应横坐标处质量时,该航空器的误差最小,因此三角形位置处的质量就是航空器估计质量。图3描述的是多高度层航空器估计质量的误差曲线,各曲线上三角形位置对应横坐标处质量就是各高度层下航空器的估计质量,图3表明选择不同高度层区间的航迹点,得到的航空器的估计质量也不同,图3中方框圈出的质量就是一个超出质量上下限的异常质量,去除该异常数值后,计算其余高度层下的估计质量的均值,就得到稳定的航空器估计质量。该结果表明本发明方法可以剔除影响质量结果的异常航迹点,通过实施多高度层的航空器质量估计,获得更加稳定的航空器估计质量。
[0107] 具体实现中,本申请提供计算机存储介质以及对应的数据处理单元,其中,该计算机存储介质能够存储计算机程序,所述计算机程序通过数据处理单元执行时可运行本发明提供的一种基于意图高度的航空器质量估计方法的发明内容以及各实施例中的部分或全部步骤。所述的存储介质可为磁碟、光盘、只读存储记忆体(read‑only memory,ROM)或随机存储记忆体(random access memory,RAM)等。
[0108] 本领域的技术人员可以清楚地了解到本发明实施例中的技术方案可借助计算机程序以及其对应的通用硬件平台的方式来实现。基于这样的理解,本发明实施例中的技术方案本质上或者说对现有技术做出贡献的部分可以以计算机程序即软件产品的形式体现出来,该计算机程序软件产品可以存储在存储介质中,包括若干指令用以使得一台包含数据处理单元的设备(可以是个人计算机,服务器单片机,MCU或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
[0109] 本发明提供了一种基于意图高度的航空器质量估计方法的思路及方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
QQ群二维码
意见反馈