2

Matlab 函数中似乎存在quad使用求积公式评估积分的错误。跑步

四边形(@(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命令似乎工作正常。这是一个已知问题还是我错过了什么?

4

2 回答 2

2

阅读quad文档,看起来使用它可能是一个好主意quadgk。在这种情况下,它给出了积分区间 (0,100) 的正确结果:

>> quadgk(@(x) (2/sqrt(2*pi))*(x.^2).*exp(-x.^2/2), 0, 100)

ans =

    1.0000
于 2013-09-21T15:05:44.860 回答
1

quad以和八度音程运行这两个表达式的quadl结果为 1。机器可能有问题。

于 2013-09-21T14:36:20.250 回答