利用【Python】【线性规划】优化工厂生产:实现智能资源配置与利润最大化的现代解决方案

目录

1. 问题背景和描述

1.1 问题背景

1.2 问题描述

2. 数学模型的建立

2.1决策变量

2.2 目标函数

2.3 约束条件

2.4 数学模型总结

3. 使用Python解决线性规划问题

3.1 导入必要的库

3.2 定义目标函数系数

3.3 定义不等式约束矩阵和向量

3.4 定义变量的边界

非负约束

变量边界在SciPy中的表示

3.5 求解线性规划问题

调用linprog函数

3.6 输出结果

获取和解释最优解

4. 结果解释和应用

4.1 结果解释

4.2 应用

4.3 验证结果

5. 扩展阅读

5.1 扩展阅读

5.2 线性规划的其他类型

​编辑

代码

结果

总结


 

ce6fbd68767d465bbe94b775b8b811db.png

731bd47804784fa2897220a90a387b28.gif

专栏:数学建模学习笔记

1. 问题背景和描述

1.1 问题背景

在现代制造业和商业运作中,资源的有效利用和利润的最大化是企业追求的重要目标。企业面临的常见问题是如何在有限的资源条件下,通过合理分配和优化使用资源,来实现利润的最大化。线性规划(Linear Programming,LP)是一种数学优化技术,能够在这些情况下发挥重要作用。它通过建立数学模型,帮助企业在众多可能的选择中找到最优解,进而指导实际操作。

假设有一家工厂,生产两种产品:产品A和产品B。每种产品的生产都需要消耗特定的资源。每天,每种资源的使用时间是有限的,这使得资源分配问题变得复杂。工厂的目标是通过合理分配资源,确定每天应该生产多少单位的产品A和产品B,以实现总利润的最大化。

具体数据如下:

  • 资源1:每天最多可用60个小时
  • 资源2:每天最多可用40个小时
  • 产品A:每单位需要资源1的2个小时和资源2的1个小时
  • 产品B:每单位需要资源1的1个小时和资源2的2个小时
  • 产品A:每单位的利润为30美元
  • 产品B:每单位的利润为20美元

通过对这些数据进行分析和建模,我们可以利用线性规划技术来制定一个优化的生产计划,确保在资源限制条件下实现利润的最大化。

1.2 问题描述

我们需要建立一个线性规划模型来描述上述问题,并使用该模型找到每天应该生产的产品A和产品B的最优数量,从而实现总利润的最大化。同时,生产过程中必须满足资源的限制条件,即不能超过每天可用的资源时间。

在这个问题中,我们的目标是构建一个数学模型,通过这个模型可以:

  1. 确定每天生产的产品A和产品B的数量。
  2. 满足资源的限制条件。
  3. 实现总利润的最大化。

2. 数学模型的建立

建立数学模型是解决线性规划问题的基础。我们将根据问题背景中的具体数据,定义决策变量,构建目标函数和约束条件。

2.1决策变量

决策变量是我们希望通过优化确定的数量。在这个问题中,决策变量是每天生产的产品A和产品B的数量。我们定义两个决策变量:

  • x1​:每天生产的产品A的单位数量
  • x2​:每天生产的产品B的单位数量

这些决策变量将用于构建目标函数和约束条件。

2.2 目标函数

目标函数是我们希望优化的表达式。在这个问题中,目标是最大化总利润。总利润可以表示为生产的产品A和产品B的利润之和。具体来说:

  • 产品A每单位的利润为30美元
  • 产品B每单位的利润为20美元

因此,总利润可以表示为: 利润=30x1​+20x2​

我们的目标是最大化总利润,因此目标函数可以表示为: 最大化 z=30x1​+20x2​

2.3 约束条件

约束条件是模型中必须满足的限制。在这个问题中,约束条件包括资源的限制和生产数量的非负性。

  1. 资源1的限制: 2*x1​+x2​≤60

  2. 资源2的限制:x1​+2*x2​≤40

  3. 非负约束:x1≥0 ,x2​≥0

这些约束条件确保生产计划不会超过可用的资源,并且生产的数量是非负的,即实际可行的。

2.4 数学模型总结

综合以上信息,我们可以建立一个完整的线性规划模型来描述这个问题。模型的形式如下:

3. 使用Python解决线性规划问题

在建立了数学模型之后,我们可以使用Python中的SciPy库来求解这个线性规划问题。SciPy库提供了许多优化算法,其中linprog函数可以用于求解线性规划问题。

3.1 导入必要的库

首先,我们需要导入必要的库:

import numpy as np
from scipy.optimize import linprog

3.2 定义目标函数系数

在求解问题时,SciPy的linprog函数默认是用于最小化问题的。因此,我们需要将最大化问题转换为最小化问题。具体来说,我们可以将目标函数的系数取负。

目标函数的系数为:c=[−30,−20]

在Python中定义目标函数系数:

c = [-30, -20]

3.3 定义不等式约束矩阵和向量

在Python中,我们可以定义约束矩阵和向量如下:

A = [[2, 1], [1, 2]]
b = [60, 40]

这里,A 是一个二维数组,表示约束条件的系数矩阵,每一行对应一个不等式约束条件,每一列对应一个决策变量。b 是一个一维数组,表示每个约束条件的右端常数项。

  • 矩阵 A:每一行代表一个约束条件,每一列代表一个决策变量。在我们的例子中,第一行 [2, 1] 表示第一个约束条件 2x1​+x2​≤60,第二行 [1, 2] 表示第二个约束条件x1​+2x2​≤40。
  • 向量 b:每个元素表示一个约束条件的右端常数项。对于我们的例子,向量 b 中的元素分别是 60 和 40,对应两个约束条件的右端值。

3.4 定义变量的边界

在实际的生产问题中,变量的取值范围通常是有限制的。对于我们的例子,每天生产的产品数量不能为负,因此我们需要设置变量的边界条件。

非负约束

线性规划问题中的非负约束是指决策变量必须是非负数,即:

x1​≥0       x2​≥0

这些非负约束条件确保了我们的生产数量是合理的(即,不能生产负数的产品)。

变量边界在SciPy中的表示

在SciPy的linprog函数中,变量的边界可以通过bounds参数来指定。每个决策变量的边界条件可以用一个元组表示,元组的第一个元素是变量的下界,第二个元素是变量的上界。如果变量没有上界,可以用None表示。

x_bounds = [(0, None), (0, None)]

这里,(0, None) 表示变量的下界是0,上界没有限制,即变量必须是非负的。

  • 变量边界定义:我们用一个列表来表示每个变量的边界。列表中的每个元素是一个元组,元组的第一个元素表示变量的下界,第二个元素表示变量的上界。在我们的例子中,我们定义了两个变量 x1​ 和 x2​,它们的边界条件都是非负的,因此我们用(0, None)来表示它们的边界。
  • 意义和应用:非负约束条件确保了我们的生产计划是现实可行的,因为生产的数量不能为负数。这在实际应用中是非常重要的,可以防止在优化过程中出现不合理的解。

3.5 求解线性规划问题

在定义了目标函数的系数、约束条件和变量的边界之后,我们可以使用SciPy的linprog函数来求解这个线性规划问题。linprog函数是SciPy库中用于求解线性规划问题的主要函数。

调用linprog函数

linprog函数的基本调用方式如下:

res = linprog(c, A_ub=A, b_ub=b, bounds=x_bounds, method='highs')

其中:

  • c:目标函数的系数向量
  • A_ub:不等式约束条件的系数矩阵
  • b_ub:不等式约束条件的右端常数项向量
  • bounds:变量的边界条件
  • method:求解方法,这里我们使用highs方法,这是SciPy推荐的高效求解方法之一
  • 目标函数的系数向量 c:之前定义的目标函数的系数向量[-30, -20]。注意这里我们将最大化问题转换为最小化问题,因此系数取负。
  • 不等式约束条件的系数矩阵 A_ub:定义了每个约束条件的系数矩阵[[2, 1], [1, 2]]
  • 不等式约束条件的右端常数项向量 b_ub:对应于不等式约束的右端常数项向量[60, 40]
  • 变量的边界条件 bounds:定义了变量的取值范围,确保变量是非负的。
  • 求解方法 method:指定使用highs方法,这是SciPy中一个高效的线性规划求解器。

3.6 输出结果

在求解完成之后,我们需要输出结果。结果对象res包含了优化的详细信息,包括最优解和最优目标函数值。

获取和解释最优解

最优解是指在满足所有约束条件的情况下,使目标函数达到最优值的决策变量值。在我们的例子中,最优解是每天应该生产的产品A和产品B的数量。

print('Optimal value:', round(res.fun * -1, ndigits=2))  # 恢复最大化问题的目标函数值
print('x:', res.x)

 这里,res.fun 是最优目标函数值,因为我们最初将目标函数系数取了负,所以需要乘以-1来恢复最大化问题的目标函数值。res.x 是最优解,即最优的 x1​ 和 x2​ 的值。

4. 结果解释和应用

在得到线性规划问题的求解结果之后,我们需要对结果进行解释和应用。最优解和最优目标函数值对于实际的生产计划具有重要的指导意义。

4.1 结果解释

  1. 最优目标函数值:1400。这个值表示在满足所有资源约束的情况下,最大化的总利润为1400美元。
  2. 最优解:每天生产20个单位的产品A和10个单位的产品B。这是指在所有约束条件下能够使总利润最大化的最优生产计划。

4.2 应用

最优解对于工厂的生产计划具有重要的指导意义。通过按照最优解安排生产,工厂可以确保资源的有效利用,最大化利润。

具体来说,工厂应该每天生产20个单位的产品A和10个单位的产品B。这将使得总利润最大化,同时不超过每天可用的资源限制。通过这种方式,工厂可以实现资源的最佳配置,提高生产效率和经济效益。

4.3 验证结果

为了验证结果的正确性,我们可以检查最优解是否满足所有约束条件:

  1. 资源1的使用情况: 2×20+1×10=40+10=50≤60   

  2. 资源2的使用情况: 1×20+2×10=20+20=40≤40   

可以看到,最优解不仅最大化了总利润,而且满足所有约束条件。这表明我们的线性规划模型和求解过程是正确的,求解结果是合理的。

5. 扩展阅读

5.1 扩展阅读

线性规划问题在实际中有很多应用,如物流运输、生产计划、资源分配等。除了SciPy,其他常用的优化库还包括PuLP和Gurobi。

  • PuLP:PuLP是一个开源的线性规划工具,可以与各种求解器结合使用。它提供了简单易用的接口,适合用于教学和简单的优化问题。
  • Gurobi:Gurobi是一款高效的商业优化软件,支持求解线性规划、整数规划和其他优化问题。它具有强大的求解能力,适用于大规模和复杂的优化问题。

5.2 线性规划的其他类型

  1. 多目标优化:多目标优化考虑多个目标函数同时进行优化。例如,一个工厂可能希望在最大化利润的同时最小化污染物排放。多目标优化可以通过加权求和法或Pareto最优解来解决。

  2. 带有等式约束的线性规划:在约束条件中包含等式约束。例如,某些资源的使用量必须精确等于特定值。等式约束可以用于表示这些严格的资源限制或平衡条件。

  3. 混合整数线性规划:决策变量不仅包括连续变量,还包括整数变量。例如,在生产计划中,某些产品的生产数量必须为整数。混合整数线性规划可以通过将部分决策变量定义为整数来解决这些问题。

代码

import numpy as np
from scipy.optimize import linprog

# 定义目标函数系数
# 我们希望最大化 30x1 + 20x2
# 在使用 linprog 时,我们需要将这个目标函数转化为最小化问题
# 因此我们取负,变为最小化 -30x1 - 20x2
c = [-30, -20]

# 定义不等式约束矩阵和向量
# 不等式约束如下:
# 2x1 + x2 <= 60
# x1 + 2x2 <= 40
# 转化为矩阵形式 A @ x <= b
A = [[2, 1], [1, 2]]
b = [60, 40]

# 定义变量的边界
# x1 >= 0
# x2 >= 0
x_bounds = [(0, None), (0, None)]

# 使用 linprog 函数求解线性规划问题
res = linprog(c, A_ub=A, b_ub=b, bounds=x_bounds, method='highs')

# 输出结果
print('Optimal value:', round(res.fun * -1, ndigits=2))  # 恢复最大化问题的目标函数值
print('x:', res.x)

# 结果解释
# 最优目标函数值
optimal_value = round(res.fun * -1, ndigits=2)
# 最优解
optimal_solution = res.x

print(f"在满足资源约束的情况下,最大化的总利润为:{optimal_value} 美元")
print(f"每天生产 {optimal_solution[0]} 个单位的产品A 和 {optimal_solution[1]} 个单位的产品B")

# 验证结果是否满足所有约束条件
# 资源1的使用情况
resource1_usage = 2 * optimal_solution[0] + optimal_solution[1]
# 资源2的使用情况
resource2_usage = optimal_solution[0] + 2 * optimal_solution[1]

print(f"资源1的使用情况:{resource1_usage} 小时(<= 60 小时)")
print(f"资源2的使用情况:{resource2_usage} 小时(<= 40 小时)")

# 验证是否满足所有约束条件
if resource1_usage <= 60 and resource2_usage <= 40:
    print("最优解满足所有约束条件。")
else:
    print("最优解不满足所有约束条件。")

结果

总结

     建立线性规划模型,以解决生产优化问题。通过定义决策变量、目标函数和约束条件,使用Python的SciPy库中的linprog函数求解模型,并验证结果的合理性。最终,确定了在资源限制条件下最大化利润的最优生产方案.

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/791084.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

【Pytorch】RNN for Image Classification

文章目录 1 RNN 的定义2 RNN 输入 input, h_03 RNN 输出 output, h_n4 多层5 小试牛刀 学习参考来自 pytorch中nn.RNN()总结RNN for Image Classification(RNN图片分类–MNIST数据集)pytorch使用-nn.RNN 1 RNN 的定义 nn.RNN(input_size, hidden_size, num_layers1, nonlinea…

游戏视频是后期配音好还是边录边配 游戏视频怎么剪辑制作才能火 视频剪辑免费软件

游戏视频后期配音是先配还是先剪&#xff1f;游戏视频后期配音没有统一的准则&#xff0c;可以先配&#xff0c;也可以后配&#xff0c;主要是根据内容而定。游戏视频剪辑在游戏玩家中十分流行&#xff0c;那么&#xff0c;游戏视频怎么剪辑制作&#xff1f;下面让我们以具体的…

apache python使用

修改httpd.conf文件。 AddHandler cgi-script .cgi .py 代码 #!自己的python.exe #-*- coding:UTF-8 -*- print ("Content-type:text/html") print () print (<html>) print (<head>) print (<meta charset"gb2312">) print (<tit…

现在国内的ddos攻击趋势怎么样?想了解现在ddos的情况该去哪看?

目前&#xff0c;国内的DDoS攻击趋势显示出以下几个特征&#xff1a; 攻击频次显著增加&#xff1a;根据《快快网络2024年DDoS攻击趋势白皮书》&#xff0c;2023年DDoS攻击活动有显著攀升&#xff0c;总攻击次数达到1246.61万次&#xff0c;比前一年增长了18.1%。 攻击强度和规…

Spring Expression表达式使用

Spring Expression 简介 Spring Expression Language&#xff08;简称 “SpEL”&#xff09;是一种功能强大的表达式语言&#xff0c;支持在运行时查询和操作对象图。其中最显著的是&#xff1a;方法调用和基本的字符串模板功能。 虽然SpEL是Spring产品组合中表达式评估的基础…

树莓派pico入坑笔记,ssd1306使用

目录 说明 后附进阶玩法&#xff1a;显示中文&#xff0c;外加简单库实现 官方模块使用 使用样例 方法说明 下面是绘图支持的方法 进阶玩法&#xff0c;显示中文 方法&#xff0c;对汉字取字模&#xff0c;然后按像素对字模进行显示 说明 circuitpython8.x使用iic接口驱…

链表 OJ(一)

移除链表元素 题目连接&#xff1a; https://leetcode.cn/problems/remove-linked-list-elements/description/ 使用双指针法&#xff0c;开始时&#xff0c;一个指针指向头节点&#xff0c;另一个指针指向头节点的下一个结点&#xff0c;然后开始遍历链表删除结点。 这里要注…

YOLOv10改进 | Conv篇 | CVPR2024最新DynamicConv替换下采样(解决低FLOPs陷阱)

一、本文介绍 本文给大家带来的改进机制是CVPR2024的最新改进机制DynamicConv其是CVPR2024的最新改进机制&#xff0c;这个论文中介绍了一个名为ParameterNet的新型设计原则&#xff0c;它旨在在大规模视觉预训练模型中增加参数数量&#xff0c;同时尽量不增加浮点运算&#x…

YOLOv10改进 | Conv篇 | 全新的SOATA轻量化下采样操作ADown(参数量下降百分之二十,附手撕结构图)

一、本文介绍 本文给大家带来的改进机制是利用2024/02/21号最新发布的YOLOv9其中提出的ADown模块来改进我们的Conv模块&#xff0c;其中YOLOv9针对于这个模块并没有介绍&#xff0c;只是在其项目文件中用到了&#xff0c;我将其整理出来用于我们的YOLOv10的项目&#xff0c;经…

Python 视频的色彩转换

这篇教学会介绍使用OpenCV 的cvtcolor() 方法&#xff0c;将视频的色彩模型从RGB 转换为灰阶、HLS、HSV...等。 因为程式中的OpenCV 会需要使用镜头或GPU&#xff0c;所以请使用本机环境( 参考&#xff1a;使用Python 虚拟环境) 或使用Anaconda Jupyter 进行实作( 参考&#x…

【TAROT学习日记】韦特体系塔罗牌学习(1)——愚者 THE FOOL 0

韦特体系塔罗牌学习&#xff08;1&#xff09;——愚者 THE FOOL 0 https://www.tarotchina.net/major-arcana0-vip/ 目录 韦特体系塔罗牌学习&#xff08;1&#xff09;——愚者 THE FOOL 0牌面分析1. 基础信息2. 图片元素 正位牌意1. 关键词/句2.爱情婚姻3. 学业事业4. 人际关…

android13 rom frameworks 蓝牙自动接收文件

总纲 android13 rom 开发总纲说明 目录 1.前言 2.源码查找 3.我们先实现第一种改法 4.实现第二种改法 5.第三种改法代码参考 6.编译测试 1.前言 我们从导航栏这里,点开这个蓝牙的接收框,弹出来的对话框,使用android studio 的layout inspector可以发现这个是 Bluetoo…

有必要找第三方软件测评公司吗?如何选择靠谱软件测评机构?

软件测试是确保软件质量的重要环节&#xff0c;而在进行软件测试时&#xff0c;是否有必要找第三方软件测评公司呢?第三方软件测评公司是指独立于软件开发公司和用户之间的中立机构&#xff0c;专门从事软件测试和测评工作。与自身开发团队或内部测试团队相比&#xff0c;选择…

大白话讲解AI大模型

大白话讲解大模型 大模型的发展重要大模型发展时间线 大模型的简单原理-训练⼤模型是如何训练并应⽤到场景中的&#xff1f;如果训练私有化模型 模型&#xff1a;model 语料库&#xff1a;用于训练模型的数据 大模型的发展 详细信息来源&#xff1a;DataLearner 2022年11月底…

JVM相关知识点汇总

JDK,JRE以及JVM的关系 我们的编译器到底干了什么事? 仅仅是将我们的 .java 文件转换成了 .class 文件,实际上就是文件格式的转换,对等信息转换。 类加载机制是什么? > **所谓类加载机制就是** > ``` > 虚拟机把Class文件加载到内存 > 并对数据进行校验,转换…

web安全及内网安全知识

本文来源无问社区&#xff08;wwlib.cn&#xff09;更多详细内容可前往观看http://www.wwlib.cn/index.php/artread/artid/7506.html Web安全 1、sql注入 Web程序中对于用户提交的参数未做过滤直接拼接到SQL语句中执行&#xff0c;导致参数中的特殊字符破坏了SQL语句原有逻…

qt 用数据画一个图,并表示出来

1.概要 想用数据绘制一个画面&#xff0c;看有相机到播放的本质是啥。 要点 // 创建一个QImage对象&#xff0c;指定图像的宽度、高度和格式 QImage image(width, height, QImage::Format_Grayscale8); // 将像素数据复制到QImage对象中 memcpy(image.bits(), pixelD…

【Linux网络】IP协议{初识/报头/分片/网段划分/子网掩码/私网公网IP/认识网络世界/路由表}

文章目录 1.入门了解2.认识报头3.认识网段4.路由跳转相关指令路由 该文诸多理解参考文章&#xff1a;好文&#xff01; 1.入门了解 用户需求&#xff1a;将我的数据可靠的跨网络从A主机送到B主机 传输层TCP&#xff1a;由各种方法&#xff08;流量控制/超时重传/滑动窗口/拥塞…

PTC可复位保险丝 vs 传统型保险丝:全面对比分析

PTC可复位保险丝&#xff0c;又称为自恢复保险丝、自恢复熔断器或PPTC保险丝&#xff0c;是一种电子保护器件。它利用材料的正温度系数效应&#xff0c;即电阻值随温度升高而显著增加的特性&#xff0c;来实现电路保护。 当电路正常工作时&#xff0c;PTC保险丝呈现低阻态&…

最新浪子授权系统网站源码 全开源免授权版本

最新浪子授权系统网站源码 全开源免授权版本 此版本没有任何授权我已经去除授权&#xff0c;随意二开无任何加密。 更新日志 1.修复不能下载 2.修复不能更新 3.修复不能删除用户 4.修复不能删除授权 5.增加代理后台管理 6.重写授权读取文件 7.修复已经知道漏洞 源码下…