等浓度Matlab程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
A1=wenti31(:,1:3);%wenti31:附件一x,y坐标值及附件二As浓度值
x=A1(:,1);
y=A1(:,2);
z=A1(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'As浓度等值线图'})
hold on
A2=wenti32(:,1:3);%wenti32:附件一x,y坐标值及附件二Cd浓度值
x=A2(:,1);
y=A2(:,2);
z=A2(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Cd浓度等值线图'})
hold on
A3=wenti33(:,1:3);%wenti33:附件一x,y坐标值及附件二Cr浓度值
x=A3(:,1);
y=A3(:,2);
z=A3(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Cr浓度等值线图'})
hold on
A4=wenti34(:,1:3);%wenti34:附件一x,y坐标值及附件二Cu浓度值
x=A4(:,1);
y=A4(:,2);
z=A4(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Cu浓度等值线图'})
hold on
A5=wenti35(:,1:3);%wenti35:附件一x,y坐标值及附件二Hg浓度值
x=A5(:,1);
y=A5(:,2);
z=A5(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Hg浓度等值线图'})
hold on
A6=wenti36(:,1:3);%wenti36:附件一x,y坐标值及附件二Ni浓度值
x=A6(:,1);
y=A6(:,2);
z=A6(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Ni浓度等值线图'})
hold on
A7=wenti37(:,1:3);%wenti37:附件一x,y坐标值及附件二Pb浓度值
x=A7(:,1);
y=A7(:,2);
z=A7(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Pb浓度等值线图'})
hold on
A8=wenti38(:,1:3);%wenti38:附件一x,y坐标值及附件二Zn浓度值
x=A8(:,1);
y=A8(:,2);
z=A8(:,3);
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)',linspace(min(y),max(y),200),'v4'); %插值
figure,contourf(X,Y,Z) %等值线图
[C,h]=contourf(X,Y,Z);
clabel(C,h)
title({'Zn浓度等值线图'})。