通过数值逼近（与符号计算相反）解决数学问题的算法。
我正在实现隐马尔可夫模型的前向算法（算法见下文）。为了防止上溢/下溢，我使用对数概率，并使用对数和表达式技巧来配合...
我正在尝试使用 SciPy 优化来找到最小化两点之间悬挂绳索势能的函数，即悬链线问题。目标函数定义为 那个...
从十进制数计算分数的算法是已知的，并且有足够的例子。就我而言，我需要找到给定小数的分数，其中 A = N/D 已优化，但是
我正在使用 scipy 的solve_bvp 求解常微分方程，以获得一系列特征态 phi_n 和特征值 omega_n。等式看起来像 y'' = -(omega_n - f(x)) ** 2 * y...
我正在尝试用 Javascript 对弹簧上的质量进行 Runge-Kutta 实现，并使用 D3 对其进行可视化。目的是将其与前向欧拉进行比较并评论差异。我的有限元工作...
我正在寻找一些四边形网格生成的示例代码（任何语言）。然而，这似乎是一项相当艰巨的任务！ 我不挑剔，我想至少对有孔的多边形进行网格划分，什么都没有......
迭代求解器比基于 LU 分解的直接求解器更稳定吗？对于基于 LU 的求解器，我们总是有 cond(A) < cond(L) * cond(U), so factorization amplifies numerical inaccuracy. So...
Burden & Faires 中的浮点问题以及高斯消元法示例
我正在阅读 Burden & Faires 第 9 版第 360 页第 6.2 节。它包含一个示例，旨在说明在没有行交换的情况下高斯消除的数值不稳定性......
有没有一种简单的方法可以判断 varchar 是否是数字？ 示例： abc123 --> 没有号码 123 --> 是的，它是一个数字
有没有一种简单的方法可以判断 varchar 是否是数字？ 示例： abc123 --> 没有号码 123 --> 是的，它是一个数字
我目前正在努力绘制 R 中的下一个函数。 (F1 / F2^m * (R1 - F1)^((1 - s) / s)) - (m * (((R1 - F1))^(1 / s) + ((R2 - F2))^ (1 / 秒)) / (F1^m + F2^m)) = 0 我想将 F2 绘制为
如果我们计算像 A=[1 2 3;4 5 6;7 8 9] 这样的秩亏矩阵的 QR 分解，R 矩阵的一些对角项将非常小。然而，当在 Ju 中进行数值计算时...
使用`solve_ivp`和`LSODA`来求解复杂的ODE
我正在尝试求解复杂的微分方程组。这些方程是刚性的，因此我需要使用一种可以处理复杂的 ODE 和刚度切换的方法。我已经降落在
我需要求解一些常微分方程 $ rac{d y}{d x} = f(x) = x^2 ln(x)$ 并继续在限制 0 之间创建数组 xpt。 <= xpt <= 2. Because I have to be careful ...
我试图用 C 级 LAPACK 替换我的 cython 代码中的 Scipy 的 eigh 例程。本质上我想计算对应于埃尔米特 m 的最大和最小特征值的特征向量...
我有这行涉及矩阵求逆的代码： X = A @ B @ np.linalg.pinv(S) A 是 n × n 矩阵，B 是 n × m 矩阵，S 是 m × m 矩阵。 m 小于 n 但通常不是
最初，我有以下等式： mgv×sin(lpha) = CdA× ho(v^2 + v_{风}^2 + 2vv_{风}cos(\phi))^(3/2) 我可以将其表示为以下非线性方程： K × v^(2/3) = v^2 + v...
所以下面是我为 Julia 准备的代码片段： 使用Printf 使用绘图 gr() # 参数定义 Ny = 21 # y 方向的点数 Aspect_ratio = 10 #
在 Python 中创建球体 S^4 上的一个简短的 d 密集点列表
我是一名数学家，这是我第一次使用 Stack Overflow，如果问题不充分或者有更好的地方可以问这个问题，我很抱歉。 我想知道是否有一个标准...
我正在计算积分值 ∫₀1 ∫₁² (r - u)r drdu 我想了解 R 中如何处理以下错误。 环顾论坛，我意识到