Matlab 之灰色预测模型GM(1,1)和GM(1,n)

1.灰色系统的定义：
灰色系统指既含有已知信息又含有未知信息的系统。
2.灰色预测模型的定义：
对灰色系统进行预测的模型。
灰色模型（Grey Model，简称GM模型）一般表达方式为GM（n,x）模型，其含义是：用n阶微分方程对x个变量建立模型。
3.灰色预测模型的目的：
通过把分散在时间轴上的离散数据看成一组连续变化的序列，采用累加和累减的方式，将灰色系统中的未知因素弱化，强化已知因素的影响程度，最后构建一个以时间为变量的连续微分方程，通过数学方法确定方程中的参数，从而实现预测目的。
4.灰色系统预测模型的特点：
无需大量数据样本，短期预测效果好，运算过程简单。
5.灰色系统预测模型的不足：
对非线性数据样本预测效果差。

常用的灰色系统预测模型主要有GM(1,1)和GM(1,n)，以下分别对这两种模型展开。
【1】.GM(1,1)模型及其matlab实现
1. GM（1,1）模型的预测原理是：对某一数据序列用累加的方式生成一组趋势明显的新数据序列，按照新的数据序列的增长趋势建立模型进行预测，然后再用累减的方法进行逆向计算，恢复原始数据序列，进而得到预测结果。
2. GM(1,1)建模过程：
(1) 设一组原始数据为，n为数据个数。对累加以便弱化随机序列的波动性和随机性,得到新的数列为：其中，
(2) 生成的邻均值等权数列其中，
(3) 根据灰色理论对建立关于t的白化形式的一阶一元微分方程GM(1,1):

其中，a，u为待解系数,分别称为发展系数和灰色作用量，a的有效区间是(-2,2)，并记a，u构成的矩阵为灰参数 ,只要求出参数a，u，就能求出,进而求出的预测值。
(4) 对累加生成数据做均值生成B与常数项向量 :

(5) 用最小二乘法求解灰参数,则
(6) 将灰参数代入，并对进行求解，得

(7) 将上述结果累减还原，即可得到预测值
(8) 利用模型进行预测:
(9) 对建立的灰色模型进行精度检验,
(9.1)残差检验：
残差：
相对误差：
(9.2)后验差检验：
均值：
方差：
残差的均值：
残差的方差：
后验差比值：
小误差概率：
(9.3) 预测精度等级对照如下：
预测精度等级
好 P>0.95 C<0.35
合格 P>0.80 C<0.45
勉强合格 P>0.70 C<0.50
不合格 P<=0.70 C>=0.65

基于Matlab实现GM(1,1)模型程序：

clear
syms a u;
c=[a,u]';%构成矩阵
A=[15 16.1 17.3 18.4 18.7 19.6 19.9 21.3 22.5];%输入数据，可以修改
Ago=cumsum(A);%原始数据一次累加,得到1-AGO序列xi(1)。
n=length(A);%原始数据个数
for k=1:(n-1)Z(k)=(Ago(k)+Ago(k+1))/2; %Z(i)为xi(1)的紧邻均值生成序列
end
Yn =A;%Yn为常数项向量
Yn(1)=[]; %从第二个数开始，即x(2),x(3)...
Yn=Yn';
E=[-Z;ones(1,n-1)]';%累加生成数据做均值
c=(E'*E)\(E'*Yn);%利用公式求出a，u
c= c';
a=c(1);%得到a的值
u=c(2);%得到u的值
F=[];
F(1)=A(1);
for k=2:(n)F(k)=(A(1)-u/a)/exp(a*(k-1))+u/a;%求出GM(1,1)模型公式
end
G=[];
G(1)=A(1);
for k=2:(n)G(k)=F(k)-F(k-1);%两者做差还原原序列，得到预测数据
end
t1=1:n;
t2=1:n;
plot(t1,A,'bo--');
hold on;
plot(t2,G,'r*-');
title('预测结果');
legend('真实值','预测值');
%后验差检验
e=A-G;
q=e/A;%相对误差
s1=var(A);
s2=var(e);
c=s2/s1;%方差比
len=length(e);
p=0;  %小误差概率
for i=1:lenif(abs(e(i))<0.6745*s1)p=p+1;end
end
p=p/len;

运行结果如下：

p=1；c=0.0148；预测等级为：好

从运行结果看，对于线性的数据使用GM(1,1)预测，其拟合效果还是不错。

【2】GM(1,n)模型及Matlab实现
1.GM（1,n）模型的预测原理：与GM(1,1)类似，不同在于输入数据变量是n个。
2. GM(1,n)模型的建模过程：
设系统有特征数据序列：
相关因素序列：

(1) 令的1-AGO序列为，其中

(2) 生成紧邻均值序列，其中

称为GM(1,n)模型。
在GM(1,n)模型中，a被称为发展系数，称为驱动系数，称为驱动项。
令和
再令，由最小二乘参数估计可得，当近似时间相应式为：
累减还原式为
差分模拟式为

基于Matlab实现GM(1,n)预测模型的程序：

A=[560823,542386,604834,591248,583031,640636,575688,689637,570790,519574,614677];
x0=[104,101.8,105.8,111.5,115.97,120.03,113.3,116.4,105.1,83.4,73.3;135.6,140.2,140.1,146.9,144,143,133.3,135.7,125.8,98.5,99.8;131.6,135.5,142.6,143.2,142.2,138.4,138.4,135,122.5,87.2,96.5;54.2,54.9,54.8,56.3,54.5,54.6,54.9,54.8,49.3,41.5,48.9];
[n,m]=size(x0);
AGO=cumsum(A);
T=1;
x1=zeros(n,m+T);for k=1:(m-1)Z(k)=(AGO(k)+AGO(k+1))/2; %Z(i)为xi(1)的紧邻均值生成序列
end
for i=1:nfor j=1:mfor k=1:jx1(i,j)=x1(i,j)+x0(i,k);%原始数据一次累加,得到xi(1)endend
end
x11=x1(:,1:m);
X=x1(:,2:m)';%截取矩阵
Yn =A;%Yn为常数项向量
Yn(1)=[]; %从第二个数开始，即x(2),x(3)...
Yn=Yn';
%Yn=A(:,2:m)';
B=[-Z',X];
C=((B'*B)\(B'*Yn))';%由公式建立GM(1,n)模型
a=C(1);
b=C(:,2:n+1);
F=[];
F(1)=A(1);
u=zeros(1,m);
for i=1:mfor j=1:nu(i)=u(i)+(b(j)*x11(j,i));end
end
for k=2:mF(k)=(A(1)-u(k-1)/a)/exp(a*(k-1))+u(k-1)/a;
end
G=[];
G(1)=A(1);
for k=2:mG(k)=F(k)-F(k-1);%两者做差还原原序列，得到预测数据
end
t1=1:m;
t2=1:m;
plot(t1,A,'bo--');
hold on;
plot(t2,G,'r*-');
title('销售预测结果');
legend('真实值','预测值');

运行结果如下：

GM(1,1)模型精度检验

Matlab 之灰色预测模型GM(1,1)和GM(1,n)相关推荐

gm11模型 MATLAB源代码,灰色预测模型 GM(1,1)- 级比检测不通过 - 平移变换常数c的确定（内含代码）...
灰色预测模型 GM(1,1) - 级比检测不通过 - 平移变换常数c的确定(内含代码) 谢谢! (哇-------没注意就5600多访问量了,激动!!!!我的初衷只是想把blog作为我的学习笔记,回首 ...
2023年数学建模国赛:灰色预测模型与MATLAB实战案例(Matlab代码)
订阅专栏后9月比赛期间会分享思路及Matlab代码灰色预测模型(Gray Forecast Model, GM)是一种基于少量.不完全的信息建立数学模型并进行预测的方法.在实际问题解决.制定发展战略 ...
利用Matlab进行灰色预测,利用matlab进行灰色预测.pdf
利用matlab进行灰色预测,灰色预测matlab程序,灰色预测模型matlab,matlab灰色预测,matlab灰色预测代码,matlab灰色预测工具箱,灰色预测的matlab程序,matlab灰 ...
python灰色预测模型步骤人口预测_超级干货：一文读懂灰色预测模型
灰色预测模型可针对数量非常少(比如仅4个),数据完整性和可靠性较低的数据序列进行有效预测,其利用微分方程来充分挖掘数据的本质,建模所需信息少,精度较高,运算简便,易于检验,也不用考虑分布规律或变化趋势 ...
灰色预测模型【GM(1,1)模型】【matlab代码】
灰色系统介绍灰色系统是由华中科技大学的邓聚龙教授于80年代初创立,该系统作为新兴的横断学科,在短短的二十年里已得到了长足的发展. 其已经成为社会,经济,科教,科技等很多领域进行预测,决策,评估,规划 ...
gm 1 n 模型matlab,灰色预测模型GM1,n模型的matlab源...
2021-02-22 15:01齐春春客户经理 function GM1_1(X0) %format long ; [m,n]=size(X0); X1=cumsum(X0); %累加 X2=[]; ...
【数学建模】灰色预测模型GM(1，1)附例题分析(MATLAB实现)
目录一.灰色预测概述 1.灰色系统.白色系统和黑色系统 2.灰色系统与其它学科的比较 3.灰色预测分类及特点二.灰色预测模型GM(1,1) 1.几个理论 1.1 累加生成数(1-AGO) 1.2 ...
代码+步骤GM（1，1）灰色预测模型-案例长江水质综合评价赛题-级比检测C的确定-matlab完整代码附送
GM(1,1)灰色预测模型-案例长江水质综合评价赛题第三题-matlab完整代码附送看到上一篇Blog在短短几天Pageviews就达到了1300多,看来大家还是比较中意建模上的笔记
推断车祸事件论文matlab,预测模型论文,关于基于优化GM(1，1)模型的水上交通事故预测相关参考文献资料-免费论文范文...
导读:本论文主要论述了预测模型论文范文相关的参考文献,对您的论文写作有参考作用. (上海海事大学商船学院中国上海 201306) [摘要]为进行水上交通事故的预测,在已有的GM(1,1)模型的基 ...

Matlab 之灰色预测模型GM(1,1)和GM(1,n)

Matlab 之灰色预测模型GM(1,1)和GM(1,n)相关推荐

最新文章

热门文章