[P, D] = eig(Z)
(Case-1)
Z = [1, j; j, 1]
Z : Symmetric=1
Z : Hermitian=0
Z : Skew-Hermitian=0
Z : Orthogonal=0
Z : Unitary=0
Z : Normal=1
P = [0.7071,0.7071; 0.7071,-0.7071]
P : Orthogonal=1
P : Unitary=1
rank(P)=2
det(P)=-1+j*-5.88785e-017
cond(P)=1
D = [1.0000+1.0000*j,0; 0,1.0000-1.0000*j]
(Case-2)
Z = [1, j; j, 2]
Z : Symmetric=1
Z : Hermitian=0
Z : Skew-Hermitian=0
Z : Orthogonal=0
Z : Unitary=0
Z : Normal=0
P = [0.7071,-0.6124+0.3536*j; 0.6124-0.3536*j,0.7071]
P : Orthogonal=0
P : Unitary=0
rank(P)=2
det(P)=0.75+j*-0.433013
cond(P)=1.73205
D = [1.5000+0.8660*j,0; 0,1.5000-0.8660*j]
(Case-3)
Z = [1, j; j, 3]
Z : Symmetric=1
Z : Hermitian=0
Z : Skew-Hermitian=0
Z : Orthogonal=0
Z : Unitary=0
Z : Normal=0
P = [0.7071,0.7071*j; -0.7071*j,0.7071]
P : Orthogonal=0
P : Unitary=0
rank(P)=1
det(P)=2.10734e-008+j*0
cond(P)=9.49063e+007
D = [2.0000,0; 0,2.0000]
Eigenvalues are degenerated and Eigenvectors are not linearly independent in Case-3.