2维高斯插值的MATLAB代码实现

以下是2维高斯插值的MATLAB代码实现：

function [zq] = gauss_interp(x,y,z,xq,yq)
% 2D Gaussian interpolation
% x, y, z: input data
% xq, yq: query points
% zq: interpolated values at query points

% calculate distances between query points and input data
dx = bsxfun(@minus,x,xq(:)');
dy = bsxfun(@minus,y,yq(:)');
dist2 = dx.^2 + dy.^2;

% calculate weights using Gaussian kernel
sigma = 1;
w = exp(-dist2/(2*sigma^2));

% normalize weights
w = bsxfun(@rdivide,w,sum(w));

% interpolate
zq = w*z(:);
zq = reshape(zq,size(xq));