一个收敛积分的九个整数,它具有大的积分极限,但不具有无限积分
用于复制的代码是下图:被积函数[x璝]:=((x+1)Abs[x]Log[(x+1)^2/Abs[x]^3])/((x+1)^2-Abs[x]^3);NIntegrate[Integrand[x],{x,-10000000,1000000
解答动态
x
是奇数,而这些贡献往往会抵消。理想情况下,我想找到一些方法让Mathematica用无穷积分范围正确地处理这个问题,这样我就不必每次需要做这样的积分时都记住用大的有限值来交换无穷积分极限可能是最好的方法,通过将负轴上的积分反映到正轴上来计算无穷远处的PV。
n积分[Integrand[x]+Integrand[-x],{x,0,infinity}](*2.813068425442558`*) @用户64494指出它是一个发散积分(在无穷远处)和方法-quot;主值gt;32,精度目标-gt;20,方法-“PrincipalValue”](*2.813065409383123(无警告)2.81306842613267712727(警告省略)*) 在-1和0处,似乎存在无穷导数,即所谓的“弱”奇点。在另外两个“奇点”处,减法对消可能会带来数值上的困难(例如,分母(1+x)^2-Abs[x]^3为零)。弱奇点往往使收敛困难,将积分分成两部分,
数值积分可以不存在任何问题零件:
f1=FullSimplify[Integrand[x],Assumptions-gt;0];f2=FullSimplify[Integrand[-x],Assumptions-gt;0];NIntegrate[f1+f2,{x,0,∞}](*2.81307*) 可以达到更高的精度(同样没有如下任何警告)
n整体[f1+f2,{x,0,1,∞},工作精度->;100](*2.81306842544297991357319331946335010348754523474644575568568457957304\5154609418329813658075049081527*)- End
免责声明:
本页内容仅代表作者本人意见,若因此产生任何纠纷由作者本人负责,概与琴岛网公司无关。本页内容仅供参考,请您根据自身实际情况谨慎操作。尤其涉及您或第三方利益等事项,请咨询专业人士处理。