问题标签 [pyomo]

我根据安装指南在我的 Ubuntu 15.10 上安装了 Pyomo，使用pip install pyomo它工作正常并返回Successfully installed pyomo PyUtilib appdirs ply nose six。（我先安装了 pip 和科学 python。）

但是，如果我现在用 测试它pyomo --help，我会得到pyomo: command not found.

我错过了什么吗？我是否需要以某种方式将 pyomo 添加到我的路径中？

我没有正确运行它吗？在快速入门中提到在 python 中导入 pyomo.core。我试过了，但它没有改变任何东西。

我是使用 Pyomo 的新手，所以如果这是一个基本问题，我提前道歉。好吧，我正在研究动力学模型，我的目标是估计动力学参数。在尝试我的复杂模型之前，我从一个“玩具模型”开始，以便更好地理解 Pyomo。

所以，我的玩具模型是一个由 3 个方程组成的简单 ODE 系统：

我的目标是估计参数 k1 和 k2。我稍微更改了本教程中的代码，如下所示：

运行此代码后，我收到以下消息：

但是，当我擦除代码中与 x3_meas 对应的所有行以及 .dat 文件中的数据时，它可以完美运行。

有谁知道是什么问题？

我的数据看起来像：

我有一个 Pyomo 模型，其形式为：

请注意，目标m.obj依赖于参数m.T。尝试运行它会给出错误：

使用一个值，例如expr=m.S[120]给出错误：

在这两种情况下，我的目标是相同的：优化S在地平线上的最大可能值。

如何创建表达这一点的抽象模型？

我有一个形式的连续时间 Pyomo 模型：

但是，当我运行上面的代码时，我收到以下错误：

看着list(m.t)，我看到，果然，122 不包括作为离散化点：

如何指定要包含在离散化中的点？

我正在使用 PyOmo 生成一个非线性模型，最终将通过 Ipopt 解决。模型如下：

运行模型会导致以下错误：

错误的第一部分来自 Ipopt，而第二部分来自 PyOmo。显然，这个问题与约束中的术语无关M.P[i]**0.75，但改变权力并不能解决问题（尽管2.0确实有效）。

我该如何解决这个问题？

你好，社区，

所以我有以下错误：

评估约束 5965 时出错：无法评估 sqrt'(0)。错误：pyomo.opt：求解器（asl）返回非零返回代码（1）错误：pyomo.opt：有关诊断信息，请参阅上面的求解器日志。

Solver (asl) 没有正常退出

实际上问题出在代码的以下部分：

有任何想法吗？？删除 sqrt 并对右侧的变量 limits_flows 进行平方，解决了这个问题，但是限制非常低（-6、-7 阶）并且问题变得不可行。

谢谢你。

编辑：

最终问题可以通过给随机初始条件来解决

我在 Windows 7 上使用 pyomo 4.3.11388 和 GLPK 求解器 4.57 在 python 中运行优化问题。同时运行 8 个调用与 multiprocessing.Pool GLPK 对 \Users\user\AppData\Local\Temp 有一个罕见的写入问题文件夹：这仅在数千次成功运行中发生一次。任何人都可以帮忙吗？谢谢！

错误：“[base]\site-packages\pyomo\opt\base\solvers.py”，575，求解器日志：GLPSOL：GLPK LP/MIP 求解器，v4.57 命令行中指定的参数：- -写 C:\Users\JGOREN~1\AppData\Local\Temp\tmp4fo_8eqx.glpk.raw --wglp C:\Users\JGOREN~1\AppData\Local\Temp\tmps7viej_q.glpk.glp --cpxlp C: \Users\JGOREN~1\AppData\Local\Temp\tmpjnqvqljx.pyomo.lp 从“C:\Users\JGOREN~1\AppData\Local\Temp\tmpjnqvqljx.pyomo.lp”读取问题数据... 25 行， 45 列，69 个非零值 223 行已读取 将问题数据写入 'C:\Users\JGOREN~1\AppData\Local\Temp\tmps7viej_q.glpk.glp'... 无法创建 'C:\Users\JGOREN ~1\AppData\Local\Temp\tmps7viej_q.glpk.glp' - 权限被拒绝无法以 GLPK 格式写入问题

我正在尝试使用 Pyomo + glpk Solver (Anaconda2 64bit Spyder) 运行一些简单的示例：

但我收到以下错误消息：

我从http://winglpk.sourceforge.net/下载了 glpk --> 解压缩 + 将部分添加到环境变量“path”。

希望有人可以帮助我-谢谢！

我有一个要解决的线性整数程序。我安装了求解器 glpk （感谢这个答案）和 pyomo。我写了这样的代码：

这产生了 file 的解决方案results.yaml。

我有很多问题要使用相同的模型但具有不同a的b、 和c值来解决。我想为、 和分配不同的值a，求解模型，获得 和 的解，并列出、、和。我阅读了文档，但示例仅将解决方案写入文件，例如.bcmodel.x[1]model.x[2]abcmodel.x[1]model.x[2]results.yaml

有什么方法可以从代码中访问解决方案值？

谢谢，

我正在使用 PyOMO 来模拟半批量反应。

考虑一个描述半间歇式反应器的 ODE 系统，其中一种反应物在t1单位时间内以给定的体积流量进料，反应一直持续到t _end，并且显然t1 < t _end。

要指定流程中的停止点，我可以使用条件规则（假设t1 = 3.5*60）：

这将产生不连续性（然后我的模型不会收敛）。我想要做的是使用一个 sigmoidal 函数，它将流转换为零而没有不连续性。

尽管我需要参考时间变量本身，但要实现 sigmoidal。

下面的 MATLAB 代码给了我想要的结果：

但是，此实现需要在函数中显式地使用时间变量。如何访问 PyOMO 规则中的时间变量？我尝试了以下方法，但出现“无法将标量组件 't_of_step' 视为数组”错误：

希望我已经清楚地描述了我所追求的。欢迎任何提示，

谢谢！ 萨尔