docplex解决一维数组和二维数组的共同作用的优化问题

最新推荐文章于 2026-06-26 20:20:32 发布

原创最新推荐文章于 2026-06-26 20:20:32 发布 · 594 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#python #docplex #cplex #最优化问题

docplex解决最优化问题--实战专栏收录该内容

7 篇文章

订阅专栏

该文展示了如何使用Python的Docplex库来构建线性规划模型。通过一维数组`x`和二维数组`y`作为决策变量，定义了它们之间的相互制约关系和约束条件，并设定优化目标，最后利用Cplex求解器找到最优解。

Python3.8

Python 是一种高级、解释型、通用的编程语言，以其简洁易读的语法而闻名，适用于广泛的应用，包括Web开发、数据分析、人工智能和自动化脚本

使用背景

使用 continuous_var_list() 定义决策变量 x
使用 continuous_var_matrix() 定义决策变量 y
x 和 y 之间的相互制约、相互限制的约束条件添加。

使用方法

#导入包
import cplex
import docplex
from docplex.mp.model import Model
from typing import List

# 输入数据
qB=[20,15,13,10,4]
cB=[2,5,7,10,11]
qQ=[2,10,15,10]
cQ=[20,18,15,13]

# 定义优化模型 
mdl = Model()

## 定义决策变量  
### continuous_var_list 定义一维决策变量 
index1: List[int]=[i for i in range(len(qQ))]
x = mdl.continuous_var_list(index1, name='X')

### continuous_var_matrix 定义一维决策变量 
index2=[i for i in range(len(qB))]
y = mdl.continuous_var_matrix(index2, index1, name='Y')

# 定义约束条件 
cst = [] 
cst += [x[j] >= 0 for j in index1] + [x[j] <= qQ[j] for j in index1] + [mdl.sum(x[j] for j in index1) <= sum(qB)] 
cst += [y[i,j] >= 0 for i in index2 for j in index1] + [mdl.sum(y[i,j] for i in index2) == x[j] for j in index1] + [mdl.sum(y[i,j] for j in index1) <= qB[i] for i in index2]
for c in cst: 
    mdl.add_constraint(c) 

# 定义优化函数
obj = mdl.sum(y[i,j]*(pow((cQ[j]-cB[i]),3)) for i in index2 for j in index1)

# 最大化 目标函数 ，并求解 
mdl.maximize(obj)

#求解模型并显示
solution = mdl.solve()
print(solution)