dwrr.zip_Help!

clc% 清屏 clear all;%删除workplace变量 close all;%关掉显示图形窗口 A=[ 0.037 0.031 0.098 0.030 0.032 0.058 0.029 0.028 0.028 0.025 0.034 0.108 0.034 0.027 0.067 0.034 0.032 0.060 0.029 0.024 0.059 0.028 0.028 0.094 0.027 0.025 0.039 0.018 0.028 0.052 0.014 0.020 0.031 0.011 0.023 0.061 0.040 0.035 0.067 0.039 0.042 0.065 0.047 0.033 0.033 0.036 0.045 0.090 0.035 0.030 0.076 0.036 0.041 0.068 0.011 0.036 0.040 0.025 0.039 0.106 0.027 0.021 0.056 0.025 0.030 0.041 0.000 0.000 0.000 0.010 0.025 0.080 0.026 0.023 0.047 0.023 0.029 0.048 0.011 0.019 0.024 0.013 0.027 0.081 0.026 0.023 0.061 0.017 0.028 0.050 0.007 0.021 0.037 0.014 0.027 0.070 0.027 0.019 0.051 0.010 0.030 0.034 0.005 0.008 0.015 0.010 0.020 0.050 0.025 0.018 0.038 0.009 0.023 0.024 0.005 0.014 0.011 0.009 0.015 0.037 0.036 0.024 0.061 0.013 0.028 0.031 0.011 0.025 0.019 0.015 0.026 0.061 0.028 0.017 0.055 0.016 0.027 0.034 0.006 0.013 0.017 0.011 0.021 0.062 0.032 0.025 0.050 0.023 0.033 0.042 0.005 0.023 0.022 0.014 0.025 0.060 0.035 0.019 0.049 0.051 0.040 0.046 0.014 0.028 0.029 0.016 0.025 0.061 0.037 0.027 0.074 0.081 0.054 0.052 0.017 0.037 0.034 0.019 0.041 0.082 0.044 0.031 0.128 0.100 0.062 0.094 0.013 0.036 0.077 0.019 0.047 0.136 0.029 0.029 0.102 0.019 0.041 0.071 0.003 0.048 0.067 0.005 0.040 0.111 0.033 0.026 0.085 0.019 0.031 0.059 0.006 0.031 0.060 0.005 0.030 0.094 0.031 0.023 0.095 0.009 0.023 0.054 0.004 0.032 0.053 0.003 0.022 0.078 0.045 0.023 0.129 0.017 0.029 0.084 0.017 0.023 0.057 0.022 0.034 0.158 0.032 0.022 0.105 0.011 0.026 0.066 0.009 0.015 0.050 0.014 0.026 0.123 0.033 0.021 0.068 0.020 0.022 0.047 0.021 0.025 0.029 0.018 0.022 0.085 0.033 0.024 0.075 0.018 0.029 0.043 0.016 0.020 0.028 0.016 0.027 0.071 0.056 0.035 0.100 0.036 0.038 0.048 0.042 0.033 0.033 0.033 0.040 0.125 0.059 0.041 0.101 0.045 0.051 0.089 0.077 0.046 0.054 0.032 0.047 0.118 0.065 0.049 0.148 0.081 0.052 0.127 0.054 0.055 0.084 0.037 0.000 0.157 0.051 0.030 0.068 0.029 0.039 0.052 0.013 0.032 0.054 0.010 0.032 0.079 0.044 0.021 0.101 0.010 0.026 0.081 0.030 0.018 0.062 0.014 0.014 0.128 0.066 0.024 0.062 0.011 0.032 0.049 0.017 0.018 0.022 0.017 0.027 0.079 0.035 0.018 0.053 0.008 0.023 0.049 0.008 0.009 0.021 0.009 0.004 0.045 0.052 0.023 0.038 0.046 0.025 0.042 0.019 0.013 0.020 0.024 0.014 0.060 0.045 0.025 0.030 0.056 0.028 0.038 0.025 0.017 0.020 0.023 0.016 0.050 0.057 0.034 0.063 0.061 0.041 0.050 0.063 0.031 0.030 0.019 0.027 0.082 0.076 0.043 0.093 0.094 0.055 0.081 0.030 0.047 0.049 0.027 0.058 0.035 0.073 0.036 0.089 0.000 0.067 0.085 0.020 0.021 0.058 0.020 0.025 0.000 0.056 0.028 0.077 0.000 0.057 0.073 0.009 0.020 0.040 0.013 0.019 0.000 0.070 0.031 0.051 0.035 0.037 0.045 0.031 0.029 0.025 0.024 0.031 0.092 0.049 0.029 0.061 0.019 0.028 0.044 0.013 0.024 0.033 0.008 0.028 0.056 0.034 0.026 0.041 0.019 0.031 0.041 0.004 0.032 0.018 0.007 0.031 0.051 0.037 0.040 0.096 0.038 0.050 0.073 0.013 0.052 0.032 0.020 0.057 0.130 0.026 0.027 0.058 0.028 0.030 0.044 0.000 0.000 0.000 0.012 0.035 0.047 0.023 0.029 0.065 0.066 0.027 0.050 0.011 0.042 0.035 0.009 0.034 0.045 0.034 0.031 0.079 0.091 0.038 0.066 0.023 0.035 0.031 0.028 0.044 0.100 0.022 0.028 0.050 0.071 0.033 0.040 0.003 0.046 0.029 0.004 0.033 0.038 0.032 0.043 0.063 0.049 0.046 0.056 0.011 0.032 0.015 0.017 0.048 0.070 0.033 0.035 0.078 0.011 0.031 0.057 0.006 0.034 0.034 0.005 0.041 0.064 0.024 0.031 0.057 0.004 0.032 0.055 0.004 0.019 0.037 0.003 0.036 0.032 0.029 0.030 0.040 0.006 0.032 0.040 0.004 0.015 0.024 0.004 0.029 0.021 0.024 0.052 0.109 0.010 0.050 0.078 0.004 0.031 0.033 0.006 0.060 0.099 0.022 0.026 0.052 0.027 0.028 0.046 0.003 0.012 0.031 0.013 0.027 0.035 0.030 0.043 0.068 0.063 0.045 0.069 0.013 0.049 0.031 0.023 0.056 0.000 0.041 0.037 0.075 0.039 0.036 0.044 0.018 0.037 0.020 0.021 0.047 0.064 0.068 0.036 0.060 0.061 0.031 0.047 0.013 0.038 0.032 0.014 0.038 0.054 0.024 0.028 0.041 0.059 0.022 0.033 0.011 0.018 0.022 0.013 0.022 0.042 0.022 0.025 0.037 0.057 0.021 0.027 0.004 0.018 0.013 0.012 0.018 0.029 ]; ca=[0.1,0.335,0.1];%中值 ca1=[0.5341 0.4492 0.5341];%δ的值 B=[0.2 0.6 1;0.67 0.67 1;0.2 0.6 1];%标准分级 a=size(A,1);%A的行数 b=size(A,2);%A的列数 MX=max(A);%A的每列的最大值 MN=min(A);%A的每列的最小值 f1=1;%B的行变量 f2=1;%B的列变量 for j=1:b%将A里面的值标准化 for i=1:a A(i,j)=(A(i,j)-MN(j))/(MX(j)-MN(j)); end end X=zeros(a,4);%建立放置每个城市每天空气质量的值矩阵 for i=1:a%计算每个城市每天的空气质量 h=1; k=1; flag0=1;%算完每个城市每天的空气质量后，跳到下一个城市（h=h+1） for j=1:b for flag=1:3 if A(i,j)<=ca(k) w=0; X(i,h)=X(i,h); break; elseif A(i,j)>ca(k)&&A(i,j)<B(f1,f2) w=1-exp(-(A(i,j)-ca(k)/ca1(k))^2); X(i,h)=X(i,h)+w*A(i,j);%加权 break; else k=k+1;%对应每个中值 f2=f2+1;%对应每个等级值 end end f1=f1+1;%对应每个属性 k=1;%每循环完一次，跳到第一个中值 f2=1;%每循环完一次，又从第一级开始 if flag0==3%判断是否算完了一个城市每天的空气质量 h=h+1; flag0=1; f1=1; else flag0=flag0+1; end end end a=size(X,1);%X的行数 b=size(X,2);%X的列数 h=1; C=zeros(1,b);%建立总体排序矩阵 for k=1:b%将所有城市总体排序 for i=1:a for j=1:b if(A(h)>A(i,j)) C(k)= C(k)+1; else end end h=h+1; end end C % Borda数