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

《数学建模大赛A题金属污染Matlab源程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、应使用excel数据源，时间紧迫用txt代替clear,clccount=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'

2、)。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)。holdonAs=ones(size(as'))*as。[X,Y]=meshgrid(x,y)。x1=lins

3、pace(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)。3/3

4、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)m

5、inzn=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)。s

6、ubplot(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]'。3/3data1=[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,lin

7、space(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:100forj=1:100ifzz(i,j)<0zz(i,j)=0。endendendxi=linspace(min(x),max(x),100)。yi=linspace(min(y),max(y),

8、100)。[xi,yi]=meshgrid(xi,yi)。zi=griddata