2011年数学建模大赛A题金属污染Matlab源程序.doc

2011年数学建模大赛A题金属污染Matlab源程序.doc

《2011年数学建模大赛A题金属污染Matlab源程序.doc》由会员分享，可在线阅读，更多相关《2011年数学建模大赛A题金属污染Matlab源程序.doc（3页珍藏版）》请在冰点文库上搜索。

应使用excel数据源，时间紧迫用txt代替

clear,clc

count=319;%count±íʾȡÑùµãµÄ¸öÊý

%fidheng=fopen（'hengzuobiao.txt','r'）;

%heng=zeros（1,count）;

%[heng,count]=fread（fidheng,inf）;

x=load（'xzuobiao.txt'）;

x=x';

y=load（'yzuobiao.txt'）;

y=y';

z=load（'zzuobiao.txt'）;

z=z';

as=load（'as.txt'）;

as=as';

cd=load（'cd.txt'）;

cd=cd';

cr=load（'cr.txt'）;

cr=cr';

cu=load（'cu.txt'）;

cu=cu';

hg=load（'hg.txt'）;

hg=hg';

ni=load（'ni.txt'）;

ni=ni';

pb=load（'pb.txt'）;

pb=pb';

zn=load（'zn.txt'）;

zn=zn';

gongneng=load（'gongnengqu.txt'）;

gongneng=gongneng';

%plot3（x,y,as）;

holdon

As=ones（size（as'））*as;

[X,Y]=meshgrid（x,y）;

x1=linspace（min（x）,max（x））;

y1=linspace（min（y）,max（y））;

[X1,Y1]=meshgrid（x1,y1）;

%As1=griddata（X,Y,As,X1,Y1）

%As1=interp2（x,y,As,X,Y）;

bianhao=1:

1:

319;

%ԭʼÊý¾Ý´¦Àí

maxas=max（as）;

minas=min（as）;

as=（as-minas）/（maxas-minas）;

maxcd=max（cd）;

mincd=min（cd）;

cd=（cd-mincd）/（maxcd-mincd）;

maxcr=max（cr）;

mincr=min（cr）;

cr=（cr-mincr）/（maxcr-mincr）;

maxcu=max（cu）;

mincu=min（cu）;

cu=（cu-mincu）/（maxcu-mincu）;

maxhg=max（hg）;

minhg=min（hg）;

hg=（hg-minhg）/（maxhg-minhg）;

maxni=max（ni）;

minni=min（ni）;

ni=（ni-minni）/（maxni-minni）;

maxpb=max（pb）;

minpb=min（pb）;

pb=（pb-minpb）/（maxpb-minpb）;

maxzn=max（zn）

minzn=min（zn）

zn=（zn-minzn）/（maxzn-minzn）;

%plot（bianhao,as,bianhao,cd,bianhao,cr,bianhao,cu,bianhao,hg,bianhao,ni,bianhao,pb,bianhao,zn）

%mesh（X1,Y1,As1）;

subplot（4,2,1）;

plot（bianhao,as）;

subplot（4,2,2）;

plot（bianhao,cd）;

subplot（4,2,3）;

plot（bianhao,cr）;

subplot（4,2,4）;

plot（bianhao,cu）;

subplot（4,2,5）;

plot（bianhao,hg）;

subplot（4,2,6）;

plot（bianhao,ni）;

subplot（4,2,7）;

plot（bianhao,pb）;

subplot（4,2,8）;

plot（bianhao,zn）;

data=[x;y;as]';

data1=[x;y;gongneng]';

x=data（:

1）;

y=data（:

2）;

as=data（:

3）;

gongz=data1（:

3）;

%scatter（x,y,z）;%É¢µãͼ

figure

[xx,yy,zz]=griddata（x,y,as,linspace（min（x）,max（x））',linspace（min（y）,max（y））,'v4'）;%²åÖµ

[xx,yy,zz]=griddata（x,y,gongz,linspace（min（x）,max（x））',linspace（min（y）,max（y））,'v4'）;%²åÖµ

daxiao=size（zz）;

zz=zz-3.6;

fori=1:

100

forj=1:

100

ifzz（i,j）<0

zz（i,j）=0;

end

end

end

xi=linspace（min（x）,max（x）,100）;

yi=linspace（min（y）,max（y）,100）;

[xi,yi]=meshgrid（xi,yi）;

zi=griddata（x,y,z,xi,yi）;%²åÖµ

gongnengi=griddata（x,y,gongneng,xi,yi）;

marker={'^','*','s','p','o'};

color={'r','b','k','c','y'};

mat={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'};

str={'µÈ¸ßÏß','Éú»îÇø','¹¤ÒµÇø','ɽÁÖÇø','½»Í¨Çø','Â̵ØÇø'};

holdon

figure

%contourf（xi,yi,zi,0:

10:

500,'b-'）;

%figure,surf（xx,yy,zz）