决策分析和应用.docx

《决策方法与应用》第一次上机练习题参看区间数多属性决策课件用Matlab编程实现：基于加权平均算子分别采用三种不同的可能度定义对方案进行排序；基于射影法对方案进行排序；基于TOPSIS法对方案进行排序；用基于欧式距离和射影的组合决策方法对方案进行排序。（1）clear;clc;r=[0.214 0.220 0.166 0.178 0.184 0.190; 0.206 0.225 0.220 0.229 0.182 0.191; 0.195 0.204 0.192 0.198 0.220 0.231; 0.181 0.190 0.195 0.205 0.185 0.195; 0.175 0.184 0.193 0.201 0.201 0.211];w=[0.4 0.4 0.4 0.4 0.2 0.2];[m,n]=size(r);%%加权规范化矩阵 z=zeros(m,2);for j=1:m for i=1:2:n z(j,1)=z(j,1)+w(1,i)*r(j,i); z(j,2)=z(j,2)+w(1,i+1)*r(j,i+1); endendz %%基于几何模型求可能定义P1=zeros(m,m);for i=1:m for j=1:m al=z(i,1);au=z(i,2);bl=z(j,1);bu=z(j,2); if au=bl P1(i,j)=0; elseif al=bu P1(i,j)=1; elseif al=blbl=bubu=au P1(i,j)=(au-0.5*(bl+bu))/(au-al); elseif bl=alal=auau=bu P1(i,j)=(0.5*(al+au)-bl)/(bu-bl); elseif al=blbl=auau=bu P1(i,j)=((au-bl)/(au-al))*((au-bl)/(bu-bl))*0.5; else bl=alal=bubu=au; P1(i,j)=(al-bl)/(bu-bl)+((bu-al)/(bu-bl))*(0.5*(bu-al)/(au-al)+(au-bu)/(au-al)); end endendP1 %%可能度判断矩阵p1=sum(P1,2)%%可能度定义2公式求可能度定义P2=zeros(m,m);for i=1:m for j=1:m al=z(i,1);au=z(i,2);bl=z(j,1);bu=z(j,2); S_overall=(bu-bl)*(au-al); % rectangle_ABCD_S if au=bl S_over=0; elseif al=bu S_over=S_overall; elseif bl=al al=bu bu=au S_under=0.5*(bu-al)^2; % triangle_tsC_S S_over=S_overall-S_under; % polygon_ABtsD_S elseif bl=al al=au au=bu S_over=0.5*((al-bl)+(au-bl))*(au-al); % trapezoid_ABCD_S elseif al=bl bl=au au=bu S_over=0.5*(au-bl)^2; % triangle_tsC_S elseif al=bl bl=bu bu=au S_over=0.5*((au-bu)+(au-bl))*(bu-bl); end P2(i,j)=S_over/S_overall; endendP2 %%可能度判断矩阵P2=sum(P2,2)%%可能度定义3公式求可能度定义P3=zeros(m,m);for i=1:m for j=1:m al=z(i,1);au=z(i,2);bl=z(j,1);bu=z(j,2); la=au-al; lb=bu-bl; A=max(au-bl,0); B=min(la+lb,A); P3(i,j)=B/(la+lb); endendP3 %%可能度判断矩阵P3=sum(P3,2)输出结果：z = 0.1888 0.1972 0.2068 0.2198 0.1988 0.2070 0.1874 0.1970 0.1874 0.1962P1 = 0.5000 0 0 0.5831 0.6296 1.0000 0.5000 0.9998 1.0000 1.0000 1.0000 0.0002 0.5000 1.0000 1.0000 0.4169 0 0 0.5000 0.5417 0.3704