2011年数学建模大赛A题金属污染Matlab源程序.doc
- 文档编号:15908194
- 上传时间:2023-07-08
- 格式:DOC
- 页数:3
- 大小:35.50KB
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)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2011 数学 建模 大赛 金属 污染 Matlab 源程序