备战国赛--散乱插值及MATLAB绘图函数

来源：互联网发布：好的算法书籍推荐编辑：程序博客网时间：2024/06/05 17:46

常见的散乱插值方法有Kriging（克里金）插值，Shepard插值等

这里使用了matlab中的工具直接进行插值并绘图，示例代码如下：

clc,clearA=xlsread('C:\Users\Administrator\Desktop\集训\8.23\3. 2011~A资料\2011A~题目\data.xls','附件1','A4:E322');C=xlsread('C:\Users\Administrator\Desktop\集训\8.23\3. 2011~A资料\2011A~题目\data.xls','附件2','A4:I322');x=A(:,2);y=A(:,3);z=A(:,4);qy=A(:,5);As=C(:,2);Cd=C(:,3);Cr=C(:,4);Cu=C(:,5);Hg=C(:,6);Ni=C(:,7);Pb=C(:,8);Zn=C(:,9);[X,Y,Z]=griddata(x,y,z,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,Z)title('采样地形图');[X,Y,AS]=griddata(x,y,As,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,AS)title('As浓度随平面地形分布图');[X,Y,CD]=griddata(x,y,Cd,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,CD)title('Cd浓度随平面地形分布图');[X,Y,CR]=griddata(x,y,Cr,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,CR)title('Cr浓度随平面地形分布图');[X,Y,CU]=griddata(x,y,Cu,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,CU)title('Cu浓度随平面地形分布图');[X,Y,NI]=griddata(x,y,Ni,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,NI)title('Ni浓度随平面地形分布图');[X,Y,HG]=griddata(x,y,Hg,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,HG)title('Hg浓度随平面地形分布图');[X,Y,PB]=griddata(x,y,Pb,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,PB)title('Pb浓度随平面地形分布图');[X,Y,ZN]=griddata(x,y,Zn,linspace(0,28654)',linspace(0,18449),'v4');figure,mesh(X,Y,ZN)title('Zn浓度随平面地形分布图');[X,Y,QY]=griddata(x,y,qy,linspace(0,28654)',linspace(0,18449),'v4');figure,surf(X,Y,QY)title('区域分布图');

目标为已知坐标、海拔、地区类型以及各重金属元素含量的319个调查点。

各二维三维绘图函数在《精通MATLABR2014b》中有简略介绍，也可以自行上网搜索。

阅读全文

0 0