MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]

MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]线性规划是数学规划中的一类最简单规划问题,常见的线性规划是一个有约束的,变量范围为有理数的线性规划。如:对于这类线性规划问题,数学理论已经较为完善,可以有多种方法求解此类问题。但写这篇文章的目的并不是为了介绍数学理论,我们这里主要讲解如果利用工具求解这一类线性规划问题。最著名,同时也是最强大的数学最优化软件是LINGO/LINDO软件包,它能够求解多种的数学规划问题,同时还提供了多

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

 

线性规划是数学规划中的一类最简单规划问题,常见的线性规划是一个有约束的,变量范围为有理数的线性规划。如:

MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]

对于这类线性规划问题,数学理论已经较为完善,可以有多种方法求解此类问题。但写这篇文章的目的并不是为了介绍数学理论,我们这里主要讲解如果利用工具求解这一类线性规划问题。

最著名,同时也是最强大的数学最优化软件是LINGO/LINDO软件包,它能够求解多种的数学规划问题,同时还提供了多种的分析能力。但LINGO软件并不容易上手,同时,应用LINGO的场合一般是大规模的线性规划问题,小小的线性规划完全可以不使用它。一个更受科研人员欢迎的数学软件是MATLAB,它以功能强大而称著,并有数学软件中的“航空母舰”之称。我们这里就是要学习使用MATLAB软件求解线性规划(含整数规划和0-1规划)问题。

为了使得不熟悉MATLAB的人员也能够使用MATLAB进行线性规划问题求解,本文将对MATALB中使用到的函数和过程以及结果进行详细的分析,最后会对每一个问题都给出一个可以完全“套用”的MATLAB程序。

我们首先从上面的线性规划问题开始,为了便于表达,将上面的式子写成矩阵形式:

MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]

于是约束就表达为了一个不等式。

求解MATLAB线性规划时,最常用的函数是linprog函数,下面来介绍一下这个函数的使用。

打开MATLAB帮助文档(PS:帮助文档的内容是最全的,只要你的英文过了专业8级),可以看到linprog函数求解的是具有如下标准形式的线性规划:

 MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]

公式中各符号的意义是自明的,在这里简单介绍下,首先MATLAB中求解的是目标函数是最小值的问题,但如果我们的目标函数是求最大值,可以通过对目标函数中每一项中乘以-1,将求最大值问题转化为求最小值问题;A,b分别为不等式约束中的系数矩阵。Aeq和beq分别为等式约束中的系数矩阵,lb,和ub分别为每个变量的上下区间;最后f为目标函数中各变量的系数矩阵。

现在,是时候动动手,使用MATLAB编写代码求解这个线性规划了。MATLAB代码如下所示:

f=[-7,-12];A=[9 4;4 5;3 10];b=[300;200;300];lb=zeros(2,1);% 生成一个2行1列的全0矩阵,很显示,上面例子中的x,y的最小值为0[x,fval]=linprog(f,A,b,[],[],lb,[])

我们来解释下linprog函数中每参数的意义,linprog中的一个原型如下:

[x,fval,exitflag] = linprog(f,A,b,Aeq,beq,lb,ub)

这7个参数的意义和上面f、A、b的意义是一样的。f为目标函数的系数矩阵,A为线性规划不等式约束的变量系数矩阵,b为不等式约束的资源数(如上面的[300;200;300]),这是一个N行1列的矩阵,N为变量的个数。Aeqbeq是相应等式约束的变量系数矩阵和资源数(很明显,上面的例子中并没有等式约束)。lbub分别为保变量的上下区间。在上面的例子中,x和y和最小值都为0但都无最大值约束。而linprog的返回值x为求得的各变量的值,这是一个向量,fval为最优化的值,一般是一个标量,exitflag意为函数的退出标志。

上面所示的代码[x,fval]=linprog(f,A,b,[],[],lb,[])中,[]代表不存在或空,因为在上面的例子中不存在等式约束,所以Aeqbeq的位置为[]。而ub也为空,是因为变量没有最大值约束。

运行上面的程序,行到结果为:

x =

 

      20.0000

      24.0000

 

 

fval =

 

      -428.0000

 

解释为:

当x=20,y=24时,可以求得最优化的值,最大值为428(因为这里的求目标最大值,但MATLAB只能求目标函数最小值,所以对目标函数进行了乘-1处理,所以也要对最后的结果乘以-1才是目标函数所求).

 

上面解决了简单的线性规划问题的求解,线性规范有两种比较特殊的情况,即整数规划和0-1整数规划。在之前(不知MATLAB几之前……),MATLAB是不能直接求解这两种规划的,bintprog函数可以用来求0-1整数规划,但求解过程比较麻烦,而且最新版的MATLAB已经遗弃了这个函数,同时提供了一个比较新的、专用于求解整数规划和0-1整数规划的函数——intlinprog。intlinprog的一个原型为:

[x,fval,exitflag]= intlinprog(f,intcon,A,b,Aeq,beq,lb,ub)

该函数的使用和linprog函数的使用十分相似,其仅仅在linprog函数的基础上多了一个参数——intcon。我们来通过下面的例子来学习该参数的意义。

 MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]

在这里例子中,变量的取值范围不再是有理数集,而是整数集。求解此规划问题的MATLAB程序如下:

f_13=[-1,-1];ic_13=[1,2];A_13=[-4,2;4,2;0,-2];b_13=[-1;11;-1];lb_13=zeros(2,1);[x_13,fval_13,flag_13]=intlinprog(f_13,ic_13,A_13,b_13,[],[],lb_13,[])

在函数intlinprog中,intcon的意义为整数约束变量的位置。如上例中,因为x1和x2都要是整数,intcon参数位置ic_13的值为[1,2]。这个位置是按照目标函数和约束条件中变量位置来排列的。如果上式中仅有x2为整数约束,那么ic_13的值应该为2。

需要说明的是,intlinprog函数在比较旧版本是不支持的(笔者使用的是MATLAB2014B),如果你发现你现在的MATLAB没有intlinprog函数,请不要吃惊,因为一直以来,MATLAB都是无法直接求解整数规划的,但今时已经不同往日了。

现在又有了一个新问题,我们解决了在MATLAB上求解一般的整数规划问题,但要是遇到0-1整数规划问题呢?到这里,我们只要转换一下思维,就可以利用MATLAB求解0-1整数规划了,这里先卖个关子,请大家看下面的例子是怎么用MATLAB求解0-1整数规划的。

MATLAB求解线性规划(含整数规划和0-1规划)问题[通俗易懂]

MATLAB程序如下:

f_12=[7 5 9 6 3];ic_12=[1,2,3,4,5];A_12=[56,20,54,42,15;1,4,1,0,0;-1,-2,0,-1,-2];b_12=[100;4;-2];lb_12=zeros(5,1);ub_12=ones(5,1);[x_12,fval_12,flag_12]=intlinprog(f_12,ic_12,A_12,b_12,[],[],lb_12,ub_12)

有木有发现在,与上面整数规划不同的地方只有一个,就是多了ub_12=ones(5,1),也就是说求解0-1整数规划只要在求解整数规划的基础上加上一个对变量最大值约束为1就行了,有木有恍然大悟的感觉???

后面两个程序并没有给出程序运行的结果,因为笔者坚信学习最好的方式就是“动手”。^_^

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/157845.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 2021 Pychram激活码_通用破解码

    2021 Pychram激活码_通用破解码,https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

  • vue刷新页面的方法_vue局部刷新页面

    vue刷新页面的方法_vue局部刷新页面业务需求/问题描述在项目中经常遇到一个问题,例如新增完表单数据和需要重新刷新页面。类似的业务还有很多。这时我们可以考虑的方式如下①(推荐)v-if刷新页面,并依赖注入(不太清楚的小伙伴可以看我之前的文章)//父组件<子组件v-if=’load’>exportdefault{ data(){ load=true }, methods:{ refresh(){ this.load=false this.$nextTick(()=>{ t

    2022年10月17日
  • Flutter 使用插件打开相册、相机「建议收藏」

    Flutter 使用插件打开相册、相机「建议收藏」需求:image_picker的使用,点击按钮底部弹出相册、拍照选择框,实现具体功能1:引入插件pubspec.yaml增加image_picker:^0.7.5+22:android添加androidx兼容gradle.properties增加android.useAndroidX=trueandroid.enableJetifier=truebuild添加依赖androidTestImplementation’androidx.tes…

  • pycharm 2021.12.13最新版激活【最新永久激活】

    (pycharm 2021.12.13最新版激活)这是一篇idea技术相关文章,由全栈君为大家提供,主要知识点是关于2021JetBrains全家桶永久激活码的内容IntelliJ2021最新激活注册码,破解教程可免费永久激活,亲测有效,下面是详细链接哦~https://javaforall.cn/100143.html2JTX0APX6F-eyJsa…

  • python量化分析系列之—使用python获取股票历史数据和实时分笔数据

    python量化分析系列之—使用python获取股票历史数据和实时分笔数据使用python获取股票历史数据和实时分笔数据

  • 查看Linux内核版本_ubuntu升级内核命令

    查看Linux内核版本_ubuntu升级内核命令想知道您的系统使用哪个Linux内核版本吗?以下是在Linux终端中检查内核版本的几种方法。当您需要了解系统上正在使用的确切Linux内核版本时,您可能会遇到这种情况。借助功能强大的Linux命令行,您可以轻松地找到答案。如何找到Linux内核版本在撰写本文时,我正在使用Ubuntu18.04。但是这些命令是通用的,可以在Fedora,Debian,CentOS,SUSELin…

    2022年10月13日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号