模糊数学323课件

第第 3 章章模糊聚类分析模糊聚类分析2.1 模糊矩阵模糊矩阵 定义定义1 设设R=(rij)mn，若，若0rij1，则称，则称R为为模模糊矩阵糊矩阵.当当rij只取只取0或或1时，称时，称R为为布尔布尔(Boole)矩矩阵阵.当模糊方阵当模糊方阵R=(rij)nn的对角线上的元素的对角线上的元素rii都都为为1时，称时，称R为为模糊自反矩阵模糊自反矩阵.定义定义2 设设A=(aij)mn,B=(bij)mn都都是模糊矩阵，是模糊矩阵，相等相等：A=B aij=bij；包含包含：AB aijbij；并并：AB=(aijbij)mn；交交：AB=(aijbij)mn；余余：Ac=(1-aij)mn.模糊矩阵的并、交、余运算性质模糊矩阵的并、交、余运算性质幂等律：幂等律：AA=A，AA=A；交换律：交换律：AB=BA，AB=BA；结合律：结合律：(AB)C=A(BC),(AB)C=A(BC)；吸收律：吸收律：A(AB)=A，A(AB)=A；分配律：分配律：(AB)C=(AC)(BC)；(AB)C=(AC)(BC)；0-10-1律：律：AO=A，AO=O；AE=E，AE=A；还原律：还原律：(Ac)c=A；对偶律：对偶律：(AB)c=AcBc,(AB)c=AcBc.模糊矩阵的合成运算与模糊方阵的幂模糊矩阵的合成运算与模糊方阵的幂 设设A=(aik)ms，B=(bkj)sn，定义模糊矩阵，定义模糊矩阵A 与与B 的合成为：的合成为：A B=(cij)mn，其中其中cij=(aikbkj)|1ks.模糊方阵的幂模糊方阵的幂 定义：若定义：若A为为 n 阶方阵，定义阶方阵，定义A2=A A，A3=A2 A，Ak=Ak-1 A.合成合成()运算的性质：运算的性质：性质性质1：(A B)C=A (B C)；性质性质2：Ak Al=Ak+l，(Am)n=Amn；性质性质3：A (BC)=(A B)(A C)；(BC)A=(B A)(C A)；性质性质4：O A=A O=O，I A=A I=A；性质性质5：AB，CD A C B D.(AB)C(A C)(B C)(AB)C (A C)(B C)模糊矩阵的转置模糊矩阵的转置 定义定义 设设A=(aij)mn,称称AT=(aijT)nm为为A的转的转置矩阵，其中置矩阵，其中aijT=aji.转置运算的性质：转置运算的性质：性质性质1：(AT)T=A；性质性质2：(AB)T=ATBT，(AB)T=ATBT；性质性质3：(A B)T=BT AT；(An)T=(AT)n；性质性质4：(Ac)T=(AT)c；性质性质5：AB AT BT.模糊矩阵的模糊矩阵的 -截矩阵截矩阵 定义定义7 设设A=(aij)mn,对任意的对任意的 0,1，称，称A=(aij()mn,为模糊矩阵为模糊矩阵A的的 -截矩阵截矩阵,其中其中 当当aij 时，时，aij()=1；当；当aij 时，时，aij()=0.显然，显然，A的的 -截矩阵为布尔矩阵截矩阵为布尔矩阵.对任意的对任意的 0,1，有，有性质性质1：AB A B；性质性质2：(AB)=A B，(AB)=A B；性质性质3：(A B)=A B；性质性质4：(AT)=(A )T.2.2 模糊关系模糊关系 与模糊子集是经典集合的推广一样，模糊关与模糊子集是经典集合的推广一样，模糊关系是普通关系的推广系是普通关系的推广.设有论域设有论域X，Y，X Y 的一个模糊子集的一个模糊子集 R 称称为从为从 X 到到 Y 的的模糊关系模糊关系.模糊子集模糊子集 R 的隶属函数为映射的隶属函数为映射R:X Y 0,1.并称隶属度并称隶属度R(x,y)为为(x,y)关于模糊关系关于模糊关系 R 的相关程度的相关程度.特别地，当特别地，当 X=Y 时，时，称之为称之为 X 上各元素之上各元素之间的间的模糊关系模糊关系.模糊关系的运算模糊关系的运算 由于由于模糊关系模糊关系 R就是就是X Y 的一个模糊子集，的一个模糊子集，因此模糊关系同样具有模糊子集因此模糊关系同样具有模糊子集的运算及性质的运算及性质.设设R，R1，R2均为从均为从 X 到到 Y 的的模糊关系模糊关系.相等相等：R1=R2 R1(x,y)=R2(x,y)；包含包含：R1 R2 R1(x,y)R2(x,y)；并并：R1R2 的隶属函数为的隶属函数为(R1R2)(x,y)=R1(x,y)R2(x,y)；交交：R1R2 的隶属函数为的隶属函数为(R1R2)(x,y)=R1(x,y)R2(x,y)；余余：Rc 的隶属函数为的隶属函数为Rc(x,y)=1-R(x,y).(R1R2)(x,y)表示表示(x,y)对模糊关系对模糊关系“R1或或者者R2”的相关程度，的相关程度，(R1R2)(x,y)表示表示(x,y)对对模糊关系模糊关系“R1且且R2”的相关程度，的相关程度，Rc(x,y)表示表示(x,y)对模糊关系对模糊关系“非非R”的相关程度的相关程度.模糊关系的矩阵表示模糊关系的矩阵表示 对于有限论域对于有限论域 X=x1,x2,xm和和Y=y1,y2,yn，则，则X 到到Y 模糊关系模糊关系R可用可用mn 阶阶模糊矩阵表示，即模糊矩阵表示，即R=(rij)mn，其中其中rij=R(xi,yj)0,1表示表示(xi,yj)关于模糊关于模糊关系关系R 的相关程度的相关程度.又若又若R为布尔矩阵时为布尔矩阵时,则关系则关系R为普通关系为普通关系,即即xi 与与 yj 之间要么有关系之间要么有关系(rij=1),要么没有关系要么没有关系(rij=0).例例 设身高论域设身高论域X=140,150,160,170,180(单位：单位：cm),体重论域体重论域Y=40,50,60,70,80(单位：单位：kg),下表给出了身高与体重的模糊关系下表给出了身高与体重的模糊关系.405060708014010.80.20.10150 0.810.80.20.1160 0.20.810.80.2170 0.10.20.810.818000.10.20.81模糊关系的合成模糊关系的合成 设设 R1 是是 X 到到 Y 的关系的关系,R2 是是 Y 到到 Z 的关系的关系,则则R1与与 R2的合成的合成 R1 R2是是 X 到到 Z 上的一个关系上的一个关系.(R1R2)(x,z)=R1(x,y)R2(y,z)|yY 当论域为有限时，模糊关系的合成化为模糊当论域为有限时，模糊关系的合成化为模糊矩阵的合成矩阵的合成.设设X=x1,x2,xm,Y=y1,y2,ys,Z=z1,z2,zn，且，且X 到到Y 的的模糊模糊关系关系R1=(aik)ms，Y 到到Z 的的模糊模糊关系关系R2=(bkj)sn，则，则X 到到Z 的的模糊模糊关系可表示为关系可表示为模糊模糊矩阵的合成：矩阵的合成：R1 R2=(cij)mn，其中其中cij=(aikbkj)|1ks.模糊关系合成运算的性质模糊关系合成运算的性质性质性质1：(A B)C=A (B C)；性质性质2：A (BC)=(A B)(A C)；(BC)A=(B A)(C A)；性质性质3：(A B)T=BT AT；性质性质4：A B，C D A C B D.注：注：(1)合成合成()运算关于运算关于()的分配律不成立的分配律不成立,即即(AB)C (A C)(B C)(2)这些性质在有限论域情况下这些性质在有限论域情况下,就是模糊矩就是模糊矩阵合成运算的性质阵合成运算的性质.设 为U到U的模糊关系，对 及 ，有：则称 为模糊传递关系。2.3 模糊等价矩阵模糊等价矩阵模糊传递关系模糊传递关系 为模糊传递关系的充要条件是：2.3 模糊等价矩阵模糊等价矩阵模糊等价关系模糊等价关系 若模糊关系若模糊关系R是是X上上各元素之间的各元素之间的模糊关系，模糊关系，且满足：且满足：(1)(1)自反性：自反性：R(x,x)=1；(2)(2)对称性：对称性：R(x,y)=R(y,x)；(3)(3)传递性：传递性：R2 R,则称则称模糊关系模糊关系R是是X上上的一个的一个模糊等价关系模糊等价关系.当论域当论域X=x1,x2,xn为有限时为有限时,X 上的上的一个一个模糊等价关系模糊等价关系R就是模糊等价矩阵就是模糊等价矩阵,即即R满足：满足：I R(rii=1)RT=R(rij=rji)R2R.R2R(rikrkj)|1kn rij).模糊等价矩阵的基本定理模糊等价矩阵的基本定理 定理定理1 若若R具有自反性具有自反性(IR)和传递性和传递性(R2R),则则 R2=R.定理定理2 若若R是模糊等价矩阵是模糊等价矩阵,则则对任意对任意 0,1，R 是等价的是等价的Boole矩阵矩阵.0,1，ABA B；(AB)=A B；(AT)=(A)T模糊相似关系模糊相似关系 若模糊关系若模糊关系 R 是是 X 上各元素之间的上各元素之间的模糊关模糊关系，且满足：系，且满足：(1)自反性：自反性：R(x,x)=1；(2)对称性：对称性：R(x,y)=R(y,x)；则称则称模糊关系模糊关系 R 是是 X 上的一个上的一个模糊相似关系模糊相似关系.当论域当论域X=x1,x2,xn为有限时，为有限时，X 上的上的一个一个模糊相似关系模糊相似关系 R 就是模糊相似矩阵，即就是模糊相似矩阵，即R满满足：足：(1)自反性：自反性：I R(rii=1)；(2)对称性：对称性：RT=R(rij=rji).模糊相似矩阵的性质模糊相似矩阵的性质 定理定理1 若若R 是模糊相似矩阵，则对任意的自是模糊相似矩阵，则对任意的自然数然数 k，Rk 也是模糊相似矩阵也是模糊相似矩阵.定理定理2 若若R 是是n阶模糊相似矩阵，则存在一阶模糊相似矩阵，则存在一个最小自然数个最小自然数 k(kn)，对于一切大于，对于一切大于k 的自然的自然数数 l，恒有，恒有Rl=Rk，即，即Rk 是模糊等价矩阵是模糊等价矩阵(R2k=Rk).此时称此时称Rk为为R的传递闭包，记作的传递闭包，记作 t(R)=Rk.上述定理表明，任一个模糊相似矩阵可诱导上述定理表明，任一个模糊相似矩阵可诱导出一个模糊等价矩阵出一个模糊等价矩阵.平方法求传递闭包平方法求传递闭包 t(R)：RR2R4R8R162.4 模糊聚类分析模糊聚类分析数据标准化数据标准化 设论域设论域X=x1,x2,xn为被分类对象为被分类对象,每每个对象又由个对象又由m个指标表示其形状个指标表示其形状:xi=xi1,xi2,xim,i=1,2,n于是于是,得到原始数据矩阵为得到原始数据矩阵为平移平移 标准差变换标准差变换其中其中平移平移 极差变换极差变换模糊相似矩阵建立方法模糊相似矩阵建立方法相似系数法相似系数法-夹角余弦法夹角余弦法相似系数法相似系数法-相关系数法相关系数法绝对值倒数法绝对值倒数法其中M选取合适的值，使0rij 1绝对值指数法绝对值指数法距离法距离法rij=1 c d(xi,xj)其中其中c为适当选取的参数为适当选取的参数.海明距离海明距离欧氏距离欧氏距离切比雪夫距离切比雪夫距离d(xi,xj)=|xik-xjk|,1kmBoole矩阵法：矩阵法：定理：设定理：设 R 是论域是论域 X=x1,x2,xn上的上的一个相似的一个相似的 Boole 矩阵，则矩阵，则 R 具有传递性具有传递性(当当R是等价是等价Boole矩阵时矩阵时)矩阵矩阵 R 在任一排列下的在任一排列下的矩阵都没有形如矩阵都没有形如的特殊子矩阵的特殊子矩阵.Boole矩阵法的步骤如下：矩阵法的步骤如下：(1)通过数据标准化以及前面的相似矩通过数据标准化以及前面的相似矩阵的建立方法，建立模糊等价矩阵。阵的建立方法，建立模糊等价矩阵。(2)求模糊等价矩阵的求模糊等价矩阵的 -截矩阵截矩阵R ；(3)若若R 在某一排列下的矩阵有形如在某一排列下的矩阵有形如的特殊子矩阵的特殊子矩阵,则将则将R 中上述特殊形式子矩阵的中上述特殊形式子矩阵的0改为改为1，直到在任一排列下，直到在任一排列下R 中不再产生上述特殊中不再产生上述特殊形式子矩阵为止形式子矩阵为止.经过4次复合可得到等价关系进行 分割，并分别取分类图例：某三个家庭，共有人口16名，各家庭成员之间有（血缘）关系。16个人各有自己的照片，且混在一起。现在要求一个不认识这三个家庭成员的人根据上述像片确定这些成员之间的相似程度。用 对该模糊等价关系进行分割因此可识别出4个性质不同的类型：1,6,8,13,16,2,5,7,11,14,3,4,9,10,12,15例：有5个地区受到近期一次地震而造成损害利用余弦幅度法，得到如下关系：经两次复合得到等价关系最佳分类的确定最佳分类的确定 在模糊聚类分析中，对于各个不同的在模糊聚类分析中，对于各个不同的 0,10,1，可得到不同的分类，从而形成，可得到不同的分类，从而形成一种动态聚类图，这对全面了解样本分类一种动态聚类图，这对全面了解样本分类情况是比较形象和直观的情况是比较形象和直观的.但在许多实际问题中，需要给出样本但在许多实际问题中，需要给出样本的一个具体分类，这就提出了如何确定最的一个具体分类，这就提出了如何确定最佳分类的问题佳分类的问题.设设X=(xij)nm为为n个元素个元素m个指标的原始数据个指标的原始数据矩阵矩阵.为总体样本的中心向量为总体样本的中心向量.对应于对应于 值的分类数为值的分类数为r，第，第 j 类的样本数为类的样本数为nj，第，第 j 类的样本标记为类的样本标记为第第 j 类样本的中心向量为类样本的中心向量为作作F-统计量：统计量：如果满足不等式如果满足不等式FF (r-1,n-r)的的F值不值不止一个，则可根据实际情况选择一个满意的分类，止一个，则可根据实际情况选择一个满意的分类，或者进一步考查差或者进一步考查差(F-F )/F 的大小，从较大的大小，从较大者中找一个满意的者中找一个满意的F值即可值即可.实际上，最佳分类的确定方法与聚类方法无实际上，最佳分类的确定方法与聚类方法无关，但是选择较好的聚类方法，可以较快地找到关，但是选择较好的聚类方法，可以较快地找到比较满意的分类比较满意的分类.模糊模糊C均值聚类均值聚类 给定数据集：给定数据集：X=（x1，x2，.，xn），其），其中每个样本包含中每个样本包含s个属性。模糊聚类就是将个属性。模糊聚类就是将X划划分为分为c个类（个类（2cn），），V=（v1,v2,.,vn）是）是c个聚类中心。在模糊划分中，每个样本个聚类中心。在模糊划分中，每个样本不能严格地划分为某一类，而是以一定的隶属不能严格地划分为某一类，而是以一定的隶属度属于某一类。其划分准则如下：度属于某一类。其划分准则如下：其中，其中，为第为第k个样本属于第个样本属于第i类的隶属度。类的隶属度。为样本为样本k与第与第i个聚类中心的距离。个聚类中心的距离。（1）选择合适的聚类中心数）选择合适的聚类中心数c和和m值。并确定初值。并确定初始隶属度矩阵。始隶属度矩阵。（2）计算聚类中心）计算聚类中心V：（3）修正隶属度矩阵：）修正隶属度矩阵：（4）对给定的）对给定的 若若 ，则算法，则算法终止，否则终止，否则 ，转向步骤（，转向步骤（2）。）。模糊模糊C均值聚类步骤均值聚类步骤 center,U,obj_fcn=fcm(date,cluster_n)参数解释：参数解释：center 最终聚类中心矩阵，每行为一个聚类最终聚类中心矩阵，每行为一个聚类中心的坐标值。中心的坐标值。U 最终的模糊分区矩阵最终的模糊分区矩阵 obj_fcn 在迭代过程中的目标函数值在迭代过程中的目标函数值 Date 要聚类的数据集合，每行为一个样本要聚类的数据集合，每行为一个样本 cluster_n 聚类数，要求聚类数，要求1。模糊模糊C均值聚类均值聚类MATLAB实现实现p经常不断地学习,你就什么都知道。你知道得越多,你就越有力量pStudyConstantly,AndYouWillKnowEverything.TheMoreYouKnow,TheMorePowerfulYouWillBe写在最后Thank You在别人的演说中思考，在自己的故事里成长Thinking In Other PeopleS Speeches，Growing Up In Your Own Story讲师：XXXXXX XX年XX月XX日