MATLAB: How to test if a matrix is unitary

conjugateinversematrixunitary

My current test for a unitary matrix is shown in the code, I'd also like to know if U = e^(iH) [i is the complex number] is coded correctly. Thanks!

 U = exp(i*H)
 Uinverse = inv(U)
 UConjTran = U'
 if UConjTran == Uinverse
     disp('U is unitary')
else
     disp('U is NOT unitary')
end

Again, e^(i*H) is not the same as exp(i*H). Check matlab's "mpower" operator.

You can find the nearest unitary matrix in each iteration using an SVD:

for t_i = 1:100
     [S,~,D] = svd(U * U);
     U       = S * D';
     norm    = trace(U * U' - eye(10))
end

Each time through the loop you replace V1 and P, so at the end of the loop you only have a single point to plot.

Try

V1 = (0:1:VCA).*2; P = V1.*IR.*35;

with no loop.