2011年全国大学生数学建模竞赛A题——城市表层土(19)

2011年全国大学生数学建模竞赛A题——城市表层土壤重金属污染分析：黄俊彬,李恪睿,陈泽君

93544035;1.463541689554989;5.161772125881541;1.491086375193536;2.080879390395252;1.490208609331122;1.333045146213944;1.740511663057965;2.25579729384667;];%由程序p2得到的各样本点综合污染指数

for a=1:9

data(:,a)=data(:,a)/data(1,a); %消除量纲并使数据交于同一点，用于灰色关联度分析 end

ans1=find(data1(:,5)==1);

newdata1=data(ans1,:); %找出1类区所属数据建立新矩阵

x01=newdata1(:,9); %以综合污染指数作为灰色关联度分析的初值序列 h1=size(newdata1,1);

for a=1:h1

d1(a,:)=abs(x01(a)-newdata1(a,1:8)); %求差序列

end

mind1=min(min(d1));

maxd1=max(max(d1));

rou=0.5; %分辨系数

for a=1:size(d1,1)

for b=1:size(d1,2)

ksi1(a,b)=(mind1+rou*maxd1)/(d1(a,b)+rou*maxd1); %计算关联系数序列

end

end

for a=1:size(ksi1,2)

r(1,a)=1/size(ksi1,1)*sum(ksi1(:,a)); %求关联度

end

ans2=find(data1(:,5)==2);

newdata2=data(ans2,:);

x02=newdata2(:,9);

h2=size(newdata2,1);

for a=1:h2

d2(a,:)=abs(x02(a)-newdata2(a,1:8));

end

mind2=min(min(d2));

maxd2=max(max(d2));

rou=0.5;

for a=1:size(d2,1)

for b=1:size(d2,2)

ksi2(a,b)=(mind2+rou*maxd2)/(d2(a,b)+rou*maxd2);

end

end

for a=1:size(ksi2,2)

r(2,a)=1/size(ksi2,1)*sum(ksi2(:,a));

end

ans3=find(data1(:,5)==3);

newdata3=data(ans3,:);