MindOpt APL:一款适合优化问题数学建模的编程语言

2023-12-09 07:36

本文主要是介绍MindOpt APL:一款适合优化问题数学建模的编程语言,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

什么是建模语言

建模语言是一种描述信息或模型的编程语言,在运筹优化领域,一般是指代数建模语言。
比如要写一个线性规划问题的建模和求解,可以采用C、Python、Java等通用编程语言来实现计算机编程(码代码),也可以换采用建模语言。
本文将以阿里达摩院研发的MindOpt建模语言(MindOpt Algebra Programming Language, MindOptAPL,简称为MAPL)来讲解。MAPL是一种高效且通用的代数建模语言,当前主要用于数学规划问题的建模,并支持调用多种求解器求解。

代数建模语言工作原理

在数学规划领域,遇到一个实际问题时候,我们需要数学建模成优化问题模型、然后编程、然后计算优化结果,得到这个实际问题的解决方案。
在这个编程过程中,可以根据选用的计算工具——优化求解器提供的通用编程语言的API来编写代码,也可以采用建模语言来编写代码。如下示例,就是一个利用MAPL建模语言来进行一个优化问题码的代码。
image.png

  • 左边是数学模型,三要素:两个变量xa和xb,目标函数是最大化一个公式,约束是最下面两行,限定取值关系。
  • 中间是用MAPL建模语言编的代码。可以看到前面4行就表达清楚了左边的数学公式。最后“option solver mindopt;”是设置计算这个问题的求解器为mindopt求解器,“solve;”是执行求解。
  • 右边就是求解器的计算结果,xa = 3,xb=5,此时目标函数最大,是1050。

为什么要用建模语言

语法更简单(代码对比)

从上面我们可以看到建模语言可以方便地进行数学建模和求解的代码。这里我们对比一下建模语言和通用的编程语言,来看看用建模语言优势。
以下面这个问题为示例:

| 线性规划模型:
max x0 + 2 * x1 + 3 * x2 + x3
s.t. (-1) * x0 + x1 + 3 * x2 + 10 * x3 <= 20
x0 - 3 * x1 + x2 = 30
x1 - 3.5 * x3 = 0
0 ≤ x0 ≤ 40
0 ≤ x1
0 ≤ x2
2 ≤ x3 ≤ 3

我们使用 MindOpt APL 建模语言 和 MindOpt 求解器的 Python APIs,分别对上面的线性规划模型建模,并求解模型。

MAPL代码:

clear model;  #清除model,多次run的时候使用
option modelname test; #运行完代码之后会自动生成.nl和.sol文件  model是存放的地址,test是文件名#--------------------------
# twoTask.mapl
var x0 >= 0;   # 声明决策变量xa |
var x1 >= 0;
var x2 >= 0;
var x3 >= 2;
maximize Reward: x0 + 2 * x1 + 3 * x2 + x3;  # 声明目标函数
subto c1: (-1) * x0 + x1 + 3 * x2 + 10 * x3 <= 20;       # 声明约束
subto c2: x0 - 3 * x1 + x2 <= 30;
subto c3: x1 - 3.5 * x3 == 0;
subto c4: x0 <= 40;
subto c5: x3 <= 3;
#--------------------------option solver mindopt;     # (可选)指定求解用的求解器,默认是MindOpt
solve;         # 求解print "-----------------Display---------------";
display;        # 展示结果
print "目标函数值 = ",x0 + 2 * x1 + 3 * x2 + x3;

Python代码:

from mindoptpy import *if __name__ == "__main__":# Step 1. Create model.model = Model("test")try:# Step 2. Input model.# Change to minimization problem.model.ModelSense = MDO.MAXIMIZE# Add variables.x = []x.append(model.addVar(0.0,         40.0, 1.0, 'C', "x0"))x.append(model.addVar(0.0, float('inf'), 2.0, 'C', "x1"))x.append(model.addVar(0.0, float('inf'), 3.0, 'C', "x2"))x.append(model.addVar(2.0,          3.0, 1.0, 'C', "x3"))# Add constraints.model.addConstr(- 1.0 * x[0] + 1.0 * x[1] + 3.0 * x[2] + 10.0 * x[3] <= 20, "c1")model.addConstr(1.0 * x[0]              - 3.0 * x[1] + x[2] <= 30, "c2")model.addConstr(1.0 * x[1]              - 3.5 * x[3] == 0, "c3")# Step 3. Solve the problem and populate optimization result.model.optimize()if model.status == MDO.OPTIMAL:print(f"Optimal objective value is: {model.objval}")print("Decision variables: ")for v in x:print(f"x[{v.VarName}] = {v.X}")else:print("No feasible solution.")except MindoptError as e:print("Received Mindopt exception.")print(" - Code          : {}".format(e.errno))print(" - Reason        : {}".format(e.message))except Exception as e:print("Received other exception.")print(" - Reason        : {}".format(e))finally:# Step 4. Free the model.model.dispose()

从上面的例子可以看到,MAPL建模语言比较简洁,没有Python运行这么多复杂的创建、添加、异常捕捉和释放的过程,就聚焦在编个模型去求解计算,更易于理解和添加。上面的例子还只是线性规划,对于非线性规划的问题,Python的API会更复杂。而采用MAPL建模语言只需要表达清楚数学公式,对于调试模型修改更方便。

支持多种求解器,换求解器的时候不用重复编程

很多人选择建模语言,最大的原因是希望切换求解器方案。因为不同品牌的求解器的求解能力不一样,遇到一个问题数学模型调整了一行公式,可能之前选择的求解器就不支持了,需要更换求解器。
此时如果选择用各家求解器的API来编程,换一个求解器,就需要重新学习对应的API,重新码代码,维护起来困难。虽然业界也有通用的 .mps 和 .nl 的优化问题数据格式,但是熟悉不同求解器的调用数据计算的方法也很耗时,或者装对应的软件也很麻烦。这个时候,建模语言的优势就很大。

比如下面是MAPL代码中,只需要换一行,就能换求解器进行计算:

option solver highs;     # 更换求解器

更多MAPL支持的求解器,可以参考上一个博客MindOpt APL,可以支持调用几十种求解器的建模语言

建模语言也支持通用编程语言的API,如Python

有很多同学喜欢Python语言,更希望用Python编程。MAPL建模语言支持Python来调用,import maplpy后就用Python的方式来编代码,能继续享受一行代码换求解器的优良属性。可一看广告流量分配:曝光和转化均衡案例中的代码对比。

常见的建模语言

市面上的建模语言有很多个,需要看各家求解器支持的建模语言,比如MindOpt求解器支持如下4种建模语言:MAPL(MindOpt APL)、AMPL、Pyomo、PuLP。
image.png
其中MindOpt APL(MAPL)就是阿里达摩院自研的建模语言,是目前中国唯一一款代数建模语言。其他建模语言的描述大家可以点击上面的链接查看:https://opt.aliyun.com/platform/docs/htmldoc/solver

MAPL的优点

MAPL国内第一款拥有自主知识产权,完全自研的国产建模语言,提供了丰富文档学习以及案例参考。并且在电力SCUC等领域问题上建模性能优秀,对标或超越已有产品。对比AMPL等建模语言,部分语法上更灵活简单,后续也会支持向量化建模等特色能力,提升建模易用性等。

灵活性:

MAPL具有非常高的灵活性,可以用来建模和求解各种类型的优化问题,包括线性规划、整数规划、非线性规划等。它支持多种数学表达式和运算符,可以方便地表示复杂的数学关系和约束条件。如下是它支持数值计算。
在这里插入图片描述

易用性

MAPL建模语言采用了类似于自然语言的语法和结构,使得用户可以很容易地理解和编写优化模型。并且支持读写csv文件,使得数据的读取和写入变得容易,数据储存也很方便。
如下将结果print输出为csv表格:
image.png

更多选择

MAPL与MindOpt Studio平台集成,可以在线上环境使用,无需下载,并且支持调用多种求解器,可直接对比结果。还支持将输出mps文件,可以在不同的计算环境和操作系统之间进行导入和导出。

这篇关于MindOpt APL:一款适合优化问题数学建模的编程语言的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


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

相关文章

Linux部署中的文件大小写问题的解决方案

《Linux部署中的文件大小写问题的解决方案》在本地开发环境(Windows/macOS)一切正常,但部署到Linux服务器后出现模块加载错误,核心原因是Linux文件系统严格区分大小写,所以本文给大... 目录问题背景解决方案配置要求问题背景在本地开发环境(Windows/MACOS)一切正常,但部署到

MySQL磁盘空间不足问题解决

《MySQL磁盘空间不足问题解决》本文介绍查看空间使用情况的方式,以及各种空间问题的原因和解决方案,文中通过示例代码介绍的非常详细,需要的朋友们下面随着小编来一起学习学习吧... 目录查看空间使用情况Binlog日志文件占用过多表上的索引太多导致空间不足大字段导致空间不足表空间碎片太多导致空间不足临时表空间

Mybatis-Plus 3.5.12 分页拦截器消失的问题及快速解决方法

《Mybatis-Plus3.5.12分页拦截器消失的问题及快速解决方法》作为Java开发者,我们都爱用Mybatis-Plus简化CRUD操作,尤其是它的分页功能,几行代码就能搞定复杂的分页查询... 目录一、问题场景:分页拦截器突然 “失踪”二、问题根源:依赖拆分惹的祸三、解决办法:添加扩展依赖四、分页

Java中InputStream重复使用问题的几种解决方案

《Java中InputStream重复使用问题的几种解决方案》在Java开发中,InputStream是用于读取字节流的类,在许多场景下,我们可能需要重复读取InputStream中的数据,这篇文章主... 目录前言1. 使用mark()和reset()方法(适用于支持标记的流)2. 将流内容缓存到字节数组

解决若依微服务框架启动报错的问题

《解决若依微服务框架启动报错的问题》Invalidboundstatement错误通常由MyBatis映射文件未正确加载或Nacos配置未读取导致,需检查XML的namespace与方法ID是否匹配,... 目录ruoyi-system模块报错报错详情nacos文件目录总结ruoyi-systnGLNYpe

深入解析Java NIO在高并发场景下的性能优化实践指南

《深入解析JavaNIO在高并发场景下的性能优化实践指南》随着互联网业务不断演进,对高并发、低延时网络服务的需求日益增长,本文将深入解析JavaNIO在高并发场景下的性能优化方法,希望对大家有所帮助... 目录简介一、技术背景与应用场景二、核心原理深入分析2.1 Selector多路复用2.2 Buffer

解决Failed to get nested archive for entry BOOT-INF/lib/xxx.jar问题

《解决FailedtogetnestedarchiveforentryBOOT-INF/lib/xxx.jar问题》解决BOOT-INF/lib/xxx.jar替换异常需确保路径正确:解... 目录Failed to get nested archive for entry BOOT-INF/lib/xxx

SpringBoot利用树形结构优化查询速度

《SpringBoot利用树形结构优化查询速度》这篇文章主要为大家详细介绍了SpringBoot利用树形结构优化查询速度,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录一个真实的性能灾难传统方案为什么这么慢N+1查询灾难性能测试数据对比核心解决方案:一次查询 + O(n)算法解决

解决hive启动时java.net.ConnectException:拒绝连接的问题

《解决hive启动时java.net.ConnectException:拒绝连接的问题》Hadoop集群连接被拒,需检查集群是否启动、关闭防火墙/SELinux、确认安全模式退出,若问题仍存,查看日志... 目录错误发生原因解决方式1.关闭防火墙2.关闭selinux3.启动集群4.检查集群是否正常启动5.

idea Maven Springboot多模块项目打包时90%的问题及解决方案

《ideaMavenSpringboot多模块项目打包时90%的问题及解决方案》:本文主要介绍ideaMavenSpringboot多模块项目打包时90%的问题及解决方案,具有很好的参考价值,... 目录1. 前言2. 问题3. 解决办法4. jar 包冲突总结1. 前言之所以写这篇文章是因为在使用Mav