1.程序源码
import numpy as np
from sympy import Matrix
import sympy
import pprint
A = np.array([[3,1,0,0],[-4,-1,0,0],[6,2,0,-1],[-2,0,1,2]])
a = Matrix(A)
P, Ja = a.jordan_form()
pprint.pprint(Ja)
pprint.pprint(P)
2.计算结果
$ python3 example.py
Matrix([
[1, 1, 0, 0],
[0, 1, 0, 0],
[0, 0, 1, 1],
[0, 0, 0, 1]])
Matrix([
[ 2, 1, 1, 0],
[-4, 0, -2, 1],
[ 6, 0, 2, 0],
[-2, 0, 0, 0]])
网友评论