Julia二维无穷限积分问题求助：结果异常及变量替换后仍不正确

问题根源分析 #

你之前的代码存在两个关键问题：

1. 无穷限积分收敛性问题：直接对[-Inf,-Inf]到[Inf,Inf]的矩形区域积分时，若被积函数f_X_opt或g_V_sub在无穷远处衰减速度不足，数值积分会因无法收敛返回NaN。
2. 未处理变量依赖的积分限：你提到内层积分限依赖x，但两次实现都把积分区域当成了x和v完全独立的矩形区域，这直接导致积分范围错误，结果自然不准确。

正确实现方案：累次积分嵌套 #

对于带变量积分限的二重积分（形式为$\int_{x=a}^{b} \left( \int_{v=v_{\text{low}}(x)}^{v_{\text{high}}(x)} f(x,v) dv \right) dx$），最直接的方法是拆成累次积分：先计算依赖x的内层积分，再对x做外层积分。

示例1：用QuadGK实现（推荐，原生支持无穷限） #

QuadGK无需手动做变量替换，直接支持无穷积分限，嵌套使用即可：

using QuadGK

function h_10(K, i)
    # 外层积分：对x积分，内层是依赖x的v积分
    outer_result, _ = quadgk(x -> begin
        # 替换成你实际的x依赖型积分限
        v_low = ...  # 比如某个关于x的表达式，如x^2
        v_high = ...  # 比如Inf或另一个x相关表达式
        # 内层积分：计算v在对应区间的积分
        inner_result, _ = quadgk(v -> f_X_opt(K, i, x) * g_V_sub(K, i, v), v_low, v_high)
        return inner_result
    end, -Inf, Inf; rtol=1e-6)  # rtol可根据精度需求调整
    return 1 - outer_result
end

res = h_10(K, i)
println(res)

示例2：用Cubature实现（自定义积分区域） #

如果坚持用Cubature，可通过被积函数过滤无效区域来实现变量积分限：

using Cubature

function h_10(K, i)
    function integrand(xv)
        x, v = xv[1], xv[2]
        # 定义x对应的v积分限，不在范围内则返回0
        v_low = ...  # 你的x依赖下限
        v_high = ...  # 你的x依赖上限
        if v < v_low || v > v_high
            return 0.0
        end
        return f_X_opt(K, i, x) * g_V_sub(K, i, v)
    end
    # 积分区域设为全局，但通过被积函数过滤无效区间
    result, _ = hcubature(integrand, [-Inf,-Inf], [Inf, Inf]; reltol=1e-6)
    return 1 - result
end

res = h_10(K, i)
println(res)

关键注意事项 #

• 被积函数收敛性：确保f_X_opt(x)*g_V_sub(v)在积分区域边界（尤其是无穷远）快速衰减，否则数值积分可能无法收敛或精度极差。
• 精度参数调整：根据需求修改rtol（相对误差）或atol（绝对误差）参数，平衡计算速度和结果准确性。
• 变量替换的正确姿势：若必须手动做变量替换，要在累次积分中分别对x和v单独处理，而非直接套用矩形区域替换逻辑。

内容的提问来源于stack exchange，提问作者KBon