刘国春

摘 要：针对四川通用航空产业发展需求，本研究应用灰色预测方法建立了四川通用航空飞行数据预测模型。针对四川通用航空各个领域飞行数据不平衡的特点，从工业、农业、其他领域分别建立了通用航空飞行数据的灰色预测模型，并进行了检验与修正。该独立模型能够较好地应用于通用航空飞行预测，并且能够准确揭示各领域的发展态势，具有一定的工程价值。

关键词：通用航空;分析模型;灰色预测;飞行数据

中图分类号：V323.1  文献标识码：A

随着国内低空空域改革的不断深化和通用航空产业环境的持续改善，通用航空产业发展呈现出如火如荼的态势，因此有必要对通用航空发展数据进行挖掘，从而了解通航产业的发展态势，便于对通用航空产业的周边布局进行目的性规划。对于通用航空的发展预测，常见的预测方法包括数据包络分析[1]（DEA）、聚类分析，蚁群优化等方法，但是此类方法需要较多的数据来源进行分析，对于较少数据的信息判断性不高。而灰色系统模型作为一种有效的预测手段，其单数列的微分模型由于较好的拟合和外推特性，广泛应用于通用航空的发展预测[2-5]，本研究由于数据来源较为单一，因此选用灰色预测方法建立四川通用航空飞行数据的分析模型。

1 四川通航飞行数据现状

本研究首先统计了2016年全国的通用航空飞行小时数据[6]，从表1中可以看到，四川的通用航空飞行小时位居全国第一，占比达到当年全国通用航空飞行小时总数764685小时的36.69%。然而分项的飞行数据却相差很大，工业位于全国第五，农业排名全国第八，而其他领域主要是为飞行培训，究其原因是由于四川拥有全国最大通用航空单位（中国民用航空飞行学院）的飞行小时数为261189小时。由于三类分项比例太过悬殊，只使用一个合计项进行预测会出现数据失真的情况，较大的数据会掩盖小数据的变化，从而造成判断失误，因此本研究针对四川通用航空飞行小时进行分项灰色预测，分别建立四川工业、农业、其他领域的通用航空飞行预测模型。

2 灰色模型检验方法描述

对系统进行灰色预测的主要步骤如下：对问题进行描述、建立数学模型、求解方程得出预测模型、模型的检验，只有检验合格才能用作预测。

2.1 灰色预测模型建立

GM（1，1）模型的建立如下：

原始数据：

x（0）=（x（0）（1），x（0）（2），x（0）（3），…，x（0）（n））（1）

一次累加生成数据

x（1）=（x（1）（1），x（1）（2），x（1）（3），…，x（1）（n））（2）

其中：

x（1）（k）=∑ki=1x（0）（i）（3）

（k=1，2，3，…，n）

建立微分方程：

dx（1）dt+Ax（1）=U（4）

其中A为发展灰数;U为内生控制灰数。

设α︿=（A，U），由最小二乘法有：

α︿=（BTB）-1BTY1（5）

其中：

B=-12x（1）（1）+x（1）（2）   1

-12x（1）（2）+x（1）（3）1

-12x（1）（n-1）+x（1）（n）1

Y1=x（0）（2）

x（0）（3）

x（0）（n）

解得（4）的解為：

x︿（1）（k+1）=x（0）（1）-AUe-ak+UA（6）

2.2 模型检验

模型的检验主要有以下三个方面：

2.2.1 关联度检验

设参考数列为：

X0=（x01，x02，x03，…x0n）

Xi=（xi1，xi2，xi3，…xin）

则有

ξi（k）=miniminkX0（k）-Xi（k）+ρmaximaxkX0（k）-Xi（k）X0（k）-Xi（k）+ρmaximaxkX0（k）-Xi（k）（7）

其中k表示时间（k=1，2，…n），n表示数列个数，ρ∈[0，1]为分辨系数，一般取ρ=0.5，式中的miniminkX0（k）-Xi（k）和maximaxkX0（k）-Xi（k）分别称为两级最小差和两级最大差，ρ越大，分辨率越大，ρ越小，分辨率越小。

数列Xi对X0关联度R表示如下：

Ri=1n∑nk=1ξi（k）（8）

根据上述计算方法，可得到生成数列x︿（0）（i）与原始数列x（0）（i）的关联系数，然后计算出关联度，根据经验，一般取ρ=0.5时关联度R>0.6就认为模型合格。

2.2.2 后验差检验

原始数列标准差：

S1= ∑X（0）（i）-X-（0）2n-1（9）

绝对误差序列标准差

S2= ∑Δ（0）（i）-Δ-（0）2n-1（10）

方差比：

C=S2S1（11）

小误差概率：

P=pΔ（0）（i）-Δ-（0）<0.6745S1

令

ei=Δ（0）（i）-Δ-（0），S0=0.6745S1

则

P=pei

在灰色预测中，将小误差概率P的值和方差比C的值作为模型等级的检验标准，具体标准见表2。

2.2.3 残差检验