问题标签 [quad]

我正在使用嵌套的 scipy.integrate.quad 调用来集成二维被积函数。被积函数由 numpy 函数组成 - 因此将输入数组传递给它比循环输入并为每个输入调用一次要高效得多 - 由于 numpy 的数组，它的速度要快约 2 个数量级。

但是....如果我只想在一个维度上集成我的被积函数 - 但在另一个维度上使用一组输入，事情就会下降 - 似乎“scipy”quadpack包无法做任何事情numpy 确实处理数组输入。有没有其他人看到这个 - 或者找到了修复它的方法 - 或者我误解了它。我从 quad 得到的错误是：

我在下面放了一个我正在尝试做的卡通版本 - 我实际上正在做的事情有一个更复杂的被积函数，但这是 gyst。

肉在顶部 - 底部正在做基准测试以表明我的观点。

这就是我想要做的：

将数值积分数组（在整个数据范围内）除以特定积分范围内的“四元”积分函数：

这给出了错误：“操作数不能与形状（n）（n，2）一起广播”，n是t循环的长度

或者换一种说法，对于循环中“t”的每一个值： （1）用trapz对一组数值数据进行积分，变成一个数组；(2) 在一定的范围内使用 quad 积分一个函数，也变成一个数组；然后 (3) 将一个数组除以另一个得到最终数组。

显然我缺少一些东西。我该怎么做？

我必须在我决定的范围内整合一个信号。

我用过这个方法：

y 是 500 微秒内 2500 个点的数组。例如，t0 是 700，t1 是 1500。0.2E-6 是信号的两个值之间的步长。

这是对的吗？

除了使用 lambda 函数将 python 2.7 中的函数集成到两个变量中的第二个之外，还有其他方法吗？例如，我想在 x 上集成这个示例函数：

这很容易通过以下方式实现：

我想知道是否有另一种方法可以做到这一点，例如使用类似的东西，我认为这更直观：

有人有替代解决方案吗？（是的，我可以将函数定义为sqrfn(x,a)并quad(sqrfn,0,10,args(2,))改为使用，但这不是重点）。

我正在尝试创建一个整数值数组，将在计算中进一步使用。问题是integrate.quad 返回（答案，错误）。我不能在其他计算中使用它，因为它不是浮点数；它是一组两个数字。

我需要使用四边形方法集成到一个循环中。我无法创建外部文件并调用它，因为积分方程中的变量会随着循环而变化，所以我为了制作简单而切断了循环。它在 quad 方法行中给我一个错误 - 不知道如何修复它 - 确切地说是 KrInitialIntegratedPart 我得到的错误消息是：

使用内联/子引用时出错（第 13 行）内联函数的输入不足。

quad 中的错误（第 68 行）y = f(x, varargin{:});

newproj4 中的错误（第 39 行）KrInitialIntegratedPart= quad('(sin(x)) (exp(-AInitial/aInitial)) (db)', 0, fe)

我一直在学习 QUADPACK，我使用 scipy.integrate.quad 来计算从 0 到无穷大的积分。它给出了非常好的结果，但现在我想知道该软件应用了哪种集成方法（QAGI、QAWF 等）。有没有办法打印技术？该软件是否应用了一些决策树？在此先感谢您的时间。问候。

我有一小段代码会产生错误消息（我认为是因为数学公式）。有人知道为什么吗？

错误日志：

谢谢！

以下示例在Scipy的参考页面上提供integration。

以下是我从IPython笔记本（在 cloud.sagemath.com 上）得到的错误：

我猜它还cloud.sagemath.com没有升级到最新版本，Scipy因此它缺少nquad. 但是，我所需要的只是集成两个变量，因此想要使用dblquad已经在云上可用的变量。

因此，我修改了最后一行以适应dblquad如下所示的语法：

但它仍然弹出错误：TypeError: 'int' object is not callable. 我的脚本中有什么错误？我在下面粘贴了整个错误消息：

编辑 1：我使用来自用户的输入Weckesser和Chen. 只是为了将来偶然发现这个问题的任何人的完整性。

我正在尝试计算具有多个变量的函数对 scipy 中的一个变量的定积分。这有点像我的代码 -

它返回此类型错误：

但是，如果我在 args 中输入一个数字，它就会起作用。但我不想这样做，因为我希望 y 保持为 y 而不是数字。有谁知道如何做到这一点？

编辑：对不起，不认为我是清楚的。我希望最终结果是 y 的函数，而 y 仍然是一个符号。