c语言gm(1 1)模型,灰色系统模型GM(1,1)的R语言实现

1、建模

##建立灰色模型GM(1，1)对应的函数

##x表示原始数据数列，k表示数据个数

gm11

{

n

x1

for(i in 1:n) ##一次累加

{

x1[i]

}

b

m

for(j in 1:m)

{

b[j+1]

}

Yn=t(t(x[2:n])) ##构造Yn矩阵

B

B[,1]

A

a

u

x2

x2[1]

for(i in 1:k-1)

{

x2[1+i]=(x[1]-u/a)*exp(-a*i)+u/a;

}

x2=c(0,x2);

y=diff(x2); ##累减生成，获得预测数据数列

y

}

##x1原始数据数列，x2是预测数据数列

x1

x2

##检验模型精度

acc

{

n

sum1=0;

for(k in 2:n-1)

{

sum1

}

s1

sum2=0;

for(k in 2:n-1)

{

sum2

}

s2

abs1

abs2

abs12

ee

ee

}

2、应用：预测2013年和2014年全国规模以上快递企业收入

##x数列是2008年-2012年全国规模以上快递企业收入数据(来源：国家邮政局，单位：亿元)

> x

> gm11(x,7)

[1] 408.4000 443.1355 585.3243 773.1370 1021.2131 1348.8894 1781.7069

> x1

> x2

> acc(x1,x2)

[1] 0.9851449

预测结果：预计2013、2014年全国规模以上快递企业收入分别为1349亿元、1782亿元

灰色绝对关联度为0.9851449，即关联度为一级，预测精度为优