·692· 智能系统学报 第11卷 发生器]为例，构建步骤如下： temp=0: 1)给定各规则2个条件对象的参数，即Ex、 for i=1:N En1、He1、Ex2、En2、He2; Ensl(1,i),Ens2(1,i)binormrand Enl, 2)由以上2组参数计算出单个规则对应的二 En2,Hel,He2): 维云发生器，MATLAB程序如下： end %二维云的合成并画出其三维图像 for i=1:N clear: q:(X1-Exl)"2/Ensl(1,i)2: cle: d=(X2-Ex2)2/Ens2(1,i)2: %输入单规则推理两个条件的参数 Y(1,i)=exp(-(1/2)*(q+d): Exl=input('pleaseinputExl='): temp=temp+Y(1,i): Ex2=input('pleaseinputEx2='): end Enl=input('pleaseinputEnl='): Y=temp/Y: En2=input('pleaseinputEn2='): 4)在计算出的u中选择最大山1和次大的山2， Hel=input('pleaseinputHel='): 它们对应的规则即为激活强度最大的两条规则， He2=input('pleaseinputHe2='): MATLAB程序图如下： N=input('pleaseinputN='): [u,X]=sot(U):%将所有规则计算出来的隶属度 temp=0: 由小到大排列 %生成二维条件云并计算(X1X2)时的隶属度 %选择隶属度最大的U和次大的U2,说明这两条 for i=1:N 规则的激活强度最大 Ensl(1,i),Ens2(1,i)binormrand Enl, ul=u(M): En2,Hel,He2): u2=u(M-1): [xI(1,i)X2(1,i)binormrand Exl,Ex2, numl=IX(M) Ensl(1,i),Ens2(1,i)): num2=IX(M-1): end fprintf('触发强度最大规则为第%i条ln',numl); for i=1:N fprintf('触发强度最大规则为第%i条\n',num2): q=(X1(1,i)-Exl)2/Ensl(1,i)2: 5)利用Y条件云由U,和U2分别计算出2组 d=(X2(1,i)-Ex2)2/Ens2(1,i)2: 云滴，选择距离最小的2个云滴，通过逆向云发生器 Y(1,i)=exp(一(1/2)*(q+d)): 计算推理恶意度的期望和方差，MATLAB程序如下： end %构造这两条规则对应的后件云发生器(Y条件云) plot3(x1,X2,Y,'.'). Y1(1),Y1(2)=ycloud(d(numl,7),d(numl,8). 生成二维隶属云如图5所示。 d(numl,9),ul,N); [Y2(1),Y2(2)ycloud (d(num2,7),d(num2, 1.0 ,08 8),d(num2,9),u2,N); 0.6 %从Y1和Y2中分别选取一点使两点距离最小， 0.41 0.2 %由ycloud可知第一个大于第二个 0 0 a=abs(Y1(1)-Y2(2)): 60 b=abs(Y2(1)一Y1(2)): 条件2参数 0 200 50100150200 条件1参数 if a>b Z1=Y2(1): 图5二维隶属云 zul=u2; Fig.5 Two dimensional cloud z2=Y1(2); 3)输入待测条件参数X,和X2,激活每一条规 zu2=ul; 则，根据每一条规则条件参数生成的二维云发生器计 else 算得到激活强度，即隶属度山；。MATLAB程序如下： z1=Y1(1): %……单规则云发生器程序…% zul=ul; function[Y onerule Exl,Enl,Hel,Ex2,En2,He2, z2=Y2(2); X1,X2,N) zu2=u2;发生器［１３］为例，构建步骤如下： １）给定各规则 ２ 个条件对象的参数，即 Ｅｘ１ 、 Ｅｎ１ 、Ｈｅ１ 、Ｅｘ２ 、Ｅｎ２ 、Ｈｅ２ ； ２）由以上 ２ 组参数计算出单个规则对应的二 维云发生器，ＭＡＴＬＡＢ 程序如下： ％二维云的合成并画出其三维图像 ｃｌｅａｒ： ｃＩｃ： ％输入单规则推理两个条件的参数 Ｅｘｌ ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＥｘｌ ＝ ’）： Ｅｘ２ ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＥｘ２ ＝ ’）： Ｅｎｌ ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＥｎｌ ＝ ’）： Ｅｎ２ ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＥｎ２ ＝ ’）： Ｈｅｌ ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＨｅｌ ＝ ’）： Ｈｅ２ ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＨｅ２ ＝ ’）： Ｎ＝ ｉｎｐｕｔ（’ｐｌｅａｓｅｉｎｐｕｔＮ＝ ’）： ｔｅｍｐ ＝Ｏ： ％生成二维条件云并计算（Ｘ１ Ｘ２）时的隶属度 ｆｏｒ ｉ ＝ １：Ｎ ［Ｅｎｓｌ（ １， ｉ）， Ｅｎｓ２ （ １， ｉ）］ ＝ ｂｉｎｏｒｍｒａｎｄ （ Ｅｎ１， Ｅｎ２，Ｈｅｌ，Ｈｅ２）： ［Ｘ１（ １，ｉ），Ｘ２ （ １，ｉ）］ ＝ ｂｉｎｏｒｍｒａｎｄ （ Ｅｘｌ，Ｅｘ２， Ｅｎｓｌ（１，ｉ），Ｅｎｓ２（１，ｉ））： ｅｎｄ ｆｏｒ ｉ ＝ １：Ｎ ｑ ＝ （Ｘ１（１，ｉ）一 Ｅｘｌ）＾２ ／ Ｅｎｓｌ（１，ｉ）＾２： ｄ ＝ （Ｘ２（１，ｉ）－Ｅｘ２）＾２ ／ Ｅｎｓ２（１，ｉ）＾２： Ｙ（１，ｉ）＝ ｅｘｐ（一（１ ／ ２）∗（ｑ＋ｄ））： ｅｎｄ ｐｌｏｔ３（ｘ１，Ｘ２，Ｙ，’．’）． 生成二维隶属云如图 ５ 所示。 图 ５ 二维隶属云 Ｆｉｇ．５ Ｔｗｏ ｄｉｍｅｎｓｉｏｎａｌ ｃｌｏｕｄ ３）输入待测条件参数 Ｘ１ 和 Ｘ２ ，激活每一条规 则，根据每一条规则条件参数生成的二维云发生器计 算得到激活强度，即隶属度 ｕｉ。 ＭＡＴＬＡＢ 程序如下： ％……………单规则云发生器程序……………％ ｆｕｎｃｔｉｏｎ［Ｙ］ ＝ ｏｎｅｒｕｌｅ（Ｅｘｌ，Ｅｎｌ，Ｈｅｌ，Ｅｘ２，Ｅｎ２，Ｈｅ２， Ｘ１，Ｘ２，Ｎ） ｔｅｍｐ ＝Ｏ： ｆｏｒ ｉ ＝ １：Ｎ ［Ｅｎｓｌ（ １， ｉ）， Ｅｎｓ２ （ １， ｉ）］ ＝ ｂｉｎｏｒｍｒａｎｄ （ Ｅｎ１， Ｅｎ２，Ｈｅｌ，Ｈｅ２）： ｅｎｄ ｆｏｒ ｉ ＝ １：Ｎ ｑ：（Ｘ１－Ｅｘｌ）＾２ ／ Ｅｎｓｌ（１，ｉ）＾２： ｄ ＝ （Ｘ２－Ｅｘ２）＾２ ／ Ｅｎｓ２（１，ｉ）＾２： Ｙ（１，ｉ）＝ ｅｘｐ（－（１ ／ ２）∗（ｑ＋ｄ））： ｔｅｍｐ ＝ ｔｅｍｐ＋Ｙ（１，ｉ）： ｅｎｄ Ｙ＝ ｔｅｍｐ ／ Ｙ： ４）在计算出的 ｕｉ 中选择最大 ｕ１ 和次大的 ｕ２ ， 它们对应的规则即为激活强度［１４］ 最大的两条规则， ＭＡＴＬＡＢ 程序图如下： ［ｕ，ＩＸ］ ＝ ｓｏｒｔ（Ｕ）：％将所有规则计算出来的隶属度 由小到大排列 ％选择隶属度最大的 Ｕｌ 和次大的 Ｕ２，说明这两条 规则的激活强度最大 ｕ１ ＝ ｕ（Ｍ）： ｕ２ ＝ ｕ（Ｍ－１）： ｎｕｍ１ ＝ ＩＸ（Ｍ）： ｎｕｍ２ ＝ ＩＸ（Ｍ－１）： ｆｐｒｉｎｔｆ（’触发强度最大规则为第％ｉ 条＼ｎ’，ｎｕｍ１）； ｆｐｒｉｎｔｆ（’触发强度最大规则为第％ｉ 条＼ｎ’，ｎｕｍ２）； ５）利用 Ｙ 条件云由 Ｕ１ 和 Ｕ２ 分别计算出 ２ 组 云滴，选择距离最小的 ２ 个云滴，通过逆向云发生器 计算推理恶意度的期望和方差，ＭＡＴＬＡＢ 程序如下： ％构造这两条规则对应的后件云发生器（Ｙ 条件云） ［Ｙ１（１），Ｙ１（２）］ ＝ ｙｃｌｏｕｄ（ｄ（ｎｕｍｌ，７），ｄ（ｎｕｍｌ，８）， ｄ（ｎｕｍｌ，９），ｕｌ，Ｎ）； ［Ｙ２（ １），Ｙ２ （ ２）］ ＝ ｙｃｌｏｕｄ （ ｄ （ ｎｕｍ２，７），ｄ （ ｎｕｍ２， ８），ｄ（ｎｕｍ２，９），ｕ２，Ｎ）； ％从 Ｙ１ 和 Ｙ２ 中分别选取一点使两点距离最小， ％由 ｙｃｌｏｕｄ 可知第一个大于第二个 ａ ＝ ａｂｓ（Ｙ１（１）－Ｙ２（２））： ｂ ＝ ａｂｓ（Ｙ２（１）一 Ｙ１（２））： ｉｆ ａ＞ｂ Ｚ１ ＝Ｙ２（１）； ｚｕｌ ＝ ｕ２； ｚ２ ＝Ｙ１（２）； ｚｕ２ ＝ ｕｌ； ｅｌｓｅ ｚ１ ＝Ｙ１（１）； ｚｕｌ ＝ ｕｌ； ｚ２ ＝Ｙ２（２）； ｚｕ２ ＝ ｕ２； ·６９２· 智 能 系 统 学 报 第 １１ 卷