Matlab矩阵论矩阵分析计算实现(四)求史密斯标准型和约当标准型
Matlab中有内置的史密斯标准型和约当标准型,所以不在用例题多做说明。
以下是代码
syms x;
A = [x*(x+1) 0 0; 0 x 0;0 0 (x+1)^2];
% 求矩阵A的Smith标准形 计算史密斯需要sI-A
S = simplify(smithForm(A));
display(S);
%
% 求矩阵A的约当标准型 约当型是原矩阵
B = [4 5 -2;-2 -2 1;-1 -1 1];
disp("求矩阵史密斯标准型,其对角元素是不变因子,其非常数的因子为初等因子");
S = simplify(smithForm(A));
disp(S);
[V,J] = jordan(B);
display(J);
运行结果
源码
https://download.csdn.net/download/Pedrotime/87355113