精品灰色系统GM11模型适用范围拓广.docx
《精品灰色系统GM11模型适用范围拓广.docx》由会员分享,可在线阅读,更多相关《精品灰色系统GM11模型适用范围拓广.docx(29页珍藏版)》请在冰点文库上搜索。
精品灰色系统GM11模型适用范围拓广
灰色系统GM(1-1)模型适用范围拓广.
1999年1月系统工程理论与实践第1期
灰色系统GM(1,1)模型适用范围拓广α
李希灿
(山东水利专科学校,山东泰安271000)
摘要 研究了灰色系统GM(1,1)模型在建模过程中由于原始数列乘以不等于零的常数对模型值及预测值的影响,得出GM(1,1)模型完全适用于负数据序列建模的结论.
关键词 灰色系统 模型 灰色参数
WideningofSuitableLimitsofGrey
SystemGM(1,1)Model
LiXican
(ShandongHydraulicEngineeringCollege,Taian271000)
Abstract Inthispaper,westudythefactthataconstantwhichmultipliesalldatainthecoarseserieswouldinfluencethevaluesofmodelandprediction.TheresultthatGM(1,1)modelissuitabletonegativedatasequenceisobtained.
Keywords greysystem;model;greyparameters
1 引言
设有时间数据序列X(0)
X(0)={x(t)t=1,2,…,n}
={x(0)
(1),x(0)
(2),…,x(0)(n)}
t
对X(0)作一次累加生成(1-AGO),令x
(1)(t)=6x(0)(k),得生成数据序列X
(1)
k=1
X
(1)={x
(1)(t)t=1,2,…,n}
={x
(1)
(1),x
(1)
(2),…,x
(1)(n)}
n
=x(0)
(1),62
x(0)(k),…,x(0)(k)
k=16k=1
利用序列X
(1)可建立如下白化方程
(1)
+
(1)
dtaX=u
式中,a,u为灰色参数.按最小二乘法求解
(a,u)T=(BTB)-1BTYN
α本文于1997年7月15日收到
©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
(1)
(2)(3)(4)
98系统工程理论与实践
-B=
1999年1月
(x
(1)
(2)+x
(1)
(1))2
(x
(1)(3)+x
(1)
(2))2
11
-
-YN=(x
(0)
1
(x2
(1)
(n)+x
(1)(n-1))
(2),x(0)(3),…,x(0)(n))T
a
a
求出a,u后,解(3)式得微分方程
x
δ
(1)(k)=
x
(0)
(1)-e-
a(k-1)
+(5)
对xδ
(1)作一次累减生成,即得xδ(0)序列
δ()δ()δ()
x0(k)=x1(k)-x1(k-1)
x
δ(0)(k)=
x
(0)
(1)-
a
(1-ea)e-
a(k-1)
(6)
因此,给定原始数据序列
(1)式,由式
(2),(3),(4),(5),(6)即可建立GM(1,1)预测模型.但GM(1,1)建模时一般要求
(1)式必须为“非负”数列.随着灰色理论研究的不断发展,GM(1,1)模型应用越来越广泛,如变形观测中,利用监测网的多期观测数据建立GM(1,1)模型进行变形预测等,当时间数据序列为高差时,就可能为负数据序列.那么负数列是否能直接用于建模呢?
为此本文加以讨论.首先导出灰色参数的显式表达式,由此对原始数据序列乘以不等于零的常数对预测结果及精度影响加以讨论,得出GM(1,1)建模完全适用于负数据序列的结论.
2 原始数据序列乘以不等于零常数对GM(1,1)模型参数及预测值的影响
设原始数据序列
(1)式乘以常数K0,K0≠0,生成新的数据序列Y(0)
Y
(0)
={K0x
(0)
(1),K0x(0)
(2),…,K0x(0)(n)}(7)
Χ{y(0)
(1),y(0)
(2),…,y(0)(n)}
对数列Y(0)类似式
(2)~(5)建模
y
δ
(1)(k)=
y
(0)
(1)-
-e
a1
a1(k-1)
+
a1
(8)
式中a1,u1为灰色参数,即
-1T
(a1,u1)T=(BTB1Y11B1)
(9)
其中
-B1=
(y
(1)
(2)+y
(1)
(1))2
(y
(1)(3)+y
(1)
(2))2
11
-
-(y
(1)(n)+y
(1)(n-1))2
1
δ
(1)(k)=K0xδ
(1)(k).首先有下面的命题成立.可以证明,y
2T
命题1 B
T1B1=K0BB
()()()
Y1=(y0
(2),y0(3),…,y0(n))T
证 记
bi=-
[x
(1)(i)+x
(1)(i+1)]2
©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第1期
灰色系统GM(1,1)模型适用范围拓广=-
6
i
x
(0)
(j)-
(0)
(i+1) i=1,2,…,n-1
j=1
2
x则
b1
1n-1
1
bbi
2
bi
BT
B=
1
b2…bn-1
2
1
i=1
6
n-i=1
11
…
1
=
b
=
6
n-1
bbin-
1
n-1
6
i=1
其中
6
n-1
n-1i
bi=-
x
(0)
(j)+
(0i=1
6
i=16
j=1
2
x)
(i+1)n
=-(n-1)x
(0)
(1)-
6
n-i+
(0)
i=2
2
x
(i)6
n-1
n-1
i
2
bi2
=
x
(0)
(j)+
2
x(0)
(i+1)i=1
6
i=1
6
j=1
n
=(n-1)(x
(0)
(1))2+2x(0)
(1)
6
n-i+
0)
)
i=2
2
x
((in
+
6
n-
i+
(x(0)(i))2+
)x(0)(i)x(0)(j)
i=2
4
2≤6
(2n-2j+1i由(10)式得
n-1
n-1
n-1
BT
B=(n-1)6bi2-b2
i
=(n-1)
(bi-λb)2
i=1
6
i=1
6
i=1
其中λb=
n-16
n-1
bi.所以当b1,b2,…,bn-1不全相等时,BTB>0.
i=1
由(13)式及(11),(12)
式,经化简整理得
n
BT
B=
6
(n+2-i)i-
n-(0)i))2
i=2
42
(x(+
6
(2n-2j+1)i-3n+3j-(i)x(0)(j)i2
x
(0)
2≤同理
n
B1T
B1=
6
(n+2-i)i-
n(y(0)(i))2
i=2
4-2
+
6
(2n-2j+1)i-3n+3j-(i)y(0)(j)2≤i2
y
(0)
由(7),(14),(15)式得BT1B1=
K20BTB
.证毕.
命题2 参数a1=a.
证 因为
x
(0)
(i)=x
(1)(i)-x
(1)(i-1),i=2,3,…,n
由(4),(10)式,
a=(BTB)-1BTY
n-1-
=
6n-1
bix
(1)
(2)-x
(1)
(1)i=1
b1b2…bn-1
(1)
BT
B
n-1
b11
x(3)-
x
(1)
(2)
-
6
i
6
n-1b2…
1
i
x
(1)
(n)-x
(1)(n-1)
i=1
i=1
©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
99
(9)
(10)
(11)
(12)
(13)
(14)
(15)
100系统工程理论与实践1999年1月
=
BB
T
n-1bi
-
66
n-1i=1
n-1
bib2i
-
i=1
[(x
(1)(n))2-2
(x
(1)
(1))2]
(16)
-
6
n-1
i=1
6
(i)(i)
2
n
x
(0)
(i)
i=2
由式(11),
a=
BB
T
T
--(i)
n
266
n
n
xx
(0)
--
(x
(0)
(1))
2
-
i=1n
6b6
i
i=1
n-1n
x
(0)
(i)
i=2
=
BB
n
(0)
2
2(n-1)xx
(0)
(x(0)
(1))2+i+
i=1
6x
i=2
(0)(0)
(1)+6
n
(n-(0)
i=2
)x(0)(i)2
(17)
=
B
T
B6
(i)
i=2
6
i=2
2
n
-i+1x(i)
同理
a1=
B1B1
T
6
y
(0)
(i)
i=2
6
n
i=2
2
-i+1y
(0)
(i)(18)
由(7)式及命题1,(17),(18)式得a1=a.证毕.
命题3 u1=K0u
证 由(16)式及(11),(12)式
u=
BB
T
T
66
n-
n
xx
(0)
(i)(i)4
i=2n
6
n-1
bi+
2
i=1
26
n-1
bi
2
i=1
6
n
x
(0)
(i)
(1)
2
-n
(x(0)
(1))2n
-i+
i=1
=
BB
(0)
(n-1)(x
(0)
(1))+2x
(0)
i=2
6
i=2
2
x
(0)
(i)
+ -
6
n
i+
(x(0)(i))2+x
(0)
i=22≤i6
(2n-2j+1)x(0)(i)x(0)(j)(i)
2
2
2x
(0)
(1)
6
n
n
(i)+
i=2
6
n
n
x
(0)
×(n-1)x
(0)
(1)+
i=2
6
n
n-i+
i=2
2
x
(0)
(i)
经整理化简后,
u=ax
(0)
(1)+
6
x
(0)T
(i)
BB
n
6
i=2
2
(x(0)(i))2+
2≤i6
n+
2
-
j+22
x
(0)
(i)x(0)(j)(19)
同理
u1=a1y
(0)
(1)+
6
y
T
(0)
(i)
B1B1
6
n
i=2
2
(y(0)(i))2+
2≤i6
n+
2
-
j+22
y
(0)
(i)y(0)(j)
(20)
由(7),(19),(20
)式及命题1,2得u1=K0u.证毕.
由(5),(8),(7)式及命题1,2,3可得
δ()δ()
y1(k)=K0x1(k)
(21)
证华.即常数乘以原始数列后的模型计算值等于常数与原始数列建模的模型计算值的乘积.
δ
(1)(k)缩小K0倍就等于以原由(21)式可知,Y(0)序列建模的模型值及精度是序列x(0)建模的K0倍,y
始序列X(0)建模的模型计算值.(21)式的实践意义是:
对于给定的时间数据序列建立GM(1,1)模型,采用不同单位建模,模型精度是相同的.
δ
(1)(k)=-xδ
(1)(k)推论 当K0=-1时,y
©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
(22)
第1期灰色系统GM(1,1)
模型适用范围拓广101
δ
(1)(k)亦为xδ
(1)(k)的相反数.因此,对由(22)式,当X(0)为正数列时,Y(0)为对应的负数列,模型计算值y
于负数据序列(序列中数据全为负数)也可以直接用于建立GM(1,1)模型.
3 应用
“非负”数据序列拓广到负数据序列,具有重要的理论与实践意义,主要表GM(1,1)模型的建模条件从现在两方面.
1)用于负数据序列建模
设给定的负数据序列Y(0)={y(0)
(1),y(0)
(2),…,y(0)(n)},正数据序列X(0)={x(0)
(1),x(0)
(2),…,(0)(0)(0)(0)
x(n)},满足y(k)=K0x(k),且K0=-1,k=1,2,…,n.则Y建模有两种方法:
δ
(1)(k)①由式(7),(8),(9)得y
δ
(1)(k),即yδ
(1)(k)=K0xδ
(1)(k)②由
(1)~(5),(21)式得y
显然,GM(1,1)建模条件拓广到负数据序列,为负数据序列直接建模提供了理论基础.
2)用于GM(1,1)模型优化
GM(1,1)模型优化就是对于给定的数据序列,寻找模型精度最优的灰色参数a,u.文献[1]证明了原始数据序列中每一数据增加同一常数a0对GM(1,1)模型值及预测值均有影响,取不同的a0进行模拟计算,
寻找模型精度最好的a0,就得到优化的GM(1,1)模型.由于GM(1,1)建模时,原来要求原始数据序列必须为非负的,所以原始数列加常数a0后仍必须满足非负条件.这就极大地限制了优化的范围.GM(1,1)建模条件由正数据序列拓广到负数据序列,就克服了这种局限性.下面以实例加以说明.
已知原始数据序列X由
(2)~(5)式
(0)
Xx
(0)
={230,226,220,218,214,208,201}
0.02189119(k-1)
δ
(1)(k)=-10455.72e-令残差q(0)(k)=x(0)(k)-xδ(0)(k),得残差序列
q
(0)
+10685.72
显然原点误差相对来看较大,为此进行模型优化,取后验误差C,绝对平均误差ϖ原点误差E0、残差E、方差Ρ、关联度G作为衡量模型精度的指标,取不同的a0模拟计算,各项指标列于表1.
表1 GM(1,1)模型优化指标表
a0
={0,-0.40,-1.50,1.30,1.99,0.58,-1.93}
数列类型正正正正正负负负
C
ϖE
1.04471.04461.09951.16451.81040.57580.81680.9944
E0Ρ
1.25561.25551.31521.39262.53980.92221.03361.2030
G
10005000-100-200-233-300-1000
0.13500.13500.14140.14970.27240.09920.11110.1293
-1.7812-1.7734-1.9287-2.123-5.24640.0297-1.2329-1.6562
0.55460.55460.55020.55030.68350.73390.60480.5587
从表1可见,a0=-233时模型精度最高.故a0=-233时得到GM(1,1)优化模型.而a0=-233时数
列为负数据,因此GM(1,1)建模条件未拓广时无法得最优化的模型.最后确定GM(1,1)模型为
δ()()
x1(k)=-30.00606e0.257818k-1+260.00606
(下转第105页)
©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第1期就业乘数新算法探讨表3 新旧方法计算结果对比
旧方法
一、各部门就业乘数 11农业 21工业 31建筑业 41运输邮电业 51商业饮食业 61其他服务部门 71居民收入 平均
二、平均就业增加量
1.14
L(I-A)-1FC
L
3
105
新方法
1
L(I-A)-
L
3
(I-A
3
)-
1
21350.870.870.740.861.12
3.991.952.051.871.992.342.112.18(I-A
3
)-1FC
0.9122.07
其中L(I-A)
-1表示各部门增加单位最终需求(如出口)等,直接和间接所增加的劳动力需求量(就业乘数);
L(I-A)-1FC表示固定资本增加一个单位直接和间接所引起的劳动力需求量,可利用它计算增加公共投资而引起
的新增劳动力就业人数.
参考文献
1 国家统计局国民经济核算司11995年度中国投入产出表1北京:
中国统计出版社,1997:
1~302 国家统计局编1中国统计年鉴19971北京:
中国统计出版社,1997:
96~973 陈锡康1现代科学管理方法基础1北京:
科学出版社4 陈锡康1中国城乡经济占用产出分析1北京:
科学出版社
(上接第101页)
4 结语
经上述讨论可知,灰色系统GM(1,1)建模条件可以拓广到负数据序列.这一拓广为负数据序列直接建
模和GM(1,1)模型优化提供了理论基础.算例说明有助于提高GM(1,1)模型精度,得到最优模型.
参考文献
1 李留藏1灰色系统GM(1,1)模型的讨论1数学的实践与认识,1993
(1):
15~222 邓聚龙1灰色系统基本方法1武汉:
华中理工大学出版社,1992
©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.