R拟合半变异函数模型，用球状模型拟合

好的，以下是用球状模型拟合半变异函数的 R 代码：

# 生成随机数据
set.seed(123)
n <- 100
x <- runif(n)
y <- runif(n)
z <- rnorm(n)

# 计算距离矩阵
dist_mat <- as.matrix(dist(cbind(x, y)))

# 计算半变异函数
gamma <- variogram(z ~ x + y, data = data.frame(x, y, z), width = max(dist_mat) / 10, cutoff = max(dist_mat) / 2)

# 拟合球状半变异函数模型
spherical_model <- fit.variogram(gamma, vgm("Sph", psill = max(gamma$gamma), range = max(dist_mat) / 2))

# 绘制拟合结果
plot(gamma, pch = 16, col = "black", main = "Spherical model")
lines(spherical_model, col = "red", lwd = 2)
legend("topleft", legend = c("Experimental", "Spherical"), col = c("black", "red"), lty = 1, lwd = 2)

上述代码中，首先生成了随机数据，并计算了距离矩阵`dist_mat`。然后，用`variogram`函数计算半变异函数`gamma`。接着，用`vgm`函数定义了球状半变异函数模型，并用`fit.variogram`函数拟合半变异函数曲线得到`spherical_model`。最后，绘制了实验结果和拟合结果的图像，并加入了图例。