我现在会若当标准型的分解了
虽然还是不明白为什么要这么做，但是我现在已经彻底会标准型的分解了。Clear["Global`*"];(*mathematica11.2,win7(64bit)Clear all variables*)
a={
{-1,0,-1,1,1,3,0},
{0,1,0,0,0,0,0},
{2,1,2,-1,-1,-6,0},
{-2,0,-1,2,1,3,0},
{0,0,0,0,1,0,0},
{0,0,0,0,0,1,0},
{-1,-1,0,1,2,4,1}
};(*定义矩阵*)
i=IdentityMatrix;(*单位矩阵*)
RowReduce]//MatrixForm
RowReduce]//MatrixForm
RowReduce]//MatrixForm
{v11,v12,v13,v14}=NullSpace](*零空间里有四个向量,证明若当块有四个*)
v21=LinearSolve
v31=LinearSolve(*求解不了,到此结束*)(*证明此向量对应的若当块是2*2*)
v22=LinearSolve(*求解不了,到此结束*)(*证明此向量对应的若当块是1*1*)
v23=LinearSolve(*求解不了,到此结束*)(*证明此向量对应的若当块是1*1*)
v24=LinearSolve
v34=LinearSolve
v44=LinearSolve(*求解不了,到此结束*)(*证明此向量对应的若当块是3*3*)
p=Transpose@{v11,v21,v12,v13,v14,v24,v34}(*得到向量矩阵，转置成列向量*)
{s,j}=JordanDecomposition(*若当分解，过度矩阵与若当矩阵*)
p==s(*两者相等*)
算得与mathematica完全一样 【转发】矩阵的极小多项式与若尔当标准型（Mathematica代码）
https://www.bilibili.com/video/BV1f4411x7DY 葡萄糖 发表于 2021-3-23 22:18
【转发】矩阵的极小多项式与若尔当标准型（Mathematica代码）
https://www.bilibili.com/video/BV1f4411x7 ...
你写的吗？
页:
[1]