助力国赛 | 第4弹 灰色预测
助力国赛 | 第4弹 灰色预测
https://pic3.zhimg.com/v2-c5dce68f8f697c33dc143d7920e90eba_b.jpg
前言上次介绍了规划问题。规划问题是一类比较重要的问题，在数学建模中出现的比率也比较高，希望各位能多去关注一下。本次将会介绍另一类问题——预测问题。在数学建模中遇到预测问题也是比较常见的事，甚至在有些赛题中还占了主导地位。灰色模型（Gray Model，又称灰色理论）在处理数据预测方面有着很大的优势，也是在数学建模上被广泛使用的一个模型，效果往往很不错，是数模参赛选手必须要掌握的模型之一。灰色系统，是由国内科学家自主提出来的一个理论，经过快40年的发展，已经有严格的理论基础了，在海内外影响十分的大。其最大的有点就是实用。用灰色模型预测的结果比较稳定，不仅仅适用于大数据量的预测，在数据量较小时（只要数据多于3个即可）预测结果依然较准确。如果各位对灰色模型有着强烈的兴趣，可以去阅读相关方面的书籍，本文只是简单介绍下灰色模型的基本概念和GM（1，1）模型的应用。
[*]模型简介
[*]GM(1,1)操作步骤
[*]实际应用
模型简介灰色系统，是由华中科技大学控制科学与工程系教授，博士生导师邓聚龙于1982年提出的，在短短的三十多年的时间里，得到了飞速的发展，应用已渗透到自然科学和社会经济等许多领域，显示出了强大生命力，具有广阔的发展前景，在海内外影响十分的深远。一些大学还开设了相关研究机构，并且招收研究生和博士生，已成为系统工程和控制领域一个比较热门的研究方向。国内在这方面有着比较深入研究的机构有南航，我想主要还是邓聚龙老先生的学生刘思峰在那里。
https://pic1.zhimg.com/80/v2-e2a451a463e0d62894ad96ff97151e9c_720w.jpg灰色系统理论认为，系统的行为现象尽管是朦胧的，不完全透明的，数据是复杂的，但是是有序的，是具有整体功能的。灰色系统的灰就是这种状态，“朦朦胧胧”，“似乎明朗但又不完全明朗”。基于数据是有序的这一点假设，需要我们对数据进行一个预处理，我们把处理后的数据序列成为生成列。进行处理的目的不是寻求它的统计规律和概率分布，而是将其变为有规律的时间序列数据，一句话概括就是——以数找数的规律，建立动态模型。灰色系统一般有两种数据处理的方式，累加和累减，一般采用累加方法。那么灰色预测就通过鉴别系统因素之间元素发展趋势的相异程度，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列，然后建立相应的微分方程，从而来预测事物的未来发展趋势的目的。灰色预测的数据就是通过生成数据的模型所得到的预测值的逆处理结果。灰色预测是以灰色模型为基础的，其中所建立的微分方程模型有很多种比如说GM（2,1），GM(1,1)等，下面主要介绍GM(1,1)。GM(1,1)操作步骤在诸多的灰色模型中，以灰色系统中单序列一阶线性微分方程模型GM(1,1)模型最为常见，其具体操作步骤如下：设有原始数据列https://pic1.zhimg.com/80/v2-60d81e564b5c296821334634407c7168_720w.jpg,https://pic1.zhimg.com/80/v2-34747986abb37f9fab9d944d1c1d1ba4_720w.jpg为数据个数。（1）数据预处理原始数据累加以便弱化随机序列的波动性和随机性，得到新数据序列：
https://pic1.zhimg.com/80/v2-7984870b6655862385e9cf5e7a222880_720w.jpg
其中https://pic4.zhimg.com/80/v2-b2bbb3c89b5bcfc3d886d1ced9b069d3_720w.jpg，https://pic4.zhimg.com/80/v2-b2bbb3c89b5bcfc3d886d1ced9b069d3_720w.jpg中各数据表示对应前几项数据的累加。https://pic3.zhimg.com/80/v2-d023270ec9e208b16423159ccd50a26a_720w.jpg（2）建立微分方程对https://pic4.zhimg.com/80/v2-b2bbb3c89b5bcfc3d886d1ced9b069d3_720w.jpg建立的一阶线性微分方程：
https://pic4.zhimg.com/80/v2-ffcd9212ee1635c85c81580b15f9c243_720w.jpg
其中，https://pic2.zhimg.com/80/v2-98bdafc648c7200b14de39c7d2c64195_720w.jpg、https://pic2.zhimg.com/80/v2-6c959e3d8450971fc6ed86e7d6a51c85_720w.jpg为待定系数，分别称为发展系数和灰色作用量；https://pic2.zhimg.com/80/v2-98bdafc648c7200b14de39c7d2c64195_720w.jpg的有效区间是https://pic1.zhimg.com/80/v2-c40d71ab0ed7c9b0e0cb722310f81c38_720w.jpg；并记https://pic2.zhimg.com/80/v2-98bdafc648c7200b14de39c7d2c64195_720w.jpg、https://pic2.zhimg.com/80/v2-6c959e3d8450971fc6ed86e7d6a51c85_720w.jpg构成的矩阵为https://pic3.zhimg.com/80/v2-19a673d610d354ae88812d3fda6586b2_720w.jpg。所以只要确定了参数https://pic2.zhimg.com/80/v2-98bdafc648c7200b14de39c7d2c64195_720w.jpg、https://pic2.zhimg.com/80/v2-6c959e3d8450971fc6ed86e7d6a51c85_720w.jpg，就可以求出https://pic4.zhimg.com/80/v2-b2bbb3c89b5bcfc3d886d1ced9b069d3_720w.jpg，进而求出https://pic4.zhimg.com/80/v2-215a285bf2d35ed971a7fc2bf480b8db_720w.jpg的未来预测值。
（3）确定参数先对累加生成数据做均值生成https://pic1.zhimg.com/80/v2-a00699ac32a0a0269b12e2167ed24744_720w.jpg与常数项向量https://pic3.zhimg.com/80/v2-5dc9aeb14cee542cc5f579365df0738a_720w.jpg，即
https://pic1.zhimg.com/80/v2-b66d2040add59ba6524ac5b611e7ba44_720w.jpg
然后用最小二乘法求解参数https://pic1.zhimg.com/80/v2-b22b7214281ee5dd95eda29acc3a3394_720w.jpg则，https://pic2.zhimg.com/80/v2-6f1cd887d9c5e52d9252df5a63ffb279_720w.jpg
将灰参数https://pic1.zhimg.com/80/v2-b22b7214281ee5dd95eda29acc3a3394_720w.jpg待人https://pic2.zhimg.com/80/v2-c0d5a0e75564cc7424a5a57da8ef57dd_720w.jpg，并对https://pic2.zhimg.com/80/v2-c0d5a0e75564cc7424a5a57da8ef57dd_720w.jpg进行求解，得：https://pic4.zhimg.com/80/v2-d2c7efdad803d5f1750ff11d0836b3b7_720w.jpg
由于https://pic1.zhimg.com/80/v2-b22b7214281ee5dd95eda29acc3a3394_720w.jpg是通过最小二乘法求出的近似值，所以https://pic3.zhimg.com/80/v2-c3c2cd03a179fe8de9a68dc68a0f8cf6_720w.jpg是一个近似表达式。最后对函数表达式https://pic3.zhimg.com/80/v2-c3c2cd03a179fe8de9a68dc68a0f8cf6_720w.jpg及https://pic3.zhimg.com/80/v2-fbb06f53c9cfedfd50a6dc6f287a9382_720w.jpg进行离散，并将二者做差以便还原https://pic4.zhimg.com/80/v2-215a285bf2d35ed971a7fc2bf480b8db_720w.jpg原序列，得到近似数据序列https://pic3.zhimg.com/80/v2-095f6318df8fcd0d642516e73cc14f36_720w.jpg如下：https://pic3.zhimg.com/80/v2-e1328d3d63dd2beb3121d6eca6dfc89a_720w.jpg（4）模型检验计算https://pic4.zhimg.com/80/v2-215a285bf2d35ed971a7fc2bf480b8db_720w.jpg与https://pic3.zhimg.com/80/v2-a270d6ecce0f211c0ed883312e170c26_720w.jpg之间的残差https://pic2.zhimg.com/80/v2-a65031c955dad334fdb88199d9fb4211_720w.jpg和相对误差https://pic4.zhimg.com/80/v2-3c8203cedd353af8313b1327876c38e7_720w.jpg：
https://pic4.zhimg.com/80/v2-d9146b9ce5e090b69c66f34ee15d8793_720w.jpg
求原始数据https://pic4.zhimg.com/80/v2-215a285bf2d35ed971a7fc2bf480b8db_720w.jpg的均值以及方差https://pic3.zhimg.com/80/v2-064f3b412f567c1b25fc3c2f8f23a35e_720w.jpg。求https://pic2.zhimg.com/80/v2-a65031c955dad334fdb88199d9fb4211_720w.jpg的平均值https://pic2.zhimg.com/80/v2-35fdb6a6f0c2c55aed89f0982726ff4d_720w.jpg以及残差的方差https://pic3.zhimg.com/80/v2-51937d845831e640e3ad23aa092fb2f6_720w.jpg。计算方差https://pic4.zhimg.com/80/v2-047e69831f95f36eb42a2b2fe82f2127_720w.jpg求小误差概率https://pic2.zhimg.com/80/v2-38b7403a79d0628f3c547f0a1177eba5_720w.jpg
灰色模型精度检验如下表https://pic1.zhimg.com/80/v2-c0083d21600845473e4786cdaf047b34_720w.jpg
需要注意的是在实际应用过程中，检验模型的精度的方法并不唯一。可以利用上述方法进行模型的检验，也可以根据https://pic1.zhimg.com/80/v2-e46a3c9da78040c8b5d5365c83881d80_720w.jpg的误差百分比并结合预测数据与实际数据之间的测试结果来确定
（5）确定参数利用模型进行预测：
https://pic2.zhimg.com/80/v2-6efc316745b00c31c647007e1ab7ff61_720w.jpg实际应用下面通过例子来介绍如何使用GM（1，1）模型。
已知一公司1999-2008年的利润为（单位：元/年）[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670]
,请预测该公司未来几年的利润。code:clearsyms a b;c=';A=;B=cumsum(A); % 原始数据累加n=length(A);for i=1:(n-1) C(i)=(B(i)+B(i+1))/2; % 生成累加矩阵end% 计算待定参数的值D=A;D(1)=[];D=D';E=[-C;ones(1,n-1)];c=inv(E*E')*E*D;c=c';a=c(1);b=c(2);% 预测后续数据F=[];F(1)=A(1);for i=2:(n+10) F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a ;endG=[];G(1)=A(1);for i=2:(n+10) G(i)=F(i)-F(i-1); %得到预测出来的数据end t1=1999:2008;t2=1999:2018;Gplot(t1,A,'o',t2,G) %原始数据与预测数据的比较legend('原始数据','预测数据');xlabel('年份')ylabel('利润')
效果：
https://pic2.zhimg.com/80/v2-5b543ea5cc07df15099527a01c3dd0cd_720w.jpg
如果用GM（1,1）来拟合数据会怎样？我们的想法是用一个确定的函数值来生成数据，并让其预测接下来的，然后两者比较一下看看效果。我们采用的函数是https://pic4.zhimg.com/80/v2-ae51ae7555eeab5016c463e0a31ec9a7_720w.jpg
code：
clearsyms a b;c=';x = 1:0.5:10;A= 1/3*x.^2 + log(x+1);B=cumsum(A); % 原始数据累加n=length(A);for i=1:(n-1) C(i)=(B(i)+B(i+1))/2; % 生成累加矩阵end% 计算待定参数的值D=A;D(1)=[];D=D';E=[-C;ones(1,n-1)];c=inv(E*E')*E*D;c=c';a=c(1);b=c(2);% 预测后续数据F=[];F(1)=A(1);for i=2:(n+4) F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a ;endG=[];G(1)=A(1);for i=2:(n+4) G(i)=F(i)-F(i-1); %得到预测出来的数据end t1=1:0.5:10;t2=1:0.5:12;t3 = 1:0.5:12;sss = 1/3*t3.^2 + log(t3+1);Gplot(t1,A,'o',t2,G,t3,sss) %原始数据与预测数据的比较legend('原始数据','预测数据','原始函数');xlabel('年份')ylabel('利润')
效果：
https://pic1.zhimg.com/80/v2-f920b5638f1fea5029261cea5d9a48a4_720w.jpg
随着自变量的增大，两者相差还是蛮大的。
页:
[1]