可能除以“零”的数值积分

计算科学 正交 准确性
2021-12-19 07:47:30

我正在尝试整合

01t2n+2exp(αr0t)dt

这是一个简单的转换

1x2nexp(αr0x)dx

使用是因为难以用数值近似不正确的积分。然而,这确实导致了评估接近零的新被积函数的问题。由于间隔的长度仅为 1,因此很容易获得适当数量的正交节点(因此可比较的可以做得非常小),但是当积分接近零时我应该考虑什么?t=1xdt

在某种程度上,我认为简单地采用是一个好主意,其中是一个小数. 但是,我应该选择什么号码?它应该是机器epsilon吗?除以机器epsilon是一个很好量化的数字吗?此外,如果除法我的机器 epsilon(或接近它)给出一个非常大的数字,那么取将变得更大。ϵ1t2n+2exp(αr0t)dtϵexp(1ϵ)

我应该如何解释这一点?有没有办法让这个函数有一个明确定义的数值积分?如果没有,集成该功能的最佳方式是什么?

2个回答

这可以通过部分积分来完成: 并继续归纳 使得

1xeax=1axeax11a1eax=eaa+eaa2=a+1a2ea
1xkeax=1axkeax1ka1xk1eax=eaa+ka1xk1eax
I(k)=eaa+kaI(k1)
I(0)=eaa

看看QUADPACK它具有在(半)无限域上集成的例程。