写一个R代码，用球状模型拟合半变异函数，参数包括步长，变程，基台，块金，块金效应

好的，这里是一个使用球状模型拟合半变异函数的 R 代码：

library(gstat)

# 定义半变异函数
h <- function(h, a, c0, c) {
  ifelse(h <= a, c0 + c * (3 * h / (2 * a) - 0.5 * (h / a)^3),
         c0 + c)
}

# 设定参数
nugget <- 0
sill <- 1
range <- 10
n <- 100
model <- vgm(psill = sill, model = "Sph", range = range, nugget = nugget)

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

# 计算距离矩阵
d <- spDists(cbind(x, y))

# 计算半变异函数值
v <- h(d, range, nugget, sill - nugget)

# 拟合球状模型
fit <- fit.variogram(vgm(psill = sill, model = "Sph", range = range), v, d)

# 绘制拟合结果
plot(fit, d, v, main = "Spherical Model Fitting")

其中，`nugget` 表示块金效应，`sill` 表示基台加块金，`range` 表示变程。`h` 函数定义了球状模型的半变异函数形式，`fit.variogram` 函数用于拟合半变异函数。最后，通过 `plot` 函数可视化拟合结果。