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

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

contourf(x1,y1,z1,10); %作出该地区综合污染指数等值线

for a=1:5

ans=find(data1(:,5)==a);

Cave(a)=sum(C(ans))/length(ans); %分别以各功能区为单位求污染物污染指数均值 end

hold on;

qu=data1(:,5);

new=[x,y,qu]; %导入水平面坐标与分区编号

b1(:,1)=find(new(:,3)==1);

b2(:,1)=find(new(:,3)==2);

b3(:,1)=find(new(:,3)==3);

b4(:,1)=find(new(:,3)==4);

b5(:,1)=find(new(:,3)==5);

plot(new(b1,1),new(b1,2),'+',new(b2,1),new(b2,2),'o',new(b3,1),new(b3,2),'*',new(b4,1),new(b4,2),'x',new(b5,1),new(b5,2),'d'); %作出该地区不同功能区的标示

三、程序p3.m: p3.m:

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

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

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

data(:,1)=[]; %删掉编号所在列，方便数据处理

data(:,9)=[5.595065696872188;2.050686879701597;3.888298563664709;19.48276529716754;16.93280375665379;22.724400402542187;3.037349192403835;278.4314856449954;326.29517448187767;2.168613213605158;5.463305948625636;3.881272052582544;2.633239515018349;9.46334805067366;5.8640039047304;11.86692132078126;9.495010113989254;5.780916469665128;4.600638095010868;18.078460014979655;2.765714198188635;75.06194642916017;6.317656062533044;3.175982075050526;1.002105755941413;13.027235962605822;4.273270498136975;2.938675903149045;31.277250875574946;17.799553790448112;7.203477918187483;4.402250215125141;5.114637582756907;7.411040169305963;6.071798817309512;30.197995035158257;4.986040063219011;4.933408478760172;6.851727650580298;9.30454412611158;38.970506299568115;13.714554002857652;7.048582202760745;7.128746796901836;12.785511357350094;3.653523595443423;3.509614909556038;2.054432017198534;18.78500460165682;4.56868763930696;3.086422991477065;3.602471019982397;3.286155698475126;15.197347372190814;2.407994127442934;3.918089639224968;3.021220225505681;3.479181017287154;4.062489944841015;2.976220965345462;39.076789250677294;1.56838923950092;1.869842022929421;2.09594617952183;4.393470573406756;3.198853554555947;2.94663449909761;3.82610699445268;2.540367896060084;1.484074279882942;2.053723421933565;1.890947323462811;3.340249275117196;1.646648526651114;2.029100519052495;0.684083342613396;2.041502463345344;0.849988961415488;2.870055943613625;1.738665489003439;1.456952671918555;3.975748198588527;1.45770399240739;6.313275792533697;2.508112770368477;3.02604549399094;1.262863159123092;1.532824738176672;1.339330145544494;5.147561