关于 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矩阵 那是个对称阵，所以行列谁优先都无所谓：
这是最后的答案：