Julia中矩阵指数的函数或算子

Question

我正在寻找一个函数或运算符来计算Julia中的第n个矩阵指数,其中n是一个非负整数.例如,在MATLAB我会写expm(A, 3)的A * A * A.我可以*多次调用该产品,但有没有办法计算矩阵指数而不这样做？

julia> A = [1 -1; 0 2]
2x2 Array{Int64,2}:
 1  -1
 0   2

julia> A * A
2x2 Array{Int64,2}:
 1  -3
 0   4

julia> A * A * A
2x2 Array{Int64,2}:
 1  -7
 0   8

Answer 1

只需使用取幂运算符即可.

julia> A = [1 -1; 0 2]
2x2 Array{Int64,2}:
 1  -1
 0   2

julia> A^2
2x2 Array{Int64,2}:
 1  -3
 0   4

julia> A^3
2x2 Array{Int64,2}:
 1  -7
 0   8

要找出它的实现位置,可以使用函数方法

julia> methods(^)
# 45 methods for generic function "^":
...

^(A::Array{T,2}, p::Integer) at linalg/dense.jl:170
^(x, p::Integer) at intfuncs.jl:108
^(A::Array{T,2}, p::Number) at linalg/dense.jl:173

有了这些信息,您可以在Julia存储库中找到它.