quadprog 的功能如下所示:solve.QP(Dmat, dvec, Amat, bvec, meq=0, factorized=FALSE)
我的目标函数是:
我的两个限制是:
现在我的目标是确定单个参数。
1)第一个论点 Dmat
:
在一般的符号是:
计算矩阵D是Dmat
:
因为D是对称的,所以解是
2)第二个参数 dvec
:目标函数的线性部分
3)第三个参数 Amat
:这是约束的左侧
4)第四个参数 bvec
:这是约束的右侧
5)第五个参数 meq
为零,因为只有不等式约束
6)第六个论点 factorized
:见问题
现在我必须将这六个参数放入函数中solve.QP()
。但在这里我确信需要进一步的“转型”......
- 当您查看包装说明时,二次部分乘以 0.5,线性部分为负
- 约束的定义也不同 - >= 符号
问题:
- 我必须如何“转换”前四个参数
solve.QP()
- 最后一个参数有什么
factorized
影响 - 如何最大化功能
solve.QP()
提前谢谢了!
编辑:我现在已经完成了我的“转换”,就像 G. Grothendieck 的链接中建议的那样——我希望我做的一切都正确?!
当我改变参数时factorized=TRUE
,我得到了 x1=2 和 x2=4 的预期解决方案。
在帮助中提到了 R^-1。但遗憾的是,我不知道这意味着什么......任何提示都表示赞赏。