如何使用Cplex（docplex）输出单纯形乘法器？

Question

我正在将python与docplex库一起使用以解决线性编程问题。我想在最终表格中获得其单纯形乘法器或系数（松弛变量的系数）。

Answer 1

您看过对偶值吗？

dual_values(cts)

Returns the dual values of a sequence of linear constraints.

Note: the model must a pure LP: no integer or binary variable, no piecewise, no SOS. The model must also be solved successfully before
calling this method.

在http://ibmdecisionoptimization.github.io/docplex-doc/mp/docplex.mp.model.html中

让我将OPL示例volsay翻译成docplex

from docplex.mp.model import Model

mdl = Model(name='volsay')
Gas = mdl.continuous_var(name='Gas')
Chloride = mdl.continuous_var(name='Cloride')

mdl.add_constraint(Gas + Chloride <= 50, 'ctMaxTotal1')
mdl.add_constraint(3 * Gas + 4 * Chloride <= 180, 'ctMaxTotal2')
mdl.add_constraint(Chloride <= 40, 'ctMaxChloride')

mdl.maximize(40 * Gas + 50 * Chloride)

mdl.solve()



for v in mdl.iter_continuous_vars():
    print(v," = ",v.solution_value)

print("ctMaxTotalDual = ",mdl.dual_values(mdl.find_matching_linear_constraints('ctMaxTotal1'))); 

在OPL CPLEX中显示我们将要编写的对偶

// --------------------------------------------------------------------------
// Licensed Materials - Property of IBM
//
// 5725-A06 5725-A29 5724-Y48 5724-Y49 5724-Y54 5724-Y55
// Copyright IBM Corporation 1998, 2013. All Rights Reserved.
//
// Note to U.S. Government Users Restricted Rights:
// Use, duplication or disclosure restricted by GSA ADP Schedule
// Contract with IBM Corp.
// --------------------------------------------------------------------------

dvar float+ Gas;
dvar float+ Chloride;


maximize
  40 * Gas + 50 * Chloride;
subject to {
  ctMaxTotal:     
    Gas + Chloride <= 50;
  ctMaxTotal2:    
    3 * Gas + 4 * Chloride <= 180;
  ctMaxChloride:  
    Chloride <= 40;
}

tuple SolutionT{ 
    float Gas; 
    float Chloride; 
};
{SolutionT} Solution = {<Gas,Chloride>};
execute{ 
    writeln(Solution);
}

execute
{
  writeln("ctMaxTotalDual = ",ctMaxTotal.dual);
}

问候