障碍因素模型程序

障碍度程序R=[0.120.520.300.06];%4¸öÒ»¼¶Ö¸±êµÄȨÖØW=[0.490.080.320.110.300.500.150.050.290.160.060.491.00];%13¸ö¶þ¼¶Ö¸±êµÄȨÖØF=zeros(1,13);fori=1:4F(1,i)=R(1,1)*W(1,i);endfori=5:8F(1,i)=R(1,2)*W(1,i);endfori=9:12F(1,i)=R(1,3)*W(1,i);endF(1,13)=R(1,4)*W(1,13);X=[30.0360.013911.88456.1638.09342.1868.573395.060.8402230.0580.002664.0039.3483.303092.79109.142474.080.4627430.0520.006460.7694.39406.531365.4364.412151.880.6011530.0700.009114.81245.4619.58251.3414.312410.420.607440.0700.006540.7942.8683.333248.37134.093267.320.4518740.0140.01121.4589.506.6765.211.65111.790.666140.0620.00444.74664.7087.20301.856.68243.680.354120.0100.0096.2268.60120.35.340.40230.540.982240.6630.00439.86134.0321.24138.617.03113.640.206230.2280.008802.36258.50150.062265.4124.7934117.940.32235];xmin=zeros(1,13);xmax=zeros(1,13);fori=1:13xmin(1,i)=X(1,i);xmax(1,i)=X(1,i);endforj=1:13fori=2:10ifX(i,j)xmax(1,j)xmax(1,j)=X(i,j);endendendendxmin;xmax;XX=zeros(10,13);forj=1:13fori=1:10XX(i,j)=(X(i,j)-xmin(1,j))/(xmax(1,j)-xmin(1,j));endendI=zeros(10,13);fori=1:10forj=1:13I(i,j)=1-XX(i,j);endendO=zeros(10,13);fori=1:10sum=0;forj=1:13sum=sum+F(1,j)*I(i,j);endforr=1:13O(i,r)=(I(i,r)*F(1,r))/sum;endendU=zeros(10,4);fori=1:10forj=1:4U(i,1)=U(i,1)+O(i,j);endforj=5:8U(i,2)=U(i,2)+O(i,j);endforj=9:12U(i,3)=U(i,3)+O(i,j);endU(i,4)=O(i,13);endOU