灰色预测GM(1,1)

简介: 灰色预测GM(1,1)

一、灰色预测模型简介

灰色预测的主要特点就是模型使用的不是原始数据列,而是生成的数据序列。其核心是灰色模型(Grey Model,GM),即对原始数据作累加生成(累加或差分,都能弱化噪声)能得到近似指数规律再进行建模的方法。

优点:

  1. 所需数据少,一般只需要4个以上数据,能解决历史数据少,序列完整性和可靠性低的问题
  2. 利用微分方程来充分挖掘系统的本质,精度高
  3. 能将无规律的原始数据进行生成得到规律性较强的生成序列
  4. 运算简便,易于检验,不考虑分布规律,不考虑变化趋势

缺点:

只适合短期预测,只适合指数增长的预测。

二、GM(1,1)灰色模型

GM(1,1),表示模型是一阶微分方程,且只含一个变量的灰色模型。

1、GM(1,1)模型预测方法

(1)原始数据(参考列)

(2)累加生成序列(Acumulated Generating Operator,1-AGO)

其中,显然,数列x(1)比数列x(0)光滑,弱化噪声。

(3)生成均值数列

(进一步弱化噪声)

其中

(4)建立灰微分方程

【1】

此方程白化为微分方程【2】

【1】的来源:x(1)是x(0)的累加,反过来,x(0)就是x(1)的差分,连续化就是导数,【1】中用z(1)而不直接用x(1),就是使数据进一步更光滑化,避免奇异值。可以利用曲线拟合来求a和b得值(就是利用最小二乘法)。用这个模型之前,最好绘制x(0),z(1)的散点图,看看大概情况。

(5)利用最小二乘估计【2】中参数a,b

【2】

的观测值序列为x(0)(k),k=2,3,…,n;

【2】中x(1)(t)的观测值为z(1)(k),k=2,3,…,n;

(5.1)

把x(0),z(1)作为【2】的因变量y和自变量x,作回归拟合,得到a,b的估计值;

(5.2)构造

根据最小二乘法,使得达到最小的u的估计值为求出a,b的估计值为方程【2】变为

(6)求解方程【3】,得到预测方程

为了求【3】的特解,引入初值确定特解。

【3】为一阶线性常系数微分方程,代入通解公式,代入定解条件,得

写成差分或者预测公式,为【4】k=0,1,2,…,n-1,…

2、GM(1,1)模型的预测步骤

(1)数据的检验与处理

设原始数列为计算数列的级比

(一个包含1的区间)

则认为原始数列适合作GM(1,1)灰色预测,否则需要对x(0)作平移:使得数列y(0)满足级比检验。

(2)建立模型

【4】

当初始数列或平移后的数列满足级比检验后,可用结论【4】进行预测,

(3)检验预测值

(3.1)残差检验

计算残差

这里,

则认为预测达到一般要求;

则认为预测达到较高要求;

(3.2)级差比偏差检验

首先用原始数据x(0)(k),x(0)(k-1)计算级比λ(k),再用发展系数a求出相应的级比偏差

若ρ(k)<0.2,认为预测达到一般要求;若ρ(k)<0.1,则认为预测达到较高要求。这里的a就是最小二乘估计的第一个参数,称为发展系数

(4)预测预报公式【5】

 


相关文章
|
4月前
|
编解码 算法
【论文速递】Remote Sensing2021 : 通过半全局匹配法的SAR立体图像DSM生成以及惩罚方程的评估
【论文速递】Remote Sensing2021 : 通过半全局匹配法的SAR立体图像DSM生成以及惩罚方程的评估
|
4月前
|
资源调度 前端开发 数据可视化
R语言GARCH模型对股市sp500收益率bootstrap、滚动估计预测VaR、拟合诊断和蒙特卡罗模拟可视化
R语言GARCH模型对股市sp500收益率bootstrap、滚动估计预测VaR、拟合诊断和蒙特卡罗模拟可视化
|
4月前
|
机器学习/深度学习 数据可视化
R语言用灰色模型 GM (1,1)、神经网络预测房价数据和可视化
R语言用灰色模型 GM (1,1)、神经网络预测房价数据和可视化
|
4月前
|
安全 Python
【视频】Python用GM(1,1)灰色模型预测模型对电力预测
【视频】Python用GM(1,1)灰色模型预测模型对电力预测
|
4月前
|
算法 计算机视觉 异构计算
yolo如何画框、如何变换目标检测框的颜色和粗细、如何运行detect脚本
yolo如何画框、如何变换目标检测框的颜色和粗细、如何运行detect脚本
653 2
|
数据挖掘 计算机视觉 索引
使用K均值聚类自动分割颜色
说明如何使用 K 均值聚类自动分割颜色。聚类是一种分离对象组的方法。K 均值聚类将每个对象视为在空间中有一个位置。它将对象划分为若干分区,使每个簇中的对象尽可能彼此靠近,并尽可能远离其他簇中的对象。您可以使用 函数将图像像素按值分成一个颜色空间内的若干个簇。此示例在 RGB 和 L*a*b* 颜色空间中执行图像的 k 均值聚类,以显示使用不同颜色空间如何改进分割结果。
159 0
|
4月前
GEE:获取sentinel2指定区域多个单景影像的值(样本点提取)
GEE:获取sentinel2指定区域多个单景影像的值(样本点提取)
124 0
|
4月前
|
机器学习/深度学习 自动驾驶 安全
【论文速递】Arxiv2019 - MultiPath:行为预测的多重概率锚点轨迹假设
【论文速递】Arxiv2019 - MultiPath:行为预测的多重概率锚点轨迹假设
119 0
|
人工智能 数据可视化
跟SCI学umap图| ggplot2 绘制umap图,坐标位置 ,颜色 ,大小还不是你说了算
跟SCI学umap图| ggplot2 绘制umap图,坐标位置 ,颜色 ,大小还不是你说了算
991 1