我正在尝试使用 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 和 google 搜索并没有帮助,但也许“矩阵”和“电源”太通用了,无法深入了解我正在寻找的内容。
我真的需要对python或 R 函数进行转码,还是我缺少一些很酷的 colt doc 站点?