灰色预测模型及应用论文

灰色系统理论的研究摘要

科学地预测尚未发生的事物是预测的根本目的和任务。无论个体还是组织，在制定和规划面向未来的策略过程中，预测都是必不可少的重要环节，它是科学决策的重要前提。在多的预测方法中，灰色预测模型自开创以来一直深受许多学者的重视，它建模不需要太多的样本不要求样本有较好的分布规律，计算量少而且有较强的适应性，灰色模型广泛运用于各种领域并得了辉煌的成就。本文详细推导GM（1,1）模型，另外对灰关联进行了进一步的改进，让改进的计算式具有唯一性和规范性。过给出的实例高校传染病发病率情况，建立了预模型，并预测了1993年传染病发病。另外对传染病发病率较高的痢疾、肝炎、疟疾三种疾病做了联度分析，发现痢疾与整个传染病关系最密切，而肝炎、疟疾与整个传染病的密切程度依次差些关键词灰色预测模型；灰关联度；灰色系统理论word档可自由复制编辑灰色系统理论研究GM(1,1)预测与关联度的拓展1引言模型按照对研究对象的了解程度可分为：黑箱模型、白箱模型、灰箱模型。黑箱模型：信息缺乏，暗，混沌。白箱模型：信息完全，明朗，纯净。灰箱模型：信息不完全，若明若暗，多种分。、究背、内究状灰色系统理论在我国提出至今已有二十几年的历史，它的应用引起了人们的广泛兴趣，不论是我国粮食发展决策中总产量预测模型是对湖北年宏观经济的发展趋势的量化分析或河南人民胜利渠的最佳灌溉决策，还是武汉汉阳火车对火车装车吨位的预测等，无一不是灰色测系统理论杰出的硕果。、外究状灰色系统理论在国际上也产生了很大的影响IBM公司求将灰色系统软件加入其为全球服务的管理软件库。目前英国、美国、德国、日本、澳大利亚、加拿大、奥地利、俄罗斯等国家、区及国际组织有许多学者从事灰色系统的研究和应用。国内外84所高校开设了灰色系课程百名博士硕士研究生运用灰色系统的思想方法开展学科研究，撰写学位论文。国际、国内200多种学术期刊发表灰色系统论文，许多会议把灰色系统列为讨论专题，SCI、、、、、等纷检索我国灰色论著。、究意邓聚龙教授提出灰色系统有着重要的意:是系统思维和系统思想在方法论上的具体体;是科学方法论上的重大进,具原创性的科学意义和深远的学术影响，对系统科学的新贡献。2灰色系统及色预测的概、色统论展况、色统论提著名学者邓聚龙教授于20世纪70年代80代初提出。word档可自由复制编辑、色统论研对灰色系统产生于控制理论的研究中。若一个系统的内部特征是完全已知的，即系统的信息是充足完全的，我们称之为白色系统。若一个系统的内部信息是一无所知，一团漆黑，只能从它同部的联系来观测研究，这种系统便是黑色系统。灰色系统介于二者之间，灰色系统的一部分信息已知的，一部分是未知的。区别白色和灰色系统的重要标志是系统各因素间是否有确定的关系。、色统论应范在工程技术、社会、经济、农业、生态、环境等各种系统中经常会遇到信息不完全的情况。比如：农业方面，农田耕作面积往往因许多非农业的因素而改变，因此很难准确计算农田产量、值，这是缺乏耕地面积信息；生物防治方面，害虫与天敌间的关系即使是明确的，但天敌与饵料、虫与害虫间的许多关系却不明确，这是缺乏生物间的关联信息；一项土建工程，尽管材料、设备施工计划、图纸是齐备的，可是还很难估计施工进度与质量，这是缺乏劳动力及技术水平的信息一般社会经济系统，除了输出的时间数据列（比如产值、产量、总收入、总支出等）外，其输入据列不明确或者缺乏，因而难以建立确定的完整的模型，这是缺乏系统信息；工程系统是客观实，有明确的“内系即统内部与系统外部，或系统本体与系统环境较清楚地明输入与输出，因此可以较方便地分析输入对输出的影响，可是社会、经济系统是抽象的对象，有明确的“内系，不是客观实体，因此就难以分析输入（投入）对输出（产出）的影响这是缺乏“模型信息用什么模型，用什么量进行观测控制等信息息不完全的情况归纳来有：元素（参数）信息不完全；结构信息不完全；关系信息（特指“内系不完全运行的行为信息不完全。一个商店可看作是一个系统，在人员、资金、损耗、销售信息完全明确的情况下，可算出该店的盈利大小、库存多少，可以判断商店的销售态势、资金的周转速度等，这样的系统是白色系。遥远的某个星球，也可以看作一个系统，虽然知道其存在，但体积多大，质量多少，距离地球多远，这些信息完全不知道，这样的系统是黑色系统。人体是一个系统，人体的一些外部参数（如身高、体温、脉搏等）是已知的，而其他一些参数如人体的穴位有多少，穴位的生物、化学、物理性能，生物的信息传递等尚未知道透彻，这样系统是灰色系统。显然，黑色、灰色、白色都是一种相对的概念。世界上没有绝对的白色系统，因为任何系统总有未确知的部分，也没有绝对的黑色系统，因为既然一无所知，也就无所谓该系统的存在了。、种确性统究法比较析项目研究对象基础集合方法依据途径手段

灰色系统贫信息不确定灰色朦胧集信息覆盖灰序列算子

概率统计随机不确定康托集映射频率统计

模糊数学认知不确定模糊集映射截集word档可自由复制编辑数据要求侧重点目标特色

任意分布内涵现实规律小样本

典型分布内涵历史统计规律大样本

隶属度可知外延认知表达凭经验表1、色系的点灰色系统理论以“部分信息已知、部分信息未知”的“样本”不确定型系统的研究对象。（）灰色数学来处理不确定量，使之量化。在数学发展史上，最早研究的是确定型的微分方程，即在拉普拉斯决定论框架内的数学。他认为一旦有了描写事物的微分方程及初值，就能确知事物任何时候的运动。随后发展了概率论与理统计，用随机变量和随机过程来研究事物的状态和运动。模糊数学则研究没有清晰界限的事物如儿童和少年之间没有确定的年龄界限加以截然划分等，它通过隶属函数来使模糊概念量化，因能用模糊数学来描述如语言、不精确推理以及若干人文科学。灰色系统理论则认为不确定量是灰，用灰色数学来处理不确定量，同样能使不确定量予以量化。，，3不确定量

量化（用确定量的方法研究）1概率论与数理统计；2、模糊数学；3、灰色数学（灰色系统理论）（）分利用已知信息寻求系统的运动规律。研究灰色系统的关键是如何使灰色系统白化、模型化、优化。灰色系统视不确定量为灰色量。提出了灰色系统建模的具体数学方法，它能利用时间序列来确定微分方程的参数。灰色预测不是把观测到的数据序列视为一个随机过程，而是看作随时间变的灰色量或灰色过程，通过累加生成和累减生成逐步使灰色量白化，从而建立相应于微分方程解模型并做出预报。这样，对某些大系统和长期预测问题，就可以发挥作用。（）色系统理论能处理贫信息系统。灰色预测模型只要求较短的观测资料即可，这和时间序列分析，多元分析等概率统计模型要求较长资料很不一样。因此，对于某些只有少量观测数据的项目来说，灰色预测是一种有用的工。、见色统型GM(1,1)模GM(1,N)模GM(0,N)模GM(2,1)模模目前，最常用、研究最多的是GM(1,1)型。word档可自由复制编辑170170(h)[157,160]、色预灰色系统分析方法是通过鉴别系统因素之间发展趋势的相似或相异程度，即进行关联度分析，并通过对原始数据的生成处理来寻求系统变动的规律。生成数据序列有较强的规律性，可以用来建立相应的微分方程模型，从而预测事物未来的发展趋势和未来状态。灰色预测是用灰色模型来行定量析的，通常分为以下几类：灰色时间序列预测。用等时距观测到的反映预测对象特征的一系列数量（如产量、销量、人口数量、存款数量、利率等）构造灰色预测模型，预测未来某一时刻的特征量，或者达到某征量的时间。畸变预测（灾变预测模预测异常值出现的时刻测异常值什么时候出现在特定时区内。波形预测，或称为拓扑预测，它是通过灰色模型预测事物未来变动的轨迹。系统预测，是对系统行为特征指标建立一族相互关联的灰色预测理论模型，在预测系统整体变化的同时，预测系统各个环节的变化。上述灰预测方法的共同特点是：允许少数据预测；允许对灰因果律事件进行预测，比:

灰因白果律事件

在粮食生产预测中，影响粮食生产的因子很多，多到无法枚举，故为灰因，然而粮食产量却是具体的，故为白果。粮食预测即为灰因白果律事件预测。

白因灰果律事件

在开发项目前景预测时，开发项目的投入是具体的，为白因，而项目的效益暂时不很清楚，为灰果。项目前景预测即为灰因白果律事件预测。c.具有可检验性，包括：建模可性的级比检验（事前检验模度检验（模型检验测的滚动检验（预测检验、本概、数概在灰色系统中，灰数（或灰色数）是指信息不完全的数，例如的高约为170cm、重大致为60kg里为（cm灰数，分别记为、。如女孩身高在157－160cm之身高的灰数为数的化默认数，~~简称白化数，则灰数为化的全体。灰数有离散灰数（属离散集）和连续灰数（属于某一区间数的运算符合集合运算规律。、色成列在灰色系统理论中，把随机变量看成灰数，即是在指定范围内变化的所有白色数的全体。对灰数的处理主要是利用苏剧处理方法寻求数据间的内在规律，通过对已知数据列中的数据尽心处而产生新的数据列，以此来研究寻找数据的规律性，这种方法称为数据的生成。数据生成的常用式有累加生成、累减生成和加权累加生成。、加成把数列各项（时刻）数据依次累加的过程称为累加生成过程。由累加生成过程所得的数列称为word档可自由复制编辑(0)(0)、关联系数ˆˆ(0)(0)、关联系数ˆˆˆ累加生成数列。设原始数列为

x

(0)

x

(0)

(1),

(0)

(0)

(n

，令x

(1)

k()

(0)

(i),k1,2,

n,i称所得到的新数列

x(1)(1)x(1)(2),x(1)(

为数列

(0)

的1次加生成数列。类似地有kx(r(k)x

(r

(i),k1,2,

r

，称为

(0)

的

r

次累加生成数列。

i、减成对于原始数据列依次做前后相邻的两个数据相减的运算过程称为累减生成过IAGO原始数据列为

x

x

(1),

(1)

(2),

(1)

(n

，令x(0)(k)(1)()(1)kk2,3,,n

称所得到的数列

(0)

为

的1次累减生成数列。注：从这里的记号也可以看到，从原始数列

(0)

，得到新数列

，再通过累减生成可以还原出原始数列。实际运用中在数列

的基础上预测出(1)，通过累减生成得到预测数列(0)。、权值成设原始数列为

x(0)x(0)(1),x(0)x(0)n

，称

(0)

(kx(0)(k)

为数列

(0)

的邻值，

(0)

(k

为后邻值，

x

(0)

()

为前邻值。对于常数

[0,1]

，令z(0)k)

(k)

(1

1),k

n由此得到的数列

称为数列

(0)

在权下邻值生成数，权也为生成系数。特别地，当生成系数

时，则称

z(0)(k)

0.5x(0)(k

x(0)(k

1),k

2,3,

为均值生成数，也称等权邻值生成数。、联设则关联系数定义为：

(k)

minmin0

式中①

为第k个

X

与

的绝对误差；②

min

为两级最小差；③

maxX

为两级最大差；ρ称为分辨率，0<ρ<1,般ρ=0.5对单位不一，初值不同的序列，在计算相关系数前应首先进行初始化，即将该序列所有数据分别除以第一个数据。、关联度word档可自由复制编辑k(0)x(2),k(0)x(2),x())，其中定义的灰导数为()(1)

n

3简单的灰色测——预测目前使用座广泛的灰色预测模型就是关于数列预测的一个变量、一阶微分的M(1,1)模型。GM(1,1)模是基于随机的原始时间序列按时间累加后所形成的新的时间序列呈现的规律可用阶线性微分方程的解来逼近。经证明，经一阶线性微分方程的解逼近所揭示的原始时间数列呈数变化规律。因此，当原始时间序列隐含着指数变化规律时，灰色模型GM(1,1)的预测将是非常成的。、GM(1,1)测型基原x

设x(0)(0)(1),x(0)(2),(0)(为原始数列，其1次累加生成数列为(1)(1)x(k)xikn(1)dk)x

(0)

i(k)(k

(1)

(k令

z

(x

(1)

,x

(

为数列

的邻值生成数列，即z

(1)

()

(1)

(

x

(1)

(于是定义GM(1,1)的灰微分方程模型为()(1)(k)即

(0)

(k)

(1)

(k)

(1)在式(1)中，

x

(0)

()

称为灰导数，

称为发展系数，

z

(1)

()

称为白化背景值，

称为灰作用量。将时刻

代入(1)式有(0)(2)az(1)(3)az(3)

(1)’

x()()(1))引入矩阵向量记号：word档可自由复制编辑(0)bˆb(0)bˆb

(2)x(0)

,

u

,

(2)(3))

数据向量

参数向量

数据矩阵于是模型可表示为二乘法求它们的估计值为

YBu.

现在问题归结为求a的用元线性回归即小

B

T

)

T

注：实际上回归分析中求估计值是用软件计算的，有标准程序求解matlab,都以。对于GM(1,1)的灰微分方程1),如果将灰数

x

(0)()的刻k2,3,,n视连续变量x(1)视为时间的数

xt)是x(0)()

对应于导数

(1)

()

背值

z(1)()对于导数xt)

。于是GM(1,1)的灰微分方程对于的白微分方程为dx

(1)

(t

(1)(t)

(2)称之为GM(1,1)的化型。式子(以值

(t(0)

的解为x

t

(0)

b(1)(ta注：的化型(2)并不是由1)直接推导出来的，仅仅是一种“借用”或“白化默认以从的白化型推导出来的结果，要在不与定义矛盾的情形下才成立。后面我们看到，对数据列有要求。令X(0)为GM(1,1)建模序列，

(0)

x

(0)

(0)

(2),...,x

(0)

())为X的1-AGO序，

(1)

x

(1)

(2),...,x

())()(0)(ikni令

(1)

为

(1)

的紧邻均值MEAN）成列Z

z

(1)

())word档可自由复制编辑ˆ(1)(0)ˆ(1)(0)z)(1)()x(1)(则GM(1,1)定义型，即GM(1,1)的灰微分方程型为x(0)()(1)(k)模型符号含义为

（3）GM

（1，

）GreyModel阶方程

个变量式中

称为发展系数，

为灰色作用量。设为估数向量，即

,b)

T

，则灰微分方程(的最小二乘估计参数列满足

B

)

Yn其中B

(2)1

Yn

(0)(3))称

为灰色微分方程x

(0)

()

(1)

()的化方程，也叫影子方程。如上所述，则有）白方程

的解也称时间响应函数为

(1)

b(t)(x(0))ta）灰微分方程

(0)()(1)(k)

的时间响应序列为b(1)(x(0))e(k)1,2,...,a）取x

(0)x

(0)

，则

(1)

b(k(0)(1))(k)na）还值

(0)

(

(1)

(

(上式即为预测方程。有关建模的问题说明如下：即

．定原序列X和不同。

中的数据不一定要全部用来建模原始数据的取舍不同得模型不同，word档可自由复制编辑(0)(i)n(0)(i)n．模的据取舍应保证建模序等时距、相连，不得有跳跃出现。．一般模数据序列应当由最的数据及其相邻数据构成，当再出现新数据时，可采用两种方法处理：一是将新信息加入原始序列中，重估参数；二是去掉原始序列中最老的一个数据，加上最新的数据，所形成的序列和原序列维数相等，再重估参数。、模型验模型的检验分为三个方面：残差检验；关联度检验；后验差检验。、差验残差大小检验，即对模型值和实际值的残差进行逐点检验。首先按模型计算

(1)

(i

，将

(i

累减生成

(0)

(i)

，最后计算原始序列

(0)(i)

与

(0)

(i)

的绝对残差序列(0)

i),in,

(i)x

(0)

i)

(0)

(i)及相对残差序列1,2,...,nii

i)(0)

%并计算平均相对残差i

i给定，当，立时，称模型为残差合格模型。、联检关联度检验，即通过考察模型值曲线和建模序列曲线的相似程度进行检验。按前面所述的关联度计算方法，计算出于便满意的。

(0)

(i与原始序列x(0)(i的联系数，然后算关联度，根据经验，关联度大、验检后验差检验，即对残差分布的统计特性进行检验。（1计算出原始序列的平均值：

(0)

(0)(i)i（2）计原始序列

X(0)的均方差：(0)(ii

(0)

（3计残差的均值：(0)

(i)i（4计算残差的均方差：word档可自由复制编辑()Si()(1),e(2),...,e()Si()(1),e(2),...,e()e令

计算方差比：计算小残差概率：S0.6745，e(i)i

S1S2Ppi)1，即Ppi0若对于给定的

当C时模为均方差比合格模型对给定的P时，称模型为小残差概率合格模型。后P>0.95>0.80>0.70<0.70

C<0.35<0.5<0.65>0.65

模型精度优合格勉强合格不合格若相对残差、关联度、后验差检验在允许的范围内，则可以用所建的模型进行预测，否则应进行残差修正。、残差型当原始数据序列

建立的模型检验不合格时，可以用残差模型来修正。如果原始序列建立的模型不够精确，可以用残差模型来提高精度。若用原始序列X

建立的模型

(ix

(0)

(1)e

可获得生成序列X(1)的测值，定义残差序列

e(0)(j)x(j)

(1)

(j)

。若取

j,in

，则对应的残差序列为：

e

(0)

(0)(0)计算其生成序列得修正模型

e(1)(k)

，并据此建立相应的GM(1,1)型：b(0)eaeword档可自由复制编辑

bb(1)(kx(0)(1)ek)(0)eaa

x(2),...,)为系统xx(2),...,xn(0)(0)(1)(1)(0)(1)(1)x(2),...,)为系统xx(2),...,xn(0)(0)(1)(1)(0)(1)(1)BBBˆ(1)12(1)(1)(1)(1)(1)(1)(0)(1)(1)(1)(1)(1)(1)(1)(1)1iiiiii其中

k()0k

为修正参数。应用此模型时要考虑：．一不是使用全残差数据来建立模型，而只是利用了部分残差。．修模型所代表是差分微分方程，其修正作用与

(k)中i的值有关。、模型如果考虑的系统由若干个相互影响的因素组成，设

(0)(0)(0)(0)1特征数据序列，而

2

(0)(0)(0)2X

N

(0)(0)N

...x(2),...,xN

(为相关因素序列。X

i

为X

i

(0)

的序iN，Z为X1x(k)((1ii

的紧邻生成序列，则称（6为灰微分方程。定义

i

...

为灰微分方程的参数列，根据最小二乘法可以得出式中

(2)x(2)(3)1...)2()

...(2)...(3)......(n)

x

(2)x

...x

(n)

T称

2

x3N

为灰微分方程）的白化方程，也称子方程。于是，我们有）白方程（）的解为1

(t)

i

(t)(0)ii1

i

(0)dtii

Nxx(0)bxt)dtword档可自由复制编辑(1)(i)ii(1)(i)iiiiaii(0)kk）当X

i

(1)

iN)变幅度很小时，可视xii

(1)

()为灰常量，这样，灰色微分方程（）的近似时间响应式为

i(1)

1N1(kx(0)bx(kb

(k

其中x

取x

。）累还原式为(0)(1)11

k)

、色系建的本路可以概括为以下几点：（）立模型常用的数据有以下几种科实验数据2验数据；3生数据；4决数据。（）列生成数据是建立灰色模型的基础数据。一般非负序列累加生成后，得到准光滑序列，对于满足光滑条件的序列，即可建微分模型。模型精度可以通过不同的灰数生成方式，数据的取舍，序列的调整、修正以及不同级别的残差模型充得到提高。灰色系统理论采用残差大小检验、关联度检验、后验差检验三种方法检验、判断模型的精度。、

灰色关联分析、色关分理及法对于两个系统之间的因素，其随时间或不同对象而变化的关联性大小的量度，称为关联度。在系统发展过程中，若两个因素变化的趋势具有一致性，即同步变化程度较高，即可谓二者关联度较高；反之，则较低。因此，灰色关联分析方法，是根据因素之间发展趋势的相似或相异程度亦即灰关联度，为衡量因素间关联程度的一种方法。灰色系统理论提出了对各子系统进行色关联度分析的概念，意图透过一定的方法，去寻求系统中各子系统（或因素）之间的数值关系因此，灰色关联度分析对于一个系统发展变化态势提供了量化的度量，非常适合动态历程分析。在实际问题中，如社会系统、经济系统、农业系统、生态系统、教育系统等都包含有许多种因素，多种因素共同作用的结果决定了系统的发展态势，人们常常希望弄清在众多因素中，哪些主要因素？哪些是次要因素？如粮食生产系统，我们希望提高粮食总产量，而影响粮食总产量的因素是多方面的，有播种面积、水利、肥料、土壤、种子、劳力、气候、耕作技术和政策环境等。为了实现少投入多产出并取得良好的经济效益、社会效益和生态效益，就必须进行系统分析。传统方法：数理统计中的回归分析、方差分析、主成分分析、因子分析等的不足之处：．求有大量数据，数据量少就难以找出统计规律；word档可自由复制编辑ii．求样本服从某个典型的概率分布求因素数据与系统特征数据之间呈现线性关系且各因素之间彼此无关，这种要求往往难以满足；．算量大，一般要靠计算机帮助；．能出现量化结果与定性分析结果不符的现象，导致系统的关和规律遭到歪曲和颠倒。关联度分析法便是对付这种情况的一副良药，它的基本思想是根据序列曲线几何形状的相似程度来判断其联系是否紧密。一般地，曲线越接近，相应序列之间的关联度就越大，反之就越小、色联术应直接应用：因素分析、方案决策、优势分析；与其他方法结合：灰色关联和聚类方法相结合、灰色关联分析和层次分析法相结合、优化方法非线性模型与灰关联分析相结合；应用新领域于全科学,如矿安全的分析与评估于环科学中如质评价、大气环境质量评价等；应用于医学诊断中；应用于油田的开发中；应用于系统水文学中。此外，在灰色关联技术的带动下，相继产生了灰色地质学、灰色育种学、灰色控制理论、灰色混沌理论、区域经济灰色系统分析、灰色价值学、灰色综防学等新兴学科。、色联计式改进设一个母因素的时间序列和

子因素时间数列分别为()x(0iiii令

i

(k)minmax

，则称

i

(k)

表示第

t

时刻子因素和母因素的关联系数，其中为分辨系数，

0.5.

x()()，imin()(k)i

，则母因素和子因素的关联度

i

max()()i计算式定义为

.1nRnk

i

()

（）式（10）关联度计算式存在如问题：①

i

与分辨系数有关，R的不唯一i

②灰联度不具有保序.即，可能.

不单函数.对某个，能i

i

j

，而对另外一个③关度不具有规范性.即对全相关的子因素与母因素的时间序列，

i

.例如，91,2,...,10),o1

显然

与

X

完全相关，但是若取

，则按上式算word档可自由复制编辑ii11iiiii11iii1

1

3

.④灰联分析不能体现负相关由关联度计算式知

i

。然而，设有两个时间序列X0

显然有

x(t)(tt10

与

X

存在典型的负相关⑤通分辨系数0.5，有灰关联度的值恒大于或等于.鉴于上述关联度计算式存在的问题要出改进的关联度计算.文献提改进的关联度计算式，步骤如下：a．作一次累减生成，即相当于在曲线上不同时点的斜(1)(k),a(1)(x(((k00iii．计算与两间数列标准差，记inx(k)x(k)，

，

.

xi

)nk

，

xi

)nk

.即

和分别表示两时间数列和的值iixc．计算时(k)，x的联系数i(t),(k))sign((1)((k))a(1)(()))0i0i

分别表示两时间数列和x0i(x(k))(1)(0ixxi

i,

的标准差.其中1,当xsignx当x0,当xd．计算与的关联度inRnk

((),(k))0i

.(11)上述改进的灰关联度式（11）具唯一性和规范性，~~且ii

当且仅当

与

i

完全相关.5传染病的问据统计，高校传染病发病率情况如下，试建立GM(1,1)预模型，并预测年传染病发病率。另外传染病发病率较高的为痢疾、肝炎、疟疾，那么哪一种疾病更危害学生的健康呢？历只有经验估计，如下统计表（表3，没有定量分析。word档可自由复制编辑(2)(3)(4)(4)(5)(1)2(6)(7)(7)(8)(9)111(2)(3)(4)(4)(5)(1)2(6)(7)(7)(8)(9)11164.53(0)(0)nb8年度19841985

1986

1987

1988

1989

199019911992发病率100.23

85.62

64.53

86.63

105.89

83.55

316.47135.93

56.56痢疾率肝炎率疟疾率

19.469.3771.34

44.198.2933.14

30.9825.817.74

45.8638.222.55

55.4725.2125.21

60.3418.574.64

271.933.4311.14

94.9336.684.32

38.2814.144.04表3、染病病的预解：设

(0)

(第1步构累加生成序列

(1)

(k100.23,185.85,250.38,337.01,442.9,526.45,842.92,978.85第2步构数据矩阵B和据向量Y2121(1)21B1(5)(1)(6)212121

1218.115910.8851

1111Y

(0)(2)85.62(0)(4)86.63(0)(5)(0)(6)(0)第3步计

aB)B

T

B

2854081.97word档可自由复制编辑0.000717670.4956944682.468986ˆˆii0.000717670.4956944682.468986ˆˆii(

B)

0.00071767

TB)

0.066654第4步得预测模型dt

0.0666544x82.4689865

(1)

(k.729e

0.0666544

(x

(0)

ba第残差检验(1)根据预测公式，计算X

(1)

()得(1))(2)累减生成

(0)

()序,

(0)

k)原始序列

(0)

((3)计算绝对残差和相对残差序列绝对残差序列：

(0)

36.8,187.82,相对残差序列：

,0.216,0.063,0.441,0.593,1.599第进行关联度检验计序列

(0)

与(0)

的绝对残差序列

k)minmax

18.7,()(k)36.8,187.82,90.44计关联系数由于只有两个序列（即一个参考序列，一个被比较序列）故不再寻求第二级最小差和最大差。)(k1,...,9,0.5)k)求得)计关联度(k)0.148kR0.148

不满足

时的检验准则

的。第后验差检验计：1x85.6264.5386.63316.47135.9356.569word档可自由复制编辑x(k)1x(k)1计X

序列的均方差：计残差的均值：

(0)119

计残差的均方差：

n

计

C

：

C

S20.964S1计小残差概率：S0S53.291ek)0.78,5.8,33.23,17.9,5.9,36.03,188.6,0.81,89.7可得小残差概率为P故P0.7表（后验检验判别参照表此模i型

x

(1)(k13380.0666544

1237.2618

勉强合格。第预测：k10,x

(0)

x

157.13即1993年的传染病发病率预测值为见附件matlab程序、、）、种染的联析建模步1：数据处理x(k)0x(k)1x(9.73,8.29,...,14.142x(k)3步2：计算关联度系数。经计算min(k)(ki((305.33i

i

(k)

18.18x((k)305.33i据前文分析，设分辨系数为0.5将相应的

(k)与(ki

的数值代入式

i

(k)

max

中，得22

0.73,0.88,0.92,0.88,0.84,0.97,0.87,0.88,11word档可自由复制编辑19202391920239步3：算出关联度

(,)i

由式

i

n

i

(k)

分别计算出痢疾肝、疟疾关于整个传染病的关联度

R,

：19R(xx)k0.8250k1R(,)9k

2

()19RR(xx)0k

3

()0.8237步4比关联度大小给出结论

23

，说明痢疾发病与整个传染病关系最密切肝炎、疟疾与整个传染病的密切程度依次差些。（见附件）所得到的结论是：通过高校传染病发病率情况