美文网首页
Jordan标准形与sympy

Jordan标准形与sympy

作者: doublell | 来源:发表于2018-09-26 22:46 被阅读0次

使用python的符号计算库sympy来计算矩阵的Jordan标准形

import numpy as np
from sympy import Matrix
import sympy
import pprint

1.求A = \left[\begin{matrix} 3 & 1 & 0 & 0 \\ -4 & -1 & 0 &0 \\ 6 & 2 & 0 & -1 \\ -2&0&1&2\\ \end{matrix}\right]Jordan标准形J_A

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)
Matrix([
[1, 1, 0, 0],
[0, 1, 0, 0],
[0, 0, 1, 1],
[0, 0, 0, 1]])

1.求B = \left[\begin{matrix} -1&4&8\\ 1&3&8\\ 0&0&-1 \end{matrix}\right]Jordan标准形J_B

B = np.array([[-1,-4,-8],[1,3,8],[0,0,-1]])
b = Matrix(B)
P,Jb = b.jordan_form()
pprint.pprint(Jb)
Matrix([
[-1, 0, 0],
[ 0, 1, 1],
[ 0, 0, 1]])

相关文章

网友评论

      本文标题:Jordan标准形与sympy

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