数学建模入门笔记(1) Python pulp库解线性规划问题

参考:Python求解线性规划——PuLP使用教程 - Only(AR) - 博客园 (cnblogs.com)

1.Define the model

model = pl.LpProblem(name="",sense=pl.LpMaximize)

name 模型的名字

sense 模型的类型(pl.LpMaximize/pl.LpMinimize)

2.Define the decision variables

用x[i]存储变量,命名为xi

x = {i: pl.LpVariable(name=f"x{i}", lowBound=,upBound=, cat=) for i in range(1,)}

name 变量名

lowBound 变量最小值(默认负无穷)

upBound 变量最大值 (默认正无穷)

cat 变量的类型(pl.Binary 0/1,pl.Integer 整数,pl.Continuous 实数(默认值))

3.Set the objective

目标函数不用括号

model += x[1]+x[2]+x[3]

4.Add constraints

约束式用括号括起来

式子可以是>=,<=,==,不能是>,<,=

model += (2*x[1]+3*x[2]+x[3]>=10)

5.Solve the optimization problem

status = model.solve()

6.Print the results

print(f"status: {pl.LpStatus[model.status]}")
print(f"objective: {model.objective.value()}")
 
for var in model.variables():
    print(f"{var.name}: {var.value()}")
 
for name, constraint in model.constraints.items():
    print(f"{name}: {constraint.value()}")
pl.LpStatus[model.status]解释
Optimal找到了最优解
Not Solved问题尚未被求解
Infeasible问题是不可行的,即不存在满足所有约束条件的解
Unbounded问题是无界的,即目标函数可以无限增加(或无限减少)而不受限制
Undefined问题是未定义的,通常表示输入数据错误或问题规格不正确

*完整代码

《数学建模算法与应用》司守奎、孙兆亮 例1.2

import pulp as pl
model = pl.LpProblem(name="problem",sense=pl.LpMaximize)
x = {i:pl.LpVariable(name=f"x{i}",lowBound=0,cat=pl.LpContinuous) for i in range(1,4)}
model += (2 * x[1] - 5 * x[2] + x[3] >= 10)
model += (x[1] + 3 * x[2] + x[3] <= 12)
model += (x[1] + x[2] + x[3] == 7)
model += 2 * x[1] + 3 * x[2] -5 * x[3]
status = model.solve()
print(f"status:{pl.LpStatus[model.status]}")
print(f"objective:{model.objective.value()}")
for var in model.variables():
    print(f"{var.name}:{var.value()}")
for name,constraint in model.constraints.items():
    print(f"{name}:{constraint.value()}")
  • 9
    点赞
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值