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

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

附录 附录

一、程序p1.m:p1.m:

clear;close all; %清空旧数据，关闭旧图像

data1=xlsread('F:\2011高教杯全国赛\cumcm2011Problems中文版\A\cumcm2011A附件_数据.xls','附件1'); %导入题目所给附件1

data2=xlsread('F:\2011高教杯全国赛\cumcm2011Problems中文版\A\cumcm2011A附件_数据.xls','附件2'); %导入题目所给附件2

x=data1(:,2); %储存水平面坐标的x值

y=data1(:,3); %储存水平面坐标的y值

z=data2(:,2:9); %储存8种不同元素的浓度值

for a=1:8

[x1,y1]=meshgrid(0:300:30000,0:200:20000); %建立网格线

z1=griddata(x,y,z(:,a),x1,y1,'v4'); %以“v4”法进行插值运算

figure;

contourf(x1,y1,z1,15); %作出8种不同元素的填充等浓度线 end

二、程序p2.m:p2.m:

clear;close all; %清空旧数据，关闭旧图像

data1=xlsread('F:\2011高教杯全国赛\cumcm2011Problems中文版\A\cumcm2011A附件_数据.xls','附件1'); %导入题目所给附件1

data2=xlsread('F:\2011高教杯全国赛\cumcm2011Problems中文版\A\cumcm2011A附件_数据.xls','附件2'); %导入题目所给附件2

x=data1(:,2); %储存水平面坐标的x值

y=data1(:,3); %储存水平面坐标的y值

s=[3.6 130/1000 31 13.2 35/1000 12.3 31 69]; %储存附件所给背景值

c=data2(:,2:9); %储存8种不同元素的浓度值

c(:,2)=c(:,2)/1000; %统一量纲，附件中Cd和Hg浓度的量纲与其他不同 c(:,5)=c(:,5)/1000;

for a=1:8

p(:,a)=c(:,a)/s(a); %污染物污染指数

end

for a=1:319

pmean(a,1)=sum(p(a,:))/8; %各采样点污染物污染指数的均值

end

for a=1:319

pmax(a,1)=max(p(a,:)); %各采样点污染物污染指数的最大值

end

C=sqrt((pmean.^2+pmax.^2)/2); %采用内梅罗指数法定义综合污染指数

[x1,y1]=meshgrid(0:300:30000,0:200:20000); %建立网格线