多元二分法
我需要一个算法来执行二维二分法，以解决2x2非线性问题。例如：两个方程f(x,y)=0和g(x,y)=0，我想同时解决它们。我非常熟悉一维二分法（以及其他数值方法）。假设我已经知道解在限制x1 ^ | C D y2 -+ o-------o | | ...
我该如何计算两个点序列之间的“差异”？
我有两个长度为n和m的序列。每个序列都是形如(x,y)的点序列，代表图像中的曲线。我需要找出这些序列有多不同（或相似），考虑到以下事实： 1. 一个序列可能比另一个序列长（即一个序列可能只有另一个序列的一半或四分之一长，但如果它们大致跟踪相同的曲线，它们就是相同的）。 2. 这些序列可能是...
scipy python模块中的odeint()函数是如何工作的？
我是一名物理学生，对于数值解 ODE 问题非常感兴趣。我通常会使用 Runge-Kutta 方法并用 C 语言编写自己的求解器。 最近我学习了 Python，并且使用 SciPy 的 odeint 函数来求解 ODE 问题。但是我很担心函数算法的工作原理，因为它不接受步长参数。那么，我该如何...
Lisp、OCaml或者什么语言用于Runge Kutta算法？
你会为解决具有以下特点的系统提出哪种语言： 一阶微分方程 复杂变量 N维空间 使用4阶Runge-Kutta或类似方法。 速度非常重要，但愿意为以下因素做出牺牲： 优雅（简洁）的代码 灵活性+可扩展性 我主要在Lisp和OCaml之间选择，但欢迎任何其他建议。 谢谢！
整数除法时最小化舍入误差
我正在尝试通过将一个整数乘积除以另一个整数乘积的比率来形成一个双精度浮点数（64位）。我希望以减少舍入误差的方式完成此操作。 我熟悉Kahan加法和减法。哪些技术适用于除法？ 分子是许多长值（成千上万）的乘积，同样是分母。我希望防止溢出和下溢。（一个应用是通过在足够数量的项后停止来估计无限...
C++库中的快速梯度下降实现？
我想运行梯度下降优化算法来最小化变量实例的成本。我的程序非常耗费计算资源，因此我正在寻找一个拥有快速GD实现的流行库。请问推荐哪个库或参考资料？
时间修正的Verlet数值积分公式
以下是翻译的结果： 在网络上有一个常用的Verlet积分公式，由Johnathan Dummer提出，称为时间校正Verlet。然而，我已经阅读了几篇论坛帖子，发现在某些条件下，人们使用此公式会得到奇怪或意外的结果。 Johnathan Dummer的公式如下： x1 = x + (x ...
我的GMRES实现有什么问题？
我正在尝试在Jupyter Notebook中实现GMRES（如果您不知道的话）, 这是一种用于解决线性方程组的迭代方法： 这是我的代码： import numpy as np def GMRes(A, b, x0, e, nmax_iter, restart=None): ...
用Python（NumPy）解决热传导方程。
我解决了金属棒的热方程，其中一端保持在100°C，另一端保持在0°C。 import numpy as np import matplotlib.pyplot as plt dt = 0.0005 dy = 0.0005 k = 10**(-4) y_max = 0.04 t_max ...
std::pow提供的分数指数近似值错误
这是我试图做的事情： double x=1.1402 double pow=1/3; std::pow(x,pow) -1; 结果为0，但我期望是0.4465。 方程是(1 + x) ^3= 1.1402，找出x的值。