通过数值逼近（与符号计算相反）解决数学问题的算法。
在 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 对其进行可视化。目的是将其与前向欧拉进行比较并评论差异。我的有限元工作...
我正在寻找一些四边形网格生成的示例代码（任何语言）。然而，这似乎是一项相当艰巨的任务！ 我不挑剔，我想至少对有孔的多边形进行网格划分，什么都没有......
迭代求解器比基于 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 --> 是的，它是一个数字