问题标签 [nonlinear-optimization]

我试图最小化以下功能：

当 x+y=1 且 0<=x,y<=1 时。要在 Rsolnp 包中使用 gosolnp，首先，我定义了 cons 以在 eqfun 参数中使用它：

然后我应用了 gosolnp 函数：

答案应该是 x = 0 和 y = 1，但是你可以在每次运行 gosolnp 时尝试，你会得到新的点，在那个点上 func 大约为 0（不完全是）。Mathematica 和 Maple 非常快地对此函数进行优化，并给出了正确的答案，即 x = 0 和 y = 1，但是 R 中的每次运行都会给出一个不正确的新解决方案。

我还在阿拉巴马州或 DEoptim 中尝试了另一个优化函数 spg()，但问题仍未解决。

所以我的问题是：

1-有什么解决方案可以最小化R中的func吗？

2- R 和 Mathematica 的精度之间有什么区别，为什么 Mathematica 可以给我确切的答案，但 R 不能？

先感谢您

我在 R 中使用nloptr，但是，我想给我的模型更多的自由，因为最好的解决方案并避免过度拟合。我在这个问题的前面已经描述了我的问题：

基本上我有一个非线性问题要解决。我有一个函数来最小化和一些非线性约束。但我不想使用找到的最佳解决方案，因为它过度拟合了样本数据并给了我极值。因此，我想找到 N 个最佳解决方案，然后选择我想要的。

所以现在我想知道是否有一种方法可以找到nloptr在迭代过程中找到的 N 个最佳解决方案。除了nloptr还有其他方法吗？

问题：无法解决非线性方程中多个已知值的未知数

(I) 是未知数，(v) 是已知值，其余是常数。我可以为 v 的一个值解决这个问题。
需要的解决方案： (1)我将如何解决 (I) 的多个 v 值？
(2) 同样，如果常量 (Rs,Io,Il,Rsh,a) 位于 8760 行和 91 列的不同数据帧中，我该如何解决 (I)？

我有一篇论文中的编程公式，想给它一个解决特定问题的工具。作者将其称为线性规划 (LP) 实例，但我不确定。配方有点像如下：

我尝试通过cplexqcp函数对其进行编程（由于二次约束，但约束不包含任何x_i^2变量）。但是我收到CPLEX Error 5002: Q in %s is not positive semi-definite error. 这是具有非凸约束的非线性规划的一个实例吗？我可以使用CPLEX或使用NLP工具来解决它吗？我是LP/NLP员工的新手（不参加任何关于他们的课程），所以非常欢迎帮助解释我的问题答案的细节。

非常感谢。

我正在寻求优化模型的拟合，该模型描述了在已知直径和物种的映射树图中描述在 0.5m^2“垃圾陷阱”网络中收集的垃圾数量。选择的模型有两个因素，即垫料生产的异速生长比例和垫料移动距离的指数衰减。

但是，我们的陷阱数据包含来自多棵树的输入（这是标题中提到的“缺失级别”）：

到目前为止，即使是模拟数据，结果也很复杂。似乎有足够的直径和距离组合，功能应该受到一定的约束。这个分析是在我抄袭的一篇文章中进行的。我还尝试了对日志（Obs.Litter）的分析，我认为这是要走的路。但是我不确定我编写日志版本的方式是否会导致您期望执行的更好。

在这一点上，我想我只是在寻找对这种类型的“隐藏过程”拟合非线性回归或模型拟合问题更有经验的人的任何类型的建议（基于代码或概念）。下面包括数据模拟和各种可能性的代码。我在使用OpenBUGS中的贝叶斯层次模型估计这些参数方面取得了更大的成功，只有信息先验。

我想找到使函数最大化的向量p（即）的值。 p(1),p(2),p(3),...A(p)

我正在使用 MATLAB 来做到这一点，我发现fsolve我认为它可以帮助我。所以我做了功能A：

然后我需要求解一个方程组，它将是：

结果p向量将解决我的问题。

如何fsolve求解这个方程组 ( p0=[1 1 1]) ？

我正在使用 scipy 中的 newton() 函数来求解适用于单个值的特定非线性方程

我正在尝试对 (I,Io,RS,a,Rsh) 的 8760 个值执行相同的操作，这些值是单独的数据帧

我收到索引错误的错误：索引超出范围
我应该得到 8760*445 的 I 值

我是使用 Python 的新手。这些天我正在尝试学习新的优化算法和python。

Python中的CMA-ES优化算法源代码可以在这里找到： CMA-ES.py

我已经安装了所有必要的 Python 包（numpy、matplotlib、winpython 等）。运行源代码提供的测试功能也很容易，例如，

所需的定制目标函数来自数据的非线性最小二乘拟合：

数据集：23x3

非线性模型：a(1)--a(5) 是参数：

我正在尝试解决 R 中的非线性方程组，但它一直给我这个错误“要替换的项目数不是替换长度的倍数”。

我的代码如下所示：

D、V、r 和 s 是数字 [1:2508] 值，我认为这就是问题所在。如果我有单个值 1x1，它可以很好地解决它，但是，如果我插入具有 2508 个值的向量，它只会计算第一个 x1 和 x2，然后会出现带有我上面写的消息的警告。

感谢您的任何帮助。

丽娜

我有一个带有约束 K3 的化学动力学模型（2-组织隔室模型）（其中 K3 是速率常数）

我已经对等离子体函数和化学动力学模型进行了建模，以绘制输出特性

我想从下面的代码中估计速率常数 k3

下面列出了所有参数的初始估计值

请建议我一种从上述非线性动力学模型代码中估计 K3 参数的方法

在上述函数中，a1、a2、a3、b1、b2、b3、td、tmax、k1、k2 的参数值将保持不变。

我想知道 K3 值如何随着时间 t 值的变化而变化，为此我想最初以 t 间隔 t=0:0.1:60 估计 K3..

非常感谢任何帮助