对于你给的矩阵, 只要

[P,J]=jordan(A);

就行了.

如果你已经通过某种途径得到了J, 只需要求P, 那么可以

n = size(A, 1);

I = eye(n);

X = null(kron(I, A)-kron(J.', I), 'r');

P = reshape(X*randi(10, size(X, 2), 1), n, n);

另外, 如果Jordan标准型J未知的话尽量小心, matlab不保证能算出真实的J, 因为这个问题对舍入误差太敏感.