栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 软件开发 > 后端开发 > Python

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

Python 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

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

转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/423677.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号