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

根据这些信息，您可以在the Julia repository中找到它。