Matlab模糊数学建模教程

Matlab模糊数学建模教程 经典数学是适应力学、天文、物理、化学这类学科的需要而发展起来的，不可能不带有这些学科固有的局限性。这些学科考察的对象，都是无生命的机械系统，大都是界限分明的清晰事物，允许人们作出非此即彼的判断，进行精确的测量，因而适于用精确方法描述和处理。而那些难以用经典数学实现定量化的学科，特别是有关生命现象、社会现象的学科，研究的对象大多是没有明确界限的模糊事物，不允许作出非此即彼的断言，不能进行精确的测量。清晰事物的有关参量可以精确测定，能够建立起精确的数学模型。模糊事物无法获得必要的精确数据，不能按精确方法建立数学模型。实践证明，对于不同质的矛盾，只有用不同质的方法才能解决。传统方法用于力学系统高度有效，但用于对人类行为起重要作用的系统，就显得太精确了，以致于很难达到甚至无法达到。 精确方法的逻辑基础是传统的二值逻辑，即要求符合非此即彼的排中律，这对于处理清晰事物是适用的。但用于处理模糊性事物时，就会产生逻辑悖论。如判断企业经济效益的好坏时，用“年利税在100万元以上者为经济效益好的企业”达，否则，便是经济效益不好的企业。根据常识，显而易见:“比经济效益好的企业年利税少1元的企业，仍是经济效益好的企业”，而不应被划为经济效益不好的企业。这样，从上面的两个结论出发，反复运用经典的二值逻辑，我们最后就会得到，“年利税为0者仍为经济效益好的企业”的悖论。类似的悖论有许多，历史上最著名的有“罗素悖论”。它们都是在用二值逻辑来处理模糊性事物时产生的。 客观实际中存在众多的模糊性事物和现象，促使人们寻求建立一种适于描述模糊事物和现象的逻辑模式。模糊集合理论便是在这种形势下应运而生的。模糊方法的逻辑基础是连续值逻辑，它是建立在[0，1]上的。如若我们把年利税在100 ，那末，相比之下，年万元以上者的属于“经济效益好”的企业的隶属度规定为1 利税少1元的企业，属于“经济效益好”的企业的隶属度就应相应减少一点，比如为0.99999，依此类推，企业的年利税每减少1元，它属于“经济效益好”的企业的隶属度就要相应减少一点。这样下去，当企业的年利税为0时，它属于“经济效益好”的企业的隶属度也就为0了，显然，模糊方法的这种处理方式，是符合于人们的认识过程的，连续值逻辑是二值逻辑的合理推广。 现代科学发展的总趋势是，从以分析为主对确定性现象的研究，进到以综合为主对不确定性现象的研究。各门科学在充分研究本领域中那些非此即彼的典型现象之后，正在扩大视域，转而研究那些亦此亦彼的非典型现象。自然科学不同 学科之间，社会科学不同学科之间，自然科学和社会科学之间，相互渗透的趋势日益加强，原来截然分明的学科界限一个个被打破，边缘科学大量涌现出来。随着科学技术的综合化、整体化，边界不分明的对象，亦即模糊性对象，以多种多样的形式普遍地、经常地出现在科学的前沿。 模糊集合理论自诞生以来，获得了长足的发展，每年全世界发表的研究论文的数量，以指数级速度增长。研究范围从开始时的模糊集合，发展为模糊数、模糊代数、模糊测度、模糊积分、模糊规划、模糊图论、模糊拓扑„„等众多的分枝。 和模糊集合理论的发展速度相比，模糊技术的应用虽稍迟一步，但也取得了令人可喜的进展。自1980年第一例应用模糊技术的产品问世以来，有关这方面的研究报告已逾7000多篇，制造出近千种模糊产品，如计算机、电饭煲、摄像机、微波炉、洗衣机、空调器等。如日本松下公司研制的智能化家用空调器，可根据 模式识别 ?2-1模式识别及识别的直接方法 在日常生活中生活中，经常需要进行各种判断、预测。如图象文字识别、故障(疾病)的诊断、矿藏情况的判断等，其特点就是在已知各种标 准类型前提下，判断识别对象属于哪个类型的问题。这样的问题就是模式识别。 一、模糊模式识别的一般步骤 模式识别的问题，在模糊数学形成之前就已经存在，传统的作法主要用统计方法或语言的方法进行识别。但在多数情况下，类型常可用模糊集表示，用模糊数学的方法进行识别是更为合理可行的，以模糊数学为基础的模式识别方法称为模糊模式识别。 模式识别主要包括三个步骤: 第一步:提取特征，首先需要从识别对象中提取与识别有关的特征，并度量这些特征，设分别为每个特征的度量值，于是每个识别对象x就对应一个向量，这一步是识别的关键，特征提取不合理，会影响识别效果。 第二步:建立标准类型的隶属函数，标准类型通常是论域的模糊集，xi是识别对象的第i个特征。 第三步:建立识别判决准则，确定某些归属原则，以判定识别对象属于哪一个标准类型。常用的判决准则有最大隶属度原则(直接法)和择近原则(间接法)两种。 二、最大的隶属度原则 若标准类型是一些表示模糊概念的模糊集，待识别对象是论域中的某一元素(个体)时，往往由于识别对象不绝对地属于某类标准类型，因而隶属度不为1，这类问题人们常常是采用称为“最大隶属度原则”的方法加以识别，这种方法(以及下面的“阈值原则”)是处理个体识别问题的，称为直接法。 最大隶属度原则:设是n个标准类型，，若 则认为x0相对隶属于Ai所代表的类型。 例1 通货膨胀识别问题 通货膨胀状态可分成五个类型:通货稳定;轻度通货膨胀;中度通货膨胀;重度通货膨胀;恶性通货膨胀.以上五个类型依次用R(非负实数域,下同)上的模糊集A1,A2,A3,A4,A5 表示,其隶属函数分别为 1, 其中对表示物价上涨x%。问时，分别相对隶属于哪种类型, 解 ， ， ， ， 由最大隶属原则，应相对隶属于A2，即当物价上涨8%时，应视为轻度通货膨胀;，应相对隶属于A5，即当物价上涨40%时，应视为恶性通货膨胀。 三、阈值原则 在使用最大隶属度原则进行识别中，还会出现以下两种情况，其一是有时待识 中每一个隶属程度都相对较低，这时说明模糊别对象x0关于模糊集 集合对元素x不能识别;其二是有时待识别对象x关于模糊集 中若干个的隶属程度都相对较高，这时还可以缩小x的识别范围，关于这两种情况有如下阈值原则。 阈值原则:是n个标准类型，]为一阈 值(置信水平)令 若则不能识别，应查找原因另作分析。 若且有，则判决x0相对地属于 例2 三角形识别问题 我们把三角形分成等腰三角形I,直角三角形R, 正三角形E,非典型三角形T,这四个标准类型,取定论域 这里A,B,C是三角形三个 现给定，，x0对上述四个标准类型的隶属度为: 由于x0关于I,R的隶属程度都相对高，故采用阈值原则，取，因 ，，按阈值原则，x0相对属于I?R,即x0可识别为等腰直角三角形。 例3 癌细胞识别 在癌细胞识别问题中细胞分成四个标准类型，即:癌细胞(M)，重度核异质细胞(N)，轻度核异质细胞(R)，正常细胞(T)。 选取表征细胞状况的七个特征: x1:核面积, x2:细胞面积, x5:核 (a正常核面积核桨比例置 D:核 上述定义中的模糊集A的隶属函数为。另两个模糊集B、1 2121212 C的隶属函数类似定义。 给定待识别细胞，设x0的核面积等七个特征值为据此可算出00012M(x0)、N(x0)、R(x0)、T(x0)，最后按最大隶属度原则识别。 例4 冬季降雪量预报 其中x1是丰镇地区若干年6、7月份气温的平均值，为方差，实际预报时取 其中x2是若干年秋季初霜日的平均值，a2是经验参数，实际预报时取x2=17(即9月17日)，a2=10(即9月10日)。 -sinx 270 取论域，“冬雪大”可以表示为论域X上的模糊集C，其隶属函数为: ?(A2(x2)?A3(x3)) 采用阈值原则，取阈值，测定当年气候因子。计算C(x)，若则预报当年冬季“多雪”，否则预报“少雪”。 用这一方法对丰镇1959~1970年间隔12年作了预报，除1965年以外均报对，历史拟合率为11/12。 ?2-2 贴近度与模式识别的间接方法 一、贴近度 表示两个模糊集接近程度的数量指标，称为贴近度，其严格的数学定义如下: 定义1 设映射 满足下列条件: ， ， (3) 若A,B满足 有 则称映射N为F(U)上的贴近度，称N(A B)为A与B的贴近度。 贴近度的具体形式较多，以下介绍几种常见的贴近度公式 (1) Hamming 贴近度 1n 或 (2)Euclid贴近度 或 (3)格贴近度 定义7 映射 ?B)c，(或?B)c]) 2 称为格贴近度，称Ng(A,B)为A与B格贴近度。其中， 称为A与B的 (称为A与B的外积) 若，则 A? 值得注意的是，这里的格贴近度是通过定义来规定的，事实上，格贴近度不满 (1)，即，但是，当时，格贴足定义1中 近度满足定义1的(1)-(3)。另外格贴近度的计算很方便，且用于表示相同类型模糊度的贴近度比较有效，所以在实际应用中也常选用格贴近度来反映模糊集接近程度。 还有许多贴近度，这里不在一一介绍。 贴近度主要用于模糊识别等具体问题，以上介绍的贴近度表示式各有优劣，具体应用时，应根据问题的实际情况，选用合适的贴近度。 在模式识别问题中，各标准类型(模式)一般是某个论域X上的模糊集，用模式识别的直接方法(最大隶属度原则、阈值原则)解决问题时，其识别对象是论域X中的元素。另有一类识别问题，其识别对象也是X上的模糊集，这类问题可以用下面的择近原则来识别判决。 择近原则:已知n个标准类型A1、A2、„、，为待识别的对象，N为F(X)上的贴近度，若 则认为B与Ai最贴近，判定B属于Ai一类。 例5 岩石类型识别 岩石按抗压强度可以分成五个标准类型:很差(A1)、差(A2)、较好(A3)、好(A4)、很好(A5)。它们都是上的模糊集，其隶属函数如下(图2-1) 图 2-1 (kg/cm2) 其它 其它 今有某种岩体，经实测得出其抗压强度为X上的模糊集B，隶属函数为(图2-3)。 图 2-3 其它 试问岩体B应属于哪一类。 计算B与的格贴近度，得: 按择近原则，B应属于A3类，即B属于“较好”类(A3类)的岩石。 例6 小麦亲本识别 在小麦杂交育种过程中，亲本选择是关键。现有五种类型的小麦亲本，它们是: A1:早熟型，A2:矮杆型，A3:大粒型， A4:高肥丰产型，A5:中肥丰产型。 判断小麦亲本类型的主要依据是以下五种性状特征: x1:抽穗期，x2:株高，x3:有效穗数， x4:主穗粒数，x5:百粒重。 第i种类型亲本的第j个特征，是模糊集Aij，这些模糊集除A11(早熟型的抽 穗期)与 A22(矮杆型的株高)外，其余都是中间型的正态分布模糊集。为简单计，将正态分布函数 展开，取前两项作它的近似值，则有 e a)2 于是Aij的隶属函数可表示为: ij 其它 ，A22的隶属函数取为偏小值型: 而A11 iiii2 为确定隶属函数中的参数值，在熟知的标准类型中，每类型选出k个新本为，分别计算各样本的第j个特征的均值xijl及方差，取 ijl 1k2 以上参数值见表(2-1) 表 2-1 现有一待识对象B，它的第j个特征Bj是中间型正态分布模糊集，隶属函数可近似表示为: 式中参数值见表(2-2) 。 其它 计算识别对象B的第j个特征与第i种标准类型对应特征Aij的格贴近度 并定义第i种标准类型Ai与识别对象B的贴近度为: 计算结果列于表(2-3) 表 2-3 表(2-3)的最后一行为B与各标准类型的贴近度。由于B与A5的贴近度最高(0.96)，故判定识别对象B为A5代表的类型，即B为中肥丰产类型的亲本。 例7 遥感土地复盖类型分类 遥感是根据不同的地物对电磁波谱有不同的响应这一原理，来识别土地复盖的类型。空间遥感的一个象元相当于地面0.45公倾地物的综合。遥感图象识别分类中，要涉及不少模糊概念，例如，“以红松为主的针叶林”就是一个没有明确界线的模糊概念。这是遥感本身的特性决定的。因此用模糊数学的方法对遥感图象进行识别分类应该是行之有效的方法。 美国爱达荷大学R.C.Heller 教授指出，国际上当以水体、沙地、森林、城镇、作物、干草作为分类单位(即标准类型)时，空间遥感的分类精度可达83.93%甚至更高。但当分类单位深入到更小的土地复盖单元时，精度就不理想了。 现在将分类单位细分阶段为以下五种标准类型: A1:公路，A2:村庄农田，A3:红松为主的针叶林， A4:阔、针混交林，A5:白桦林。 对于多波段遥感技术，假设采用p个波段，则每一地物对应一个p维数据向量 。1975年1月22日美国发射LandSat-2，提供了MSS-4，5，6，7这四个波段的数据，故有。取论域 其中x1,x2,x3,x4分别为象元对应于MSS-4，5，6，7各波段的光谱强度。于是五种标准类型可表为X上的模糊集。 由于各波段光谱强度是正态分布模糊集，故第i个标准类型的(j+3)波段光谱强度的隶属函数为: 定义第i种标准类型Ai为: 因而 其中aij为若干个第i种类型第(j+3)个波段光谱强度的均值，为方差， -4) 东北凉水林场的这些参数值见表(2 表 2-4 设B为识别对象，定义Ai与B的贴近度为: (1) 4 其中 Ng(Aij,Bj)= 1 ?Bj)c] (2) 2 表 2-5 按及(A? c ccc 1 2 )2 ] (3-2 6) (这里aj与是Bj的均值与方差)。 现有东北凉水林场空间遥感象元(待识别对象)五个，按(1)与(2)计算它们与五个标准类型的贴近度，计算结果在表(2-5)按择近原则进行识别判决，准确率100%。 例8 雷达识别 现有n个雷达类，每个雷达类可用发射频率、脉冲重复频率、脉冲宽度等特征来刻画，假设共有j个特征，第i类雷达的第j个特征可以取nij个值。由于保密的需要及信号环境的日益复杂，这些特征及其取值都带有一定的模糊性。设第i类雷达的k个特征为 类雷达的第j个特征取值为 m其隶属函数为 ， 中间型柯西分布，即 设X为待识别对象，它的k个特征为的第j个特征Xj的隶属函数也取中间型柯西分布: 采用格贴近度，令 m 则dij为识别对象X的第j个特征与i类雷达第j个特征贴近程度的度量。 一般情况可令 (di是各dij的加权平均值，权系数aj表示j个特征的重要性程度)di可作为识别对象X与第i类雷达总贴近的度量。根据di的大小可判定X属于何类雷达，但是，由于权系数aj的确定有一定的模糊性，Aij及Xj的隶属函数的确定带有一定的主观性，从而导致贴近度dij有一定的模糊性。因此对aj及dij进行模糊化处理，设 mm 这里Aj，Dij都是模糊数(见第五章)，取。 c 令 Di的隶属函数为 则Di为识别对象X与第i类雷达的贴近程度的模糊测度。 ，令 为得到X所属雷达类别的确切判决，类似于阈值法则，给定水平值 若 且i0唯一，则判定X为i0类雷达; 若 且，则判定X为i1类雷达。 用上述方法(将权系数及贴近度模糊化)，经上千次仿真试验，比传统的贴近度及线性加弘平均法，误判率有所下降。 第三章 模糊规划 ?3-1 模糊极值 一、有界函数的模糊极值 设 (R为实数集) 是有界函数，求函数f(x)的普通极值问题是求x使 满足上式的x为f(x)在X上的最大值点，f(x)为最大值，最大值点不一 定唯一. 设f(x)的一切最大值点的集合为 称Mf为f(x)的优越集.当时，函数在x处取到最大值f(x)，x使f(x)达到最优.当时，f(x)虽不是最大值，但对不同的x，f(x)与最大值的差异有所不同,也就是说，对于不属于Mf的x，它们的“优越性”程度有所不同，为了反映X中各点不同的优越程度，将优越集Mf模糊化，并利用它将极值模糊化. 定义1设是有界函数，定义Mf的隶属函数为 称Mf为f的无条件模糊优越集称f(Mf)的f的无条件模糊极大值.这里 ，它的求属函数按扩张原理为 约定 注 (1)当为f(x)的极大点，即时，当为f(x)的极小点，即时充分必要条件是 (2)当时，当 时，当 时， 因此,f(Mf)(y)反映了在模糊意义下,y对f的模糊数大值的求属程度. 例1 设， 定义则 并且 于是 又 故 f的无条件模糊极小集mf定义为的无条件极大集，显然有 且有， 所有极小集mf是极大集Mf的余集. 二、模糊约束下有界函数的模糊极值 设:是有界函数，，考虑f在C约束下的最大值问题，这是一个模糊规划问题，求解这个问题意味着既要最大限度地满足约束，又要最大限度地达到理想目标，为此定义如下: 定义2 设目标函数是有界函数，是模糊约束，令 这里的Mf是定义1中f的无条件模糊优越集，称D为f在C约束下的条件模糊优越集，称f(D)为f在C约束下的条件模糊极大值.它们的求属函数分别为: 求解目标函数f(x)在模糊约束C下的条件极大值有如下三个步骤: (1)求无条件模糊优越集Mf (2)求条件模糊优越集求条件最佳决策，即选择x,使 就是所求的条件极大点，就是在模糊约束C下的条件极大值. 例2采区巷道布置是矿井开拓中的重要 (2)采煤机械化程度高; (3)采区生产系统十分完善; (4)可靠性好; (5)煤炭损失率低; (6)巷道掘进费用尽可能低. 上述问题,实际上就是一个模糊约束下的条件极值问题，我们可以把(1),(5)作为模糊约束，而把(6)作为目标函数. 设某矿井的采区巷道布置有六种方案可供选择，即X=,x1(方案?), x2 (方案?), x3(方案?), x4(方案?), x5(方案?), x6(方案?),. 经过对六种方案进行审议，评价后，将其结果列于表1 将表 1 中的语言真值(评价结果)转化为各模糊约束集 的隶属度转化的对应关系如下: 对C1, C2, C3, C4而言，对应关系为: 对 C5而言，对应关系为 将表1中的巷道掘进费用目标函数f用公式 计算出，因此得表2 其值语言与隶属函数转换表2 计算模糊判决集D为 按列求最小) 由maxD( 根据最大求属度原则，方案四最优 例3 在某种食品中投放某种调味剂，每公斤食品中的含量设为x克，对顾客爱 好作调查统计，得爱好函数为 对于使爱好函数值越大的x值，所制产品越畅销，因而收益越大，但是由于成本核算等等原因，对x值需要进行限制，这种限制集合的边界是模糊的，即x的约束条件为一模糊集A，其隶属函数为 试确定合理的剂量x，使得在接受约束的条件下，获得最优收益. 解 这是一个规划问题，分三步进行. (1) 求无条件模糊优越集Mf，由于 ， 220 令，得又当时，，时，，因而 ，因此 (2) 求条件模糊优越集Af 其中x满足方程 (3) 选择x，使 A， 即A对目标f的可能度为45.93%，而要实现这种可能性，应选择调味剂的最佳剂量为2.085克. 需要说明的是，在本例中如果将约束条件确切化，以A的核[0,1]为约束，这是一个普通规划问题，所得结论是选择最佳剂量为1克.从约束条件看，已是100%遵守，但所能达到的最高目标相对整个目标函数来说是很低的，由，说明相对整个目标来说，其优越程度仅达24.6%.如果把条件放松为模糊约束条件A，且适当降低A(x)的水平，却可以获得较好的目标值.如例中的结果，当时，从接受约束条件来看虽仅达45.9%，但目标函数的优越程度也升到了45.9%，从而提高了整体优化水平.由于在实际问题中，约束条件往往不是绝对的，有一定的伸缩性，模糊规划的思想就是利用这点灵活性，兼顾目标函数与约束条件综合地选 择最优方案. 例4 植物的种植密度与产量有密切的关系.已知某种杉树的种植密度与产量V的关系如下: V表示每公顷土地产出木材的体积.这里表示每公顷土地上种植的棵数， 现有一片杉树森林，其密度不均匀，估计大约是三千”.试估计该森林每公顷木材最高产量. 解 设C表示“大约是三千”这一模糊，C的隶属函数为 估计木材产量的问题，就是求在C的约束下函数f的模糊条件极大值.为此先求有界函数f的无条件模糊优越集.因，，所 以 f在约束条件C下的条件模糊优越集为: 条件模糊极值为f(Cf)，其隶属函数为: ，即满足条件 为求条件最佳决策 的 注意到Mf的隶属函数曲线是单调降的，而C是正态分布模糊集，f在约束C下的模糊最佳决策(即模糊条件极大点)，是方程 103 的两个根当中的较小者，解之得 由Cf(2130)可知，时，接受约束的程度为46.9%，同时，相对于整体目标函数，优越程度也是46.9%. 由米3)可知，该森林每公顷木材最高产量估计为334.74(米3). ?3-2 模糊线性规划 一、普通线性规划 普通线性规划的一般形式为 目标函数 约束条件 矩阵表达形式 其中 线性规划问题的标准形式 二、模糊线性规划 在实际问题中，有时线性规划的约束条件带有模糊性，这就是解谓的模糊线性 3-1) 规划，其模型为 ( 这是表示一种弹性约束，可读作“近似小于等于”.“近似小于等于”是一个模糊概念，可以用一个模糊集来表示它 表示第i个约束的左 这一事实，当 ~边表达式，模糊集Di表示 时，完 全接受约束，应有适当选择一个伸缩系数di，约定当时，不认为 ，这时应有当 时，Di(x)应从1下降到0，表示约束程度降低.为了简 单可行，Di(x)规定如下: n设 ，对每一个约束，相应地有X中一个 模糊渠Di与之对应，它的隶属函数为 其中di是适当选择的常数，叫做伸缩指标，，这样一来,我们将弹性约束转化成模糊约束,再令就将全部 约束条件转化成一个模糊约束. 当时,D退化为普通约束集D，模糊约束条件中退化为 模糊线性规划的模型简记为 ~ (3-2) 约束的弹性必然导致目标的弹性，为将目标函数模糊化，先求解普通线性规划问题: 满足 (3-3) 以及 满足 (3-4) 其中称为(3-2)的伸缩指标向量. 设Z0是(3-32)的最优值，Z1是(3-4)的最优值.Z0所满足的约束条件为，对应的模糊约束若适当降低模糊约束的隶属度D(x)，可以相应提高目标函数值Z,Z1所满足的约束条件已放到最宽，对应的模糊约束D(x)也接近于0.于是目标函数的弹性可表示为为此构造模糊目标集 其隶属函数为 其中 时，，要提高目 由模糊目标的上述隶属函数可知，当 标函数值使之大于Z0.就必须降低D(x).为了兼顾目标与约束，可采用模糊决策 ，最佳决策为，满足 为 若令则有 于是求最佳决策的问题，就转化为求普通线性规划问题: 0 即 求解上述普通规划问题，可得 最佳决策 n) 目标函数值 jxj. 例5:求解模糊线性规划问题 3-5) ( -6) 取伸缩指标 解 (一)解普通线性规划 最优解最优值 (二)解普通线性规划 最优解 (三) 解普通线性规划 解 这个线性规划采用大M法 原线性规划改写为 充分大的正数) ? 从而(3-4)的最优值 例6某企业根据市场信息及自身生产能力，准备开发甲、乙两种系列产品.甲种 .据测算，系列产品最多大约能生产400套，乙种系列产品最多大约能生产250套甲种产品每套成本3万元，每套获纯利润7万元;乙种系列产品每套成本2万元，每套获纯利润3万元.生产甲、乙两种系列产品的资金总投入大约不能超过1500万元.在上述条件下，如何安排两种系列产品的生产，才能使企业获利最大, 解 设甲种系列产品生产x1套，乙种系列产品生产x2套，则 目标 约束(3-7) x 设约束条件(1)、(2)、(3)的伸缩系数分别取为(元)，(套)， (套).为将目标函数模糊化，解经典线性规划问题 使 (4) 用单纯形法求解，得，， 再解经典线性规划问题 max 解得 ，， 于是 将z0、d0、d1、d2、d3代入(3-5)，将原问题经为经典线性规划问题: 使 上述线性规划问题最优解为，，因此安 排甲种系列产品403套、乙种系列产品159套(取整数)时，能获得最大利润，最大利润为: 万元 对比经典线性规划问题(4)，利润提高43.75万元，这是因为甲种系列产品403 套比400套多3套;乙种系列产品生产159套比150套多9套，这 是在伸缩指标允许范围模糊数学及其应用，四川大学出版社，1994. 2、张跃、邹寿平、宿芬，模糊数学方法及其应用，煤炭工业出版社，1992. 3、李安贵、张志宏、段凤英、模糊数学及其应用，冶金工业出版社，1994. 4、杨雄、李崇文，模糊数学和它的应用，天津科技出版社，1993. 5、唐林炜、樊铭渠伪模糊度的公理化定义及余滋伪模糊度，系统理论与 (8). 实践，1999 6、唐林炜、高国成、樊铭渠、无解的模糊关系方程的最优近似解，模糊系统与数学，1999(3).