数学建模第四次优秀论文研读报告

特等奖论文

1、用到知识

clear

load shuju0.txt load shuju1.txt

n1=load('1.txt') %生活区横纵坐标数据 n2=load('2.txt') %工业区横纵坐标数据 n3=load('3.txt') %山区横纵坐标数据 n4=load('4.txt') %交通区横纵坐标数据 n5=load('5.txt') %公园绿地区横纵坐标数据 plot(n1(:,1),n1(:,2),'*') %生活区位置图 hold on

plot(n2(:,1),n2(:,2),'.') %工业区位置图 hold on

plot(n3(:,1),n3(:,2),'^') %山区位置图 hold on

plot(n4(:,1),n4(:,2),'+') %交通区位置图

hold on

plot(n5(:,1),n5(:,2),'o') %公园绿地区位置图

legend('生活区','工业区','山区','交通区','公园绿地区')

hold on

cx=shuju0(:,2); %x轴坐标 cy=shuju0(:,3); %y轴坐标

cz=shuju1(:,5); %所有重金属浓度数据（As Cd Cr Cu Hg Ni Pb Zn） x=0:73:28654; %插值横坐标 y1=0:73:18449;

y=y1'; %插值纵坐标

z1=griddata(cx,cy,cz,x,y,'cubic') %插值后的重金属浓度值 z=log2(z1/(1.5*35)); %计算插值后的地累积指数 [x,y]=meshgrid(x,y); %自动生成网格

[c,h]=contour(x,y,z) %插值后的地累积指数等高线图 c是等高线矩阵，h是等高线句柄 set(h,'ShowText','on') %把h值写上去

21.8x 1041.61.4-2-42--2-200-22-241-202--4002-40020200.60.40.202200-200.5-400020802204-2120-20444-22024400040240200.84-22460-22220-22200008222201.522.5x 100000-2-4-41.2-20-2生活区0工业区山区交通区0公园绿地区-20-02-42-60020000-20-40-20-2-20-44-20-20-222-22-24-40cx=shuju0(:,2); %x轴坐标 cy=shuju0(:,3); %y轴坐标

cz=shuju1(:,5); %所有重金属浓度数据（As Cd Cr Cu Hg Ni Pb Zn） x=0:73:28654; %插值横坐标 y1=0:73:18449; y=y1'; %插值纵坐标 %z=300.*x+350.*y-385000;

z=griddata(cx,cy,cz,x,y,'cubic') %插值后的重金属浓度值 mesh(x,y,z) h=shuju0(:,4)

xlabel('x'),ylabel('y'); title('浓度');

02-20-20462200-2-2-26202234

浓度20000150001000050000-500021.5x 104310.5y001x2x 104

2、区域污染负荷指数计算