通过数值逼近(与符号计算相反)解决数学问题的算法。
我正在尝试在Python中实现一个简单的MD模拟(我是新手),我将使用LJ电位和力方程与Verlet方法一起使用。 我昨天发布了另一个有关同一代码的问题,终于感谢您取得了进展! 这是我的代码:
为了确保它可以在最简单的系统上使用它,即由距离分隔为4Å的2个原子。
如何有效地确定在Skyfield的卫星通行证期间是否越过某些方位角标题 我试图确定Skyfield即将到来的卫星通行证是否会越过一定方位角,例如北0.0度。现在,我在整个p上“手动”迭代1秒步骤...
(背景故事是我有一个跟踪平台,由于电缆包装的0.0-360.0,它的方位范围有限方位角旅行极限。)
在 20 世纪 70 年代的经典著作中,Richard Hamming 建议将浮点计算设置为更容易发生下溢而不是溢出。 例如,如果您想计算 1 - 1/(exp(x)...
TL;DR:问题是关于乘法精度 我必须将矩阵 A (100x8000)、B (8000x27) 和 C (27x1) 相乘。 由于矩阵 B 和 C 是常数,而 A 是变量,所以我更喜欢计算...
时间序列分散的步数近似,使用 BIC 改变每 K 个步数的均值
首先,所使用的合成数据是按以下方式生成的: 将 sympy 导入为 sp 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt 随机导入 导入数学 np.随机.种子(2) n_样本 = 180
我使用前向差分进行数值微分的 k 阶导数代码不起作用。它给了我 -38 而不是 -2.8,我找不到问题 #k阶导数 #f'(x) = 1/h * (yn+1 ...
积分为 其中变量遵循以下分布: 因此,积分变为: 现在,我在 python 中的代码是将表达式与 \(W_{1:3}\) 集成为 \(1,2,3\) 尊重...
我有一个名为 trap_1D 的 Python 函数,它使用复合梯形规则计算单变量函数 f 在区间 [a,b] 上的(近似)积分。现在我想评估一个
我尝试使用牛顿法近似 255 / sqrt(x) 以避免在架构上使用除法或 sqrt() 操作,因为这些操作的成本很高。 我的推导是: y = 255 / 秒...
我这样做是为了使用射击方法解决BVP功能。 我不明白为什么当我的步长 h = 0.1 时,限制 x_end 超出了我在函数solve下的 while 循环中设置的值...
我正在编写代码,我有一个函数可以用满足条件的值来计算不满足条件的值,但是我在管理方面遇到了很多麻烦...
我正在编写代码,我有一个函数可以用满足条件的值来计算不满足条件的值,但是我在管理方面遇到了很多麻烦...
我正在实现隐马尔可夫模型的前向算法(算法见下文)。为了防止上溢/下溢,我使用对数概率,并使用对数和表达式技巧来配合...
我正在尝试使用 SciPy 优化来找到最小化两点之间悬挂绳索势能的函数,即悬链线问题。目标函数定义为 那个...
从十进制数计算分数的算法是已知的,并且有足够的例子。就我而言,我需要找到给定小数的分数,其中 A = N/D 已优化,但是
我正在使用 scipy 的solve_bvp 求解常微分方程,以获得一系列特征态 phi_n 和特征值 omega_n。等式看起来像 y'' = -(omega_n - f(x)) ** 2 * y...
我正在尝试用 Javascript 对弹簧上的质量进行 Runge-Kutta 实现,并使用 D3 对其进行可视化。目的是将其与前向欧拉进行比较并评论差异。我的有限元工作...
我正在寻找一些四边形网格生成的示例代码(任何语言)。然而,这似乎是一项相当艰巨的任务! 我不挑剔,我想至少对有孔的多边形进行网格划分,什么都没有......