我正在尝试编写有关辛普森定律的程序。我正在尝试做的是用作错误,如图所示: . 在我写的代码中,Ih 是我的 f1,Ih/2 是我的 f2。如果错误没有发生,那么步骤就会减半。 但是我得到这个错误 我的代码是 我没有预料到任何错误,并且基本上期望能够正确计算。 ...
我正在尝试编写有关辛普森定律的程序。我正在尝试做的是用作错误,如图所示: . 在我写的代码中,Ih 是我的 f1,Ih/2 是我的 f2。如果错误没有发生,那么步骤就会减半。 但是我得到这个错误 我的代码是 我没有预料到任何错误,并且基本上期望能够正确计算。 ...
使用内置提供的辛普森积分方法 Scipy,以下正 y 值的数值积分计算为负积分值。这让我非常惊讶。 我是这方面的初学者。 任何帮助,将不胜感激。 我期望辛普森的积分值是正值,但我得到的是负值。 ...
我想从 Matlab Hermite 函数的输出中创建一个函数(例如,如果我们有 Hermite 函数[8 0 -12 0]的输出,它将是8x^3 - 12x多项式),然后使用Simpson's集成该函数3/8 规则。 我已经在 Matlab 中创建了一个函数,该函数使用此规则集成任何函数,并且我 ...
我有一个 N x N 复杂的 NumPy 数组 U_0,我必须对它进行操作: 首先,如何有效地增加零数组? 我可以简单地将它复制到 np.zeros((2N-1, 2N-1)) 但也许你们知道更好的方法。 感谢 Alexander Riedel 用 numpy.pad 的解决方案回答了这个问题 ...
我在 Python 中实施了辛普森规则,并测试了三次多项式的误差。 辛普森法则的误差界限表明它应该等于 0,但事实并非如此。 有人知道为什么吗? 我希望每个 $n$ 的误差都等于 0,但只有较小的 $n$ 才等于 0。 ...
我创建了一个 function,它可以使用辛普森规则求解二重积分,但它仅在积分极限为常量时有效 但它没有通过这些测试,并且会给出错误:TypeError: unsupported operand type(s) for *: 'function' and 'float'。 但我不确定如何修复它 ...
我有三列数据(时间(t)、推力(T)、距离(s)),我想从中计算完成的工作。 包含数据的文本 (logdata.txt) 文件结构如下: 我遇到的示例仅显示了在为 f(x) 设置自己的 function 时如何使用 scipy.integrate 函数,但是在这种情况下,我有原始数据(因此没有函数) ...
我有一个从文本文件中绘制数据的代码。 我已经能够得到每条数据曲线下的积分(虽然我不知道它们中的任何一个的函数,所以我只是使用了integral = s = simps(y, x) )。 两条曲线水平重叠,我想找到它们重叠的区域。 我遇到了一些麻烦,因为它们水平重叠而不是垂直重叠。 我想找到的区域在这 ...
我正在尝试整合下面给出的加速度信号: 对于时间点: 我正在尝试使用来自 scipy 的 simps/simpson package 将 f2 集成到 t 但我收到此错误: 如果我使用如下列表理解创建加速数据点列表: 并使用辛普森规则整合上面的列表,我只得到一个值: 这是 我想在 t 周期内使用 S ...
我已经实现了这个even()和simpson() ,但无法获得所需的 output。 其次得到错误([lambda x:1/x, 1, 11, 6]]) ZeroDivisionError: division by zero ,无法理解这个错误 预期 Output 我正在尝试这样 ...
任何人都可以帮助解决这个问题, 所以,我有数学表达式(如下图所示),b 的积分极限设置在 0 到 1 的范围内,因为 b 值的变化,那么 c 的值(积分结果)将是每个b都不一样。 另外,我只能使用辛普森法则来解决积分问题。 这个任务的主要目标是得到 f(c,b) 的图在此处输入图像描述 ...
我是并行化范式的新手。 我已经有了以串行形式工作的算法,但我无法对其进行并行化。 我四处打听,有些人告诉我,我的程序编写方式无法并行化。 这是我的串行代码,任何帮助或建议将不胜感激。 ...
我一直在使用 np.trapz 通过 arrays 进行集成,直到现在还没有遇到任何问题。 我得到了一个面积明显小于 1 的分布,因为它的最大值为 0.16,分布的宽度大约为 6,但似乎返回分布下方的面积 > 60。 这是我的代码: ...
以下源代码是辛普森积分规则的实现。 C# 源代码。 Output: . Python 源代码。 Output: C# 程序的 output 不同于 Python 的 output。 这是为什么? ...
所以我必须为辛普森的双积分编写一个算法,这样我才能以更快的方式找到答案。 我有一个指南,显示了编写这个程序的步骤。 在遵循本指南并在 netbeans 中运行它之后,我发现从程序中出来的值并没有真正接近真正的答案。 这是我在 java 中编写的代码: 我尝试了几次调试程序,但我还没有找到我遇到这 ...
我一直在尝试编写一个函数来使用复合辛普森规则来近似积分的值。 我所做的是将函数的每个值乘以 2 或 4(和 h/3)与pow(2,i%2 + 1)并减去边缘,因为它们的权重仅为 1。 起初,我认为它工作得很好,但是,当我将它与我的梯形方法函数进行比较时,它更不准确,这不应该是这种情况。 这是 ...
我正在尝试以这种方式从 f(x) function 获取数组: 请注意,x 是 integer,z 和 change 是在我的代码中建立的浮点数。 接下来我想要使用辛普森的规则,使用 scipy 中的 simps function。 我试过这个: 它说有一个错误我该如何解决这个问题? ...
我正在尝试使用 for 循环在 python 中编写辛普森规则,但我不断收到断言错误并且无法找出原因。 F: 样本: ...
我不确定我是否在这里很厚,但我不习惯 function 手柄等。 我试图找到我拥有的一些原始加速度计数据的积分,一位朋友建议我使用辛普森一家可能是最好的方法。 所以我发现 MATLAB 具有“quad()” function 但它接受 function 句柄作为参数。 如果我有两个向量,一个用于时 ...
我已经编写了一种方法,该方法通过合成Simpson规则近似定积分。 对于“简单”函数,例如e^(-x^2) , simpson函数起作用。 但是,对于更复杂的函数f(x) 其中generator(θ, plotsol)是一个以百分比表示的缺陷θ的函数,一个布尔值plots ...