关于 python—cvxopt.solvers.qp 这个玩意儿算不对的问题
参考了这两篇博客:
https://blog.csdn.net/qq_45669448/article/details/104678910
https://blog.csdn.net/QW_sunny/article/details/79793889
计算这么一个破玩意儿
不就是个小小的不等式约束么?? 瞧不起我???
然后他就算不对(
from cvxopt import solvers, matrix
P = matrix([[1., 0., 0.],
[0., 1., 0.],
[0., 0., 0.]])
q = matrix([0.,
0.,
0.])
G = matrix([[-3., -3., -1.],
[-4., -3., -1.],
[ 1., 1., 1.]])
h = matrix([-1.,
-1.,
-1])
solvers.options['show_progress'] = True
sol = solvers.qp(P,q,G=G,h=h)
print(sol['x'])
print(sol['primal objective'])
结果是:
显然不对,DEBUG半天,就离谱…
上边儿那两篇博客,不都是那么做的嘛?? 为啥按这个就是不行
最后在我的好友校长的其启发下,发现这个cvxopt.matrix
是个列优先的矩阵…
所以,这里转置一下就OK了
G = matrix([[-3., -3., -1.],
[-4., -3., -1.],
[ 1., 1., 1.]]).T # 看见没,这里加个转置
那你要问我,为啥人家上边的博客没啥问题
你去看看他们G和P矩阵 那是个对称阵,所以行列谁优先都无所谓:
这是最后的答案: