我现在正在阅读一篇批评一种方法的论文,因为它使用梯形规则,而不是“像高斯求积这样的更高级的求积规则”......
高斯求积法则要求,给定,计算一堆权重和节点。在知道这些权重和节点之后,积分可以计算为一个简单的和.
好的。
但站在脚上,没人知道什么和是。他们需要计算。这需要时间。
另一方面,您可以直接应用梯形规则。如果使用统一网格,则可以立即计算节点和权重,因此唯一真正的计算开销是实际评估总和。
那么我在这里错过了什么?为什么这里还有辩论?我们不是总是更喜欢梯形规则,因为权重和节点更容易计算吗?
举个例子,在 R 中,我可以计算所需的节点和权重使用函数 的点高斯求积gauss.quad。和,gauss.quad(4096)在我的计算机上大约需要一整秒的时间来计算。
相比之下,如果我们要在均匀网格上计算梯形规则所需的评估点,则需要 0.01毫秒!.
那么我错过了什么?当梯形规则的准确度略低且速度快一百万倍时,“高级正交规则”有什么大不了的?