中国电信新建网站备案管理系统 录完信息,苏州市网站制作,发表文章的平台有哪些,飞机查询网站开发的创新点2023年第十三届MathorCup高校数学建模挑战赛
D题 航空安全风险分析和飞行技术评估问题
原题再现 飞行安全是民航运输业赖以生存和发展的基础。随着我国民航业的快速发展#xff0c;针对飞行安全问题的研究显得越来越重要。2022 年 3 月 21 日#xff0c;“3.21”空难的发生…2023年第十三届MathorCup高校数学建模挑战赛
D题 航空安全风险分析和飞行技术评估问题
原题再现 飞行安全是民航运输业赖以生存和发展的基础。随着我国民航业的快速发展针对飞行安全问题的研究显得越来越重要。2022 年 3 月 21 日“3.21”空难的发生终结了中国民航安全飞行 1 亿零 59 万飞行小时的历史最好安全记录。严重飞行事故的发生不仅会给航空公司带来巨大的经济损失更会对乘客造成极大的生命威胁。因而需要聚焦飞行安全问题强化航空安全研究综合利用现有数据强化科学管理通过有针对性、系统性的管控手段有效提升从业人员的素质监测和预警风险进而降低飞行事故的发生几率。 航空安全大数据主要包括快速存取记录器(Quick Access RecorderQAR)数据该数据主要记录飞机在飞行过程中的各项飞行参数在飞行品质监控Flight Operational Quality AssuranceFOQA中QAR 中超出人为设定限制值的数据记为超限数据。除此之外在实际研究过程中还会涉及到飞行中的舱音数据等。本问题主要涉及的是 QAR 数据QAR 数据相对比较规范。 在飞行品质监控具体研究和应用方面目前我国民航业内的研究主要分为两个方面一是针对超限事件的研究、分析和应用二是对非超限数据的统计分析和应用。对于超限事件的研究一般是通过规定飞行参数的集中区域设置超限阈值将超出阈值部分的飞行记录找出来进行重点分析防范潜在隐患造成严重飞行事故。目前此类分析是飞行品质监控工作的主体较好地保证了现阶段的安全工作其不足之处在于缺少对超限原因的分析。由于超限并非全部是人为因素引发例如许多是由于特殊环境条件造成的甚至有可能是飞机本身的设计、制造因素所致因此仅通过单纯的超限分析很难识别出来如果仅基于超限事件对飞行机组进行管理很容易误入歧途。QAR 超限可用于航空安全管理和飞行训练的数据支持。目前并不倾向于仅以少量的 QAR 超限数据为依据开展飞行训练工作因此飞行品质监控工作逐渐衍生出另外一种倾向性即通过挖掘 QAR 全航段数据开展分析形成特定人员的飞行品质记录。基于不同飞行机组、飞行航线、机场、特定飞行条件下的飞行记录通过对数据进行建模、分析计算评估风险倾向性开展有针对性的安全管理排查安全隐患改进安全绩效。目前类似研究主要是大规模读取飞行数据并进行存储和分析形成飞行品质服务平台为风险评估和趋势分析提供数据基础。G 值是飞机飞行过程中过载情况的直接反应在着陆安全分析中G 值通常是描述落地瞬间安全性的重要指标。着陆瞬间 G 值指的是飞机接地瞬间前 2 秒和后5 秒数据的最大 G 值。 基于以上背景请你们团队解决以下问题 问题 1有些 QAR 数据存在错误需要对数据进行预处理去伪存真以减少错误数据对研究分析带来的影响。请你们的队伍对附件 1 的数据质量开展可靠性研究提取与飞行安全相关的部分关键数据项并对其重要程度进行分析。 问题 2飞机在从起飞到着陆的整个飞行过程中通过一系列的飞行操纵确保飞行安全这些操纵主要包括横滚操纵、俯仰操纵等。目前国内航空公司通过超限监控飞行操纵动作这种监控方法虽然能够快速分辨出飞机的状态偏差但是只能告诉安全管理人员发生了什么而不能立刻得出发生这种偏差的原因。为此可以通过操纵杆的过程变化情况来分析产生这种偏差的原因。根据附件 1请你们对飞行操纵进行合理量化描述。下图为 3 次着陆过程中的杆位变化曲线其中红色曲线描述了一次重着陆着陆 G 值超过给定限制值过程该重着陆主要是由于飞行机组在低空有一次不当松杆操纵所致红色曲线中的接地前 5 秒有一个明显下凸这就是需要进行量化描述的一次松杆操纵。 问题 3导致不同超限发生的原因各不相同有时是特定机场容易出现特定的超限有时是特定的天气容易出现特定的超限有时是特定的飞行员容易出现特定的超限。请研究附件 2 的数据对超限的不同情况进行分析研究不同超限的基本特征如分析飞机在哪些航线或者在哪些机场容易出现何种超限等。 问题 4飞机运行数据的研究一般分为两大类一类是通过航线运行安全检查Line Operations Safety AuditLOSA获取的飞行员的运行表现另外一类是根据相关学者建议基于飞行参数开展飞行技术评估。根据附件 3请你们建立数学模型探讨一种基于飞行参数的飞行技术评估方法分析飞行员的飞行技术数据表中的“不同资质”代表飞行员的不同技术级别。 问题 5随着技术的进步未来在民航客机上安装实时传输的 QAR 数据记录系统已成为可能这种“实时飞行数据”技术可以在接近实时的情况下把航班飞行数据传输到地面分析系统极大地提高风险识别能力和预防水平。假设飞行数据已能实现陆空实时传输如果你是该航空公司的安全管理人员请建立航空公司实时自动化预警机制预防可能的安全事故发生结合附件 1 的数据给出仿真结果。
整体求解过程概述(摘要) 本题主要是开发一个基于航空安全风险分析及飞行技术评估的预警模型等问题展开深入研究。通过熵权法、随机森林等数学模型使用 Python、Matlab、Spsspro 等软件对操作进行量化描述超限不同情况进行统计分析最后基于参数评估飞行技术建立一个自动化预警模型。从而降低我国航空事故的发生。 对于问题一首先使用拉伊达法则和箱线图对附件 1 数据进行预处理再用 Python对于附件 1 中异常的 QAR 数据进行去除之后用 Python 可视化[1]对数据进行可靠性研究最后用熵权法得出来每个影响因素的具体权重其中影响最大的因素是飞机重量其余前五项关键数据项分别是风速、无线电高度、风向、右发油门杆位置、左发油门杆位置。 对于问题二首先用线性回归分析 F 值验证结果后可以显着拒绝总回归系数为0 的原假设P0.05。表明之间存在着线性关系并且通过 R²值分析模型拟合情况同时对 VIF 值进行分析模型出现共线性后 X 呈现显著性结合 B 值探究 X、Y 之间关系并确定模型公式。最终得出结论CAP CLM 1 POSN 和 CAP WHL 1 POSN 的和作为表征量在其和较大时说明操作不平稳反之平稳。 对于问题三首先对于附件 2 中数据进行超限分析得到不同情况下的超限类型下滑道偏差过高或过低、俯仰角率过高或过低等随后利用频数分析算法和智能分析对于超限数据进行频数分析使用 spsspro 对于飞行超限频数所占比例进行绘图。得出频数分析图表最终可以得出飞机最容易出现的超限类型为50 英尺至接地距离远在飞行阶段为 6 时目的机场为 68 时最容易出现超限。 对于问题四首先对附件 3 数据进行预处理去除无关数据以及空白数据等。随后基于飞行参数评估飞行技术建立随机森林、逻辑回归、LGBM 等数学模型运用机器学习算法来衡量训练集准确度和测试集准确度以及平均绝对误差进行模型评估从两种方法来评估飞行员的飞行技术训练分数、交叉验证。从模型最终结果来看LGBM 模型对于数据的预测精度良好。 对于问题五首先进行挖掘学习行为特征的时间序列设计了一种基于时间记忆增强模型的预警机制挖掘和预测决策和干预机制。随后选取三种临近最优算法进行对比实验经过充分的数据训练和测试计算了图中三个典型数据子集的精度。随着预测戳数的增加ILA GAN 的优势更加明显。
模型假设 结合本题实际为确保模型求解准确性和合理性提出以下几点假设 1. 假设本题中给出所有附件中的数据都真是可信。 2. 假设本题考虑到因素外其他未考虑的因素都不影响本题的计算。 3. 假设本题中的样本互不影响、互不关联、相互独立各指标之间互不干扰。
问题分析 针对问题一的分析 需要我们对 QAR 数据预处理去伪取真减少错误数据对研究的影响和误差根据附件 1 数据进行可靠性研究并且运行数据分析。初步数据处理是指在分类或分组之前对收集到的数据进行必要的审查、筛选和整理。本题首先我们利用拉依达和箱线图进行数据预处理去除已知 QAR 数据里的错误数据进行更深程度的分析使用熵权法和 python 的可视化[2]计算出权重对各指标的权重进行分析得到权重可视化图形最终得出结论权重值所在越高重要程度越高反之越低。 针对问题二的分析 在整个飞行过程中从起飞到着陆飞行安全都受到飞行控制主要包括横滚操纵、俯仰操纵等的影响航空公司只能监测飞行操纵动作从而告诉管理者发生了什么却不能立刻得出这种偏差的原因无法立即确定。因此我们必须根据附件 1 对飞行控制进行合理描述通过操纵杆的变化来分析偏差的原因。本题我们采用线性回归最小二乘法进行分析 F 的值分析 F 值能否明显的拒绝总体回归系数为 0 的原始假设P0.05如呈现显著性则表示 F 值与其存在线性关系并且经过 R²值分析模型拟合状况并且对VIF 的值开展分析VIF 大于 5 或者 10如若模型呈现出共线性则使用逐步回归法。开展 X 显著性的分析P0.05若呈现显著性则使用其分析 X 对 Y 的影响关系。结合回归因子 B比较分析 X 对 Y 的影响。确定模型模式在应用线性回归之前可以通过正态性检验或残差数据处理等统计方法对数据进行验证和纯化。 针对问题三的分析 本题要求我们根据附件 2对于不同的超限情况进行分析分析不同超限时的特征。首先我们对附件 2 数据进行超限分析得出不同超限情况的不同超限类型之后利用频数分析算法对不同机场和不同超限类型进行频数分析得到频数分析图和分析表最终得出不同地点和时间出现的不同超限特征。 针对问题四的分析 本题要求我们基于飞行参数进行飞行评估并根据附件 3 建立数学模型使用飞行技术评估方法分析飞行员的飞行技术水平。并且我们基于对给出的附件 3 数据进行预处理去除了较多无关数据以及附件 3 中空缺数据。随之我们建立了随机森林、逻辑回归、LGBM 回归等数学模型运用机器学习算法来衡量训练集准确度和测试集准确度以及平均绝对误差进行模型评估从飞行员训练分数和交叉验证的评估方法来分析飞行员的飞行技术。 针对问题五的分析 对于问题五的分析我们根据问题一的飞行安全相关的部分关键数据项问题二的重着陆分析问题三的不同超限情况的分析问题四的飞行技术评估设计了一种基于时间记忆增强模型的预警机制。通过大量的数据训练和测试该方法对学习行为特征的分析有效具有较强的有效性和可靠性。其次挖掘和预测决策和干预机制。本研究分析了三种密度的特征计算特征之间的相关性将强相关系数≥0.85的特征聚类。当一个聚类中的特征数量超过总特征的60%时该数据子集被标记为典型否则被定义为非典型。选取三种临近最优算法进行对比实验分别是:(1)将原始行为与递归神经网络(OB RNN)相结合(2)基于总频率行为和递归神经网络(TF RNN)的预警算法(3)基于LSTM-Autoencoder和注意机制(LA AM)的预警算法。
模型的建立与求解整体论文缩略图 全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可
程序代码(代码和文档not free)
import numpy
import pandas
from spsspro.algorithm import statistical_model_analysis
#生成案例数据
data_x1 pandas.DataFrame({
A: numpy.random.random(size100),
B: numpy.random.random(size100)
})
data_x2 pandas.DataFrame({C: numpy.random.choice([1, 2, 3],
size100)})
data_y pandas.Series(datanumpy.random.choice([1, 2], size100), nam
eY)
#线性回归输入参数详细可以光标放置函数括号内按 shifttab 查看输出结果参考
spsspro 模板分析报告
result statistical_model_analysis.linear_regression(data_ydata_y, dat
a_x1data_x1, data_x2data_x2)
print(result)import pandas as pd
datapd.read_excel(附件/附件 3 飞行参数测量数据/附件 3飞行参数测量数据.xlsx)
data
for i in data.columns:
if data[i].isnull().sum()0.15*2369:
print(i,data[i].isnull().sum())
del data[i]
data.shape
data.dropna().shape
data.dropna(inplaceTrue)
tt2[]
for i in data.columns:
if data[i].dtypes not in [int64,float64]:
print(i,data[i].dtypes)
else:
tt2.append(i)
tt[机型, 落地主操控, V2_Method, Vref_Method, RoD_Method, MACH_Method, Turnoff_GS]
newdatapd.concat([data[tt2],pd.get_dummies(data[tt])],axis1)
newdata
机器学习代码
import pandas as pd
from numpy import *
import numpy as np
from sklearn.neural_network import MLPClassifier
from sklearn.tree import DecisionTreeClassifier
from sklearn.ensemble import RandomForestClassifier
from sklearn.linear_model import LogisticRegression
import xgboost as xgb
import lightgbm as lgb
import numpy as np
import matplotlib.pyplot as plt
from sklearn.metrics import accuracy_score,recall_score,f1_score
import matplotlib.pyplot as plt
from sklearn.metrics import mean_absolute_error
import warnings
warnings.filterwarnings(ignore)
from sklearn.model_selection import train_test_split
def mape(actual, pred):
actual, pred np.array(actual), np.array(pred)
return np.mean(np.abs((actual - pred) / actual))
XnewX
Ydata[落地主操控人员资质]
tr_x,te_x,tr_y,te_ytrain_test_split(X ,Y,test_size0.1,random_state
5)
modelMLPClassifier(hidden_layer_sizes10,max_iter1000).fit(tr_x,tr
_y)
y_predmodel.predict(te_x)
print(\n 逻辑回归:)
logreg LogisticRegression(C1,solverliblinear,multi_class auto)
logreg.fit(tr_x, tr_y)
score logreg.score(tr_x, tr_y)
score2 logreg.score(te_x, te_y)
print(训练集准确度:{:.3f}.format(logreg.score(tr_x,tr_y)))
print(测试集准确度:{:.3f}.format(logreg.score(te_x,te_y)))
y_pred logreg.predict(te_x)
print(平均绝对误差:,mean_absolute_error(te_y, y_pred))
print(ACC,accuracy_score(te_y,y_pred))
print(REC,recall_score(te_y,y_pred,averagemicro))
print(F-score,f1_score(te_y,y_pred,averagemicro))
print(\n 决策树:)
treeDecisionTreeClassifier(max_depth50,random_state0)
tree.fit(tr_x,tr_y)
y_pred tree.predict(te_x)
print(训练集准确度:{:.3f}.format(tree.score(tr_x,tr_y)))
print(测试集准确度:{:.3f}.format(tree.score(te_x,te_y)))
print(平均绝对误差:,mean_absolute_error(te_y, y_pred))
print(ACC,accuracy_score(te_y,y_pred))
print(REC,recall_score(te_y,y_pred,averagemicro))
print(F-score,f1_score(te_y,y_pred,averagemicro))
print(\n 随机森林:)
rfRandomForestClassifier(max_depth20,n_estimators1000,random_state0)
rf.fit(tr_x,tr_y)
print(训练集准确度:{:.3f}.format(rf.score(tr_x,tr_y)))
print(测试集准确度:{:.3f}.format(rf.score(te_x,te_y)))
y_pred rf.predict(te_x)
print(平均绝对误差:,mean_absolute_error(te_y, y_pred))
print(ACC,accuracy_score(te_y,y_pred))
print(REC,recall_score(te_y,y_pred,averagemicro))
print(F-score,f1_score(te_y,y_pred,averagemicro))
print(\nLGBM:)
lgb_modellgb.LGBMClassifier()
lgb_model.fit(tr_x,tr_y)
print(训练集准确度:{:.3f}.format(lgb_model.score(tr_x,tr_y)))
print(测试集准确度:{:.3f}.format(lgb_model.score(te_x,te_y)))
y_pred lgb_model.predict(te_x)
print(平均绝对误差:,mean_absolute_error(te_y, y_pred))
print(ACC,accuracy_score(te_y,y_pred))
print(REC,recall_score(te_y,y_pred,averagemicro))
print(F-score,f1_score(te_y,y_pred,averagemicro))
# 学习曲线函数
def plot_learning_curve(estimator, title, X, y, ylimNone, cvNone, n_jobs1, train_sizesnp.linspace(.1, 1.0, 5)):
plt.figure()
plt.title(title)
if ylim is not None:
plt.ylim(*ylim)
plt.xlabel(game num)
plt.ylabel(score)
train_sizes, train_scores, test_scores learning_curve(
estimator, X, y, cvcv, n_jobsn_jobs, train_sizestrain_sizes)
train_scores_mean np.mean(train_scores, axis1)
train_scores_std np.std(train_scores, axis1)
test_scores_mean np.mean(test_scores, axis1)
test_scores_std np.std(test_scores, axis1)
plt.grid()
plt.fill_between(train_sizes, train_scores_mean - train_scores_std, train_scores_mean train_scores_std, alpha0.1, colorr)
plt.fill_between(train_sizes, test_scores_mean - test_scores_std, test_scores_mean test_scores_std, alpha0.1, colorg)
plt.plot(train_sizes, train_scores_mean, o-, colorr, labelTraining score)
plt.plot(train_sizes, test_scores_mean, o-, colorg, labelCross-validation score)
plt.savefig(./Q4/%s.jpg%title)
plt.legend(locbest)
plt.show()
return plt
cv ShuffleSplit(n_splits5, test_size0.2, random_state0)
plot_learning_curve(logreg, 逻辑回归, tr_x, tr_y, ylimNone, cvcv, n_jobs1)
plot_learning_curve(tree, 决策树, tr_x, tr_y, ylimNone, cvNone, n_jobs1)
plot_learning_curve(rf, 随机森林, tr_x, tr_y, ylimNone, cvNone, n_jobs1)
plot_learning_curve(lgb_model, LGBM, tr_x, tr_y, ylimNone, cvNone, n_jobs1)想要模型和程序的请见下方“ 只会建模 QQ名片” 点击QQ名片即可
