美文网首页
python计算Jordan标准形

python计算Jordan标准形

作者: 一路向后 | 来源:发表于2020-10-17 13:59 被阅读0次

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

相关文章

网友评论

      本文标题:python计算Jordan标准形

      本文链接:https://www.haomeiwen.com/subject/hkaxmktx.html