获取一个三维轮廓的最小表面解决方案
我有一组定义3D轮廓的3D点。 我想要做的是获得与此轮廓相对应的最小曲面表示形式 (详见维基百科上的最小曲面)。基本上，这需要解决一个非线性偏微分方程。 在Matlab中，使用pdenonlin函数(请参阅Matlab的文档)几乎是很简单的。这里可以找到一个使用它来解决最小曲面问题的示例：Un...
如何将ODE系统与FEM系统结合
我有一个动态模型，它被设置为一组ODE（僵硬的）。目前，我使用CVODE（来自Assimulo python包中的SUNDIALS软件包）来解决这个问题，一切顺利。 现在，我想将一个新的3D散热器（具有温度相关的热参数）添加到问题中。我的想法是，不必从头开始编写所有3D热方程的方程式，而是使...
如何在函数定义中使用Maxima的diff？
我想在另一个函数中使用一个函数的导数。在Maxima中应该如何实现？ 例如：f(x) := 2*x^4; g(x) := diff(f(x),x)-8; 现在g(x)的结果如预期所示为8x^3-8，但是g(0)会出错，因为diff(f(0),0)没有意义。那么我该如何正确定义g呢？
有没有适用于普通微分方程（ODE）求解器的C++库？
更具体地说，我对8阶Dormand-Prince嵌入式方法感兴趣，它基于龙格-库塔法和刚性方程。 我使用Numerical Recipes 3，但我经常在编译其库时遇到麻烦。我想知道是否有其他替代方案。
Python中的数值ODE求解
我该如何用Python语言数值求解ODE？ 考虑以下方程： \ddot{u}(\phi) = -u + \sqrt{u} 以下是需要满足的条件： u(0) = 1.49907 并且 \dot{u}(0) = 0 有了这个限制条件 0 <= \phi <= ...
使用辛积分器的帮助
我正在尝试开发物理模拟，并希望实现一个四阶辛积分方法。问题在于我的数学肯定有问题，因为使用辛积分器时我的模拟根本不工作（相比之下，对于模拟而言，使用四阶龙格-库塔积分器的表现还算可以）。我已经在Google上搜索了很长时间，但我只能找到与此相关的科学文章。我试图改编这些文章中使用的方法，但没有...
索引错误：索引10超出了大小为10的轴0的范围。
我正在为x轴网格、x向量和时间网格设定数值，但我已经为x（位置）设置了一个数组，它应该只在0到20之间，而t（时间）将从0到1000，因此可以解决热传导方程。但是每当我想要设置步骤数为10时，就会出现错误： "Traceback (most recent call last): File "...
在Java中的微分方程
我正在尝试用Java创建一个简单的SIR流行病模型的仿真程序。 基本上，SIR由三个微分方程组成： S'(t) = - l(t) * S(t) I'(t) = l(t) * S(t) - g(t) * I(t) R'(t) = g(t) * I(t) 其中，S - 易感人群，I - 感染者...
使用Python通过Scipy和Numpy将数据拟合到ODE系统
我在通过Scipy和Numpy将我的MATLAB代码转换为Python时遇到了一些困难。我不知道如何找到最佳参数值（k0和k1），以便将我的ODE系统拟合到十个观测数据点上。我目前有一个关于k0和k1的初始猜测。在MATLAB中，我可以使用称为“fminsearch”的函数。该函数接受ODE系...