我有
形式的二次优化问题min w.r.t.到x 1/2 x' x + q' x S.t. Gx< = h
我有一个相当大的问题(几百万点和约束),虽然cvxopt的默认求解器被证明是有效的,但我很好奇在SCS中实现它应该更快(并且不使用CVXPY)接口)。
在进行文献检索(主要是Boyd&#39> 凸面分析)之后,对SOCP形式的重新制定应该产生
min w.r.t.在Z
中受到Az + s = b,s的s的zc =(1 q)'和z =(t x)'其中t是标量,K是与我的原始约束(Gx \ leq h)和二次锥Q = {(t,x)|相关的线性锥的笛卡尔积。 t> = || x ||}
但是,我该如何定义A和b?
我想象A = [[0 G],[2,1,..,1]]和b =(h 0)之类的东西。 然而,随着''在python scs.solve cone字典参数中将选项设置为[1]我无法使其工作? A的最后一行的预期语法是什么? (也就是说,假设我的数学重新设计是正确的......)
感谢您的帮助!