将矩阵的数字快速提高到序列的幂?

Raising a matrix's digits to the power of a sequence fast?

提问人:matehorvath 提问时间:11/14/2023 最后编辑:user20650matehorvath 更新时间:11/14/2023 访问量:37

问:

有一个矩阵叫做 A,它是一个 n x k 矩阵。它的数字形式很复杂。我正在尝试将矩阵 A 的每一行提升到 rownum(th) 指数。

n <- 1000
k <- 10
x <- complex(real = runif(n*k), imaginary = runif(n*k))
A <- matrix(x, ncol = k, nrow = n)

Exponent_Matrix <- matrix(1:n, ncol = k, nrow = n)

Result <- A^Exponent_Matrix

library(microbenchmark)

microbenchmark(
  Result <- A^Exponent_Matrix,
  times = 1000
)
#0.000487secs

这是一个有效的解决方案,但它可以以更有效的方式完成吗?此外,矩阵 (n) 的长度可能会有所不同,但列数 (k) 通常在 1 到 20 之间。

R 矩阵 指数

评论


答: 暂无答案