Python 线性规划 求解

案例分析

一家医药制造公司生产6种产品，其在资源方面的需求、单位利润和需求量的情况如表1-1所示。对每种产品的需求也是有限的，每月生产的产品不能超过需求的数量。已知公司可提供的劳动力工时数为5000h，可提供的原料为2000lb（1lb=0.453592kg）。公司想知道各生产每种产品多少磅才可以使得利润最大化？

消耗系数 产品1 产品2 产品3 产品4 产品5 产品6
劳动力/h 1 2 3 4 5 6
原料/lb 3.4 2.8 1.9 0.5 0.6 0.1
单位利润/元 1.1 3.2 4.2 5.6 5.2 6
需求量/lb 560 460 760 300 1003 1043

现有 5000 2000

代码

#导入包
from scipy import optimize
import numpy as np

#确定c,A,b,Aeq,beq
c = np.array([1.1,3.2,4.2,5.6,6.5,2.6])
A = np.array([[1,2,3,4,5,6],[3.4,2.8,9,0.5,0.6,0.1],[1,0,0,0,0,0],[0,1,0,0,0,0],[0,0,1,0,0,0],[0,0,0,0,1,0],[0,0,0,0,0,1]])
b = np.array([5000,2000,560,460,760,1003,1043])
# Aeq = np.array([[1,1,-2,1]])
# beq = np.array([10])

#求解
res = optimize.linprog(-c,A,b)
print(res)

结果

重点关注的就是fun和最后一行了，fun是整个结果，最后一行是每个x的结果。为什么第一行是负的呢？原来这个函数其实是求最小值的，那么求最大值，怎么办呢？很简单，仔细观察的人应该发现，之前的函数里面，我写的是-c，而不是c。那么这个函数的出来的结果其实就是-c的最小值，但很明显这恰恰是c最大值的相反数。那么答案就是7183.9999了，以上。