用R语言编写一下OBJ=((n-1) C_P+C_m (∑_i^n▒∫_0^(T_i)▒〖α^(i-1) (t)dt〗)+C_f)/(∑_i^n▒T_i )
时间: 2023-12-07 19:05:23 浏览: 21
可以使用以下代码来计算OBJ:
```
n <- 10 # 替换为实际的n值
C_P <- 100 # 替换为实际的C_P值
C_m <- 50 # 替换为实际的C_m值
C_f <- 200 # 替换为实际的C_f值
lambda <- function(t) {
# 替换为实际的λ(t)函数
return(0.1)
}
T <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10) # 替换为实际的T值
OBJ <- ((n - 1) * C_P + C_m * sum(sapply(1:n, function(i) integrate(function(t) (0.1^(i-1)) * lambda(t), 0, T[i])$value)) + C_f) / sum(T)
OBJ # 输出OBJ的值
```
请注意,上述代码假定lambda(t)是一个固定的函数,可以通过替换函数体来进行修改。