这个问题跟我上一个问题有关GSL integration behaves strange
由于缩放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
效果很好。
提前致谢, NIKHIL