3.灰色预测模型（Gray Forecast Model）

一、引言

1、什么是灰色预测法？

（1）灰色系统、白色系统和黑色系统
**白色系统（白箱系统）**是指一个系统的内部特征是完全已知的，既系统信息是完全充分的。
**黑色系统（黑箱系统）**是一个系统的内部信息对外界来说是一无所知的，只能通过它与外界的联系来加以观测研究。
灰色系统介于白色和黑色之间，灰色系统内的一部分信息是已知的，另一部分信息是未知的，系统内各因素间有不确定的关系。
一般地说，社会系统、经济系统、生态系统都是灰色系统。例如物价系统，导致物价上涨的因素很多，但已知的却不多，因此对物价这一灰色系统的预测可以用灰色预测方法。

（2）灰色系统理论
灰色系统理论认为对既含有已知信息又含有未知或非确定信息的系统进行预测，就是对在一定方位内变化的、与时间有关的灰色过程的预测。尽管过程中所显示的现象是随机的、杂乱无章的，但毕竟是有序的、有界的，因此这一数据集合具备潜在的规律，灰色预测就是利用这种规律建立灰色模型对灰色系统进行预测。

（3）灰色预测
灰色预测通过鉴别系统因素之间发展趋势的相异程度，即进行关联分析，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列，然后建立相应的微分方程模型，从而预测事物未来发展趋势的状况。其用等时距观测到的反应预测对象特征的一系列数量值构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。

（4）关联度
一般表达方式为GM（n,x）模型，其含义是：用n阶微分方程对x个变量建立模型

（5）关联分析
关联分析实际上是动态过程发展态势的量化比较分析。所谓发展态势比较，也就是系统各时期有关统计数据的集合关系的比较。例如，某地区1977~1983年总收入与养猪、养兔收入资料见下表格。

二、 灰色生成数列

灰色系统理论认为，尽管客观表象复杂，但总是有整体功能的，因此必然蕴含某种内在规律。关键在于如何选择适当的方式去挖掘和利用它。
灰色系统是通过对原始数据的整理来寻求其变化规律的，这是一种就数据寻求数据的现实规律的途径，也就是灰色序列的生产。
一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。数据生成的常用方式有累加生成、累减生成和加权累加生成，将灰色系统中的未知因素弱化，强化已知因素的影响程度，最后构建一个以 时间 为变量的连续微分方程，通过数学方法确定方程中的参数，从而实现预测目的。

（1）累加生成（AGO）
设原始序列为
令
那么

为数列 x⁽⁰⁾ 的1次累加生成数列。
类似的有：

称为 x⁽⁰⁾的r次累加生成数列。

累加生成的特点：
一般经济数列都是非负数列。累加生成能使任意非负数列、摆动的与非摆动的，转化为非减的、递增的。

（2）累减生成（IAGO）
如果原始数列为

令

称x⁽⁰⁾为为x⁽¹⁾的1次累减生成数列。
可以看出，通过累加数列的到的新数列，可以通过累减生成还原出原始数列。

（3）加权邻值生成
如果原始数列为

称任意一对相邻元素x⁽⁰⁾(k-1)，x⁽⁰⁾(k)互为邻值。
对于常数α∈[0,1]，令
由此得到的数列称为邻值生成数，权α 也称为生成系数。特别地，当α=0.5 时，则称该数列为均值生成数，也称为等权邻值生成数。

三.、灰色模型GM(1,1)

灰色系统理论是基于关联空间、光滑离散函数等概念定义灰导数与灰微分方程，进而用离散数据列建立微分方程形式的动态模型，即灰色模型是利用离散随机数经过生成变为随机性被显著削弱而且较有规律的生成数，建立起的微分方程形式的模型，这样便于对其变化过程进行研究和描述。
G表示grey（灰色），M表示model（模型）
定义x⁽¹⁾的灰导数为：
令z⁽¹⁾(k)为数列x⁽¹⁾的邻值生成数列，即

于是定义GM(1,1,)的灰微分方程模型为：

或

其中x⁽⁰⁾(k)称为灰导数，α称为发展系数，z⁽¹⁾(k)称为白化背景值，b称为灰作用量。
将时刻k=2,3,…,n代入上式有

引入矩阵向量记号：

于是GM(1,1)模型可表示为Y=Bu
那么现在的问题就是求a和b的值，我们可以用一元线性回归，也就是最小二乘法求它们的估计值：

GM(1,1)的白化型：
对于GM(1,1)的灰微分方程，如果将时刻k=2,3,…,n视为连续变量t，则之前的x⁽¹⁾视为时间t函数，于是灰导数x⁽⁰⁾(k)变为连续函数的导数dx⁽¹⁾/dt，白化背景值z⁽¹⁾(k)对应于导数x⁽¹⁾(t)，于是GM(1,1)的灰微分方程对应于的白微分方程为：

四、GM（1，1）灰色预测的步骤

（1）数据的检验与处理
为了保证GM（1，1）建模方法的可行性，需要对已知数据做必要的检验处理。设原始数据列为

首先计算数列的级比：

如果所有的级比都落在可容覆盖区间

内，则数列x⁽⁰⁾可以建立GM（1，1）模型且可以进行灰色预测。否则，对数据做适当的变换处理，如平移变换：

取c使得数据列的级比都落在可容覆盖内。

（2）建立GM(1,1)模型
不妨设
满足上面的要求，以它为数据列建立GM（1，1）模型：
用回归分析求得a,b的估计值，于是相应的白化模型为：

解为
于是得到预测值：

从而相应地得到预测值：

其中k=1，2，……，n-1

五、检验预测值

（1）残差检验
计算相对残差：

如果对所有的|ε(k)|<0.1，则认为到达较高的要求；否则，若对所有的|ε(k)|<0.2，则认为达到一般要求。

（2）级比偏差值检验
计算：

如果对所有的|ρ(k)|<0.1，则认为达到较高的要求；否则，对于所有的 |ρ(k)|<0.2，则认为达到一般要求。

六、. 计算实例

七、六中计算实例的Python代码

import pandas as pd
import numpy as np

def