GM(1,1)灰色预测模型计算说明

提示:本指南只针对GM(1,1)模型特有的内容进行描述,关于软件的通用操作流程(如文件上传、参数设置、结果下载等),请查看:通用计算说明

1. 方法概述

GM(1,1)(Grey Model(1,1))中文一般称作GM(1,1)灰色预测模型,它是一种基于灰色系统理论的预测方法,特别适用于具有较少样本量、信息不完全确定的时间序列预测。

GM(1,1)通过对原始时间序列进行一次累加生成(1-AGO),构建灰色微分方程,求解发展系数和灰色作用量,进而对未来进行预测。该方法在经济、工程、管理等领域具有广泛应用。

2. 适用条件与数据要求

2.1 适用条件

  • 小样本预测:适合样本量在4-20个之间的时间序列预测
  • 灰色信息:适用于信息不完全、数据有限的系统
  • 指数增长规律:原始数据序列应具有一定的指数增长或衰减规律
  • 非负序列:原始数据应为非负数序列

2.2 级比检验要求

GM(1,1)模型对数据有严格的级比检验要求:

级比检验标准:

  • 级比值:$$\lambda(k)=\frac{x^{(0)}(k-1)}{x^{(0)}(k)}, \quad k=2,3, \cdots, n$$
  • 检验区间:$$(e^{-\frac{2}{n+1}}, e^{\frac{2}{n+1}})$$
  • 所有级比值都必须落在检验区间内

不满足级比检验的处理:

  • 如果级比检验不通过,系统会自动进行数据平移处理
  • 平移量 $C$:整体都加上 $C$ 之后满足级比检验要求
  • 预测完成后再减去平移量得到最终预测值

3. 数据文件格式

GM(1,1)需要上传包含时间序列数据的Excel文件。

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/gm11示例数据.png” width=“130” alt=“GM11数据格式”> </div>

文件格式说明:

  • 第一列:列名为【数据】,要求为时间序列数据
  • 只需要数值数据,不需要时间列
  • 数据类型:非负数值
  • 数据量:建议4-20个数据点
  • 数据应该具有一定的增长或变化趋势
  • ==数据必须按时间顺序排列,不能有缺失值==

4. 参数设置

GM(1,1)包含以下计算和绘图参数:

4.1 基本计算参数

  • 预测长度:设置需要预测的期数(默认10期)
  • 结果保留小数位:设置计算结果的小数位数(1-10位)

4.2 时间设置

  • 时间类型:年份或月份
  • 起始年份:历史数据的起始年份
  • 起始月份:当选择月份类型时的起始月份

4.3 绘图设置

GM(1,1)提供丰富的图表自定义选项:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/绘图设置.png” width=“1200” alt=“绘图设置”> </div>

基本设置:

  • 绘图风格:默认、经典、ggplot、bmh、深色背景、灰度
  • 字体选择:宋体、黑体、微软雅黑、楷体、仿宋
  • 字体大小:8-20号字体

坐标轴设置:

  • x轴标签:自定义x轴标签文字
  • y轴标签:自定义y轴标签文字

数据点样式:

  • 原始值散点样式:星形(*)、圆形(o)、三角形(^)、正方形(s)等
  • 预测值散点样式:与原始值相同的样式选项
  • 散点大小:1-20的大小设置

颜色和线型:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/颜色选择器.png” width=“350” alt=“颜色选择器”> </div>

  • 原始值颜色:支持颜色选择器
  • 预测值颜色:支持颜色选择器
  • 线型选择:实线(-)、虚线(–)、点划线(-·)

5. 计算结果

5.1 计算结果预览

计算结果页面可以预览到以下内容,直观看到部分结果,更多结果可以下载结果文件到本地查看。

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/计算结果.png” width=“1200” alt=“GM11计算结果”> </div>

  • 级比检验结果:包括原始值、级比值列表,级比检验区间,是否需要平移,平移量
  • 模型参数:发展系数 $a$,灰色作用量 $b$,保留10位小数,供精确验算
  • 相对误差、级比偏差、后验差比和小误差概率检验:检验指标数值和达到要求的程度描述,如达到较高要求等。

除了计算结果页面可以预览部分结果,在页面下方还可以下载结果文件,结果文件的Excel里包含以下工作表:

5.2 预测结果

主要的预测结果表格,包含:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/预测结果.png” width=“570” alt=“预测结果表”> </div>

  • 时间:根据时间类型显示的时间序列
  • 原始值:历史观测数据
  • 预测值:包括拟合值和未来预测值
  • 残差:原始值与预测值的差值
  • 相对误差(%):残差的相对百分比(取了绝对值)
  • 级比偏差:级比偏差值

5.2 级比检验

级比检验详细信息:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/级比检验.png” width=“280” alt=“级比检验结果”> </div>

  • 原始值:原始时间序列数据
  • 级比值:相邻数据的比值
  • 平移后值:如需平移处理,显示平移后的数值
  • 平移后级比值:平移后数据的级比值

5.3 建模过程量

展示GM(1,1)建模的中间计算过程:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/建模过程量.png” width=“450” alt=“建模过程量”> </div>

  • 原始序列:$X^{(0)}$ 序列
  • 平移后序列:如有平移,显示平移后序列
  • 累加序列:$X^{(1)}$ 序列(一次累加生成序列)
  • 紧邻均值序列:$Z^{(1)}$ 序列,用于构建灰色微分方程

5.4 模型参数

模型的关键参数和精度指标:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/模型参数.png” width=“440” alt=“模型参数”> </div>

  • 级比检验区间:计算得出的检验区间
  • 是否需要平移:数据预处理结果
  • 平移量:如有平移,显示平移量
  • 发展系数 $a$:灰色微分方程参数
  • 灰色作用量 $b$:灰色微分方程参数
  • 精度评价指标:平均相对误差、平均级比偏差、后验差比C、小误差概率P

5.5 预测图表

系统自动生成预测图表,直观展示:

<div style=“text-align: center;”> <img src=“/guide/figures/gm11/预测图.png” width=“570” alt=“GM11预测图表”> </div>

  • 蓝色线(原始值):历史观测数据
  • 红色线(预测值):包括拟合值和未来预测值
  • 时间轴:根据设置的时间类型显示
  • 图例:区分原始值和预测值

提醒:原始值与预测值对比曲线,可以在绘图设置页面右侧预览的图片上右键-图片另存为,也可以在结果文件的【预测图】sheet表里在图片上右键另存为,即可保存到本地使用。

6. 结果解读

6.1 精度评价指标

指标计算公式

平均相对误差: $$\omega^{(0)}(k) = \frac{|x^{(0)}(k) - \hat{x}^{(0)}(k)|}{|x^{(0)}(k)|}$$

平均级比偏差: $$\rho(k)= \left|1-\left(\frac{1-0.5 a}{1+0.5a}\right)\lambda(k)\right|, \quad k=2,3 \cdots,n$$

后验差比C: $$C = \frac{S_2}{S_1}$$

小误差概率p: $$p = P\left{\left|{\varepsilon}^{(0)}(i) - \bar{\varepsilon}^{(0)}\right| < 0.6745S_1\right}$$

相对误差类指标评价标准

指标名称 非常优秀 优秀 良好 一般
平均相对误差 < 5% 5%-10% 10%-20% > 20%
平均级比偏差 < 0.05 0.05-0.1 0.1-0.2 > 0.2

统计检验类指标评价标准

指标名称 合格 勉强 不合格
后验差比C ≤ 0.35 0.35-0.50 0.50-0.65 > 0.65
小误差概率p ≥ 0.95 0.80-0.95 0.70-0.80 < 0.70

说明:

  • 平均相对误差:小于10%为优秀
  • 平均级比偏差:小于0.1为优秀
  • 后验差比C:小于0.35为优秀
  • 小误差概率p:大于0.95为优秀

6.2 模型参数意义

发展系数 $a$: 反映系统内在的​​发展趋势或速率​​,是序列增长或衰减的固有特性的量化体现。

  • $a > 0$:系统处于衰减状态
  • $a < 0$:系统处于增长状态
  • $|a|$的大小​​ 反映了这种增长或衰减的​​速度​​。$|a|$越大,说明趋势越强劲,变化越快;$|a|$越小,则变化越平缓

灰色作用量 $b$: 代表所有​​外部因素​​(已知和未知)对系统产生的​​综合影响​​的量化,其​​具体数值的绝对值大小​​反映了外部作用的强度

6.3 预测可靠性判断

模型适用性评判:

  1. 级比检验必须通过(软件自动处理)
  2. 平均相对误差 $< 20%$
  3. 后验差比 $C < 0.65$
  4. 小误差概率 $p > 0.70$

7. 常见问题

7.1 数据要求问题

:数据量太少或太多会怎样?

:GM(1,1)最适合4-20个数据点。数据太少(<4个)会导致建模不稳定;数据太多(>20个)可能违背灰色系统的"贫信息"假设,建议考虑其他预测方法。

7.2 级比检验问题

:什么是级比检验,为什么重要?

:级比检验是GM(1,1)的前提条件,确保数据适合建立灰色微分方程。如果检验不通过,软件会自动进行平移处理。系统会在结果中明确显示是否进行了平移处理。

7.3 精度评价问题

:如何判断预测结果是否可靠?

:主要看四个指标:

  • 平均相对误差 $< 20%$(越小越好)
  • 平均级比偏差 $< 0.2$(越小越好)
  • 后验差比 $C < 0.65$(越小越好)
  • 小误差概率 $p > 0.70$(越大越好)

7.4 预测长度问题

:可以预测多长时间?

:理论上可以预测任意长度,但GM(1,1)属于短期预测方法,建议预测长度不超过原始数据量的1/2。预测时间越长,不确定性越大。

7.5 数据趋势问题

:什么样的数据适合用GM(1,1)?

:GM(1,1)最适合具有单调增长或衰减趋势的数据,对于波动剧烈或周期性明显的数据,预测效果可能不理想。如果数据波动较大,可以考虑使用GM(1,1)-马尔可夫修正模型。

8. 扩展应用

8.1 GM(1,1)-马尔可夫修正

对于预测精度要求更高的场景,平台还提供GM(1,1)-马尔可夫修正模型,通过马尔可夫链修正GM(1,1)的预测误差。

8.2 多变量灰色预测

对于多变量系统预测,可以考虑使用GM(1,N)模型(需要其他专业软件实现)。

8.3 应用领域

  • 经济预测:GDP、销售额、股价等
  • 工程预测:设备寿命、故障率、性能衰减等
  • 环境预测:污染物浓度、气候变化等
  • 人口预测:人口增长、老龄化趋势等