如何利用Julia的IntervalArithmetic.jl包求解向量函数的雅可比矩阵？

步骤1：安装并导入包 #

首先在Julia REPL中安装IntervalArithmetic.jl包：

using Pkg
Pkg.add("IntervalArithmetic")

导入包以使用其功能：

using IntervalArithmetic

步骤2：定义向量函数 #

编写接受区间向量输入的目标向量函数，示例如下（可替换为你实际需要计算的函数）：

function f(x::Vector{Interval})
    x1, x2, x3 = x
    return [
        x1^2 + x2*x3,  # f1(x,y,z)
        sin(x1) + x2^3, # f2(x,y,z)
        exp(x3) - x1*x2 # f3(x,y,z)
    ]
end

步骤3：实现雅可比矩阵计算逻辑 #

通过固定其他变量、对单个变量求区间导数的方式，构建雅可比矩阵：

function jacobian(f::Function, x::Vector{Interval})
    n_vars = length(x)
    n_outputs = length(f(x))
    J = Matrix{Interval}(undef, n_outputs, n_vars)
    
    for i in 1:n_outputs
        for j in 1:n_vars
            # 构造单变量函数：仅让第j个变量变化，其余固定为输入区间
            single_var_func(t) = f([k == j ? t : x[k] for k in 1:n_vars])[i]
            # 计算该单变量函数在对应区间上的导数区间
            J[i,j] = derivative(single_var_func, x[j])
        end
    end
    
    return J
end

步骤4：测试计算流程 #

定义输入区间向量，调用jacobian函数得到区间雅可比矩阵：

# 示例输入区间：x∈[1,2], y∈[0,1], z∈[-1,0]
input_intervals = [Interval(1,2), Interval(0,1), Interval(-1,0)]
result_jacobian = jacobian(f, input_intervals)

# 打印结果
println("区间雅可比矩阵：")
display(result_jacobian)

说明 #

上述方法生成的是区间雅可比矩阵，每个元素对应偏导数在输入区间内所有可能值的包络区间，能严格保证输入区间内任意点的偏导数都落在对应区间中，适合不确定性分析、数值验证等场景。

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