《关于python cvxopt.solvers.qp 算不对这件事》 cvxopt.matrix是个列优先的矩阵

关于 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矩阵 那是个对称阵,所以行列谁优先都无所谓:在这里插入图片描述
这是最后的答案:
在这里插入图片描述

评论 7
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值