补充1 设计区域网络配置的优化模型(工厂选址模型)——基于java和cplex的实现
这一部分是第二部分的供应链网络的设计(二)供应链的网络设计第四部分设计区域网络配置的优化模型。
目录
一、所需输入
二、有能力约束的工厂选址模型
1.利用excel进行规划求解
2.利用java和cplex进行规划求解
三、考虑税收、关税和顾客要求
一、所需输入
使用网络优化模型来进行区域网络配置时,管理者需要获得各个区域的以下信息:需求、期望的响应水平、兴建设施的固定成本、劳动力和材料的变动成本、库存持有成本、每两个地区之间的运输成本、产品售价、税收和关税、设施产能。结合这些成本及税收和关税信息,就可以了解在每个地区兴建设施的年固定成本,以及从一个地区到另一个地区提供单位产品的可变成本。
MoonLight公司是一家经营不锈钢制造的公司,其产品销往全球各个地方。接下来,就借MoonLight公司为例说明区域网络的设计。假如你是该公司负责供应链的副总裁,当前有两种备选方案,你会选择哪一种方案呢?
第一种方案是在每个地区都兴建一个工厂。这种方案的好处是可以降低运输成本,而且有助于避免从其他地区进口产品可能被征收的关税。这种方案的缺点是工厂的规模是按满足当地需求来规划的,有可能无法充分利用规模经济。
第二种方案将工厂集中建在少数几个地区,这样可以提高规模经济水平,但会增加运输成本和关税。
作为总裁的你,必须同时你还要考虑对这些可量化因素进行权衡,同时还必须考虑竞争环境、政治风险等一些不可量化的因素。面对这样的方案选择,当然不能拍脑袋就决定,必须进行一番细致的考察和计算,接下来就正式开始介绍工厂选址模型。
第一步是收集可用于量化模型的数据。对于MoonLight公司来说,你决定要全球需求划分为5个区域:北美、南美、欧洲、非洲和亚洲。所收集的数据如下图所示。上述几个区域的年需求都在B9 :F9中。单元格B4:F8中的数据标识在一个区域组织生产来满足另一个区域需求的可变生产、库存和运输成本(包括关税和税收)。所有成本都是以千美元为单位。这个阶段收集的数据是一个总量数据。
MoonLight公司的成本(以千美元计)和需求数据(以百万单位计)
在每个设施,与设施、运输和库存相关的成本既有固定成本也有可变成本。设施、运输和库存成本通常存在规模经济效应,边际成本随着设施生产数量的增加而降低。在这里讨论时,我们假定所有的可变成本与生产量或运输量程线性正相关。
MoonLight公司对于每个选址的工厂规模也有两个方案。低产能的工厂每年生产1000万单位产品,高产能工厂每年可生产2000万单位产品,分别如单元格H4:H8和J4:J8所示。由于高产能工厂具有一定的规模经济,所以其固定成本不到低产能工厂固定成本的两倍,如单元格G4:G8和I4:I8所示。所有的固定成本都按年计算。你想知道哪种网络方案的成本最低。为了回答这个问题接下来将要讨论这种情况的、有能力约束的工厂选址模型。
二、有能力约束的工厂选址模型
有能力约束的工厂选址网络优化模型需要以下输入:
在MoonLight公司的例子中,n=10,因为有5个潜在工厂位置和2种可供选择的工厂产能。m=5,因为有5个市场。年需求Dj见B9:F9单元格。工厂潜在产能Ki见单元格H4:H8和J4:J8。单元格G4:G8和I4:I8中数据标识年固定成本fi。变动成本cij如表格B4:F8所示。
供应链小组的目标时确定一个可以使税后利润最大化的网络设计。为了简化起见,假设所有需求都必须被满足且不考虑所得税。因此,模型仅关注使满足全球需求的成本最小化。不过,也可以对模型稍加修改,将利润和税收纳入考虑之中。相关决策变量定义如下:
则问题可以表述为以下混合整数规划问题:
约束条件为:
目标函数是使简历和运营该网络的总成本(固定成本+可变成本)最小化。式(1)的约束条件是指每个区域的市场需求都必须被满足。式(2)的约束条件是指每个工厂的供应不能超出其产能。式(3)的约束条件是指一个工厂或是开工,或是关闭。函数的姐将表明哪些工厂将开工及其产能,以及各区域的需求在这些工厂之间如何分配。
1.利用excel进行规划求解
该模型可利用Excel中的规划求解工具来进行求解。由于篇幅原因,这部分内容不过介绍,有兴趣的朋友可以自己上网查询。
上图便是EXCEL得到的结果,分别在南美、亚洲和非洲建设一个高产能的工厂。南美的工厂满足北美的需求。欧洲的需求由亚洲和非洲的工厂来满足。这种方案下,最低成本为23751千美元。
2.利用java和cplex进行规划求解
Cplex是IBM公司开发的一款商业版的优化引擎,该引擎专门用于求解大规模的线性规划(LP)、二次规划(QP)、带约束的二次规划(QCQP)、二阶锥规划(SOCP)等四类基本问题,以及相应的混合整数规划(MIP)问题,并且与众多优化软件及语言兼容(与C++,JAVA,EXCEL,Matlab等都有接口)。
java引入cplex在这不多讲述了,不了解的直接在论坛里搜就可以了,废话不多说,直接上代码。
public class InputData {//产品需求double[] demands={12,8,14,16,7};//低产能固定成本double[] lowFixedCost={6000,4500,6500,4100,4000};//高产能固定成本double[] highFixedCost={9000,6750,9750,6150,6000};//低产能double[] lowProductivity = {10,10,10,10,10};//高产能double[] highProductivity = {20,20,20,20,20};//成本double[][] cost = {{81,92,101,130,115},{117,77,108,98,100},{102,105,95,119,111},{115,125,90,59,74},{142,100,103,105,71}};int nodeNumber = demands.length;
}
import ilog.concert.IloException;
import ilog.concert.IloNumExpr;
import ilog.concert.IloNumVar;
import ilog.concert.IloNumVarType;
import ilog.cplex.IloCplex;public class LocationDemo {InputData data;public LocationDemo(InputData data) {this.data = data;}//定义cplex内部对象IloCplex model;//定义变量public IloNumVar[][] y;public IloNumVar[][] x;//求解函数public void solve() throws IloException{if(model.solve()==false){System.out.println("模型不可解");} else {System.out.println(model);System.out.println("目标值:"+model.getObjValue());System.out.println("=========选择变量===========");for(int i=0;i<data.nodeNumber;i++){for(int j=0;j<2;j++) {if (model.getValue(y[i][j]) != 0) {System.out.println("变量值y[" + (i + 1) + "]["+(j+1)+"]:" + model.getValue(y[i][j]));}}}System.out.println("=========分配关系===========");for(int i=0;i<data.nodeNumber;i++){for(int j=0;j<data.nodeNumber;j++) {if (model.getValue(x[i][j]) != 0) {System.out.println("变量值x[" + (i + 1) + "]["+(j+1)+"]:" + model.getValue(x[i][j]));}}}}}//根据数学模型建立求解模型public void BuildModel() throws IloException{//初始化modelmodel = new IloCplex();model.setOut(null);//定义决策变量x和y的数据类型及取值范围y = new IloNumVar[data.nodeNumber][2];x = new IloNumVar[data.nodeNumber][data.nodeNumber];for (int i = 0; i < data.nodeNumber; i++) {for(int j=0;j<data.nodeNumber;j++) {x[i][j] = model.numVar(0,20, IloNumVarType.Int, "x[" + i + "]["+j+"]");}}for (int i = 0; i < data.nodeNumber; i++) {for(int j=0;j<2;j++) {y[i][j] = model.numVar(0,1, IloNumVarType.Bool, "y[" + i + "]["+j+"]");}}//设置目标函数IloNumExpr obj = model.numExpr();for(int i=0;i<data.nodeNumber;i++){for (int j = 0; j < data.nodeNumber; j++) {obj = model.sum(obj,model.prod(data.cost[i][j],x[i][j]));}}for (int i = 0; i < data.nodeNumber; i++) {obj = model.sum(obj,model.prod(data.lowFixedCost[i],y[i][0]));}for (int i = 0; i < data.nodeNumber; i++) {obj = model.sum(obj,model.prod(data.highFixedCost[i],y[i][1]));}model.addMinimize(obj);//设置约束条件1for(int j=0;j<data.nodeNumber;j++) {IloNumExpr expr1 = model.numExpr();for (int i = 0; i < data.nodeNumber; i++) {expr1 = model.sum(expr1, x[i][j]);}model.addEq(expr1,data.demands[j]);}//设置约束条件2for (int i = 0; i < data.nodeNumber; i++) {IloNumExpr expr2 = model.numExpr();for(int j=0;j<data.nodeNumber;j++){expr2=model.sum(expr2,x[i][j]);}model.addLe(expr2,model.sum(model.prod(y[i][0],data.lowProductivity[i]),model.prod(y[i][1],data.highProductivity[i])));}}
}
public class test {public static void main(String[] args) throws Exception {InputData data =new InputData();LocationDemo lp=new LocationDemo(data);lp.BuildModel();lp.solve();}
}
结果和excel运行结果相同,下图为运行结果:
三、考虑税收、关税和顾客要求
所有网络设计模型的构建应确保最终的供应链网络在满足顾客服务要求的同时,实现税后利润最大化。以下所介绍的方法可用于将任何成本最小化的网络设计模型转换为利润最大化模型。如果rj是在市场j卖出一单位产品的收入,则有能力约束的工厂选址模型的目标函数可以修改为:
该目标函数是使企业的利润最大化。当使用一个利润最大化的目标函数时,管理者应当将上面的(1)式修改为:
上面的约束条件更为适合,因为它允许网络设计者识别处哪些需求可以获利,而满足哪些需求会使企业造成损失,用上面的式子为利润最大化的工厂选址模型将进服务于能够获利的那部分需求。除非有其他约束条件,否则这可能导致一些市场中的部分需求被放弃。因为服务这些需求企业无利可图。下面是改进后的java代码。
public class ImproveData {//产品需求double[] demands={12,8,14,16,7};//低产能固定成本double[] lowFixedCost={6000,4500,6500,4100,4000};//高产能固定成本double[] highFixedCost={9000,6750,9750,6150,6000};//低产能double[] lowProductivity = {10,10,10,10,10};//高产能double[] highProductivity = {20,20,20,20,20};//成本double[][] cost = {{81,92,101,130,115},{117,77,108,98,100},{102,105,95,119,111},{115,125,90,59,74},{142,100,103,105,71}};//单位产品的收入double[] income = {550,525,555,530,520};int nodeNumber = demands.length;
}
import ilog.concert.IloException;
import ilog.concert.IloNumExpr;
import ilog.concert.IloNumVar;
import ilog.concert.IloNumVarType;
import ilog.cplex.IloCplex;public class ImproveDemo {ImproveData data;public ImproveDemo(ImproveData data) {this.data = data;}//定义cplex内部对象IloCplex model;//定义变量public IloNumVar[][] y;public IloNumVar[][] x;//求解函数public void solve() throws IloException{if(model.solve()==false){System.out.println("模型不可解");} else {System.out.println(model);System.out.println("目标值:"+model.getObjValue());System.out.println("=========选择变量===========");for(int i=0;i<data.nodeNumber;i++){for(int j=0;j<2;j++) {if (model.getValue(y[i][j]) != 0) {System.out.println("变量值y[" + (i + 1) + "]["+(j+1)+"]:" + model.getValue(y[i][j]));}}}System.out.println("=========分配关系===========");for(int i=0;i<data.nodeNumber;i++){for(int j=0;j<data.nodeNumber;j++) {if (model.getValue(x[i][j]) != 0) {System.out.println("变量值x[" + (i + 1) + "]["+(j+1)+"]:" + model.getValue(x[i][j]));}}}}}//根据数学模型建立求解模型public void BuildModel() throws IloException{//初始化modelmodel = new IloCplex();model.setOut(null);//定义决策变量x和y的数据类型及取值范围y = new IloNumVar[data.nodeNumber][2];x = new IloNumVar[data.nodeNumber][data.nodeNumber];for (int i = 0; i < data.nodeNumber; i++) {for(int j=0;j<data.nodeNumber;j++) {x[i][j] = model.numVar(0,20, IloNumVarType.Int, "x[" + i + "]["+j+"]");}}for (int i = 0; i < data.nodeNumber; i++) {for(int j=0;j<2;j++) {y[i][j] = model.numVar(0,1, IloNumVarType.Bool, "y[" + i + "]["+j+"]");}}//修改后的目标函数IloNumExpr obj = model.numExpr();for (int i = 0; i < data.nodeNumber; i++) {for (int j = 0; j < data.nodeNumber; j++) {obj = model.sum(obj,model.prod(data.income[i],x[j][i]));}}for(int i=0;i<data.nodeNumber;i++){for (int j = 0; j < data.nodeNumber; j++) {obj = model.diff(obj,model.prod(data.cost[i][j],x[i][j]));}}for (int i = 0; i < data.nodeNumber; i++) {obj = model.diff(obj,model.prod(data.lowFixedCost[i],y[i][0]));}for (int i = 0; i < data.nodeNumber; i++) {obj = model.diff(obj,model.prod(data.highFixedCost[i],y[i][1]));}model.addMaximize(obj);//修改后的设置约束条件1for(int j=0;j<data.nodeNumber;j++) {IloNumExpr expr1 = model.numExpr();for (int i = 0; i < data.nodeNumber; i++) {expr1 = model.sum(expr1, x[i][j]);}model.addLe(expr1,data.demands[j]);}//设置约束条件2for (int i = 0; i < data.nodeNumber; i++) {IloNumExpr expr2 = model.numExpr();for(int j=0;j<data.nodeNumber;j++){expr2=model.sum(expr2,x[i][j]);}model.addLe(expr2,model.sum(model.prod(y[i][0],data.lowProductivity[i]),model.prod(y[i][1],data.highProductivity[i])));}}
}
public class test {public static void main(String[] args) throws Exception {ImproveData data =new ImproveData();ImproveDemo lp=new ImproveDemo(data);lp.BuildModel();lp.solve();}
}
运行结果:
顾客的偏好和要求可以通过期望的响应时间、对运输方式或运输承运人的选择等方式表现出来。例如,假设在工厂选址i和市场j之间存在两种运输方式:方式1是海运,方式2是空运。那么,可以对工厂选址模型做以下修正:定义两个独立的决策变量x1ij和x2ij,分别标识采用方式1和方式2从选址i运送到市场j的数量。只有当使用某种运输方式所花费的时间小于所期望的响应时间时,才允许利用这种运输方式进行装运。例如,如果采用方式1(海运)从地址i运送产品到市场j所需时间大于顾客所能接受的时间,那么可以在工厂选址模型中将决策变量x1ij赋值为0。存在多个运输承运人是可以对模型做类似修改。改进的代码就不放在这里了。
补充1 设计区域网络配置的优化模型(工厂选址模型)——基于java和cplex的实现相关推荐
- 补充3 需求分配和工厂选址模型(工厂和仓库同时选址)——基于java和cplex的实现
这一部分是第二部分的供应链网络的设计(二)供应链的网络设计第六部分在区域内识别潜在地点的模型. 目录 一.所需输入 二.将需求分配给现有生产设施的需求分配模型 三.生产设施选址模型 四.更复杂的.有能 ...
- java网络爬虫论文_毕业设计(论文)-基于JAVA的网络爬虫的设计与实现.doc
您所在位置:网站首页 > 海量文档  > 计算机 > Java 毕业设计(论文)-基于JAVA的网络爬虫的设计与实现. ...
- 玩转“网上邻居”之网络配置(二)
接上篇 二.混合系统网络 目前实际应用中单一系统的对等网络已不是主流,主要是由于这种网络性能较低,缺乏很好看安全及网络管理能力.目前在各企事业单位中普遍应用的还是基于多系统的域结构混合网络.但因目前主 ...
- 基于JAVA读书网络社区设计计算机毕业设计源码+系统+lw文档+部署
基于JAVA读书网络社区设计计算机毕业设计源码+系统+lw文档+部署 基于JAVA读书网络社区设计计算机毕业设计源码+系统+lw文档+部署 本源码技术栈: 项目架构:B/S架构 开发语言:Java语言 ...
- Solaris11之网络配置(更新:2021-10-21)
目录 Solaris11网络配置新特性 手动配置网络 创建静态IP地址 创建IP接口 创建IP地址 配置默认网关 配置DNS服务器 配置主机名 放开root远程登录 Solaris11网络配置新特性 ...
- windows网络配置
前言 Windows 操作系统是两大主流之一的网络操作系统,因此相关的基础知识和对应的网络管理.配置指令在历年的考试中经常出现.接下来我来简单的介绍windows的网络配置. 文章目录 前言 域模型与 ...
- 基于Java的飞机大战游戏的设计与实现论文
源码下载 http://www.byamd.xyz/hui-zong-1/ 摘 要 现如今,随着智能手机的兴起与普及,加上4G(the 4th Generation mobile communicat ...
- 基于Java的飞机大战游戏的设计与实现(含源文件)
欢迎添加微信互相交流学习哦! 项目源码:https://gitee.com/oklongmm/biye 基于Java的飞机大战游戏的设计与实现 摘 要 现如今,随着智能手机的兴起与普及,加上4G( ...
- 基于Java的飞机大战游戏的设计与实现(论文+源码)
毕业设计(论文)任务书 第1页 毕业设计(论文)题目:基于Java的飞机大战游戏的设计与实现毕业设计(论文)要求及原始数据(资料):1.综述基于Java的飞机大战游戏的现状:2.深入了解Java的设计 ...
最新文章
- php 格式化评论量函数,深入剖析PHP中printf()函数格式化使用
- 字符串 mysql索引,长字符串的MySQL索引
- 数据挖掘——我们能从股市数据得出什么,以及一些算法
- b g opencv读入的图片 r_OpenCV提取显示一张图片(或者视频)的R,G,B颜色分量
- 使用单例模式加载properties文件
- foursquare nyc数据集_炫酷的python地理数据可视化
- NIO客户端主要创建过程
- 电脑如何测网速_物联网卡的网速到底怎么样呢
- 未给员工足额缴纳公积金!董明珠曾豪言:每人一套房不需要公积金
- 使用FlexBox和Json实现类似ComboBox(类似Google的输入提示和自动)功能-基于JQuery-ASP.NET...
- NodeJs和ReactJs单元测试工具——Jest
- WDS+ADK+MDT部署Windows pxe
- LTP学习笔记(基本安装+执行)
- 《从零开始的 RPG 游戏制作教程》第十二期:刷怪与防守(上)
- 计算机应用基础7次作业答案,北京中医药大学远程教育“计算机应用基础”第7次作业(14页)-原创力文档...
- 查看文章影响因子的插件_查询文献可实时显示影响因子与分区排名的2个强大浏览器插件...
- 算法图解 第7章 狄克斯特拉算法
- 【物联网服务NodeJs-5天学习】第三天实战篇③ ——基于MQTT的环境温度检测
- JPA操作数据库中的脏数据更新作用
- 跨界融合引尴尬,智能珠宝能向前走远?
热门文章
- ubuntu编译安装libtool
- 真无线蓝牙耳机哪个延迟最低?2023延迟超低的蓝牙耳机推荐
- vue中table表格导出为图片格式
- 谈谈面试题之什么是面向对象?谈谈你对面向对象的理解?
- php 7编译 phar,解决PHP 7编译安装错误:cannot stat ‘phar.phar’: No such file or directory...
- 根据日期查询当前万年历信息
- 查询android手机端包名,查看android手机原生app包名的方法
- 微信小游戏入门实战《小猫钓鱼》
- Romberg求积分算法
- python3 生成条形码_python3转换code128条形码的方法