我想明白为什么要这样打电话
G = const * quadv(fun,a,b,tol)
返回的值不同于
lenB = length(b)
for 1 = 1:lenB
G(i) = const .* quadl(fun,a,b,tol)
end
以及如何实现两次调用返回相同的值？
编辑: 我也想为数组更快地运行 quadl。所以我会提出上面的问题:
如果可能的话，如何也为数组调用 quadl(类似于上面的 quadv 调用)，从而通过计算提高性能？
最佳答案
Matlab的quadl使用自适应Lobatto求积； quadv 使用自适应辛普森规则。
您得到的答案的质量可能取决于您承担的职能。您的测试乐趣是什么样的？
我假设您将相同的函数、限制和公差传递给两个调用。这将意味着差异小于公差。
这两种方法不同——高斯求积与辛普森法则不同:
关于matlab - (MATLAB) 了解 quadl 与 quadv，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/9383066/