%--------------------------------------------------------------------------
%??
%??%mohu.m
%本程序為模糊綜合評判的程序，人為給定權重
%程序運行，直接輸入命令mohu即可
%本程序可以進行擴展
%對于進行二級模糊綜合評判，可編制mohufun.m函數來實現
%如[yy1qdhqdh1]=mohufun（RLMwXXyy]來實現運行兩次這個函數，并編制相應的m文件，其中放至兩個這樣的函數即可實
%現
%模糊評判的結果說明：
%1:對于所有等級的隸屬度之和為1。
%2：輸出結果隨某個變量的增大而qdh1的結果呈現出單調遞增或遞減
%3:最后的等級出[0.2?0?0?0.8]這種情況是正常的，和可拓評判中的對于等級的距離的概念是不一樣的。
%只表現有變量是屬于0.2這個等級的。
%--------------------------------------------------------------------------

clc
clear?all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%此部分由計算時，人為給定部分，即相關求解條件
%要求計算前輸入部分
R=4;??%影響因素的個數
L=4;??%評判的等級個數
M=29;?%待評價洞段的個數
w=[0.1?0.7?0.1?0.1];%各指標的權重
%可拓評判的經典域，每天列構成一組，共有L組，由左到右依此為1，2，3，4
%評判結果為1，則表示屬于最左邊的經典域，為4則屬于最右邊的經典域
%XX為隸屬函數的范圍，具體可參見論文基于粗集的模糊綜合評判，一個有四個隸屬函數
%前兩個數，確定是一個下降段的隸屬函數，分成三段
%以后依此4上數為一組，確定一個梯形的隸屬函數，分成五段
%最后兩個是單獨一組，為一個上升段的隸屬函數分成三段
%下面的每行有12個數，其分成4個等級，分辦界點數個2?4?4?2?之和為12
%也就是構成四個等級范圍
%XX也是計算的核心數據
xx=[70	90	70	90	110	130	110	130	170	190	170	190
???45	35	45	35	30	20	30	20	17	10	17	10
??1.5	2.5	1.5	2.5	3	4	3	4	4.5	5.5	4.5	5.5
0.25	0.35	0.25	0.35	0.45	0.55	0.45	0.55	0.65	0.75	0.65	0.75];

%pp1為個指標為越大，等級越高，還是越大，等級越小的指標
%0為指標值越大，等級越大
%1為指標值越大，等級越小
%這一點要嚴格執行，XX（i:）從小到大為0XX（i:）從大到小為1
%若將以上的XX第行的順序顛倒，則為pp1=[1?0?1?1?]
pp1=[0?1?0?0];
%各待評價洞段參數
yy=[200	8	6	0.8
200	9.5	6	0.8
200	11	6	0.8
200	12.5	6	0.8
200	14	6	0.8
200	15.5	6	0.8
200	17	6	0.8
200	18.5	6	0.8
200	20	6	0.8
200	21.5	6	0.8
200	23	6	0.8
200	24.5	6	0.8
200	26	6	0.8
200	27.5	6	0.8
200	29	6	0.8
200	30.5	6	0.8
200	32	6	0.8
200	33.5	6	0.8
200	35	6	0.8
200	36.5	6	0.8
200	38	6	0.8
200	39.5	6	0.8
200	41	6	0.8
200	42.5	6	0.8
200	44	6	0.8
200	45.5	6	0.8
200	47	6	0.8
200	48.5	6	0.8
200	50	6	0.8
];
%下升段和下降段的隸屬函數形式
%注意隸屬函數的形式和論文中的隸屬函數形式也要一樣，否則要更改程序中的基些部分
%這里的a指區間的中間點，b指區間的寬度
f1=inline（‘0.5-0.5*sin（（x-a）*pi/b）‘‘a‘‘b‘‘x‘）;
f2=inline（‘0.5+0.5*sin（（x-a）*pi/b）‘‘a‘‘b‘‘x‘）;




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%
%求每個指標對于第一個等級的隸屬度
for?i=1:M
????for?k=1:R
????????j=1;
???if?pp1（k）==0?????%先是求越大，等級越高的指標????
if?yy（ik）????yy1（ikj）=1;
elseif?yy（ik）<=xx（kj+1）
????yy1（ikj）=f1（0.5*（xx（kj）+xx（kj+1））abs（xx（kj）-xx（kj+1））yy（ik））?;
else?yy1（ikj）=0;
end
%再是求值越小，等級越高的指標
%下在和一個相比大于和小于要互換，函數f1和f2也要互換
???elseif?yy（ik）>xx（kj）
????yy1（ikj）=1;
elseif?yy（ik）>=xx（kj+1）
????yy1（ikj）=f2（0.5*（xx（kj）+xx（kj+1））abs（xx（kj）-xx（kj+1））yy（ik））;??
else?yy1（ikj）=0;
end

????????end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%求每個指標對于最后一個等級的隸屬度
?for?i=1:M
????for?k=1:R
????????j=L;
????????if?pp1（k）==0???????%先是求越大，等級越高的指標?
if?yy（ik）????yy1（ikj）=0;
elseif?yy（ik）<=xx（kj*4-4）
????yy1（ikj）=f2（0.5*（xx（kj*4-5）+xx（kj*4-4））abs（xx（kj*4-5）-xx（kj*4-4））yy（ik））;??
else?yy1（ikj）=1;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????5766??2008-03-01?21:54??mohu.m

-----------?---------??----------?-----??----

?????????????????5984????????????????????2