Clojure中的分数矩阵幂

Question

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

在python中，我们有来自scipy的fractional_matrix_power为：

>>> 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和谷歌搜索并没有帮助，但也许"矩阵"和"权力"太笼统了，无法深入到我正在寻找的内容。

我真的需要对 python 或 R 函数进行转码，还是缺少一些很酷的柯尔特文档站点？

Answer 1

你可能已经看过 Clojure 的 core.matrix 库。 它专注于基本的线性代数，并且没有(尚未！(您所描述的矩阵指数函数。

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

您还应该研究尼安德特人项目，该项目专注于使用 GPU 进行矩阵运算。