我正在使用 Verlet 积分（更具体地说，Verlet 积分维基百科页面中“非恒定时差”部分中的最后一个方程）与 C++ 和 Eigen 库一起创建数值 N 体模拟。 我试过测量模拟的动量和能量守...

我正在寻找一种算法来生成随机无幂矩阵。 （我在这里有一个关于检查矩阵是否为无能的类似问题：确定矩阵是否为无能的数值算法) 我已经尝试了以下方法： 排名 1 矩阵 生成两个随机向量 v，w，其中 ，...

我有一个目标函数，我想使用 MLE 来查找模型的参数。这是我到目前为止的代码。此外，我的returnarr是股票的日志返回数组 以下是我的 returnarr 的前 30 个值，仅供参考： retu...

我对二维机制的动态模拟很感兴趣，有一些特定于我的领域的特点。我曾尝试过使用 simscape multibody，但我发现它不适合应用程序：强制 3D 会带来不必要的复杂性和计算时间，我无法修改应用程...

我正在为我的数值方法科目做一项任务，他们要求使用 Runge-Kutta 和 Taylor 方法来估计微分系统（Chua 电路）的解，其场由下式确定： F(x,y,z)=\[2\*(y-0.5\*x...

对于一个学校项目，我正在 Python 中实现一种有限差分方法，以数值求解以下反应扩散偏微分方程组： 到目前为止，我已经实现了以下目标： import numpy as np import mat...

我正在尝试评估（使用数值方法求解超越方程）当我改变初始猜测 M0 的值时，迭代次数会变化多少。我想我不应该期望它为如此不同的 M0 值返回相同数量的迭代，所以我认为这段代码可能不会使用我希望它使用的所...

我想在scipy函数scipy.optimize.curve_fit（）背后获得一些直觉。下面，我写了一个简单的版本，试图最小化最小二乘函数的总和（权重都是 1），模型是 sin（omega*t）。我...

我正在研究一个用 JAX 编写的微分方程求解器。我遇到的一个常见工作流程是这样的： import jax.numpy as jnp from jax import jit # Function t...

def hermite_interpolation(x, y, yp, xi): n = len(x) - 1 result = 0.0 for j in range(n + 1): ter...