问题标签 [simpsons-rule]

我不确定我在这里是否很厚，但我不习惯使用手柄等。

我试图找到我拥有的一些原始加速度计数据的积分，一位朋友建议我使用辛普森一家可能是最好的方法。所以我发现 MATLAB 有“quad()”函数，但它接受函数句柄作为参数。

如果我有两个向量，一个用于时间，一个用于加速，我如何将它们创建为函数句柄以发送到 quad()？

谢谢

我试图在我的逻辑中实现这个问题，以使用 epsilon = 1e-8 计算辛普森的方法自适应实现。以下是解释：

以下是我的辛普森方法代码：

以下是我的 f2(x) 函数：

**结果只有 1 次迭代**：

和预期结果类似这样的多次迭代：

我想用一些理解来探索代码，因为这个问题已经被要求提供其他值作为 n 但不是数组

我已经实现了以下逻辑，并针对不同的问题（数组范围）提出了这个问题。我得到了输出，但它没有通过 for 循环进行迭代，因为我已经给出了frange(start, stop, range)

解释

代码：

调用上述函数的另一个函数

以下是结果

它应该给我更多的迭代

任何人都可以帮助我迭代 for 循环以获得更多结果

我正在尝试使用 for 循环在 python 中编写辛普森规则，但我不断收到断言错误并且无法找出原因。

F：

样本：

我正在尝试以这种方式从 f(x) 函数中获取数组：

请注意，x 是一个整数，z 和 change 是我的代码中建立的浮点数。

我想要的下一件事是使用 scipy 中的 simps 函数使用 simpson 规则。我试过这个：

它说有一个错误我该如何解决这个问题？

我一直在尝试编写一个函数来使用复合辛普森规则来近似积分的值。

我所做的是将函数的每个值乘以 2 或 4（和 h/3）pow(2,i%2 + 1)并减去边缘，因为这些边缘的权重应该只有 1。

起初，我认为它工作得很好，但是，当我将它与我的梯形方法函数进行比较时，它更加不准确，这不应该是这种情况。

这是我之前写的代码的一个更简单的版本，它有同样的问题，我认为如果我稍微清理一下，问题就会消失，但是唉。从另一篇文章中，我了解到类型和我正在对它们执行的操作会导致精度损失，但我只是没有看到。

编辑：

为了完整起见，我将 e^x 从 1 运行到零

所以我必须为辛普森的双积分编写一个算法，这样我才能以更快的方式找到答案。我有一个指南，显示了编写这个程序的步骤。在遵循本指南并在 netbeans 中运行它之后，我发现从程序中出来的值并没有真正接近真正的答案。这是我用java编写的代码：

我尝试了几次调试程序，但我还没有找到我遇到这个错误的原因。如果您在我的代码中发现任何错误，请告诉我是否可以修复它。对于给定的示例，我得到的值是 0.9069281684027777，而不是正确的值 0.7838542。谢谢您的帮助。您还可以查看我遵循的指南来创建此程序。

以下源代码是辛普森积分规则的实现。

C#源代码。

输出：

.

Python 源代码。

输出：

C# 程序的输出不同于 Python 的输出。

这是为什么？

一段时间以来，我一直在使用 np.trapz 对数组进行集成，并且到目前为止还没有遇到任何问题。我得到了一个面积明显小于 1 的分布，因为它的最大值为 0.16，分布的宽度大约为 6，但似乎返回分布下方的面积 > 60。

这是我的代码：

我正在尝试实现这部分代码。最后一行出现错误“解析时出现意外的 EOF”。我尝试了各种方法，但似乎找不到问题所在 (T^T)

（我正在尝试实现 simpsons 规则。f 是一个函数；[a,b] 积分的区间和 N 是 [a,b] 中的部分数，也许这有点有用的上下文）

先感谢您！