通过数值逼近（与符号计算相反）解决数学问题的算法。
我正在运行梯度下降来找到非线性方程组的根，我想知道你如何检测该方法是否陷入局部最小值，因为我相信随着设置...
Python 中的 Hermite 插值，无需使用 scipy 库
def hermite_interpolation(x, y, yp, xi): n = len(x) - 1 结果 = 0.0 对于范围 (n + 1) 内的 j： 项 = y[j] 对于范围 (n + 1) 内的 i： 如果 i != j 且 x[j] != x[i]： ...
我正在实施 Dekker 的数值方法来求根，但我对更新间隔有点困惑。维基百科上提到了这一点 选择新的对位点使得 f(a_k+1) 和 f(b_k...
我正在使用 Verlet 积分（更具体地说，是 Verlet 积分维基百科页面“非常量时间差”部分中的最后一个方程），以及 C++ 和 Eigen 库...
我正在尝试实现 Dekker 的求根方法，但我似乎无法找出算法似乎无法正确收敛到解决方案并且有时会突然停止的问题......
如何在没有库的情况下数值求解隐式函数（方程）a*y^2 - b*x*y + c = 0？
我想知道如何用 x 来求解 y a*y**2 - b*x*y + c = 0 使用数值方法，例如 Newton-Raphson、Bisection 等，无需库，除非用于定义...
我有一个Matlab代码 模型 = createpde(); GeometryFromEdges（模型，@lshapeg）； applyBoundaryCondition(模型,"狄利克雷", ... “边缘”，1：模型.Geom...
我有一个函数 f(x1,…,xN)，其中变量 xi 仅采用二进制值，xi = 0 或 1。输入空间是以汉明距离为度量的空间。 我想在数值上优化这个函数...
我在 np.gradient 中遇到了相对平滑的数据的高方差。 假设我们要计算 dx/dt： 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt x = np.array([13.11149679, 13.2141427 ,...
我在 np.gradient 中遇到了相对平滑的数据的高方差。 假设我们要计算 dx/dt： 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt x = np.array([13.11149679, 13.2141427 ,...
我正在开发一个用 JAX 编写的微分方程求解器。我遇到的常见工作流程是这样的： 将 jax.numpy 导入为 jnp 从 jax 导入 jit # 要集成的函数。 @jit 定义
我编写了一个 Gauss Seidel 函数，它接受一个系数矩阵、一个保存常量的向量数组和一个整数 n，其中 n-1 是矩阵的长度。它返回一个数组，其中包含...
我想获得 scipy 函数 scipy.optimize.curve_fit() 背后的一些直觉。下面，我写了一个简单的版本，试图最小化最小二乘函数之和（权重均为1），
在我目前正在阅读的书中，确定性数值方法（例如基于梯度的：牛顿方法，无梯度：Nelder-Mead）和算法方法（例如...
我正在尝试求解包括空气阻力在内的自由落体的常微分方程。 因此我将 ODE 定义为： def f(v, g, k, m): 返回 g - k/m * v**2 我认为这应该代表系统
我有两个 3D 点，例如 a = (100, 100, 10) 和 b = (0, 100, 60)，并且想通过这些点拟合一条线。 我知道，3D 直线方程可以有不同的形状： 矢量形式： ...
尊敬的各位专家、各位老同事： 我想使用提供的代码创建下图。然而，这似乎花费了很长的时间，但没有明显的结果。 这是我的
k=3;m=123 欧拉C_x=[A] eulerC_v=[0] eulerC_t=[0] c=0 对于 np.arange(t_euler,t_max,dt_euler) 中的 z： a=-(k/m)*欧拉C_x[c] eulerC_v.append(eulerC_v[c]+a*dt_euler) eulerC_x.append(eulerC_x[c]+
我需要一个服从麦克斯韦分布的 10,000 个随机数。 对于正态分布，我知道我必须使用 Box-Muller 变换，但我的问题是为什么 正常分布...
我正在寻找一个R代码，它可以帮助我在以数字形式给出pdf时找到分位数。也就是说，假设我的数据是 x = c(0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65...