灰色系统理论及其应用

灰色系统理论及其应用第一章灰色系统的概念与基本原理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）用灰色数学来处理不确定量，使之量化。

在数学发展史上，最早研究的是确定型的微分方程，即在拉普拉斯决定论框架内的数学。

他认为一旦有了描写事物的微分方程及初值，就能确知事物任何时候的运动。

随后发展了概率论与数理统计，用随机变量和随机过程来研究事物的状态和运动。

模糊数学则研究没有清晰界限的事物，如儿童和少年之间没有确定的年龄界限加以截然划分等，它通过隶属函数来使模糊概念量化，因此能用模糊数学来描述如语言、不精确推理以及若干人文科学。

灰色系统理论则认为不确定量是灰数，用灰色数学来处理不确定量，同样能使不确定量予以量化。

1，2，3不确定量量化（用确定量的方法研究）1、概率论与数理统计；2、模糊数学；3、灰色数学（灰色系统理论）（2）充分利用已知信息寻求系统的运动规律。

研究灰色系统的关键是如何使灰色系统白化、模型化、优化。

灰色系统视不确定量为灰色量。

提出了灰色系统建模的具体数学方法，它能利用时间序列来确定微分方程的参数。

灰色预测不是把观测到的数据序列视为一个随机过程，而是看作随时间变化的灰色量或灰色过程，通过累加生成和累减生成逐步使灰色量白化，从而建立相应于微分方程解的模型并做出预报。

这样，对某些大系统和长期预测问题，就可以发挥作用。

（3）灰色系统理论能处理贫信息系统。

灰色预测模型只要求较短的观测资料即可，这和时间序列分析，多元分析等概率统计模型要求较长资料很不一样。

因此，对于某些只有少量观测数据的项目来说，灰色预测是一种有用的工具。

1.6灰数灰数是灰色系统理论的基本“单元“或”细胞“。

我们把只知道大概范围而不知道其确切值的数称为灰数。

在应用中，灰数实际上指在某一个区间或某个一般的数集内取值的不确定数。

通常用记号“⊗”表示灰数。

灰数有以下几类：1. 仅有下界的灰数。

有下界而无上界的灰数记为⊗∈[,]a -∞，其中a 是灰数⊗的下确界，是确定的数，我们称[,]a -∞为⊗的取数域，简称⊗的灰域。

2. 仅有上界的灰数。

有上界而无下界的灰数记为⊗∈[,]a --∞ ，其中a --是灰数⊗的上确界，是确定的数。

3. 区间灰数。

既有下界又有上界的灰数称为区间灰数，记为⊗∈[,]a a ----4.连续灰数与离散灰数。

5. 黑数与白数。

当⊗∈[,]-∞+∞，称⊗为黑数；当⊗∈[,]a a ----且a a ----=时，称⊗为白数。

6. 本征灰数与非本征灰数。

本征灰数是指不能或暂时还不能找到一个白数作为其“代表”的灰数，比如一般的事前预测值，宇宙的总能量等。

非本征灰数是指凭先验信息或某种手段，可以找到一个白数作为其代表的灰数。

我们称此白数为相应灰数的白化值。

第二章 序列算子与灰色序列生成灰色系统理论的主要任务之一，是根据社会，经济，生态等系统的行为特征数据，寻找不同系统变量之间或某些系统变量自身的数学关系和变化规律。

灰色系统理论认为任何随机过程都是在一定幅值范围和一定时区内变化的灰色量，并把随机过程看成灰色过程。

灰色系统理论是通过对原始数据的挖掘，整理来寻求其变化规律的，这是一种就数据寻找数据的现实规律的途径，我们称为灰色序列生成。

灰色系统理论认为，尽管客观系统表象复杂，数据离乱，但它总是有整体功能的，因此必然蕴含某种内在规律。

关键在于如何选择适当的方式去挖掘它和利用它。

一切灰色序列都能通过某种生成弱化其随机性，显现其规律性。

例如考虑4个数据，记为)4(),3(),2(),1()0()0()0()0(X X X X ，其数据见下表：序号1 2 3 4 符号)1()0(X )2()0(X )3()0(X )4()0(X 数据1 2 1.5 4将上表数据作图得123451234X Y上图表明原始数据)0(X 没有明显的规律性，其发展态势是摆动的。

如果将原始数据作累加生成，记第K 个累加生成为)()1(K X ,并且1)1()1()0()1(==X X321)2()1()2()0()0()1(=+=+=X X X5.45.121)3()2()1()3()0()0()0()1(=++=++=X X X X5.735.121)4()3()2()1()4()0()0()0()0()1(=+++=+++=X X X X X得到数据如下表所示 序号1 2 3 4 符号)1()1(X )2()1(X )3()1(X )4()1(X 数据 1 3 4.5 7.5123456781234X Y上图表明生成数列X (1)是单调递增数列。

2.1冲击扰动系统与序列算子定义2.1.1 设0000((1),(2),,())X x x x n = 为系统真实行为序列，而观察到的系统行为数据序列为000012((1),(2),,())((1),(2),,())n X x x x n x x x n X εεεε==+++=+其中，12(,)n εεεε=为冲击扰动项（干扰项）。

X 称为冲击扰动序列。

所以本章我们的讨论围绕：由X →X 0展开（扰动还原真实）2.2缓冲算子公理定义 2.2.1 设系统行为数据序列为((1),(2),,())X x x x n =，1. 若2,3,,()(1)0k n x k x k ∀=-->，则称X 为单调增长序列；2. 若1中不等号反过来成立，则称X 为单调衰减序列；3. 若,{2,3,},()(1)0,()(1)0k k n x k x k x k x k '''∃∈-->--<有，则称X 为随机振荡序列。