java - Clojure 中的分数矩阵幂

Question

我正在尝试在 Clojure 中使用转换矩阵。将年度债券评级转换矩阵转换为季度债券评级时，我需要方阵的 0.25 次幂。

在Python中，我们有fractional_matrix_power from scipy如:

>>> from scipy.linalg import fractional_matrix_power
>>> a = np.array([[1.0, 3.0], [1.0, 4.0]])
>>> b = fractional_matrix_power(a, 0.5)
>>> b
array([[ 0.75592895,  1.13389342],
       [ 0.37796447,  1.88982237]])

正在搜索Incanter和 Parallel Colt我还没有找到任何东西。费力地浏览 javadocs 和谷歌搜索并没有帮助，但也许“矩阵”和“幂”太笼统，无法深入了解我正在寻找的内容。

我确实需要 transcode a python或 R 函数，或者我缺少一些很酷的 colt 文档网站吗？

Answer 1

您可能已经看过 the core.matrix library for Clojure 。它专注于基本线性代数，并且(还!)没有您所描述的矩阵指数函数。

但是，如果您确实将上述算法从 Python 转码为 Clojure，我相信这将是 core.matrix 中受欢迎的补充。

您还应该查看 the Neanderthal project ，重点是使用 GPU 进行矩阵运算。