灰色预测模型【行业内容】

灰色预测模型,马文丽 201620154,主要内容,灰色预测模型,2,一灰色预测的概念 ；,二 灰色生成数列；,四 灰色预测计算实例。,三 灰色模型GM；,一、灰色预测的概念,（1）灰色系统、白色系统和黑色系统,白色系统是指一个系统的内部特征是完全已知的，即系统的信息是完全充分的。 黑色系统是指一个系统的内部信息对外界 来说是一无所知的，只能通过它与外界的联系来加以观测研究。 灰色系统内的一部分信息是已知的，另一部分信息是未知 的，系统内各因素间有不确定的关系。,灰色预测法是一种对含有不确定因素的系 统进行预测的方法。 灰色预测是对既含有已知信息又含有不确定信息的系统进行预则，就是对在一定范围内 变化的、与时间有关的灰色过程进行预测。,（2）灰色预测法,一、灰色预测的概念,灰色预测通过鉴别系统因素之间发展趋势的相异程度，即进行关联分析，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列, 然后建立相应的微分方程模型，从而预测事物未来发展趋势的状况。 灰色预测法用等时距观测到的反映预测对象特征的一系列数量值构造灰色预测模型，预测未来某一时刻的特征量，或达到某一特征量的时间。,一、灰色预测的概念,（3）灰色预测的四种常见类型, 数列预测 对某现象随时间的顺延而发生的变化所做的预测定义为数列预测。例如对消费物价指数的预测，需要确定两个变量，一个是消费物价指数的水平,另一个是这一水平所发生的时间。 灾变预测 对灾害或异常突变可能发生的时间预测称为灾变预测。例如对地震时间的预测。,一、灰色预测的概念,系统预测 对系统中众多变量间相互协调关系的发展变化所进行的预测称为系统预测。例如市场中替代商品、相互关联商品销售量互相制约的预测。 拓扑预测 将原始数据作曲线，在曲线上按定值寻找该定值发生的所有时点，并以该定值为框架构成时点数列，然后建立模型预测未来该定值所发生的时点。,一、灰色预测的概念,二、灰色生成数列,灰色系统理论认为，尽管客观表象复杂，但总是有整体功能的，因此必然蕴含某种内在规律。关键在于如何选择适当的方式去挖掘和利用它。灰色系统是通过对原始数据的整理来寻求其变化规律的，这是一种就数据寻求数据的现实规律的途径，即为灰色序列的生成。一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。数据生成的常用方式有累加生成、累减生成和加权累加生成。,（1）累加生成 把数列各项（时刻）数据依次累加的过程称为累加生 成过程（AGO ）。由累加生成过程所得的数列称为 累加生成数列。设原始数列为 令 称所得到的新数列为数列 的1次累加生成数列。类似 地有 称为 的r次累加生成数列。,累加生成计算示例,例：x (0)=(x (0) (k) k=1，2，3，4，5) =x(0)(1)，x(0)(2)，x(0)(3)，x(0)(4)，x(0)(5) =(3.2，3.3，3.4，3.6，3.8) 求 x(1)(k) 解：,累加生成的特点,一般经济数列都是非负数列。累加生成 能使任意非负数列、摆动的与非摆动的，转 化为非减的、递增的。 原始数列作图 1AGO作图,某市的汽车销售量,递增的规律,原始数列作图,1AGO作图,有明显的指数关系的规律,某钢厂产量,某地区作物产量,s型变化规律,（2）累减生成 对于原始数据列依次做前后相邻的两个数据相减的 运算过程称为累减生成过程IAGO。如果原始数据列 为 令 称所得到的数列 为 的1次累减生成数列。 注：从这里的记号也可以看到，从原始数列 ，得 到新数列 ，再通过累减生成可以还原出原始数列 。实际运用中在数列 的基础上预测出 ，通过 累减生成得到预测数列 。,累减生成计算示例,（3）加权邻值生成 设原始数列为 称 为数列 的邻值。 为后邻值， 为前邻值，对于常 数 ，令 由此得到的数列 称为数列 在权 下的邻值生 成数，权 也称为生成系数。 特别地，当生成系数 时，则称 为均值生成数，也称等权邻值生成数。,灰色系统理论是基于关联空间、光滑离散函数等概念定义灰导数与灰微分方程，进而用离散数据列建立微分方程形式的动态模型，即灰色模型是利用离散随机数经过生成变为随机性被显著削弱而且较有规律的生成数，建立起的微分方程形式的模型，这样便于对其变化过程进行研究和描述。 G表示grey（灰色），M表示model（模型）,灰色模型(1,1),设 为原始数列，其1次累 加生成数列为 ，其中 定义 的灰导数为 令 为数列 的邻值生成数列，即 于是定义GM（1，1）的灰微分方程模型为,即或 （1） 在式（1）中， 称为灰导数，a称为发展系数， 称为白化背景值，b称为灰作用量。 将时刻表 代入（1）式有 引入矩阵向量记号： 于是GM（1，1）模型可表示为,现在问题归结为求a,b在值。用一元线性回归，即最小二乘法求它们的估计值为 注：实际上回归分析中求估计值是用软件计算的，有标准程序求解，如matlab等。,GM（1，1）的白化型 对于GM（1，1）的灰微分方程（1），如果将灰导数 的时刻 视为连续变量t，则 视为时间t函数, 于是 对应于导数 白化背景值 对应于 。于是GM（1，1）的灰微 分方程对应于的白微分方程为 （2） 称之为GM(1，1)的白化型。,三、GM（1，1）灰色预测的步骤,1.数据的检验与处理 为了保证GM（1，1）建模方法的可行性，需要对已知数据做必要的检验 处理。 设原始数据列为了 ，计算数列的级比 如果所有的级比都落在可容覆盖区间 内，则数据列 可以建立GM（1，1）模型且可以进行灰 色预测。否则，对数据做适当的变换处理，如平移变换： 取C使得数据列 的级比都落在可容覆盖内。,2. 建立GM（1，1）模型 不妨设 满足上面的要 求，以它为数据列建立GM（1，1）模型 用回归分析求得a,b的估计值，于是相应的白化模型为 解为 （3） 于是得到预测值 从而相应地得到预测值：,3. 检验预测值 （1）残差检验：计算相对残差 如果对所有的 ，则认为达到较高的要求：否则，若 对所有的 ，则认为达到一般要求。 （2）级比偏差值检验：计算 如果对所有的 ，则认为达到较高的要求；否则 若对所有的 ，则认为达到一般要求。,四、灰色预测计算实例,例 北方某城市19861992 年道路交通噪声平均声级数据见表6 表6 市近年来交通噪声数据dB(A),序号 年份 eq L 1 1986 71.1 2 1987 72.4 3 1988 72.4 4 1989 72.1 5 1990 71.4 6 1991 72.0 7 1992 71.6,第一步: 级比检验 建立交通噪声平均声级数据时间序列如下：,=(71.1, 72.4, 72.4, 72.1, 71.4, 72.0, 71.6),（1）求级比 (k),=(0.982,1,1.0042,1.0098, 0.9917,1.0056) （2）级比判断 由于所有的 (k)0.982,1.0098，k = 2,3,7，故可以用x(0)作满意的GM（1，1）建模。 第二步: GM（1，1）建模 （1）对原始数据 作一次累加，即 = (71.1,143.5,215.9,288,359.4,431.4,503） （2）构造数据矩阵B 及数据向量Y,（3）计算u 于是得到a = 0.0023，b = 72.6573。 （4）建立模型,求解得,（5）求生成数列值 及模型还原值 ： 令k = 1,2,3,4,5,6,由上面的时间响应函数可算得 ，其中取 由 取k = 2,3,4,7，得 (71.1, 72.4, 72.2, 72.1, 71.9, 71.7, 71.6),第三步: 模型检验 模型的各种检验指标值的计算结果见表 7. 表7 GM(1,1)模型检验表 序号 年份 原始值 模型值 残差 相对误差 级比偏差 1 1986 71.1 71.1 0 0 0 2 1987 72.4 72.4 -0.0057 0.01% 0.0023 3 1988 72.4 72.2 0.1638 0.23% 0.0203 4 1989 72.1 72.1 0.0329 0.05% -0.0018 5 1990 71.4 71.9 -0.4984 0.7% -0.0074 6 1991 72.0 71.7 0.2699 0.37% 0.0107 1992 71.6 71.6 0.0378 0.05% -0.0032 经验证，该模型的精度较高，可进行预测和预报。,