问题标签 [simpsons-rule]

我是并行范式的新手。我已经让算法以串行形式工作，但我无法并行化它。

我四处询问，有人告诉我，我的程序编写方式不能并行化。

这是我的序列号，任何帮助或建议将不胜感激。

我正在尝试在 Matlab 中使用辛普森一家 3/8 规则求解数值积分，但我遇到了一个我不明白的错误。

图片中给出了数值积分（公式 26）

这也是 Matlab 代码：

任何人都可以帮助解决这个问题，

所以，我有数学表达式（如下图所示），b的积分限制设置在0到1的范围内，因为b值的变化，那么c的值（积分结果）将是每个 b 不同。

另外，我只能使用辛普森规则来解决集成问题。这个任务的主要目标是得到 f(c,b) 的图

在此处输入图像描述

我有一个数字问题。我需要计算两个函数 f1 和 f2 之间的积分。我已经有一个功能的辛普森规则。我需要帮助如何找到两个函数之间的截距并计算积分。

我正在尝试使用 R 中的数值积分来近似以下积分：

,

其中函数mu由以下公式定义：

为此，我将复合辛普森规则实现为 R 中的一个函数，它以函数 ( integrand )、积分区间 ( [a,b] ) 和所需子区间数 ( n ) 作为参数。

我已经在各种不同的数学函数上测试了我的代码，它似乎工作得很好。但是，当我尝试逼近图中所示的积分时，逼近变得很大。

我的方法是首先根据其复合辛普森近似值将内部积分定义为 R 中t的函数。然后，再次使用复合辛普森规则，以便通过将内部近似视为函数来计算外部积分被整合。

这样做时，内部近似值在自行计算时是正确的，正如预期的那样，但是整个表达式的近似值变得太大，我似乎无法弄清楚为什么。

我正在将近似值与 Maple 给出的近似值进行比较；使用t=20自己计算的内部表达式应该给出0.8157191，整个表达式应该是12.837。R 正确计算0.8157191，但给出整个表达式的32.9285 。

我尝试使用许多不同的数学函数进行简化，并使函数独立于 R 中的t，但似乎都导致相同的错误。所以，总而言之，我的问题是，为什么只有外部积分被错误地近似？

我将非常感谢任何提示或指针 - 我在此处包含了说明问题的代码：

我已经实现了这个even()& simpson()，但无法获得所需的输出。其次得到错误([lambda x:1/x, 1, 11, 6]]) ZeroDivisionError: division by zero，无法理解这个错误

我正在尝试这样的预期输出

我试图创建一个使用辛普森方法解决多个积分的函数，我的想法是一个一个地积分变量，直到我只剩下一个积分，但我的递归深度超出了我无法弄清楚为什么（我怀疑那是因为我使用了太多的 lambda 函数）。

我正在尝试整合下面给出的加速度信号：

对于时间点：

我正在尝试使用 scipy 中的 simps/simpson 包将 f2 集成到 t 上，但出现此错误：

如果我使用如下列表理解创建加速数据点列表：

并使用辛普森规则整合上面的列表，我只得到一个值：

这是

我想在 t 期间使用 Simpson 积分获取速度值的列表或数组，并将其与 t 进行对比，但我尝试的一种方法是给我一个错误，而另一种方法是给一个奇异值。

我有一个从文本文件中绘制数据的代码。我已经能够得到每条数据曲线下的积分（虽然我不知道其中任何一个的函数，所以我只使用了integral = s = simps(y, x)）。其中两条曲线水平重叠，我想找到它们重叠的区域。我遇到了一些麻烦，因为它们水平重叠而不是垂直重叠。我想找到的区域在这张图上的红线和蓝线之间。bx数据位于 numpy 数组中，例如by蓝线 x 和 y 以及rx红线ryx 和 y。我没有代表它们的功能，但它们在那些数组中。很难知道如何进行，因为它们水平重叠而不是垂直重叠。

我非常坚持这一点，任何帮助将不胜感激！

更新：JohanC提供的代码在一种情况下有效，但在另一种情况下无效。非工作图最初看起来像这样（我只想要那个中间重叠区域的区域），但是当它也不起作用时，我尝试限制 x 范围以希望帮助代码缩小到我想要的区域就像要计算的重叠积分（这是具有较小范围的图的来源）。我目前使用的代码是JohanC提供的代码，如下所示：

文本 bx、by、rx 和 ry 是如下值：

我不确定为什么此功能不适用于一种情况，但会完美地适用于另一种情况。任何帮助将不胜感激！

所以基本上我们必须为一个班级项目完成这个，到目前为止我有以下内容：

对于数字部分。但是我不知道如何用任何通用函数来表示它，然后生成一个对数图，显示当面板数量 N 变化时，数值计算结果和分析得出的结果之间的分数误差。有人可以告诉我如何去做吗？