求jordan标准型matlab代码

以下是求Jordan标准型的MATLAB代码：

function [J, V] = jordan(A)
[V, D] = eig(A);
J = zeros(size(A));
tol = 1e-10;
for i = 1:size(A,1)
if abs(D(i,i) - round(D(i,i))) < tol
J(i,i) = round(D(i,i));
else
J(i,i) = D(i,i);
end
for j = i+1:size(A,1)
if abs(D(j,j) - J(i,i)) < tol
J(i,j) = 1;
elseif abs(D(j,j) - conj(J(i,i))) < tol
J(i,j) = 1i;
end
end
end
end

其中，A是要求Jordan标准型的方阵，J是Jordan标准型，V是相应的特征向量矩阵。tol参数为一个非常小的误差容限，用于处理复数相等的情况。