matlab 二维IDW代码

以下是MATLAB中进行二维IDW插值的示例代码：

% 样本点
x = [1 2 3 4 5];
y = [1 2 3 4 5];
z = [2 5 6 5 1];

% 插值点
xi = linspace(1, 5, 20);
yi = linspace(1, 5, 20);

[XI, YI] = meshgrid(xi, yi);

% 计算距离并求权重
r = hypot(bsxfun(@minus, x', XI(:)), bsxfun(@minus, y', YI(:)));
w = 1 ./ r;

% 归一化权重
w = bsxfun(@rdivide, w, sum(w, 1));

% 计算插值结果
ZI = reshape(w * z', length(yi), length(xi));

% 绘制插值结果
contourf(XI, YI, ZI, 20, 'LineStyle', 'none');
colorbar;

在这个示例中，我们先定义了一组样本点(x, y, z)，然后定义了一组插值点(xi, yi)，并通过meshgrid函数生成了插值网格。接下来，我们计算了每个插值点与样本点之间的距离，并求出了距离的倒数作为权重。为了避免插值结果受到距离远近的影响，我们对权重进行了归一化处理。最后，我们将权重和样本点的值相乘并进行累加，得到了插值结果(ZI)，并通过contourf函数绘制了插值结果的等高线图。