问题标签 [numerical-analysis]

我正在做一个项目，发现自己需要一个函数来至少获得 W(x) 的值的近似值，即Lambert W 函数，其中 x 可以是任何实数。我在 Java 中工作。当我搜索时，我在 java 中找不到 W 的任何实现。如果需要，我愿意自己编写实现代码，但我不确定现在将如何完成。任何朝着正确方向的推动都将不胜感激。

Monniaux 的优秀文章“验证浮点运算的陷阱”给出了意外数值行为的示例。然而，大多数示例依赖于扩展精度 FPU 或扩展精度运算（例如乘法累加）。除了带符号的零和涉及 NaN 的比较之外，当扩展精度硬件被禁用时，还有任何令人惊讶的数值行为的好例子吗？

使用双类型我制作了三次样条插值算法。这项工作看起来很成功，但是当计算出非常小的值时，存在大约 6% 的相对误差。

双数据类型是否足以进行准确的科学数值分析？

是否存在一个简单的、类似备忘单的文档来编译 R 中数学计算的最佳实践？有没有人列出他们的最佳实践？例如，它将包括以下项目：

1. 对于较大的数值向量，应该计算x而不是计算。这加快了计算速度。x^2x*x
2. 要求解系统 $Ax = b$，永远不要求解 $A^{-1}$ 和左乘 $b$。存在低阶算法（例如，高斯消除）

我确实在这里找到了一个不错的数值分析备忘单。但我正在寻找更快、更脏、更特定于 R 的东西。

我想计算向量，

s = AB u,

其中 s 和 u 是 N 维复数向量，A 是 N×M 复数矩阵，B 是 M×N 复数矩阵。当 A、B 和 u 的元素表示为浮点数时，以下两种方式中哪一种精度更高（有效位数更高）？

(1) 先计算 B u。

首先做矩阵向量乘法，

y = B u

然后，另一个矩阵向量乘法

s = A y

(2) 先计算AB。

首先做矩阵-矩阵乘法，

C = AB

然后，矩阵向量乘法

s = C u

有没有已知的一般规则？

顺便说一句，我知道方法（1）比方法（2）更有效。

我首先要说这感觉应该很容易......但这对我来说并不完全明显。我正在尝试使用 BrentOptimizer 来查找函数的局部最小值和最大值。我对这些的周期性有一个想法，我觉得我应该能够使用 BracketFinder 将最优值括起来，然后将其发送到 BrentOptimizer。

她是文档：http ://commons.apache.org/math/api-2.2/org/apache/commons/math/optimization/univariate/BracketFinder.html

因此，对于一个简单的情况，请考虑：

我们知道 Pi/2 有一个最大值，3Pi/2 有一个最小值。如果我从零开始并沿着函数移动，我将如何在 Pi/2 处寻找那个根？它实际上归结为构造函数参数和初始点。是否有任何最佳实践（假设您对函数的形状有所了解）可以用来以合理的方式设置这些参数？

谢谢

我正在使用卡尔曼滤波器进行一些数值估计和校正，并希望更好地估计我的 Q 和 R 参数，最好是动态的。 http://en.wikipedia.org/wiki/Kalman_filter#Estimation_of_the_noise_covariances_Qk_and_Rk

那篇文章提到GNU Octave是目前从数据中确定这些参数的最佳方式： http ://en.wikipedia.org/wiki/GNU_Octave#C.2B.2B_integration 不幸的是，它是为 Matlab 编写的，据说有一个 C++ 实现. 我在 C++ 方面非常薄弱，甚至不知道如何导入 C++ 库并在 XCode 中正确链接它。迄今为止，我所有的 C++ 库都包含在第 3 方 Objective-C 类中。

有没有人在 iPhone 上将 C++ 实现用于科学计算或工程应用程序？我将不胜感激有关如何使用 Objective-C 进行此类分析的任何指针或教程。

附加关键字：从数据中估计协方差

自协方差最小二乘 (ALS) 技术

噪声协方差

谢谢！

http://en.wikipedia.org/wiki/Binary_GCD_algorithm

这个 Wikipedia 条目有一个非常不令人满意的含义：二进制 GCD 算法曾经比标准 Euclid 算法效率高出 60%，但直到 1998 年，Knuth 得出结论，他的同时代算法的效率只有 15%。电脑。

又过了 15 年……今天这两种算法如何与硬件的进步相叠加？

二进制 GCD 是否继续在低级语言中优于欧几里得算法，但由于其在 Java 等高级语言中的复杂性而落后？还是现代计算中的差异没有实际意义？

我为什么在乎你可能会问？我今天碰巧不得不处理其中的 1000 亿个 :) 这是为生活在计算时代干杯（可怜的欧几里德）。

我到处看了看，但似乎我看到的例子只允许数字 0-9

我正在写一个勾股定理程序。我希望手机 (Windows Phone 7) 检查文本框中是否有任何alpha (AZ, az)、符号 (@,%) 或除数字以外的任何内容。如果没有，那么它将继续计算。我想检查一下，这样以后就不会出现错误了。

这基本上是我想要它做的一个糟糕的伪代码

txtOne-->任何字母？--否-->任何符号--否-->继续...

我实际上更喜欢一个命令来检查字符串是否完全是一个数字。

提前致谢！

测试用例 ：

35000

-> 数字的标准化科学记数法是 3.5 * 10E4

-> 工程符号为 35 * 10E3

执行此操作的简单算法将继续将数字除以 10，直到我们获得所需的符号。然而，这意味着算法将是 O（零的数量）。我们能做得更好吗？