optimization - 动态优化

Question

我正在尝试设置一个问题，使随着时间（10 年）的消费（u）最大化，这取决于财富的变化（x _t+1 - x _t = -c _t），其中 x ₀ =1000。我试图解决它，但无法在初始参数处评估函数。此外，最终条件是：x ₁₀ >= 0 和 u ₁₀ >= 0。这是我正在使用的代码：

library('optimx')

# Parameters
t <- seq(0,10) # time
desc <- 0.1 # tax
beta <- 1/(1+desc) # parameter of function
alpha <- 0.8 # function
xo = 1000 # wealth in t=0

# Variables
x <- rep(xo, length(t)) # wealth in time

# Function
utilidade_total <- function(u, npar=TRUE, print=TRUE){
  for (i in lenght(t)){
    obj[i]=(beta^i)*(u[i]^alpha)
    x[i]=x[i-1]-u[i-1]
  }
  result <- sum(obj)
  return(result)
}

maxuti <- optimx(par=u,fn= utilidade_total,
                 control = list(maximize=TRUE))

umax <- coef(maxuti)

for (i in 2:lenght(t)){
   print(obj[i]=(beta^i)*(u^alpha))
    print(stock[i]=x[i-1]-umax[i-1])
}

# Result
for (i in 2:lenght(t)){
  obj[i]=(beta^i)*(u^alpha)
  x[i]=x[i-1]-u[i-1]
}

plot(t, x, type='b')
plot(t, u, type='b')