用r语言写一下三次指数平滑法的代码

下面是用R语言实现三次指数平滑法的代码：

# 三次指数平滑法函数
# x 为原始数据，alpha为平滑系数，beta和gamma为趋势系数，n为预测期数
triple_exponential_smoothing <- function(x, alpha, beta, gamma, n){
  # 初始化
  m <- length(x)
  y <- rep(0, m+n)
  s <- rep(0, m+n)
  b <- rep(0, m+n)
  y[1] <- x[1]
  s[1] <- x[1]
  b[1] <- x[2] - x[1]
  
  # 三次指数平滑法公式
  for (i in 2:(m+n)){
    if (i <= m){
      y[i] <- x[i]
    } else {
      y[i] <- alpha * (x[i-m] - b[i-m]) + (1 - alpha) * (y[i-1] + b[i-1])
    }
    s[i] <- gamma * (y[i] - b[i-1]) + (1 - gamma) * s[i-1]
    b[i] <- beta * (s[i] - s[i-1]) + (1 - beta) * b[i-1]
  }
  
  # 返回预测值
  return(y[(m+1):(m+n)])
}

# 示例
x <- c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100)
n <- 5
alpha <- 0.2
beta <- 0.3
gamma <- 0.4
y_pred <- triple_exponential_smoothing(x, alpha, beta, gamma, n)
print(y_pred)

其中，x为原始数据，alpha为平滑系数，beta和gamma为趋势系数，n为预测期数。函数输出预测期的预测值。