matlab求约旦标准型的代码

求解约旦标准型（Jordan Canonical Form）的 MATLAB 代码如下：

% 输入矩阵 A
A = [2 1 0; 0 2 0; 0 0 3];

% 求特征值和特征向量
[V, D] = eig(A);

% 初始化 Jordan 标准型矩阵 J 和变换矩阵 P
J = zeros(size(A));
P = zeros(size(A));

% 遍历每个特征值
for i = 1:size(A, 1)
    % 计算特征值的重数
    m = sum(D(i,i) == diag(D));
    
    % 初始化 Jordan 方块
    J_block = zeros(m, m);
    
    % 计算 Jordan 方块
    for j = 1:m
        J_block(j, j) = D(i,i);
        if j < m
            J_block(j, j+1) = 1;
        end
    end
    
    % 将 Jordan 方块插入到 Jordan 标准型矩阵中
    J(sum(m(1:i-1))+1:sum(m(1:i)), sum(m(1:i-1))+1:sum(m(1:i))) = J_block;
    
    % 计算变换矩阵
    P_block = zeros(size(A));
    for j = 1:m
        P_block(:, j) = V(:, i);
    end
    
    % 将变换矩阵插入到 P 中
    P(:, sum(m(1:i-1))+1:sum(m(1:i))) = P_block;
end

% 输出结果
disp('Jordan Canonical Form:');
disp(J);
disp('Transformation Matrix:');
disp(P);

其中，输入矩阵 `A` 可以根据实际情况自行修改。输出结果包括 Jordan 标准型矩阵和变换矩阵。