如何解决积分的奇怪结果？ [复制]

Question

我想计算积分(@(v)v^46*exp(-v^2),0,inf)。从数学上讲，它是可以计算出来的。但是，当我在 matlab 中运行时，它会给出 NaN。是不是电脑精度的问题？如果是，我该如何解决？我还尝试将 inf 更改为一个很大的不，比如从 10000,100000,1000000,... 值增加然后下降到接近 0 的非常小的不。

在极端情况下，我需要计算积分(@(v)(v-a)^200*exp(-v^2),0,inf)，所以我很欣赏基于最后一个积分的任何校正。

Answer 1

根本问题是不可能通过替换浮点数inf 来积分到无穷大。 integral 函数是纯数字的。

q = integral(fun,xmin,xmax) 对函数 fun 进行数值积分 使用全局自适应正交和默认误差的 xmin 到 xmax 公差。

您将需要使用符号工具箱，它处理符号而不是浮点值as per this answer。