Matlabquad
函数中似乎存在一个使用正交公式评估积分的错误。跑步
四(@(x)(2 / sqrt(2 * pi))*(x。^ 2)。* exp(-x。^ 2/2),0,10)
给出1.0000
正确的答案,但将上限说成100,即
四(@(x)(2 / sqrt(2 * pi))*(x。^ 2)。* exp(-x。^ 2/2),0,100)
退还3.4715e-8
。似乎也是如此quadl
。integral
但是,该命令似乎可以正常工作。这是一个已知问题还是我缺少什么?
阅读quad
的文档似乎可以quadgk
改用它。在这种情况下,它以积分间隔(0,100)给出正确的结果:
>> quadgk(@(x) (2/sqrt(2*pi))*(x.^2).*exp(-x.^2/2), 0, 100)
ans =
1.0000
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句