模糊综合评判matlab源程序(共6页)

1、精选优质文档-倾情为你奉上模糊综合评判matlab源程序2009-02-09 10:161.原理模糊综合评判方法即将评价目标看成是由多种因素组成的模糊集合（称为因素集u），再设定这些因素所能选取的评审等级，组成评语的模糊集合（称为评判集v），分别求出各单一因素对各个评审等级的归属程度（称为模糊矩阵），然后根据各个因素在评价目标中的权重分配，通过计算（称为模糊矩阵合成），求出评价的定量解值。它是应用模糊变换原理和最大隶属度原则，对各因素作综合评价的。其原理表示为：B=E×R式中E=e1,e2,. ,ei,.,em,为模糊向量或称模糊变换器，是评价因素集X= x1,x2, ,xi, ,x

2、m 的权重分配。式中R为评价因素X= x1,x2, .,xi, .,xm 与评判集U= u1,u2, ,ui, ,un 构成的模糊关系矩阵。    |R1| |r11 r12 . r1m|R=|R2|= |r21 r22 . r2m|    | |   |                |       

3、60;        |Rn |     |rn2     rn2 . rnm|其中n为评价集合的个数，m为评价因素(或评判指标)的个数。2.程序算法下面是采用环境中的常用的超标加权法计算权重，使用“线性降半阶”函数计算隶属度的matlab程序，供各位参考。clear;clc;a=48. 26. 819. 334. 1032. 262. 2374.15.84 6. 981. 756. 991. 82. 2535.70.7225 29.8525 8

4、95.5325 294.5875 1059.1875 462.9525 2724.5; %a为评价集标准值d=43.49 28.05 737.98 391.12 1025.66 7.41 21342.61 1.82 920.75 636.41 1005.45 74.09 233028.66 8.51 774.99 322.78 1013.96 93.84 20013.01 1.95 897.53 614.44 889.87 123.27 220025.85 5.72 759.59 302.03 1001.96 76.15 19722 1.95 1161.68 1003.73 1077.06 1

5、10.1 30043.81 1.09 820.51 396.61 1004.74 37.04 20195.21 2.92 814.08 419.8 1005.8 31.49 20184.41 2.8 824.65 449.06 998.36 38.28 20473.01 1.58 1220.54 956.14 1244.75 3.91 30716.01 2.43 1791.61 2338.17 1278.08 30.87 43621.2 2.67 1160.54 821.29 1100.82 85.41 29426.61 7.3 865.57 389.31 1065.27 46.51 2244

6、9.82 3.77 1240.77 939.71 1165.24 177.19 324817.64 6.44 884.2 473.49 930.29 218.95 241726.65 7.9 839.5 474.71 941.99 184.18 236325.25 4.74 808.33 486.31 881.01 191.6 221725.45 6.93 852.01 478.37 966.45 182.54 233935.27 18.48 785.11 331.32 979.57 4.04 22452.81 2.31 1601.02 2533.55 486.73 2.47 38014.21

7、 4.86 1815.52 2584.68 963.61 0 439915.23 5.35 813.85 787.16 688.79 205.18 209367.01 36.65 864.23 357.76 1035.8 426.31 260984.65 34.24 892.72 381.19 1060 466.64 273115.43 21.52 898.68 88.47 1061.46 414.48 2483115.8 27 926.5 350.93 1079.49 544.38 3075; %b为待判样品值m,n=size(a);x,y=size(d);TheResultMoHu=;%f

8、or s=1:x    for p=1:n        b(p)=d(s,p);    end    for i=1:n %计算每一个列的平均值        ColAverage(i)=0;        for j=1:m      

9、      ColAverage(i)=ColAverage(i)+a(j,i);        end        ColAverage(i)=ColAverage(i)/m;    end    %    ETotal = 0; %超标加权法计算权重    for j = 1: n&

10、#160;        ETotal = ETotal + (b(j) / ColAverage(j);    end    for i = 1: n        EResult(i) = (b(i) / ColAverage(i) / ETotal; %EResult为计算结果    end    %  

11、0; %使用隶属函数，做预处理对每一列进行排序    SortedMatrix=a;    for j=1:n        for i=1:m            for k=i:m              

12、60; if SortedMatrix(i,j)>SortedMatrix(k,j)                    tmp=SortedMatrix(i,j);                    Sor

13、tedMatrix(i,j)=SortedMatrix(k,j);                    SortedMatrix(k,j)=tmp;                end      

14、      end        end    end    %    %计算单因素隶属度    c=SortedMatrix;    for j = 1 : n              

15、 for i = 1 : m                    for k = 1 : m                        if a(i, j) = c(k, j

16、)                             if k = 1                      

17、           if b(j) < c(k, j)                                     LSDResu

18、lt(j, i) = 1;                                end                &

19、#160;               if b(j) >= c(k, j) & b(j) < c(k + 1, j)                           

20、0;         LSDResult(j, i) = (c(k + 1, j) - b(j) / (c(k + 1, j) - c(k, j);                               

21、 end                                if b(j) >= c(k + 1, j)               

22、;                      LSDResult(j, i) = 0;                         &#

23、160;      end                            end                

24、            if k > 1 & k < m                                if b(j) < c(k - 1

25、, j)                                     LSDResult(j, i) = 0;          

26、0;                     end                             

27、   if b(j) >= c(k - 1, j) & b(j) < c(k, j)                                     LSDResult(j, i) = (b(j) -

28、c(k - 1, j) / (c(k, j) - c(k - 1, j);                                end            &

29、#160;                   if b(j) >= c(k, j) & b(j) < c(k + 1, j)                       

30、;             LSDResult(j, i) = (c(k + 1, j) - b(j) / (c(k + 1, j) - c(k, j);                           &

31、#160;    end                                if b(j) >= c(k + 1, j)           

32、                          LSDResult(j, i) = 0;                     

33、60;          end                            end            &

34、#160;               if k = m                                 if b(j) < c

35、(k - 1, j)                                     LSDResult(j, i) = 0;         

36、0;                      end                            

37、    if b(j) >= c(k - 1, j) & b(j) < c(k, j)                                     LSDResult(j, i) = (b

38、(j) - c(k - 1, j) / (c(k, j) - c(k - 1, j);                                end