SymPy可以计算任意精度:
from sympy import exp, N, Sfrom sympy.matrices import Matrixdata = [[S("-800.21"),S("-600.00")],[S("-600.00"),S("-1000.48")]]m = Matrix(data)ex = m.applyfunc(exp).applyfunc(lambda x:N(x, 100))vecs = ex.eigenvects()print vecs[0][0] # eigen valueprint vecs[1][0] # eigen valueprint vecs[0][2] # eigen vectprint vecs[1][2] # eigen vect输出:
-2.650396553004310816338679447269582701529092549943247237903254759946483528035516341807463648841185335e-2612.650396553004310816338679447269582701529092549943247237903254759946483528035516341807466621962539464e-261[[-0.9999999999999999999999999999999999999999999999999999999999999999999999999999999999999994391176386872][ 1]][[1.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000560882361313][ 1]]
您可以将N(x,100)中的100更改为其他精度,但是,当我尝试1000时,本征vect的计算失败。



