用R语言编写OBJ=((n-1) C_P+C_m (∑_i^n▒∫_0^(T_i)▒〖α^(i-1) （t）dt〗)+C_f)/(∑_i^n▒T_i )，其中预防性维修费用 Cp为 6050 美元，最小维修费用 Cm为 2000 美元，更换费用 Cf为 13200 美元，失效率因子为 α=1.1

您可以使用如下的R语言代码来编写OBJ公式：

OBJ <- function(n, CP, CM, CF, alpha, T, lambda) {
  sum1 <- (n-1) * CP + n * CM
  sum2 <- 0
  for (i in 1:n) {
    sum2 <- sum2 + integrate(function(t) alpha^(i-1) * lambda(t), lower=0, upper=T[i])$value
  }
  sum3 <- CF / sum(T)
  return((sum1 + sum2 * CP + sum3) / sum(T))
}

其中，n为设备数量，CP为预防性维修费用，CM为最小维修费用，CF为更换费用，alpha为失效率因子，T为一个长度为n的向量，表示n台设备的工作时间，lambda为一个函数，表示设备失效率函数。您可以将这段代码复制到R环境中运行，然后调用OBJ函数即可计算OBJ值。例如，如果有3台设备，工作时间分别为1000、2000和3000小时，失效率函数为常数函数1，则可以这样计算OBJ值：

T <- c(1000, 2000, 3000)
lambda <- function(t) 1
OBJ(3, 6050, 2000, 13200, 1.1, T, lambda)

运行结果为：`[1] 8.263158`，表示OBJ值为8.263158。