1

在将 NLOPTR 的 ISRES 算法应用于具有不等式约束的非线性问题时，我陷入了困境。我是这样制定的：

library(nloptr)
fn <- function(x) {
(x[1]-10)^2 + 5*(x[2]-12)^2 + x[3]^4 + 3*(x[4]-11)^2 + 10*x[5]^6 + 7*x[6]^2 + x[7]^4 - 4*x[6]*x[7] - 10*x[6] - 8*x[7]
}

hin <- function(x) {
h <- numeric(4)
h[1] <- 127 - 2*x[1]^2 - 3*x[2]^4 - x[3] - 4*x[4]^2 - 5*x[5]
h[2] <- 282 - 7*x[1] - 3*x[2] - 10*x[3]^2 - x[4] + x[5]
h[3] <- 196 - 23*x[1] - x[2]^2 - 6*x[6]^2 + 8*x[7]
h[4] <- -4*x[1]^2 - x[2]^2 + 3*x[1]*x[2] -2*x[3]^2 - 5*x[6] +11*x[7]
return(h)
}

x0 <- c(1, 2, 0, 4, 0, 1, 1)

isres(x0 = x0, fn = fn, hin = hin)

我收到一条消息，上面写着“匹配错误（hin）：缺少参数“表”，没有默认值”

我想我在不等式约束下做的不对。你能告诉我如何解决这个问题吗？非常感谢！

4

1 回答 1

3

这不是一个真正的答案。但我目前无法发表评论。

我认为问题不在于您的限制，而在于功能isres。如果您查看函数的代码（只需isres在控制台中输入），开始部分如下所示：

function (x0, fn, lower, upper, hin = NULL, heq = NULL, maxeval = 10000, 
          pop.size = 20 * (length(x0) + 1), xtol_rel = 1e-06, nl.info = FALSE, 
          ...) 
{
  opts <- list()
  opts$maxeval <- maxeval
  opts$xtol_rel <- xtol_rel
  opts$population <- pop.size
  opts$algorithm <- "NLOPT_GN_ISRES"
  fun <- match.fun(fn)
  fn <- function(x) fun(x, ...)
  if (!is.null(hin)) {
    .hin <- match(hin)
    hin <- function(x) (-1) * .hin(x)
  } 
...

现在，请注意这一行.hin <- match(hin)。match 函数是基础函数R，应该接受两个参数，即：x和在table执行时x %in% table（查看文档以获取更多信息）。但是，不提供表格给功能。 ?matchRisres

就目前而言，您有 3 个选项来解决您的问题：

  • (I)nloptr直接与算法一起使用"NLOPT_GN_ISRES"。
  • (II) 包的管理员告诉我，这个 bug 可能已经知道了。因此，如果您使用github上的开发包而不是 CRAN 版本，您可能会很好。注意这个包没有在我的电脑上编译
  • （三）等到bug修复。这应该很快，因为他们现在知道了这个问题。
于 2016-07-19T09:02:34.080 回答