灰色系统建模 灰色系统理论在建模中的应用:灰色系统理论在建模中被广泛用来处理数据。与插值拟合相比，利用灰色模型处理数据不仅对数据没有很强的限制，而且精度更高，计算更简便。常用的灰色系统生成方式有: 累加生成,累减生成,均值生成,级比生成等,下面对这几种生成做简单介绍. 累加生成:
(0)(0)(0)(0)(0)
(1)(1)(1)(1)(1)(1)(0)
[(1),(2),,()],,[(1),(2),,()],:xxxxxnxxxxxnxx令为原始序列,记生成数为如果与之间满足如下关系
(1)(0)1()();1,2,,(21)kixkxikn
,1AGO一次累加生成则称为记为 累加生成在灰色系统理论中有着非常重要的地位,它能使任意非负数列,摆动的或非摆动的,转化为非减的的,递增的数列. 累减生成:累减生成,即对数列求相邻两数据的差,累减生成是累加生成的逆运算,常简记为IAGO(Inver se Accumulated Generating Operation), 累减生成可将累加生成还原为非生成数列,在建模过程中用来获得增量信息,其运算符号为∆.
()()(),,:rr
ixrxi令为次生成数列对作次累减生成记为
其基本关系式为
(0)()()(1)()(0)()(0)()(2)()(1)()(1)()()()(1)()(1)()[()]()[()][()][(1)][()][()][(1)](25)[()][()][(1)]rrrrrrrr
iririrxkxkxkxkxkxkxkxkxkxkxk
(0)(1)(),(0)0,;(0)1(0)11.ikkikki式中为次累减即无累减为1次累减,即与时刻两个零次累减量求差,为次累减,即与时刻两个次累减量求差
(25):从式还可得到以下关系 (1)()(0)()(0)()()()1(1)(1)11(1)[()][()][(1)]()(1)(26)()()()rrrrrkkrriirxkxkxkxkxkxixixk
(2)()(1)()(1)()(1)(1)1(2)(2)11(2)[()][()][(1)]()(1)(27)()()()rrrrrkkrriirxkxkxkxkxkxixixk
:同理可得 ()()()[()]()(28)irrixkxk
()()(0)[()]()(29)rrxkxk (29),,.,,.:1,rrr从式可以看出对次生成数列作次累减即还原为非生成数列事实上累加中包含着累减累减中包含着累加比如时有 1(1)(0)(0)(0)11(1)(0)()()()()(1)()(210)kkiixkxixixkxkxk
(0)(1)(1)()()(1)xkxkxk
进一步有(1)()()()()(1)(211)rrrxkxkxk
.均值生成分为邻均值生成与非邻均值生成两种
级比生成是一种常用的填补序列端点空穴的方法.对数列端点值的生成,我们无法采用均值生成填补空缺,只能采用级比生级比生成.成是级比级比生(k成在建模中可以获得较好的灰)与光滑比(k)生成指数律.的总称. (0)(0)(0)(0)[(1),(2),,()],(),(),XxxxnKk设序列为原始序列
称为级比为光滑比其表达式为 (0)(0)(0)(1)()()/(1)()()/(1)(212)kxkxkkxkxk
(0)(0)(0)(0)(0)(0)(0)
[(1),(2),,(1),()],(1)(1),()(),(1)()Xxxnnxnxnxxn设为端点是空穴的序列若用右邻的级比生成用的左邻级比生成则称和为级比生成 GM(1.1)模型建模机理 GM(1.1)原理步骤 原始数列：
(0)(0)(0)(0)(1),(2),,()Xxxxn
对(0)X进行一次累加，得到新数列： (1)(1)(1)(1)(1),(2),,()Xxxxn
(1)0,()()kixkxi其中
于是(0)()xk的GM（1.1）白化形式的微分方程为： (1)(1)(216)dxaxudt
其中，a,u为待定系数，将（2-16）式离散化，即得： (1)(1)(1)((1))((1))(217)xkazxku
其中，(1)(1)((1))xk为(1)x在（k+1）时刻的背景值 因为：(1)(1)(1)(1)(0)((1))(1)()(1)(218)xkxkxkxk (1)(1)(1)1(1)((1)())(219)2zkxkxk
将（2-18），（2-19）式代入（2-17）式，得 (0)(1)(1)1(1)[(()(1))](220)2xkaxkxku
将（2-20） (1)(1)(0)
(1)(1)(0)
(0)(1)(1)
1((1)(2))12(2)1((2)(3))1(3)(221)2()1((1)())12xxxxxxxnxnxn
(1)(1)(0)
(1)(1)(0)
(0)(1)(1)
1((1)(2))12(2)1((2)(3))1(3),2()1((1)())12xxxxxxYBxnxnxn
令
[]Tau为待辨识参数向量，则（2-21）可写成：
(222)YB 参数向量可用最小二乘法求取，即 1ˆˆˆ[,]()TTTauBBBY
把求取的参数代入（2-16）式，并求出其离散解
ˆ(1)(1)
ˆˆˆ(1)[(1)](224)
ˆˆakuuxkxeaa
(0)(1)(1)ˆˆ(1)
ˆˆˆ(1)(1)()
ˆ(1)[(1)](225)ˆaakxkxkxkuexea
还原到原始数据得 (224),(225)(1.1),(1.1).GMGM式称为模型的时间相应函数模型它是模型灰色预测的具体计算公式
（GM1.1）模型的精度检验 模型选定之后,一定要经过检验才能判定其是否合理,只有通过检验的模型才能用来作预测,灰色模型的精度检验一般有三种方法:相对误差大小检验法,关联度检验法和后验差检验法.下面对这三种方法做个简单介绍. 1 级比检验： 为了保证建模方法的可行性，需要对已知数据列做必要的检验处理。 原始数列：
(0)(0)(0)(0)(1),(2),,()Xxxxn
级比表达式为： 2 后验差检验法 (0)(0)2212
ˆ(1.1)(231),(232),,GMXXESS设按建模法所求出的如所示残
差如所示原始序列及残差序列的方差分别为和则
2(0)21122211[()]1[()](235)nknkSxkxnSeken
(0)1111,(),()nnkkxxkeeknn其中
计算后验差比为： 21/(236)CSS 计算小误差概率： 1()0.6745(237)pPekeS
1212
,.CpCCSSSSC指标和是后验差检验的两个重要指标.指标越小越好越小表示大而越小大表示原始数据方差
大,即原始数据离散程度大.小表示残方差小,即残差离散程度小.小就表明尽管原始数据很离散,而模型所得计算值与实际值之差并不太离散.
1,,0.6745,,,.ppCp指标越大越好越大表明残差与残差平均值之差小于给定值的点较多即拟合值(或预测值)分布比较均匀.按两个指标可综合评定预测模型的精度模型的精度由后验差和小误差概率共同刻划.一般地,将模型的精度分为四级,见表2-1
21表精度检验等级参照表 模型精度等级 均方差比值C 小误差概率p 1级（好） C<=0.35 0.95<=p 2级（合格） 0.353级（勉强） 0.54级（不合格） 0.65
,,MaxpC模型的精度级别的级别于是的级别
3 序列光滑度的理论分析
*[(1),(2),,()],:[(1),(2),,()],,()0.5()0.5(1),1,2,,;,(1),:3.1XxxxnZXZzzznzkzkzkknXxnX设序列是的均值生成序列其中是某一可导函数的代表序列我们将删去后所得的序列仍记为定若满足义
11**11(1),()();(2)max()()max()()kiknknkxkxixkxkxkzk
当充分大时
,(1),(2).X则称为光滑序列称为序列光滑条件