0

quadprog 的功能如下所示:solve.QP(Dmat, dvec, Amat, bvec, meq=0, factorized=FALSE)

我的目标函数是:

我的两个限制是:

现在我的目标是确定单个参数。

1)第一个论点 Dmat

在一般的符号是:

计算矩阵DDmat

因为D是对称的,所以解是

2)第二个参数 dvec:目标函数的线性部分

3)第三个参数 Amat:这是约束的左侧

4)第四个参数 bvec:这是约束的右侧

5)第五个参数 meq为零,因为只有不等式约束

6)第六个论点 factorized:见问题

现在我必须将这六个参数放入函数中solve.QP()。但在这里我确信需要进一步的“转型”......

  • 当您查看包装说明时,二次部分乘以 0.5,线性部分为负

  • 约束的定义也不同 - >= 符号

问题

  1. 我必须如何“转换”前四个参数solve.QP()
  2. 最后一个参数有什么factorized影响
  3. 如何最大化功能solve.QP()

提前谢谢了!

编辑:我现在已经完成了我的“转换”,就像 G. Grothendieck 的链接中建议的那样——我希望我做的一切都正确?!

在此处输入图像描述

当我改变参数时factorized=TRUE,我得到了 x1=2 和 x2=4 的预期解决方案。

在此处输入图像描述

在帮助中提到了 R^-1。但遗憾的是,我不知道这意味着什么......任何提示都表示赞赏。

在此处输入图像描述

4

0 回答 0