Matlab中的符号积分

Question

我只是想在 matlab 中求解一个方程，但它确实给了我一个警告。 “警告无法找到显式积分” 经过一些搜索，人们建议用 Double() 方法包装我的方程，但即使这样，我还是收到了错误的答案（ans=0），我在mathematica/maple 上计算为 4.62240566。

我认为它与浮点有关，但我不确定如何修复它

syms t
int( (t^2+100)^(-1/2)*exp(-10^-3*(t^2 + 100)^(1/2)),t , 1, Inf) 
   Warning: Explicit integral could not be found. 

Answer 1

从http://www.mathworks.nl/help/symbolic/int.html 来看，您需要的代码是：

int(1/(exp(1000*(t^2 + 100)^(1/2))*(t^2 + 100)^(1/2)), t , 1,Inf)

或者你可以使用

t == 1..inf

并像示例中那样用vpa 包装它，而不是用double。