1

这个问题是对我之前的问题的跟进,这里GSL 集成行为很奇怪

由于x -> (1-t)/t在无限积分方法中缩放导致不想要的答案,gsl_integration_qagi我现在求助于在有限支持上使用积分。我现在要做的是:

我有一个(离散S的)实数系列,我将(每个)与给定的指数函数卷积

    exp(-t/T) ... T = decay constant

我选择对集成的支持,(min(S) - 10*T, max(S) + 10*T)以便涵盖该功能的大部分“重大”贡献。

使用此支持进行集成gsl_integration_qag需要几秒钟,而 gsl_integration_qagi`(无限支持)几乎不需要几毫秒,但会产生错误的结果。有人知道原因吗?

gsl_integration_qag如果卷积是高斯而不是指数,则同样有效。

在此先感谢,尼基尔

4

0 回答 0