python计算Jordan标准形

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]])