基于PSO-SVR的飞行员工作负荷预测

杨琪，黄磊，陆中，张子文，韩冰; YANG Qi; HUANG Lei; LU Zhong; ZHANG Ziwen; HAN Bing

2025年6月22日 5:15 星期日

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

基于PSO‑SVR的飞行员工作负荷预测 PDF

- ORCID：
杨琪 ¹
- ORCID：
黄磊 ²
- ORCID：
陆中 ¹
✉
- ORCID：
张子文 ¹
- ORCID：
韩冰 ¹

1. 南京航空航天大学民航学院，南京 211106； 2. 中国直升机设计研究所总体气动室，天津 300450

中图分类号： V323

最近更新：2021-12-21

DOI：10.16356/j.1005-2615.2021.06.014

摘要

飞行员工作负荷是影响飞机运行安全的重要因素。开展飞行员工作负荷预测是适航审定过程中，验证驾驶舱设计是否符合适航规章的重要手段。本文针对某型民用飞机设计了模拟飞行试验，用于采集飞行员生理指标数据和国家航空航天局任务负荷指数（National Aeronautics and Space Administration task lood index， NASA‑TLX）量表评价数据。以飞行员生理指标数据为输入，NASA‑TLX量表主观评价数据为输出，建立了基于粒子群算法优化的支持向量回归机（Particle swarm optimization‑support vector regression， PSO‑SVR）模型的飞行员工作负荷预测模型。对本文建立的PSO‑SVR模型与默认参数的支持向量回归机（Support vector regression， SVR）模型的预测精度进行了对比，针对4个不同场景，预测精度分别提高了7.5%、9.5%、7%和5.8%，结果表明基于PSO‑SVR的预测模型得到的飞行员工作负荷预测值精度更高。

关键词

飞行安全; 国家航空航天局任务负荷指数; 粒子群优化; 支持向量回归机; 工作负荷预测

飞行员工作负荷的大小与飞行安全紧密相关。同一时间段内的飞行事故，由飞行员差错导致的致命事故占比58%，机械故障占比17%^［

1］。飞行员出现差错的原因为执行飞行任务时工作负荷强度大，应激水平高。因此，飞行员执行飞行任务时的工作负荷已成为民用航空器适航审定中的一个重要指标。通过测量飞行员在执行任务时的工作负荷，可以分析驾驶舱的设计是否合理，是否满足适航规章的要求。对飞行员工作负荷进行预测分析是针对25.1302“飞行机组使用的安装系统和设备”、25.1523“最小飞行机组”等适航条款开展符合性验证的重要内容^{［参考文献 2

百度学术}2］，是保证飞机设计满足安全性要求的重要手段。

飞行员工作负荷评价是对飞行员进行工作负荷预测的基础，针对飞行员工作负荷评价国内外已有学者做过一定研究。飞行员工作负荷评价方法主要分为基于生理指标参数的客观评价和基于量表的主观评价。针对工作负荷客观评价，Wilson等^［

3］在真实飞行中，采用10名飞行员在目视飞行仪表飞行状态下，发现视觉要求高的飞行场景下（工作负荷大），飞行员眨眼率下降；Veltman^{［参考文献 4

百度学术}4］验证了多种生理指标（心率、心率变异性、呼吸、血压和眨眼）在实际飞行中的可用性；Wilson等^{［参考文献 3

百度学术}3］在采用多生理指标研究飞行员脑力负荷时发现，脑电波（Electroencephalogram， EEG）数据中的α波和δ波随脑力负荷的变化而变化；针对工作负荷主观评价，常用方法主要有美国国家航空航天局任务负荷指数（National Aeronautics and Space Administration task load index， NASA‑TLX）量表、库珀‑哈珀（Cooper‑Harper， CH）量表、主观工作负荷评价技术（Subjective workload assessment techniques， SWAT）量表等^{［参考文献 5

百度学术}5］；辛甜等^{［参考文献 6

百度学术}6］采用NASA‑TLX量表对飞行员工作负荷进行分析，结果表明NASA‑TLX方法在细化工作负荷来源具有一定优越性；李琳君等^{［参考文献 7

百度学术}7］利用SWAT量表和贝德福德（Bedford）量表对直升机飞行员工作负荷进行评价。

建立飞行员生理指标与主观评价值之间的映射关系是对工作负荷进行准确评价的关键，人工神经网络广泛应用于该领域。刘树强建立了基于误差反向传播（Back propagation， BP）神经网络的机组工作量评估模型，并采用了遗传算法优化了BP模型的权值和阈值^［

8］；李金波等以飞行员心理努力、任务主观难度、注视时间、注视次数、主任务正确率和主任务反应时6项指标为输入，以认知负荷的3个等级为输出，建立了基于BP神经络和自组织神经网络（Self‑organization‑map， SOM）的认知负荷综合评估模型^{［参考文献 9

百度学术}9］。Mazaeva等对30个EEG指标通过主成分降维后，采用自组织神经网络预测了信息加工要求不同的任务负荷的变化^{［参考文献 10

百度学术}10］。

从上述文献可以看出，目前对飞行员工作负荷评价的研究主要分为基于生理指标和基于主观量表两方面。而生理指标只能间接推断工作负荷大小，无法求得生理指标与工作负荷之间的函数关系，主观量表的评价结果受个体差异性的影响较大。在建立生理评价和主观评价两者之间关系的研究中，目前的研究主要集中在通过神经网络建立两者之间的关系，而在神经网络建模时，容易陷入局部最优，并且在网格训练中往往需要大量数据作为训练样本，并且在预测模型精度优化方面的研究较少。

针对以上所述各类方法的不足，本文提出了一种以粒子群算法优化的支持向量回归机（Particle swarm optimization‑support vector regression， PSO‑SVR）为核心的飞行员工作负荷评价模型，综合考虑飞行员的客观生理指标数据和主观评价数据，可对飞行员工作负荷进行准确评价及预测。通过PSO优化SVR的核函数参数和惩罚因子系数。将优化后的参数代入SVR模型中，以飞行员生理指标数据作为网络输入，主观评价数据作为网络输出，建立生理指标与主观评价值之间的非线性映射关系。通过主观评价实际值与预测值之间的对比，验证了所建模型的准确性。

本文首次将支持向量回归模型应用在飞行员工作负荷评估上，为工作负荷评估提出一种新方法。同时为了进一步提高工作负荷评估的精度，引入了粒子群优化算法。基于本文所提方法，在后续需要对飞行员进行工作负荷评价时，只需获取其客观生理数据，有效缩短工作负荷评估周期，同时无需飞行员进行主观评价，避免了由于飞行员主观评价带来的不确定性，提高工作负荷评价的准确性。

1 模拟飞行试验设计

基于某型运输类飞机模拟驾驶舱，构建模拟飞行试验。基于模拟飞行试验，采集飞行员生理指标数据和主观评价数据。

1.1　模拟飞行人员选取

选取民航学院45名男飞行员参与模拟飞行试验，年龄在20~25岁之间，视力正常、身心健康，飞行技术过硬。

1.2　模拟飞行阶段确定

参考某型运输类飞机飞行手册，对整机运行所有阶段进行模拟，各阶段划分如图1所示。

图1 模拟飞行阶段

Fig.1 Simulated flight phase

各阶段描述如下。

滑出：飞机从停放位置移动到跑道；

起飞滑跑：释放刹车到离地之前；

起飞：离地35 英尺；

起飞爬升：离地35~1 500 英尺；

爬升：离地1 500~20 000 英尺；

巡航：飞机水平，离地高度大于20 000 英尺；

下降：离地大于20 000 英尺的任何高度到1 500 英尺；

等待：离地1 500~20 000 英尺的任何高度；

进近：离地1 500~100 英尺；

着陆：离地100 英尺到接触地面；

着陆滑跑：接触地面到速度达到滑行速度；

滑入：跑道至机场航站楼。

各阶段的持续时间分别为：滑出5 min，起飞滑跑26.9 s，起飞5.5 s，起飞爬升1 min，爬升6 min，巡航45 min，下降12 min，等待根据实际情况，进近3.5 min，着陆10 s，着陆滑跑20 s，滑入5 min。

1.3　模拟飞行状态选取

为了增加试验数据量，辅助SVR建模分析，通过改变飞机飞行状态，增加飞行员的工作负荷，设计正常模拟飞行状态和非正常模拟状态，非正常模拟状态包括主飞行显示器（Primary flight display， PFD）显示失效、单发失效、PFD显示和单发同时失效。

1.4　模拟飞行场景设计

模拟飞行场景由模拟飞行阶段和模拟飞行状态组合而成，设计如表1所示的模拟飞行场景。

表1 模拟飞行场景

Table 1 Simulated flight scene

场景编号	飞行阶段	模拟飞行状态
1	所有阶段	正常飞行
2	巡航	PFD显示失效
3	巡航	单发失效
4	巡航	PFD、单发同时失效

1.5　试验数据采集

本飞行试验主要采集飞行员的生理指标数据以及主观评价数据，生理指标数据与主观评价数据均为多次测量后的平均值。

飞行员在进行信息处理（工作负荷发生变化）时，中枢神经系统会活动，与之相关的生理指标也会有所变化，因此可以通过测量生理指标的变化进而评价工作负荷^［

11］。生理指标评价方法中所使用的指标按照涉及的生理器官功能分为3大类：与心脏活动相关的生理指标、与眼睛活动相关的指标以及与大脑活动相关的生理指标。飞行员执行飞行任务需要对界面进行监控，而在视觉监控任务中，被试人员的瞳仁、眨眼率变化显著^{［参考文献 12

百度学术}12］；刘树强通过相关性分析对飞行员心率、眨眼率、脑电功率与工作负荷主观评价值的相关性进行分析，结果表明所选生理指标数据与主观评价数据存在相关性^{［参考文献 8

百度学术}8］；Kilseop等通过实验也证明了眼动、脑电功率和心率与工作负荷存在相关性，可以在一定程度上反映工作负荷大小^{［参考文献 13

百度学术}13］。因此，本文选取脑电功率、心率和眨眼率3项生理指标作为PSO‑SVR模型的输入。脑电功率通过KT88‑2400数字脑电仪进行采集分析；心率通过MD100A心电检测仪进行采集；眨眼率通过ASL 504眼部跟踪仪进行采集。

目前常用的工作负荷主观评价法是CH量表法、SWAT量表法和NASA‑TLX量表法。SWAT量表通过3个维度对工作负荷进行评估，NASA‑TLX量表通过6个维度对工作负荷进行评估，相对于SWAT量表更为全面，评估值更准确。CH量表通过飞行员操纵飞机的主观感受，主要用于评估飞机的操纵性是否良好。因此在选用主观评估量表时采用NASA‑TLX量表获取飞行员的工作负荷主观评估值^［

14］。

2 基于PSO‑SVR的工作负荷预测模型

粒子群算法通过优化SVR模型中核函数的两个重要参数：惩罚系数c和核函数参数g，从而提高算法的预测精度。本节包含构建基于SVR算法的飞行员生理指标数据与主观评价数据之间的非线性映射关系以及基于PSO的SVR模型参数优化两部分。

2.1　面向飞行员工作负荷预测的SVR模型

支持向量回归算法是基于统计学原理对样本数据进行预测，以建立样本输入与输出之间的非线性映射关系。本文采用非线性回归ε‑SVR对飞行员工作负荷预测进行研究。

样本数据集 $\{(x_{1}, y_{1}), (x_{2}, y_{2}), \dots, (x_{30}, y_{30})\}$ ， $x_{i} = (x_{i 1}, x_{i 2}, x_{i 3}), i = 1 ~ 30$ ，表示第i组样本的生理指标数据；x_i₁为第i组样本的脑电功率值；x_i₂为第i组样本的心率；x_i₃为第i组样本的眨眼率；y_i为第i组样本的NASA‑TLX评分值。

本文通过径向基核函数将样本数据集中的 $x$ 映射到高维特征空间并进行线性回归，即有

K (x_{i}, x_{j}) = e x p (- g {‖x_{i} - x_{j}‖}^{2})

（1）

式中g为核函数参数。

SVR算法结构如图2所示，回归函数为

r (x) = \sum_{i = 1}^{30} y_{i} \cdot α_{i}^{*} \cdot K (x_{i}, x) + b

（2）

式中： $α_{i}^{*}$ 为（0，c）之间的常数，c为惩罚系数；b为偏置常数。

图2 工作负荷预测的SVR算法结构图

Fig.2 SVR algorithm structure diagram for workload prediction

为求c和b，引入损失函数有

l (x_{i}) = \{\begin{matrix} 0 |y_{i} - r (x_{i})| \leq ε \\ |y_{i} - r (x_{i})| - ε |y_{i} - r (x_{i})| > ε \end{matrix}

（3）

式中ε为不敏感系数，本文取值为0.01。

所有样本点与回归函数 $r (x)$ 之间距离都小于ε，结合结构分险最小（Structure risk minimization， SRM）准则，回归问题等价于二次凸优化问题，即

m i n \frac{1}{2} (y_{i} \cdot α_{i}^{*})^{2} + c \sum_{i = 1}^{30} l (x_{i})

（4）

s . t . \{\begin{array}{l} r (x_{i}) - y_{i} \leq ε + ξ_{i} \\ y_{i} - r (x_{i}) \leq ε + ξ_{i}^{*} \\ ξ_{i}, ξ_{i}^{*} \geq 0 \end{array}

式中 $ξ_{i}$ 、 $ξ_{i}^{*}$ 为松弛因子。

通过引入拉格朗日函数，结合对偶公式可求得 $α_{i}^{*}$ ，结合KKT条件可求得偏置常数b，即可建立生理指标数据与主观评价数据之间的非线性映射关系。

2.2　基于粒子群算法的SVR参数优化

粒子群优化算法是一种启发式优化算法，通过一组候选解（粒子）借助速度和位置更新公式在解空间内移动以解决优化问题。粒子的移动由自身经验和最优个体的经验进行调整，最终找到最佳位置^［

15］。

目前对SVR参数进行优化的方法主要有试凑法、网格搜索法和随机算法等。试凑法依赖于分析人员的经验，主观性较强；网格搜索法需要遍历网格内的所有参数点，操作复杂；随机算法中使用较多的是遗传算法，遗传算法在处理高维样本问题时，算法的收敛速度、精度有一定的限制。粒子群算法相对于其他算法而言参数设置较少、收敛速度较快且全局搜索能力较强，且模型简单易于实现，通用性较强^［

16］。因此，本文采取粒子群算法对SVR参数进行全局优化搜索。

本文中粒子的位置代表待优化参数c和g的值，粒子的速度代表c和g的变化率。设种群中粒子总数为n，第k代种群中各粒子的位置可以表示为

θ^{(k)} = [\begin{matrix} \begin{array}{l} c_{1}^{(k)} \\ c_{2}^{(k)} \\ ⋮ \\ c_{n}^{(k)} \end{array} & \begin{array}{l} g_{1}^{(k)} \\ g_{2}^{(k)} \\ ⋮ \\ g_{n}^{(k)} \end{array} \end{matrix}]

（5）

第k代群中各粒子的速度可以表示为

V_{θ}^{(k)} = [\begin{matrix} \begin{array}{l} v_{c 1}^{(k)} \\ v_{c 2}^{(k)} \\ ⋮ \\ v_{c n}^{(k)} \end{array} & \begin{array}{l} v_{g 1}^{(k)} \\ v_{g 2}^{(k)} \\ ⋮ \\ v_{g n}^{(k)} \end{array} \end{matrix}]

（6）

对惩罚系数c和核函数参数g进行优化，其流程如图3所示。

图3 基于PSO的SVR参数优化流程

Fig.3 SVR parameter optimization process based on PSO

（1）初始化参数c和g的值及其变化率

惩罚系数c和核函数参数g初始化为

\{\begin{matrix} c_{i}^{(0)} = c^{(m i n)} + (c^{(m a x)} - c^{(m i n)}) \times r \\ g_{i}^{(0)} = g^{(m i n)} + (g^{(m a x)} - g^{(m i n)}) \times r \end{matrix}

（7）

式中： $c^{(m a x)}$ 、 $c^{(m i n)}$ 、 $g^{(m a x)}$ 、 $g^{(m i n)}$ 分别为参数c和g的最大值和最小值，r为（0，1）之间的随机数，i为粒子序号。

惩罚系数c和核函数参数g变化率为

\{\begin{matrix} v_{c i}^{(0)} = (c^{(m a x)} - c^{(m i n)}) + 2 (c^{(m a x)} - c^{(m i n)}) \times r \\ v_{g i}^{(0)} = (c^{(m a x)} - c^{(m i n)}) + 2 (c^{(m a x)} - c^{(m i n)}) \times r \end{matrix}

（8）

（2）计算参数c和g对应的初始适应度值

选用K折交叉验证法样本数据集y_i的均方误差（Mean square error， MSE）作为适应度函数有

f_{i} (θ) = \frac{1}{30} \sum_{i = 1}^{30} [y_{i} (θ) - {\hat{y}}_{i} {(θ)]}^{2}

（9）

式中： $y_{i} (θ)$ 为第i位飞行员的NASA‑TLX评分值； ${\hat{y}}_{i} (θ)$ 为第i位飞行员的NASA‑TLX预测值。

（3）初始化每一组c和g的最优值、全局最优c和g的值

将c和g初始最优值记为

p_{b e s t}^{(0)} = θ^{(0)}

（10）

c和g全局最优初始值 $g_{b e s t}^{(0)}$ 通过寻找 $f (θ^{(0)})$ 中最小值确定。

（4）更新参数c和g的变化率及其值

第k代种群中c和g的变化率更新为

\begin{array}{l} V_{θ}^{(k)} = w^{(k)} \times V_{θ}^{(k)} + C_{1} \times r \times (p_{b e s t}^{(k - 1)} - θ^{(k - 1)}) + \\ C_{2} \times r \times (g_{b e s t}^{(k - 1)} - θ^{k - 1}) \end{array}

（11）

式中：r为（0，1）之间的随机数； $w^{k}$ 为第k代种群中变化率的惯性权重，即^［

17］

w^{(k)} = 0.5 \times \frac{K - k}{K} + 0.4

（12）

式中：K为最大迭代次数；k为当前迭代次数。

为了防止变化率超出范围，对更新后的变化率进行限制有

\{\begin{matrix} v_{i} = (1 - {(\frac{k}{K})}^{α}) \times v_{m a x} v_{i} > (1 - {(\frac{k}{K})}^{α}) \times v_{m a x} \\ v_{i} = (1 - {(\frac{k}{K})}^{α}) \times v_{m i n} v_{i} < (1 - {(\frac{k}{K})}^{α}) \times v_{m i n} \end{matrix}

（13）

式中 $v_{m a x}$ 、 $v_{m i n}$ 为参数c或g变化率的最值。

动态学习因子C₁、C₂^［

18］

\{\begin{matrix} C_{1} = r_{1} + \frac{r_{2} \times k}{K} \\ C_{2} = r_{3} - \frac{r_{4} \times k}{K} \end{matrix}

（14）

式中r₁、r₂、r₃、r₄为常数^［

19］。

第 $k$ 代种群中参数c和g的值更新为

θ^{(k)} = θ^{(k - 1)} + V_{θ}^{(k)}

（15）

超出范围的参数处理如下

\{\begin{matrix} δ_{i}^{(k)} = δ_{i}^{(m a x)} δ_{i}^{(k)} > δ_{i}^{(m a x)} \\ δ_{i}^{(k)} = δ_{i}^{(m i n)} δ_{i}^{(k)} < δ_{i}^{(m i n)} \end{matrix}

（16）

式中： $δ_{i}^{(m a x)}$ 为待优化参数的最大值； $δ_{i}^{(m i n)}$ 为待优化参数的最小值。

（5）更新每组最优的c和g、全局最优的c和g

根据参数c和g当前适应度值更新最优c和g的值为

p_{b e s t}^{(k)} = \{\begin{matrix} θ^{(k)} & f (θ^{(k)}) < f (p_{b e s t}^{(k - 1)}) \\ p_{b e s t}^{(k - 1)} & 其 他 \end{matrix}

（17）

根据每组最优c和g的适应值更新全局最优c和g的值为

g_{b e s t}^{(k)} = \{\begin{matrix} p_{b e s t}^{(k)} & f (p_{b e s t}^{(k)}) < f (p_{b e s t}^{(k - 1)}) \\ g_{b e s t}^{(k - 1)} & 其 他 \end{matrix}

（18）

（6）判断是否满足终止条件

当迭代次数达到指定值时，或者适应度函数值收敛于指定值时停止迭代，否则返回步骤（4）。

2.3　PSO‑SVR预测模型误差分析

根据预测结果相对实际结果的误差百分比对预测结果进行误差分析。相对误差百分比

e r r o r = \frac{|a_{0} (m) - t_{0} (m)|}{t_{0} (m)} \times 100 %

（19）

式中： $a_{0} (m)$ 表示第m组样本数据的预测值； $t_{0} (m)$ 表示第m组样本数据的实际值。

平均相对误差百分比

E = \frac{1}{M} \sum_{m = 1}^{M} (\frac{a_{0} (m) - t_{0} (m)}{t_{0} (m)} \times 100 %)

（20）

式中m=1，2，…，M表示第m个样本数据。

3 实例分析

以某型号民用飞机为例，对该型飞机飞行员的工作负荷进行预测分析。基于本文设计的模拟飞行试验，获取飞行员的脑电功率、心率变化量、眨眼率以及NASA‑TLX量表评分值，其中脑电功率为绝对功率值，心率变化量和眨眼率以“次/min”为单位，脑电功率和NASA‑TLX评分值量纲为一。场景1~4分别获取45组试验数据（见附录1~4）。取每个场景前30组作为训练样本，后15组作为预测样本。

PSO算法基本参数中的种群大小和最大迭代次数是通过多次仿真取预测精度最大时确定，r₁、r₂、r₃、r₄、c^（max）、c^（min）、 $v_{c}^{(m a x)}$ 、 $v_{c}^{(m i n)}$ 、g^（max）、g^（min）、 $v_{g}^{(m a x)}$ 、 $v_{c}^{(m i n)}$ 等参数的取值通过参考相关文献结合实例给出。种群大小为20，最大迭代次数为200，学习因子C₁、C₂计算公式中r₁=1，r₂=0.5，r₃=6，r₄=2^［

20‑21］，参数c的取值边界为c^（max）=100，c^（min）=0.1，参数c的速度边界为

v_{g}^{(m a x)}

=20，

v_{c}^{(m i n)}

=-20，参数g的取值边界为g^（max）=1 000，g^（min）=0.01，参数g的速度边界为

v_{g}^{(m a x)}

=200、

v_{g}^{(m i n)}

=-200^{［参考文献 22

百度学术}22］。

为了比较优化后的预测精度，采用默认参数的SVR预测模型和本文建立的PSO‑SVR模型分别进行训练，NASA‑TLX预测结果如表2所示。将SVR模型与PSO‑SVR模型的NASA‑TLX预测值与实际值之间进行对比如图4所示。

表2 SVR模型与PSO‑SVR模型NASA‑TLX预测结果

Table 2 NASA‑TLX prediction results of SVR model and PSO‑SVR model

编号	场景1			场景2			场景3			场景4
编号	实际值	SVR	PSO‑SVR	实际值	SVR	PSO‑SVR	实际值	SVR	PSO‑SVR	实际值	SVR	PSO‑SVR
1	54	45	56	57	64	60	72	65	69	80	71	75
2	45	47	48	60	66	57	59	50	56	76	70	73
3	47	53	45	55	45	53	63	54	61	78	83	79
4	43	49	40	50	41	48	67	59	64	84	76	80
5	55	50	51	57	49	56	64	70	62	90	81	86
6	48	41	52	59	66	57	60	51	57	72	66	69
7	55	46	49	50	58	53	70	64	64	80	69	76
8	45	52	49	54	61	52	65	60	61	74	65	76
9	48	40	51	49	56	46	59	52	64	78	70	75
10	44	51	47	46	53	51	63	54	60	79	72	83
11	58	50	55	48	57	50	68	59	64	83	74	81
12	53	45	49	57	64	53	70	62	73	92	86	89
13	51	43	53	45	53	42	73	61	70	85	76	80
14	46	50	46	49	58	47	63	54	59	81	71	76
15	52	43	56	53	59	57	69	63	66	76	69	74

图4 SVR/PSO‑SVR模型的NASA‑TLX预测值与实际值对比

Fig.4 Comparison of NASA‑TLX prediction results of SVR/PSO‑SVR model with actual results

对表2中各样本进行误差分析，由式（19）得SVR/PSO‑SVR预测模型的样本相对误差如表3所示。

表3 SVR/PSO‑SVR预测模型相对误差

Table 3 Relative error of SVR/PSO‑SVR prediction models ( % )

编号	场景1		场景2		场景3		场景4
编号	SVR	PSO‑SVR	SVR	PSO‑SVR	SVR	PSO‑SVR	SVR	PSO‑SVR
1	16.7	3.7	12.3	5.3	9.7	4.2	11.3	6.3
2	4.4	6.7	10.0	5.0	15.3	5.1	7.9	3.9
3	12.8	4.3	18.2	3.6	14.3	3.2	6.4	1.3
4	14.0	7.0	18.0	4.0	12.0	4.5	9.5	4.8
5	9.1	7.3	14.0	1.8	9.4	3.1	10.0	4.4
6	14.6	8.3	11.9	3.4	15.0	5.0	8.3	4.2
7	16.4	10.9	16.0	6.0	8.6	8.6	13.8	5.0
8	15.6	8.9	13.0	3.7	7.7	6.2	12.2	2.7
9	16.7	6.3	14.3	6.1	11.9	8.5	10.3	3.8
10	16.0	6.8	15.2	10.9	14.3	4.8	8.9	5.1
11	13.8	5.2	18.8	4.2	13.3	5.9	10.8	2.4
12	15.1	7.6	12.3	7.0	11.4	4.3	6.5	3.3
13	15.7	3.9	17.8	6.7	16.4	4.1	10.6	5.9
14	8.7	0.0	18.4	4.1	14.3	6.3	12.3	6.2
15	17.3	7.7	11.3	7.5	8.7	4.3	9.2	2.6

根据表3中SVR/PSO‑SVR误差，由式（20）可得场景1~4的SVR预测模型样本平均相对误差分别为13.8%、14.8%、12.2%和9.9%，PSO‑SVR预测模型样本平均相对误差分别为6.3%、5.3%、5.2%和4.1%，将表3中SVR误差与PSO‑SVR误差以折线图的形式进行对比如图5所示。

图5 SVR/PSO‑SVR模型相对误差对比

Fig.5 Relative error comparison of SVR/PSO‑SVR model

由图5可知，除个别样本出现优化后误差增大的情况，绝大多数样本经过PSO优化后预测误差都减小，与SVR相比，PSO‑SVR的预测精度更高，表明利用粒子群算法对SVR参数进行优化的方法非常有效。

为了验证上述模型的有效性与正确性，设置区别于场景1~4的4组不同的飞行场景记为场景5~8，获取飞行员的工作负荷主观评价值与客观生理评价数据。利用上述评估模型预测出新设置的4组不同飞行场景下飞行员工作负荷主观评价值，并与实际获取到的值进行对比，分析其误差是否在可接受的范围之内，4组场景的误差对比结果的如图6所示。从图6的误差对比可以看出，针对不同场景优化后的预测模型能够将误差控制在10%以内。

图6 新增场景下的SVR/PSO‑SVR模型相对误差

Fig.6 Relative error of SVR/PSO‑SVR models in new scenes

4 结　　论

本文对飞行员工作负荷评价方法进行了研究，提出了基于PSO‑SVR的飞行员工作负荷预测方法，通过对比验证说明了本文方法的合理性。本文将支持向量回归模型应用在飞行员工作负荷预测上，为飞行员工作负荷评估提供了一种新方法。为了进一步提高工作负荷预测的精度，本文提出了基于PSO‑SVR的飞行员工作负荷预测模型，对SVR模型进行改进，通过与本文所提的SVR模型的对比以及与基于BP神经网络的机组工作量评估模型的对比^［

8］，PSO‑SVR工作负荷预测模型具有更高的预测精度，表明本文所提方法是有效的。

实例分析中多组场景的预测结果表明，本文所提模型能够将飞行员工作负荷预测误差控制在可接受的范围之内。表明本文所提的PSO‑SVR工作负荷预测方法对不同飞行场景的适用性。利用本文所提方法，在进行飞行员工作负荷适航符合性验证时，只需测量飞行员生理指标数据而不用多次对飞行员进行工作负荷主观评价，既能够避免飞行员的主观不确定性，提高符合性验证的准确性，同时提高了符合性验证的效率。

本文所提PSO‑SVR模型在预测精度上的确要优于SVR模型，但是从不同场景的预测精度上看，仍然存在一定误差，原因在于建模所用数据存在的局限性，数据量和建模所用的工作负荷主观评估值的准确性是主要原因，后续需要通过增加数据量以及增加建模所用输出数据（工作负荷主观评估值）的准确性提高模型的预测精度。

附录

场景1~4试验样本数据，如附表1~4所示。

附表1 场景1样本数据

Annexed table 1　Sample data for scene 1

编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值	编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值
1	10.30	12.0	19.3	45.0	24	12.40	14.0	16.6	51.0
2	12.00	13.6	17.0	50.0	25	13.40	15.0	15.2	53.7
3	10.60	12.3	18.9	46.0	26	12.40	14.0	16.6	51.0
4	13.80	15.3	14.7	55.0	27	9.90	11.6	19.8	44.0
5	13.10	14.6	15.7	53.0	28	9.60	11.3	20.3	43.0
6	11.30	13.0	18.0	48.0	29	10.30	12.0	19.3	45.0
7	13.10	14.0	15.7	52.0	30	9.60	11.3	20.3	43.0
8	12.70	14.5	16.1	53.0	31	13.40	15.0	15.2	54.0
9	13.80	15.3	14.7	55.0	32	10.30	12.0	19.3	45.0
10	13.00	15.3	14.7	54.0	33	11.00	12.6	18.4	47.0
11	9.90	11.6	19.8	44.0	34	9.61	11.3	20.3	43.0
12	11.00	12.6	18.4	47.0	35	13.80	15.3	14.7	55.0
13	12.00	13.6	17.0	50.0	36	11.30	13.0	18.0	48.0
14	13.40	15.0	15.2	54.0	37	13.80	15.3	14.7	55.0
15	11.70	13.3	17.5	49.0	38	10.30	12.0	19.3	45.0
16	11.70	14.0	17.5	52.7	39	11.30	13.0	18.0	48.0
17	11.30	13.0	18.0	48.0	40	9.96	11.6	19.8	44.0
18	11.70	13.3	17.5	49.0	41	14.80	16.3	13.4	58.0
19	13.40	15.0	15.2	54.0	42	13.10	14.6	15.7	53.0
20	9.90	11.6	19.8	44.0	43	12.40	14.0	16.6	51.0
21	11.00	12.6	18.4	47.0	44	10.60	12.3	18.9	46.0
22	12.00	13.0	17.0	50.0	45	12.70	14.3	16.1	52.0
23	13.40	15.0	15.2	54.0

附表2 场景2样本数据

Annexed table 2　Sample data for scene 2

编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值	编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值
1	19.30	19.2	15.1	57	24	18.30	18.5	15.7	58
2	19.60	19.5	14.9	60	25	22.80	21.9	13.0	53
3	19.10	19.1	15.2	54	26	17.80	18.1	16.0	60
4	20.00	19.8	14.6	62	27	22.60	21.7	13.1	57
5	23.90	22.7	12.3	58	28	21.30	20.7	13.9	60
6	21.50	20.9	13.8	59	29	17.80	18.1	16.0	54
7	17.30	17.8	16.3	54	30	22.80	21.8	13.0	62
8	23.90	22.7	12.3	60	31	18.00	18.2	15.9	57
9	20.60	20.2	14.3	53	32	11.10	13.1	20.0	60
10	18.20	18.4	15.8	60	33	12.60	14.2	19.1	55
11	24.40	23.1	12.0	58	34	9.55	11.9	21.0	50
12	17.20	17.7	16.3	53	35	18.70	18.8	15.4	57
13	18.00	18.3	15.9	60	36	13.40	14.8	18.7	59
14	24.00	22.8	12.3	57	37	18.70	18.8	15.4	50
15	20.80	20.4	14.2	60	38	11.10	13.1	20.0	54
16	23.00	22.1	12.8	54	39	13.40	14.8	18.7	49
17	18.00	18.3	15.9	62	40	10.30	12.5	20.5	46
18	20.30	20.0	14.5	58	41	21.00	20.5	14.1	48
19	22.70	21.8	13.1	59	42	17.20	17.7	16.4	57
20	18.00	18.3	15.9	54	43	15.70	16.5	17.3	45
21	24.80	23.4	11.8	60	44	11.80	13.6	19.6	49
22	21.40	20.8	13.8	53	45	16.40	17.1	16.8	53
23	22.40	21.6	13.2	60

附表3 场景3样本数据

Annexed table 3　Sample data for scene 3

编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值	编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值
1	29.9	27.7	13.4	70	24	29.7	27.5	13.5	69
2	23.6	22.7	15.8	61	25	26.8	25.2	14.6	66
3	28.9	26.9	13.8	68	26	23.0	22.2	16.1	61
4	23.4	22.5	15.9	61	27	23.6	22.7	15.8	61
5	23.7	22.8	15.8	61	28	29.7	27.5	13.5	69
6	24.5	23.4	15.5	63	29	26.7	25.2	14.6	65
7	25.7	24.3	15.0	64	30	26.2	24.8	14.8	65
8	31.5	28.9	12.8	72	31	31.7	29.1	12.7	72
9	26.7	25.2	14.6	65	32	21.8	21.3	16.5	59
10	24.6	23.5	15.4	63	33	24.9	23.7	15.3	63
11	29.4	27.3	13.6	69	34	27.9	26.1	14.2	67
12	23.3	22.4	15.9	61	35	25.6	24.3	15.0	64
13	28.7	26.7	13.9	68	36	22.6	21.9	16.2	60
14	26.6	25.1	14.6	65	37	30.2	27.9	13.3	70
15	25.2	24.0	15.2	64	38	26.4	24.9	14.7	65
16	27.5	25.8	14.3	66	39	21.8	21.3	16.5	59
17	31.2	28.7	12.9	71	40	24.9	23.7	15.3	63
18	24.9	23.8	15.3	63	41	28.7	26.7	13.9	68
19	30.8	28.4	13.1	71	42	30.2	27.9	13.3	70
20	27.3	25.6	14.4	66	43	32.5	29.7	12.4	73
21	22.8	22.1	16.1	60	44	24.9	23.7	15.3	63
22	30.1	27.8	13.3	70	45	29.4	27.3	13.6	69
23	30.4	28.1	13.2	70

附表4 场景4样本数据

Annexed table 4　Sample data for scene 4

编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值	编号	脑电功率	心率变化量	眨眼率	NASA‑TLX值
1	47.7	34.8	11.70	84	24	47.3	34.5	11.80	83
2	50.6	36.6	10.80	86	25	45.3	33.2	12.40	81
3	42.1	31.2	13.40	78	26	52.9	38.1	10.00	88
4	45.8	33.5	12.30	82	27	53.7	38.6	9.70	89
5	52.7	38.0	10.10	88	28	44.8	32.9	12.60	81
6	52.8	38.0	10.00	85	29	45.7	33.5	12.30	82
7	43.5	32.1	13.00	80	30	50.0	36.3	10.90	86
8	53.6	38.5	9.82	89	31	43.8	32.2	12.90	80
9	41.7	30.9	13.60	78	32	39.5	29.5	14.30	76
10	43.5	32.1	13.30	80	33	41.6	30.9	13.60	78
11	42.4	31.4	13.40	79	34	48.1	35.0	11.60	84
12	48.8	35.5	11.30	85	35	54.5	39.1	9.52	90
13	55.3	39.7	9.26	91	36	35.2	26.7	15.60	72
14	47.6	34.7	11.70	84	37	43.8	32.2	12.90	80
15	47.2	34.5	11.80	83	38	37.3	28.1	14.90	74
16	51.3	37.1	10.50	87	39	41.6	30.9	13.60	78
17	54.4	39.1	9.50	90	40	42.7	31.6	13.30	79
18	55.5	39.8	9.20	91	41	47.0	34.3	11.90	83
19	51.7	37.3	10.40	87	42	56.6	40.5	8.85	92
20	47.4	34.6	11.70	83	43	49.1	35.7	11.20	85
21	55.1	39.5	9.30	91	44	44.8	32.9	12.60	81
22	43.8	32.3	12.90	80	45	39.5	29.5	14.30	76
23	47.2	34.5	11.80	83

参考文献

刘选民，李凡. 大中型运输机飞行事故［M］. 北京：航空工业出版社， 2012： 2‑3. [百度学术]

LIU Xuanmin， LI Fan. Flight accidents of large and medium transport aircraft［M］. Beijing： Aviation Industry Press， 2012： 2‑3. [百度学术]

FAA. Installed systems and equipment for use by the flight crew： AC 25.1302‑1［S］. Washington： FAA， 2013. [百度学术]

Wilson G F. An analysis of mental workload in pilots during flight using multiple psychophysiological measures［J］. The International Journal of Aviation Psychology， 2002， 12（1）： 3‑18. [百度学术]

Veltman J A. A comparative study of psychophysiological reactions during simulator and real flight［J］. The International Journal of Aviation Psychology， 2002， 12（1）： 33-48. [百度学术]

Mansikka H， Virtanen K， Harris D. Comparison of NASA‑TLX scale， modified Cooper‑Harper scale and mean inter‑beat interval as measures of pilot mental workload during simulated flight tasks［J］. Ergonomics， 2019， 62（2）： 246-254. [百度学术]

辛甜，罗志燕，蒋利添，等. 飞机着陆阶段飞行员工作负荷分析［J］. 人类工效学， 2019， 25（2）： 62‑65. [百度学术]

XIN Tian， LUO Zhiyan， JIANG Litian， et al. Analysis of pilot's workload during the landing phase ［J］. Chinese Journal of Ergonomics. 2019， 25（2）： 62‑65. [百度学术]

李琳君，李向南.飞行负荷评价技术与方法［J］. 直升机技术， 2019（2）： 37‑41，46. [百度学术]

LI Linjun， LI Xiangnan. The technology and methods of flight workload evaluation［J］. Helicopter Technique， 2019（2）： 37‑41，46. [百度学术]

刘树强. 民机驾驶舱最小机组工作量评估方法研究［D］. 南京：南京航空航天大学， 2012. [百度学术]

LIU Shuqiang. Research on assessment method of min crew workload application to civil aircraft cockpit［D］. Nanjing： Nanjing University of Aeronautics and Astronautics， 2012. [百度学术]

李金波，许百华. 人机交互过程中认知负荷的综合测评方法［J］. 心理学报， 2009， 41（1）： 35‑43. [百度学术]

LI Jinbo， XU Baihua. Synthetic assessment of cognitive load in human‑machine interaction process［J］. Acta Psychologica Sinica， 2009， 41（1）： 35‑43. [百度学术]

Mazaeva N， Ntuen C. Self‑organizing map （SOM） model for mental workload classification［J］ IEEE Transaction on Systems， Man， and cyber nentics‑part A： Systems and Humans， 2001， 31（6）： 1822‑1825. [百度学术]

卫宗敏，郝红勋，徐其志，等. 飞行员脑力负荷测量指标和评价方法研究进展［J］. 科学技术与工程， 2019， 19（24）： 1‑8. [百度学术]

WEI Zongmin， HAO hongxun， XU Qizhi， et al. Review of research on pilots' mental workload evaluation indices and measurement methods［J］. Science Technology and Engineering， 2019， 19（24）： 1‑8. [百度学术]

Noel J B， Bauer K W， Lanning J W. Improving pilot mental workload classification through feature exploitation and combination： A feasibility study［J］. Computers & Operations Research， 2005， 32（10）： 2713‑2730. [百度学术]

Kilseop R， Myung R. Evaluation of mental workload with a combined measure based on physiological indices during a dual task of tracking and mental arithmetic［J］. International Journal of Industrial Ergonomics， 2005， 35（11）： 991‑1009. [百度学术]

完颜笑如，庄达民. 飞行员脑力负荷测量与应用［M］. 北京：科学出版社， 2014： 15， 49‑52. [百度学术]

WANYAN Xiaoru， ZHUANG Damin. Pilot mental workload measurement and application ［M］. Beijing： Science Press， 2014： 15， 49‑52. [百度学术]

李丽，牛奔. 粒子群优化算法［M］. 北京：冶金工业出版社， 2009. [百度学术]

LI Li， NIU Ben. Particle swarm optimization algorithm ［M］. Beijing： Metallurgical Industry Press， 2009. [百度学术]

王建国，张文兴. 支持向量机建模及其智能优化［M］.北京：清华大学出版社， 2015. [百度学术]

WANG Jianguo， ZHANG Wenxing. Support vector machine modeling and intelligent optimization［M］. Beijing： Tsinghua University Press， 2015. [百度学术]

郭文忠，陈国龙. 离散粒子群优化算法及其应用［M］.北京：清华大学出版社， 2012. [百度学术]

GUO Wenzhong， CHEN Guolong. Discrete particle swarm optimization algorithm and its application［M］. Beijing： Tsinghua University Press， 2012. [百度学术]

康晓凡. 基于优化支持向量回归机的充电站负荷预测研究［D］. 西安：西安科技大学， 2019. [百度学术]

KANG Xiaofan. Research on load forecasting method of charging station based on support vector regression machine［D］. Xi'an： Xi'an University of Science and Technology， 2019. [百度学术]

张文成，周穗华，吴志东，等. 改进策略的粒子群算法及其实验测试［J］. 计算机工程与设计， 2012， 33（10）： 3945‑3949. [百度学术]

ZHANG Wencheng， ZHOU Huihua， WU Zhidong， et al. Improved particle swarm optimization and experimentation［J］. Computer Engineering and Design， 2012， 33（10）： 3945‑3949. [百度学术]

邱春荣. 基于PSO‑SVR的瓦斯浓度预测方法研究［D］. 徐州：中国矿业大学， 2017. [百度学术]

QIU Chunrong. Study on gas concentration prediction methods based on the PSO‑SVR［D］. Xuzhou： China University of Mining and Technology， 2017. [百度学术]

王瑞，李芯蕊，马双斌. 基于PSO‑SVR的网络态势预测模型［J］. 信息安全研究， 2018， 4（8）： 734‑738. [百度学术]

WANG Rui， LI Xinrui， MA Shuangbin. Network situation prediction model based on PSO‑SVR［J］. Journal of Information Security Research， 2018， 4（8）： 734‑738. [百度学术]

张冉. 基于PSO‑SVR的渲染时间预估研究［D］. 镇江：江苏科技大学， 2016. [百度学术]

ZHANG Ran. The research of time prediction for rendering based on PSO‑SVR［D］. Zhenjiang： Jiangsu University of Science and Technology， 2016. [百度学术]

请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器，其他浏览器不建议使用!

基于PSO‑SVR的飞行员工作负荷预测 PDF

摘要

关键词

1 模拟飞行试验设计

1.1 模拟飞行人员选取

1.2 模拟飞行阶段确定

1.3 模拟飞行状态选取

1.4 模拟飞行场景设计

1.5 试验数据采集