案例分析
一家医药制造公司生产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了,以上。