- 博客(4)
- 收藏
- 关注
RSS订阅原创 gurobi不可解溯源——大M法中M取值问题
由于代码中多次用到了M，每处变量的量级不同，都更改工作量有点大，故其它地方都先暂时改小为1e5。再次求解，可解次数增加了！跑模型的时候用了分段线性化和用户均衡条件，都是用大M法线性化的，结果总是不可解。输出冲突约束发现都是分段线性化那部分的，但乍一看又很莫名其妙，看不出有哪几个约束冲突。总结：冲突约束物理意义不明显且涉及到M值时，考虑M设置不当的原因，将M改为相关变量的10-100倍，不能太大或太小。看到这次的冲突中出现了用户均衡部分的内容，故技重施。这里的变量最大值约为几十，故M用100，再次求解。
2023-12-15 22:18:08
388
原创 KKT条件下算节点边际电价为负解决方法
今天在学用KKT条件，算了电价导出，发现每个节点的电价都变-300了，给我弄迷糊了，怎么我用电还得电网给我贴钱（虽然上网查说确实会有这种情况？但总之在我这个模型里不太正常）。查了目标函数值，发现比较大，而且有几个变量取值都跑到边界上了，隐约感觉是边界条件限制得太紧了。free的边界意思是即使这些变量上下限都无穷，模型也不可解，我们可以排除它们。找到这几个>=的，发现这些都是用KKT时的对偶变量，那应该是说这些对偶变量的下限限制了模型的求解。因此我们去把这几个变量的下限改宽一点就好啦，解决！
2023-11-16 20:51:24
162
1
原创 pandas中concat函数下合并列重复
对第一大行第1个方阵列索引命名为['q1', ' ', ' ', ...]（方阵阶数为n，则希望第一个索引为q1，之后的n-1个都为空格，后同），对第一大行第2个方阵列索引命名为['q2', ' ', ' ', ...]；解决：对列索引进行更详细的命名，如第一大行第1个方阵列索引命名为['q1, 1', 'q1, 2 ', ' q1, 3', ...]，第一大行第2个方阵列索引命名为['q2, 1', 'q2, 2 ', ' q2, 3', ...]，这样就可以对应着合并了，没有数据的地方用NaN补充。
2023-11-13 20:46:17
354
原创 gurobi求解时间过长
结果发现这样会导致运行起来特别慢。即使是可以求解到最优解的情况，也会拖慢运行速度，并且得到的解精度很低。发现确实加上后预处理的内容减少了，导致后续计算内容增多（13→131，10→3316，40→3358）推测原因：查找了官方文档中关于这个参数的介绍，大概是这个参数的设置会影响预处理，导致计算量增大。解决：在确定模型求解情况以后及时删去该行代码，能不用就尽量不用orz。
2023-11-12 21:25:25
1352
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人