例如,我计算矩阵A的倒数,

import numpy as np

A = np.diag([1, 2, 3])
A_inv = np.linalg.pinv(A)
print(A_inv)

我有,

[[ 1.          0.          0.        ]
 [ 0.          0.5         0.        ]
 [ 0.          0.          0.33333333]]

但是我想要这个

[[ 1.          0.          0. ]
 [ 0.          1/2         0. ]
 [ 0.          0.          1/3]]

我尝试了np.set_printoptions

import fractions
np.set_printoptions(formatter={'all':lambda x: str(fractions.Fraction(x))})
print(A_inv)

但是我明白了

[[1 0 0]
 [0 1/2 0]
 [0 0 6004799503160661/18014398509481984]]

如何在NumPy中将小数转换为分数？

解决方法:

这是一个浮点问题-请记住,在Pythons表示中,2/3并不完全是2/3.

Fraction类有一个内置方法limit_denominator()来解决此问题：

import fractions
np.set_printoptions(formatter={'all':lambda x: str(fractions.Fraction(x).limit_denominator())})
print(A_inv)

给出所需的答案：

[[1 0 0]
 [0 1/2 0]
 [0 0 1/3]]