灰色系统理论及其应用

灰色系统理论及其应用第一章灰色系统的概念与基本原理1.1灰色系统理论的产生和发展动态1982年，北荷兰出版公司出版的《系统与控制通讯》杂志刊载了我国学者邓聚龙教授的第一篇灰色系统理论论文”灰色系统的控制问题”，同年，《华中工学院学报》发表邓聚龙教授的第一篇中文论文《灰色控制系统》，这两篇论文的发表标志着灰色系统这一学科诞生1985灰色系统研究会成立，灰色系统相关研究发展迅速。1989海洋出版社出版英文版《灰色系统论文集》，同年，英文版国际刊物《灰色系统》杂志正式创刊。目前，国际、国内300多种期刊发表灰色系统论文，许多国际会议把灰色系统列为讨论专题。国际著名检索已检索我国学者的灰色系统论著3000多次。灰色系统理论已应用范围已拓展到工业、农业、社会、经济、能源、地质、石油等众多科学领域，成功地解决了生产、生活和科学研究中的大量实际问题，取得了显著成果。1.2几种不确定方法的比较概率统计，模糊数学和灰色系统理论是三种最常用的不确定系统研究方法。其研究对象都具有某种不确定性，是它们共同的特点。也正是研究对象在不确定性上的区别，才派生了这三种各具特色的不确定学科。模糊数学着重研究“认识不确定”问题，其研究对象具有“内涵明确，外延不明确”的特点。比如“年轻人”内涵明确，但要你划定一个确定的范围，在这个范围内是年轻人，范围外不是年轻人，则很难办到了。概率统计研究的是“随机不确定”现象，考察具有多种可能发生的结果之“随机不确定”现象中每一种结果发生的可能性大小。要求大样本，并服从某种典型分布。灰色系统理论着重研究概率统计，模糊数学难以解决的“小样本，贫信息”不确定性问题，着重研究“外延明确，内涵不明确”的对象。如到2050年，中国要将总人口控制在15亿到16亿之间，这“15亿到16亿之间“是一个灰概念，其外延很清楚，但要知道具体数值，则不清楚。1.3灰色系统理论的基本概念定义1.3.1信息完全明确的系统称为白色系统。定义1.3.2信息未知的系统称为黑色系统。定义1.3.3部分信息明确，部分不明确的系统称为灰色系统。1.4灰色系统理论的基本原理公理1(差异信息原理)“差异“是信息，凡信息必有差异。公理2（解的非唯一性原理）信息不完全，不确定的解是非唯一的。公理3（最少信息原理）灰色系统理论的特点是充分开发利用已占有的“最少信息“。公理4(认知根据原理)信息是认知的根据。公理5（新信息优先原理）新信息对认知的作用大于老信息。公理6（灰性不灭原理）：信息不完全是绝对的1.5灰色系统理论的主要内容灰色系统理论经过20多年的发展，现在已经基本建立起一门新兴学科的结构体系。其主要内容包括以灰色代数系统，灰色方程、灰色矩阵等为基础的理论体系。以灰色序列生成为基础的方法体系，以灰色关联空间为依托的分析体系。以灰色模型（GM）为核心的模型体系，以系统分析，评估，建模，预测，决策，控制，优化为主体的技术体系。1.6灰数灰数是灰色系统理论的基本“单元“或”细胞“。我们把只知道大概范围而不知道其确切值的数称为灰数。在应用中，灰数实际上指在某一个区间或某个一般的数集内取值的不确定数。通常用记号“”表示灰数。灰数有以下几类：1.仅有下界的灰数。有下界而无上界的灰数记为∈[,]a，其中a是灰数的下确界，是确定的数，我们称[,]a为的取数域，简称的灰域。2.仅有上界的灰数。有上界而无下界的灰数记为∈[,]a，其中a是灰数的上确界，是确定的数。3.区间灰数。既有下界又有上界的灰数称为区间灰数，记为∈[,]aa4.连续灰数与离散灰数。5.黑数与白数。当∈[,]，称为黑数；当∈[,]aa且aa时，称为白数。6.本征灰数与非本征灰数。本征灰数是指不能或暂时还不能找到一个白数作为其“代表”的灰数，比如一般的事前预测值，宇宙的总能量等。非本征灰数是指凭先验信息或某种手段，可以找到一个白数作为其代表的灰数。我们称此白数为相应灰数的白化值。第二章序列算子与灰色序列生成灰色系统理论的主要任务之一，是根据社会，经济，生态等系统的行为特征数据，寻找不同系统变量之间或某些系统变量自身的数学关系和变化规律。灰色系统理论认为任何随机过程都是在一定幅值范围和一定时区内变化的灰色量，并把随机过程看成灰色过程。灰色系统理论是通过对原始数据的挖掘，整理来寻求其变化规律的，这是一种就数据寻找数据的现实规律的途径，我们称为灰色序列生成。灰色系统理论认为，尽管客观系统表象复杂，数据离乱，但它总是有整体功能的，因此必然蕴含某种内在规律。关键在于如何选择适当的方式去挖掘它和利用它。一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。例如考虑4个数据，记为)4(),3(),2(),1()0()0()0()0(XXXX，其数据见下表：序号1234符号)1()0(X)2()0(X)3()0(X)4()0(X数据121.54将上表数据作图得0123451234XY上图表明原始数据)0(X没有明显的规律性，其发展态势是摆动的。如果将原始数据作累加生成，记第K个累加生成为)()1(KX,并且1)1()1()0()1(XX321)2()1()2()0()0()1(XXX5.45.121)3()2()1()3()0()0()0()1(XXXX5.735.121)4()3()2()1()4()0()0()0()0()1(XXXXX得到数据如下表所示序号1234符号)1()1(X)2()1(X)3()1(X)4()1(X数据134.57.50123456781234XY上图表明生成数列X(1)是单调递增数列。2.1冲击扰动系统与序列算子定义2.1.1设0000((1),(2),,())Xxxxn为系统真实行为序列，而观察到的系统行为数据序列为000012((1),(2),,())((1),(2),,())nXxxxnxxxnX其中，12(,)n为冲击扰动项（干扰项）。X称为冲击扰动序列。所以本章我们的讨论围绕：由XX0展开（扰动还原真实）2.2缓冲算子公理定义2.2.1设系统行为数据序列为((1),(2),,())Xxxxn，1.若2,3,,()(1)0knxkxk，则称X为单调增长序列；2.若1中不等号反过来成立，则称X为单调衰减序列；3.若,{2,3,},()(1)0,()(1)0kknxkxkxkxk有，则称X为随机振荡序列。4.设max()|12,3,,,()|12,3,,Mxkknmxkkn，，，则称M-m为序列X的振幅定义2.2.2设((1),(2),,())Xxxxn为系统行为数据系列，D为作用于X的算子，X经过算子D作用后所得序列记为((1),(2),,())XDxdxdxnd称D为序列算子，称XD为一阶算子作用序列。序列算子的作用可以多次，相应的，若123,,DDD都是序列算子，我们称12DD为二阶算子，并称12121212((1),(2),,())XDDxddxddxndd为二阶算子作用序列，同理，123DDD为三阶序列算子……定义2.2.3称下述三公理为缓冲算子三公理，满足缓冲算子三公理的序列算子D称为缓冲算子，一阶，二阶，三阶……缓冲算子作用序列称为一阶，二阶，三阶……缓冲序列。公理1（不动点公理）设((1),(2),,())Xxxxn为系统行为数据系列，D为序列算子，则D满足()()xndxn。不动点公理限定在序列算子作用下，系统行为数据序列的数据()xn保持不变。根据定性分析的结论，亦可使()xn以前的若干个数据在序列算子作用下保持不变。例如，令()()()()xjdxjxidxi且,1,2,1,1,,.jkikkn其中公理2.（信息充分利用公理）系统行为数据序列X中的每一个数据(),1,2,xkk，都要充分地参与算子的作用全过程公理3（解析化、规范化公理）任意的(),(1,2,)xkdk，皆可由一个统一的(1),(2),,()xxxn的初等解析式表达。定义2.2.4设X为原始数据序列，D为缓冲算子，当X分别为增长序列，衰减序列或振荡序列时：1.若缓冲序列XD比原始序列X的增长速度（或衰减速度）减缓或振幅减小，则称缓冲算子D为弱化算子。2.若缓冲序列XD比原始序列X的增长速度（或衰减速度）加快或振幅增大，则称缓冲算子D为强化算子。2.3实用缓冲算子的构造定理2.3.1设原始数据序列((1),(2),,())Xxxxn令缓冲序列((1),(2),,())XDxdxdxnd其中1()[()(1)()]1xkdxkxkxnnk；k=1,2,……，n,则当X为增长序列，衰减序列或振荡序列时，D为弱化算子，并称为平均弱化缓冲算子（AWBO）证明：直接利用(),(1,2,)xkdk的定义，可知定理成立。推论2.3.1对于定理1中定义的弱化算子D，令2222((1),(2),,())XDXDDxdxdxnd21()[()(1)()],1,21xkdxkdxkdxndknnk，则2D对于增长序列，衰减序列或振荡序列时，皆为二阶弱化算子。定理2.3.2设原始序列和其缓冲算子序列分别为((1),(2),,())Xxxxn((1),(2),,())XDxdxdxnd其中(1)(2)(1)()(),1,2,121xxxkkxkxkdknk()()xndxn则当X为增长序列，衰减序列或振荡序列时，D为强化算子。推论2.3.2设D为定理2中定义的强化算子，令2222((1),(2),,())XDXDDxdxdxnd，其中2()()()xndxndxn，2(1)(2)(1)()(),1,2,121xdxdxkdkxkdxkdknk则2D对于增长序列，衰减序列或振荡序列皆为二阶强化算子。定理2.3.3原始数据序列和其缓冲算子序列分别为((1),(2),,())Xxxxn((1),(2),,())XDxdxdxnd其中()(1)(1)()(),1,2,()(1)/2kxkkxknxnxkdknnknk，则当X为增长序列，衰减序列或振荡序列时，D为弱化算子，并称D为加权平均弱化缓冲算子（WAWBO）定理2.3.4设((1),(2),,())Xxxxn为非负的系统行为数据序列，令((1),(2),,())XDxdxdxnd其中1111()[()(1)()][()],1,2,nnknkikxkdxkxkxnxikn。则当X为增长序列，衰减序列或振荡序列时，D为弱化缓冲算子，并称D为几何平均弱化缓冲算子（GAWBO）定理2.3.5设((1),(2),,())Xxxxn为系统行为数据序列，各时点的权重向量为12(,)n，则((1),(2),,())XDxdxdxnd其中11()(1)()(),1,2,kknkknxkxkxnxkdkn。则当XD皆为弱化缓冲算子，并称D为加权平均弱化缓冲算子（WAWBO）。定理2.3.6设((1),(2),,())Xxxxn，各时点的权重向量为12(,)n0，令((1),(2),,())XDxdxdxnd其中11111()[()(1)()][()],1,2,kknkknkknnikxkdxkxkxnxikn则当XD为弱缓冲算子，并称D为加权几何平均弱化缓冲算子（WGAWBO）。定理2.3.7设((1),(2),,())Xxxxn为系统行为数据序列，令((1),(2),,())XDxdxdxnd其中2(1)()(),1,2,()(1)()nkxkxkdknxkxkxn。则当X为增长序列，衰减序列或振荡序列时，D为强化缓冲算子，并称D为平均强化缓冲算子（ASBO）定理2.3.8设((1),(2),,())X