模糊数学基本知识new

1、一模糊数学的基础知识1模糊集、隶属函数及模糊集的运算。普通集合A，对，有或。如果要进一步描述一个人属于年轻人的程度大小时，仅用特征函数就不够了。模糊集理论将普通集合的特征函数的值域推广到0，1闭区间内，取值的函数以度量这种程度的大小，这个函数（记为）称为集合的隶属函数。即对于每一个元素，有0,1内的一个数与之对应。（1）模糊子集的定义：射给定论域U，U到0，1上的任一映射：都确定了U上的一个模糊集合，简称为模糊子集。称为元素属于模糊集的隶属度。映射所表示的函数称为隶属函数。例如：设论域U=0，100，U上的老年人这个集合就是模糊集合：若在集合U上定义了一个隶属函数，则称为模糊集。（2）模糊集合

2、的表示：，称为元素属于模糊集的隶属度；则模糊集可以表示为：。或 ，（3）模糊集合的运算：，并集：，交集：，补集：，包含：，2模糊集的截集已知U上模糊子集对，则称为模糊集的-截集；称为模糊集的-强截集；称为、的置信水平或阀值。二模糊数学的基本定理1 模糊截积：已知U上模糊子集对，也是U上模糊集，其隶属函数为：；称为为与的模糊截积。2 分解定理1：已知模糊子集，则推论1：对3 分解定理2：已知模糊子集，则推论2：对 三模糊关系与模糊聚类1 模糊关系与模糊关系的合成（1） 模糊关系普通集合的经典关系，模糊关系：从U到V 上的一个模糊关系：，表示具有的关系程度，。（满足01）称为U到V 上的一个模糊关

3、系的模糊矩阵。（2）设和=为两个模糊矩阵，令 ，1，2，,=1,2,。则称矩阵=为模糊矩阵与的褶积，记为 ,其中“”和“”的含义为 显然，两个模糊矩阵的褶积仍为模糊矩阵2. 模糊等价矩阵及其矩阵设方阵为以模糊矩阵，若满足 = 则称为模糊等价矩阵。模糊等价矩阵可以反映模糊分类关系的传递性，即描述诸如“甲像乙，乙像丙，则甲像丙”这样的关系。设为一个模糊等价阵，01为一个给定的数，令 则称矩阵为的截阵例如，为一个模糊等价阵，取0.4<,则=若取，则=2 模糊聚类：模糊划分的概念最早由Ruspini提出，利用这一概念人们提出了多种聚类方法，比较典型的有：基于相似性关系和模糊关系的方法(包括聚合法

4、和分裂法)，基于模糊等价关系的传递闭包方法、基于模糊图论最大树方法，以及基于数据集的凸分解、动态规划和难以辨识关系等方法. 然而由于上述方法不适用于大数据量情况，难以满足实时性要求高的场合，因此其实际的应用不够广泛，故在该方面的研究也就逐步减少了. 实际中受到普遍欢迎的是基于目标函数的方法，该方法设计简单、解决问题的范围广，最终还可以转化为优化问题而借助经典数学的非线性规划理论求解，并易于计算机实现. 因此，随着计算机的应用和发展，该类方法成为聚类研究的热点.（1）模糊聚类的基本概念模糊聚类目标函数的演化模糊聚类方法 模糊聚类法和一般的聚类方法相似，先将数据进行标准化，计算变量间相似矩阵或样品

5、间的距离矩阵，将其元素压缩到0与1之间形成模糊相似矩阵，进一步改造为模糊等价矩阵，最后取不同的标准，得到不同的截阵，从而就可以得到不同的类。具体步骤如下：第一步：数据标准化1.数据矩阵设论域为被分类的对象，每个对象又由个指标表示其性状： （）于是得到原始数据矩阵为2.数据标准化 在实际问题中，不同的数据一般有不同的量纲。为了使有不同的量纲的量也能进行比较，通常需要对数据作适当的变换。但是，即使这样得到的数据也不一定在区间0，1上。因此，这里所说的数据标准化，就是要根据模糊矩阵的要求，将数据压缩到区间0,1上。通常需要作如下变换：（1）平移·标准差变换： （）其中。经过变化后，每个变量

6、的均值为0，标准差为1，且消除了量纲的影响。但是，这样得到的还不一定在区间0,1上。（2）平移·级差变换 （）显然有，而且也消除了量纲的影响。第二步：标定（建立模糊相似矩阵）设论依照传统的方法确定相似系数，建立模糊相似矩阵，与的相似程度。可根据问题的性质，选取下列公式之一计算1. 数量积法其中显然，若中出现负值，也可采用下面的方法将压缩在0,1上令，则。当然也可用上述的平移·级差变换。2 夹角余弦法若将变量的个观测值与变量的相应个观测值看成维空间中的两个向量，正好时这两个向量夹角的余弦。3.相关系数法从统计角度看，两个随机变量的相关系数是描述这两个变量关联性（线性关系）强弱

7、的一个很有用的特征数字。因此，用任意两个变量的个观测值对其相关系数的估计可作为两个变量关联性的一种度量，其定义为，其中（1，2，,）见（，1，2，, ，）。（1）其实就是的样本相关矩阵中的各元素。4指数相似系数法，其中，而需要注意的是，相关系数法与指数相似系数法中的统计指标的内容是不同的。5最大最小法6算术平均最小法 7几何平均最小法（上述5，6，7三种方法均要求，否则也要做适当变换）8绝对值减数法适当选取，使得0。9绝对值倒数法其中适当选取，使得0。10绝对值指数法11距离法 其中为适当选取的参数，它使得0，经常采用的距离有（1）绝对距离（）.（2）欧式距离：（）（3）Chebishov距离

8、：（）.12主观评分法：请有实际经验者直接对与的相似程度评分，作为的值。上述方法究竟选哪一种，需要根据问题的性质及应用方便来选择。第三步：进行模糊聚类1基于模糊等价矩阵聚类方法一般来说。上述模糊矩阵是一个模糊相似矩阵，不一定具有等价性，即不一定是模糊等价矩阵。这可以通过模糊矩阵的褶积将其转化为模糊等价阵，具体方法如下： 计算,直到满足这时模糊矩阵便是一个模糊等价矩阵。记。将按由大到小的顺序排列，从1开始，沿着由大到小的次序依次取，求的相应的截阵，其中元素为1的表示将其对应的两个变量（或样品）归为一类，随着的变小，其合并的类越来越多，最终当时，将全部变量（或样品）归为一个大类。按值画出聚类的谱系

9、图2直接聚类法所谓直接聚类法是指：在建立模糊相似矩阵之后，不去求传递闭包，直接从相似矩阵出发，求得聚类图。其步骤如下：（1）取1（最大值），对每个作相似类:=|,即将满足的与放在一类，构成相似类。相似类与等价类的不同之处是，不同的相似类可能有公共元素，即可出现=,=,.此时只要将有公共元素的相似类合并，即可得1水平上的等价分类。（2）取为次大值，从中直接找出相似程度为的元素对（，）（即），相应的将对应于1的等价分类中所在类与所在类合并，将所有这些情况合并后，即得对应的等价分类。（3）取为第三大值，从中直接找出相似程度为的元素对（，）（即），类似的将对应于的等价分类中所在类与所在类合并，将所有这

10、些情况合并后，即得对应的等价分类。(4)依次类推，直到合并到成为一类为止。直接聚类法与传递闭包法所得的结果是一致的，直接聚类法要明显简单一些，下面再介绍直接聚类法的图形化方法，即最大树法。所谓最大树法，就是画出以被分类元素为顶点，以相似矩阵的元素为权重的一棵最大的树，取定，去掉权重低于的枝，得到一个不连通的图，各个连通的分支便构成了在水平上的分类。下面介绍求最大树的Kruskal法设，先画出所有顶点从模糊相似矩阵中按从大到小的顺序依次画枝，并标上权重，要求不产生圈，直到所有顶点连通为止，这就得到一棵最大树（最大树可以不唯一）。 上述两个聚类方法各有优劣，使用传递闭包法分类，当矩阵阶数较高时，手

11、工计算量大，但在计算机上还是容易实现的，因此，人们还是乐于使用它。当矩阵阶数不高时，直接聚类法比较直观，也便于操作，适合推广使用。最佳阙值的确定在模糊聚类分析中，对于各个不同的，可得到不同的分类，从而形成一种动态聚类图，这对全面了解样本的分类情况是比较形象和直观的。但许多实际问题需要选择某个阙值的问题。现介绍下面两种方法。1. 按照实际需要，在动态聚类图中，调整的值以得到适当的分类，而不需要事先准确地估计好样本应分为几类。当然，也可由具有丰富经验的专家结合专业知识来确定阙值，从而得出在水平上的等价分类。2. 用统计量确定最佳值 设论域为样本空间（样本总数为），而每个样本有个特征（即由试验或观察

12、得到的个数据）；（）（）。于是，得到原始数据矩阵，如下表所示样本指 标 1 2 · · · · · · · · · · 其中，（），称为总体样本的中心向量。设对应于值的分类数为，第类的样本数为，第类的样本记为：，第类的聚类中心为向量（，, ），其中，为第个特征向量的平均值： （） 作统一量其中为与的距离，为第类样本与中心的距离，称式（*）为统一量。它的分子表征类与类之间的距离，分母表征类样本间的距离。因此，值越大，说明分类越合理，对应统一值最大的阙值为最佳值。（二）模型实例分析 例： 设某地区设

13、置有11个雨量站，其分布图见图51，10年来各雨量站所测得的年降雨量列入表51中。现因经费问题，希望撤销几个雨量站，问撤销那些雨量站，而不会太多的减少降雨信息？ 图1表1年序号12763241594132922583113031752433202251287349344310454285451402307470319243329056347950222122032041123242462322432812673102733152853273525291311502388330410352267603290292646615822417816420350232024027835072583274

14、324013613813014134021994218453365357452384420482228360316252915827141030828341020117943034218510324406235520442520358343251282371应该撤销那些雨量站，涉及雨量站的分布，地形，地貌，人员，设备等众多因素。我们仅考虑尽可能地减少降雨信息问题。一个自然的想法是就10年来各雨量站所获得的降雨信息之间的相似性，对全部雨量站进行分类，撤去“同类”（所获降雨信息十分相似）的雨量站中“多余”的站。问题求解 假设为使问题简化，特作如下假设（1） 每个观测站具有同等规模及仪器设备；（2）

15、 每个观测站的经费开支均等；具有相同的被裁可能性。分析：对上述撤销观测站的问题用基于模糊等价矩阵的模糊聚类方法进行分析，原始数据如上。求解步骤：1利用相关系数法，构造模糊相似关系矩阵，其中 其中，1，2，,11。 ，1，2，,11。用C语言编程计算出模糊相似关系矩阵，具体程序如下#include<stdio.h>#include<math.h>double r1111; double x11;void main() int i,j,k; double fenzi=0,fenmu1=0,fenmu2=0,fenmu=0;int year1011=276,324,159,4

16、13,292,258,311,303,175,243,320,251,287,349,344,310,454,285,451,402,307,470,192,433,290,563,479,502,221,220,320,411,232,246,232,243,281,267,310,273,315,285,327,352,291,311,502,388,330,410,352,267,603,290,292,466,158,224,178,164,203,502,320,240,278,350,258,327,432,401,361,381,301,413,402,199,421,453,3

17、65,357,452,384,420,482,228,360,316,252,158,271,410,308,283,410,201,179,430,342,185,324,406,235,520,442,520,358,343,251,282,371;for(i=0;i<11;i+) for(k=0;k<10;k+) xi=xi+yearki; xi=xi/10;for(i=0;i<11;i+)for(j=0;j<11;j+) for(k=0;k<10;k+) fenzi=fenzi+fabs(yearki-xi)*(yearkj-xj); fenmu1=fen

18、mu1+(yearki-xi)*(yearki-xi); fenmu2=fenmu2+(yearkj-xj)*(yearkj-xj); fenmu=sqrt(fenmu1)*sqrt(fenmu2);rij=fenzi/fenmu; fenmu=fenmu1=fenmu2=fenzi=0; for(i=0;i<11;i+) for(j=0;j<11;j+) printf("%6.3f",rij); printf("n");getchar();得到模糊相似矩阵1.000 0.839 0.528 0.844 0.828 0.702 0.995 0

19、.671 0.431 0.573 0.712 0.839 1.000 0.542 0.996 0.989 0.899 0.855 0.510 0.475 0.617 0.572 0.528 0.542 1.000 0.562 0.585 0.697 0.571 0.551 0.962 0.642 0.568 0.844 0.996 0.562 1.000 0.992 0.908 0.861 0.542 0.499 0.639 0.607 0.828 0.989 0.585 0.992 1.000 0.922 0.843 0.526 0.512 0.686 0.584 0.702 0.899 0.697 0.908 0.922 1.000 0.726 0.455 0.667 0.596 0.511 0.995 0.855 0.571 0.861 0.843 0.726 1.000 0.676 0.489 0.587 0.719 0.671 0.510 0.551 0.542 0.526 0.455 0.676 1.0