求教Matlab使用过Kriging插值函数的大虾

%克里金插值 Kriging interpolation

load data1

theta = [10 10]; lob = [1e-1 1e-1]; upb = [20 20];

[dmodel, perf] =dacefit(S, Y, @regpoly0, @corrgauss, theta, lob, upb)

X = gridsamp([0 0;100 100], 40);

[YX MSE] = predictor(X, dmodel);

X1 = reshape(X(:,1),40,40); X2 = reshape(X(:,2),40,40);

YX = reshape(YX, size(X1));

figure(1), mesh(X1, X2, YX)

hold on,

plot3(S(:,1),S(:,2),Y,'k', 'MarkerSize',10)

hold off

怎么用matlab实现kriging插值算法

第一步：% 已知某天海面2:00，8：00，14：00，20：00的

% 温度分别为 2210645， 2337419，2507742，2296129 ，（间隔为6小时）

% （1）采用线性内插、多项式内插和样条内插方法插值到一小时间隔，并比较不同；

% （2）对上述一小时内插结果进行调和分析；

% （3）对上述海表温度一小时内插结果进行Fourier分析。

clc; clear all; close all;

三维或两维半成图都牵扯到数据的网格化问题，而这却是matlab成图计算的最大缺陷。

对于数据的网格化，matlab通常对于有显性表达式的函数（z＝f（x，y））有很好的支持，但对于应用中的实际数据，其网格化的手段就捉襟见肘了（只有四种方法'linear'，'cubic'， 'nearest'， 'v4'），因为实际数据是没有函数表达式的。这使得matlab在处理实际问题时大受限制，很不方便；

数据的网格化，在数学领域也可以写一本厚厚的著作了，它牵扯到数据的插值，抽值，函数计算等许多方面，在地球科学领域广泛采用Golden Software公司的Surfer软件，其数据的网格化常用克里金方法，通常用surfer中处理完的网格化数据在matlab中成图；你的问题也是matlab初学者的一个典型问题，看似简单实际牵扯信号处理和数学方面的基础较多，数学是编程的基础。

形式是这样的：

[Vp, Ep]=krig3d(xp,yp,zp,xn,yn,zn,var,model_para,hmsg) performs 3-D kriging

% xn, yn, zn - coordinates of the input data

代表输入的变量 xn, yn, zn 是输入数据的空间坐标

xp, yp, zp 是输出kriging点的坐标

以上就是关于求教Matlab使用过Kriging插值函数的大虾全部的内容，包括:求教Matlab使用过Kriging插值函数的大虾、怎么用matlab实现kriging插值算法、关于matlab三维绘图的问题等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！