Problem in assigning proper eigevalues in a for-loop

5 visualizaciones (últimos 30 días)
Dear all, I encounter a problem in assigning and plotting the eigenvalues, the final graphs shown 3 separate cosine bands but the locations get some serious problems, causing many discontinuities. I have no idea how to fix it, but I think the problem goes to the fact that matlab does not know how to allocate the eigenvalues correctly for each k in the loop. Here is my code:
theta1= 0.4;
theta2= 0.2;
theta3= 0.4;
theta4= 0.5;
theta5= 0.4;
theta6= 0.3;
theta7= 0.8;
theta8= 0.5;
GellMann1 = [0 1 0; 1 0 0; 0 0 0];
GellMann2 = [0 -1i 0; 1i 0 0; 0 0 0];
GellMann3 = [1 0 0; 0 -1 0; 0 0 0];
GellMann4 = [0 0 1; 0 0 0; 1 0 0];
GellMann5 = [0 0 -1i; 0 0 0; 1i 0 0];
GellMann6 = [0 0 0; 0 0 1; 0 1 0];
GellMann7 = [0 0 0; 0 0 -1i; 0 1i 0];
GellMann8 = [1/sqrt(3) 0 0; 0 1/sqrt(3) 0; 0 0 -2/sqrt(3) ];
normtheta= sqrt(theta1^2 + theta2^2 + theta3^2 + theta4^2 + theta5^2 + theta6^2 + theta6^2 + theta8^2);
PhaseMatrix = theta1*GellMann1 + theta2*GellMann2 + theta3*GellMann3 + theta4*GellMann4 + theta5*GellMann5 + theta6*GellMann6+ theta7*GellMann7 + theta8*GellMann8 ;
J=0.8;
omega_1 = 0.08;
omega_2 = 0.1;
counter_Eigenvalues1_vector =1;
counter_Eigenvalues2_vector =1;
counter_Eigenvalues3_vector =1;
for k = -2*pi : 4*pi/500 : 2*pi
fprintf('The corresponding k is %d \n', k);
Hamiltonian= J*(expm(1i*(k*eye(3) - PhaseMatrix)) + expm(-1i*(k*eye(3) - PhaseMatrix)));
eig_Correct = eig(Hamiltonian);
e1correct = eig_Correct(1);
e2correct = eig_Correct(2);
e3correct = eig_Correct(3);
%%%Save the three eigenvalues into an array%%%%
Eigenvalues1_vector(1, counter_Eigenvalues1_vector) = (e1correct);
counter_Eigenvalues1_vector = counter_Eigenvalues1_vector +1;
Eigenvalues2_vector(1, counter_Eigenvalues2_vector) = (e2correct);
counter_Eigenvalues2_vector = counter_Eigenvalues2_vector +1;
Eigenvalues3_vector(1, counter_Eigenvalues3_vector) = (e3correct);
counter_Eigenvalues3_vector = counter_Eigenvalues3_vector +1;
end
The corresponding k is -6.283185e+00 The corresponding k is -6.258053e+00 The corresponding k is -6.232920e+00 The corresponding k is -6.207787e+00 The corresponding k is -6.182654e+00 The corresponding k is -6.157522e+00 The corresponding k is -6.132389e+00 The corresponding k is -6.107256e+00 The corresponding k is -6.082123e+00 The corresponding k is -6.056991e+00 The corresponding k is -6.031858e+00 The corresponding k is -6.006725e+00 The corresponding k is -5.981592e+00 The corresponding k is -5.956460e+00 The corresponding k is -5.931327e+00 The corresponding k is -5.906194e+00 The corresponding k is -5.881061e+00 The corresponding k is -5.855929e+00 The corresponding k is -5.830796e+00 The corresponding k is -5.805663e+00 The corresponding k is -5.780530e+00 The corresponding k is -5.755398e+00 The corresponding k is -5.730265e+00 The corresponding k is -5.705132e+00 The corresponding k is -5.680000e+00 The corresponding k is -5.654867e+00 The corresponding k is -5.629734e+00 The corresponding k is -5.604601e+00 The corresponding k is -5.579469e+00 The corresponding k is -5.554336e+00 The corresponding k is -5.529203e+00 The corresponding k is -5.504070e+00 The corresponding k is -5.478938e+00 The corresponding k is -5.453805e+00 The corresponding k is -5.428672e+00 The corresponding k is -5.403539e+00 The corresponding k is -5.378407e+00 The corresponding k is -5.353274e+00 The corresponding k is -5.328141e+00 The corresponding k is -5.303008e+00 The corresponding k is -5.277876e+00 The corresponding k is -5.252743e+00 The corresponding k is -5.227610e+00 The corresponding k is -5.202477e+00 The corresponding k is -5.177345e+00 The corresponding k is -5.152212e+00 The corresponding k is -5.127079e+00 The corresponding k is -5.101946e+00 The corresponding k is -5.076814e+00 The corresponding k is -5.051681e+00 The corresponding k is -5.026548e+00 The corresponding k is -5.001416e+00 The corresponding k is -4.976283e+00 The corresponding k is -4.951150e+00 The corresponding k is -4.926017e+00 The corresponding k is -4.900885e+00 The corresponding k is -4.875752e+00 The corresponding k is -4.850619e+00 The corresponding k is -4.825486e+00 The corresponding k is -4.800354e+00 The corresponding k is -4.775221e+00 The corresponding k is -4.750088e+00 The corresponding k is -4.724955e+00 The corresponding k is -4.699823e+00 The corresponding k is -4.674690e+00 The corresponding k is -4.649557e+00 The corresponding k is -4.624424e+00 The corresponding k is -4.599292e+00 The corresponding k is -4.574159e+00 The corresponding k is -4.549026e+00 The corresponding k is -4.523893e+00 The corresponding k is -4.498761e+00 The corresponding k is -4.473628e+00 The corresponding k is -4.448495e+00 The corresponding k is -4.423362e+00 The corresponding k is -4.398230e+00 The corresponding k is -4.373097e+00 The corresponding k is -4.347964e+00 The corresponding k is -4.322831e+00 The corresponding k is -4.297699e+00 The corresponding k is -4.272566e+00 The corresponding k is -4.247433e+00 The corresponding k is -4.222301e+00 The corresponding k is -4.197168e+00 The corresponding k is -4.172035e+00 The corresponding k is -4.146902e+00 The corresponding k is -4.121770e+00 The corresponding k is -4.096637e+00 The corresponding k is -4.071504e+00 The corresponding k is -4.046371e+00 The corresponding k is -4.021239e+00 The corresponding k is -3.996106e+00 The corresponding k is -3.970973e+00 The corresponding k is -3.945840e+00 The corresponding k is -3.920708e+00 The corresponding k is -3.895575e+00 The corresponding k is -3.870442e+00 The corresponding k is -3.845309e+00 The corresponding k is -3.820177e+00 The corresponding k is -3.795044e+00 The corresponding k is -3.769911e+00 The corresponding k is -3.744778e+00 The corresponding k is -3.719646e+00 The corresponding k is -3.694513e+00 The corresponding k is -3.669380e+00 The corresponding k is -3.644247e+00 The corresponding k is -3.619115e+00 The corresponding k is -3.593982e+00 The corresponding k is -3.568849e+00 The corresponding k is -3.543717e+00 The corresponding k is -3.518584e+00 The corresponding k is -3.493451e+00 The corresponding k is -3.468318e+00 The corresponding k is -3.443186e+00 The corresponding k is -3.418053e+00 The corresponding k is -3.392920e+00 The corresponding k is -3.367787e+00 The corresponding k is -3.342655e+00 The corresponding k is -3.317522e+00 The corresponding k is -3.292389e+00 The corresponding k is -3.267256e+00 The corresponding k is -3.242124e+00 The corresponding k is -3.216991e+00 The corresponding k is -3.191858e+00 The corresponding k is -3.166725e+00 The corresponding k is -3.141593e+00 The corresponding k is -3.116460e+00 The corresponding k is -3.091327e+00 The corresponding k is -3.066194e+00 The corresponding k is -3.041062e+00 The corresponding k is -3.015929e+00 The corresponding k is -2.990796e+00 The corresponding k is -2.965663e+00 The corresponding k is -2.940531e+00 The corresponding k is -2.915398e+00 The corresponding k is -2.890265e+00 The corresponding k is -2.865133e+00 The corresponding k is -2.840000e+00 The corresponding k is -2.814867e+00 The corresponding k is -2.789734e+00 The corresponding k is -2.764602e+00 The corresponding k is -2.739469e+00 The corresponding k is -2.714336e+00 The corresponding k is -2.689203e+00 The corresponding k is -2.664071e+00 The corresponding k is -2.638938e+00 The corresponding k is -2.613805e+00 The corresponding k is -2.588672e+00 The corresponding k is -2.563540e+00 The corresponding k is -2.538407e+00 The corresponding k is -2.513274e+00 The corresponding k is -2.488141e+00 The corresponding k is -2.463009e+00 The corresponding k is -2.437876e+00 The corresponding k is -2.412743e+00 The corresponding k is -2.387610e+00 The corresponding k is -2.362478e+00 The corresponding k is -2.337345e+00 The corresponding k is -2.312212e+00 The corresponding k is -2.287079e+00 The corresponding k is -2.261947e+00 The corresponding k is -2.236814e+00 The corresponding k is -2.211681e+00 The corresponding k is -2.186548e+00 The corresponding k is -2.161416e+00 The corresponding k is -2.136283e+00 The corresponding k is -2.111150e+00 The corresponding k is -2.086018e+00 The corresponding k is -2.060885e+00 The corresponding k is -2.035752e+00 The corresponding k is -2.010619e+00 The corresponding k is -1.985487e+00 The corresponding k is -1.960354e+00 The corresponding k is -1.935221e+00 The corresponding k is -1.910088e+00 The corresponding k is -1.884956e+00 The corresponding k is -1.859823e+00 The corresponding k is -1.834690e+00 The corresponding k is -1.809557e+00 The corresponding k is -1.784425e+00 The corresponding k is -1.759292e+00 The corresponding k is -1.734159e+00 The corresponding k is -1.709026e+00 The corresponding k is -1.683894e+00 The corresponding k is -1.658761e+00 The corresponding k is -1.633628e+00 The corresponding k is -1.608495e+00 The corresponding k is -1.583363e+00 The corresponding k is -1.558230e+00 The corresponding k is -1.533097e+00 The corresponding k is -1.507964e+00 The corresponding k is -1.482832e+00 The corresponding k is -1.457699e+00 The corresponding k is -1.432566e+00 The corresponding k is -1.407434e+00 The corresponding k is -1.382301e+00 The corresponding k is -1.357168e+00 The corresponding k is -1.332035e+00 The corresponding k is -1.306903e+00 The corresponding k is -1.281770e+00 The corresponding k is -1.256637e+00 The corresponding k is -1.231504e+00 The corresponding k is -1.206372e+00 The corresponding k is -1.181239e+00 The corresponding k is -1.156106e+00 The corresponding k is -1.130973e+00 The corresponding k is -1.105841e+00 The corresponding k is -1.080708e+00 The corresponding k is -1.055575e+00 The corresponding k is -1.030442e+00 The corresponding k is -1.005310e+00 The corresponding k is -9.801769e-01 The corresponding k is -9.550442e-01 The corresponding k is -9.299114e-01 The corresponding k is -9.047787e-01 The corresponding k is -8.796459e-01 The corresponding k is -8.545132e-01 The corresponding k is -8.293805e-01 The corresponding k is -8.042477e-01 The corresponding k is -7.791150e-01 The corresponding k is -7.539822e-01 The corresponding k is -7.288495e-01 The corresponding k is -7.037168e-01 The corresponding k is -6.785840e-01 The corresponding k is -6.534513e-01 The corresponding k is -6.283185e-01 The corresponding k is -6.031858e-01 The corresponding k is -5.780530e-01 The corresponding k is -5.529203e-01 The corresponding k is -5.277876e-01 The corresponding k is -5.026548e-01 The corresponding k is -4.775221e-01 The corresponding k is -4.523893e-01 The corresponding k is -4.272566e-01 The corresponding k is -4.021239e-01 The corresponding k is -3.769911e-01 The corresponding k is -3.518584e-01 The corresponding k is -3.267256e-01 The corresponding k is -3.015929e-01 The corresponding k is -2.764602e-01 The corresponding k is -2.513274e-01 The corresponding k is -2.261947e-01 The corresponding k is -2.010619e-01 The corresponding k is -1.759292e-01 The corresponding k is -1.507964e-01 The corresponding k is -1.256637e-01 The corresponding k is -1.005310e-01 The corresponding k is -7.539822e-02 The corresponding k is -5.026548e-02 The corresponding k is -2.513274e-02 The corresponding k is 0 The corresponding k is 2.513274e-02 The corresponding k is 5.026548e-02 The corresponding k is 7.539822e-02 The corresponding k is 1.005310e-01 The corresponding k is 1.256637e-01 The corresponding k is 1.507964e-01 The corresponding k is 1.759292e-01 The corresponding k is 2.010619e-01 The corresponding k is 2.261947e-01 The corresponding k is 2.513274e-01 The corresponding k is 2.764602e-01 The corresponding k is 3.015929e-01 The corresponding k is 3.267256e-01 The corresponding k is 3.518584e-01 The corresponding k is 3.769911e-01 The corresponding k is 4.021239e-01 The corresponding k is 4.272566e-01 The corresponding k is 4.523893e-01 The corresponding k is 4.775221e-01 The corresponding k is 5.026548e-01 The corresponding k is 5.277876e-01 The corresponding k is 5.529203e-01 The corresponding k is 5.780530e-01 The corresponding k is 6.031858e-01 The corresponding k is 6.283185e-01 The corresponding k is 6.534513e-01 The corresponding k is 6.785840e-01 The corresponding k is 7.037168e-01 The corresponding k is 7.288495e-01 The corresponding k is 7.539822e-01 The corresponding k is 7.791150e-01 The corresponding k is 8.042477e-01 The corresponding k is 8.293805e-01 The corresponding k is 8.545132e-01 The corresponding k is 8.796459e-01 The corresponding k is 9.047787e-01 The corresponding k is 9.299114e-01 The corresponding k is 9.550442e-01 The corresponding k is 9.801769e-01 The corresponding k is 1.005310e+00 The corresponding k is 1.030442e+00 The corresponding k is 1.055575e+00 The corresponding k is 1.080708e+00 The corresponding k is 1.105841e+00 The corresponding k is 1.130973e+00 The corresponding k is 1.156106e+00 The corresponding k is 1.181239e+00 The corresponding k is 1.206372e+00 The corresponding k is 1.231504e+00 The corresponding k is 1.256637e+00 The corresponding k is 1.281770e+00 The corresponding k is 1.306903e+00 The corresponding k is 1.332035e+00 The corresponding k is 1.357168e+00 The corresponding k is 1.382301e+00 The corresponding k is 1.407434e+00 The corresponding k is 1.432566e+00 The corresponding k is 1.457699e+00 The corresponding k is 1.482832e+00 The corresponding k is 1.507964e+00 The corresponding k is 1.533097e+00 The corresponding k is 1.558230e+00 The corresponding k is 1.583363e+00 The corresponding k is 1.608495e+00 The corresponding k is 1.633628e+00 The corresponding k is 1.658761e+00 The corresponding k is 1.683894e+00 The corresponding k is 1.709026e+00 The corresponding k is 1.734159e+00 The corresponding k is 1.759292e+00 The corresponding k is 1.784425e+00 The corresponding k is 1.809557e+00 The corresponding k is 1.834690e+00 The corresponding k is 1.859823e+00 The corresponding k is 1.884956e+00 The corresponding k is 1.910088e+00 The corresponding k is 1.935221e+00 The corresponding k is 1.960354e+00 The corresponding k is 1.985487e+00 The corresponding k is 2.010619e+00 The corresponding k is 2.035752e+00 The corresponding k is 2.060885e+00 The corresponding k is 2.086018e+00 The corresponding k is 2.111150e+00 The corresponding k is 2.136283e+00 The corresponding k is 2.161416e+00 The corresponding k is 2.186548e+00 The corresponding k is 2.211681e+00 The corresponding k is 2.236814e+00 The corresponding k is 2.261947e+00 The corresponding k is 2.287079e+00 The corresponding k is 2.312212e+00 The corresponding k is 2.337345e+00 The corresponding k is 2.362478e+00 The corresponding k is 2.387610e+00 The corresponding k is 2.412743e+00 The corresponding k is 2.437876e+00 The corresponding k is 2.463009e+00 The corresponding k is 2.488141e+00 The corresponding k is 2.513274e+00 The corresponding k is 2.538407e+00 The corresponding k is 2.563540e+00 The corresponding k is 2.588672e+00 The corresponding k is 2.613805e+00 The corresponding k is 2.638938e+00 The corresponding k is 2.664071e+00 The corresponding k is 2.689203e+00 The corresponding k is 2.714336e+00 The corresponding k is 2.739469e+00 The corresponding k is 2.764602e+00 The corresponding k is 2.789734e+00 The corresponding k is 2.814867e+00 The corresponding k is 2.840000e+00 The corresponding k is 2.865133e+00 The corresponding k is 2.890265e+00 The corresponding k is 2.915398e+00 The corresponding k is 2.940531e+00 The corresponding k is 2.965663e+00 The corresponding k is 2.990796e+00 The corresponding k is 3.015929e+00 The corresponding k is 3.041062e+00 The corresponding k is 3.066194e+00 The corresponding k is 3.091327e+00 The corresponding k is 3.116460e+00 The corresponding k is 3.141593e+00 The corresponding k is 3.166725e+00 The corresponding k is 3.191858e+00 The corresponding k is 3.216991e+00 The corresponding k is 3.242124e+00 The corresponding k is 3.267256e+00 The corresponding k is 3.292389e+00 The corresponding k is 3.317522e+00 The corresponding k is 3.342655e+00 The corresponding k is 3.367787e+00 The corresponding k is 3.392920e+00 The corresponding k is 3.418053e+00 The corresponding k is 3.443186e+00 The corresponding k is 3.468318e+00 The corresponding k is 3.493451e+00 The corresponding k is 3.518584e+00 The corresponding k is 3.543717e+00 The corresponding k is 3.568849e+00 The corresponding k is 3.593982e+00 The corresponding k is 3.619115e+00 The corresponding k is 3.644247e+00 The corresponding k is 3.669380e+00 The corresponding k is 3.694513e+00 The corresponding k is 3.719646e+00 The corresponding k is 3.744778e+00 The corresponding k is 3.769911e+00 The corresponding k is 3.795044e+00 The corresponding k is 3.820177e+00 The corresponding k is 3.845309e+00 The corresponding k is 3.870442e+00 The corresponding k is 3.895575e+00 The corresponding k is 3.920708e+00 The corresponding k is 3.945840e+00 The corresponding k is 3.970973e+00 The corresponding k is 3.996106e+00 The corresponding k is 4.021239e+00 The corresponding k is 4.046371e+00 The corresponding k is 4.071504e+00 The corresponding k is 4.096637e+00 The corresponding k is 4.121770e+00 The corresponding k is 4.146902e+00 The corresponding k is 4.172035e+00 The corresponding k is 4.197168e+00 The corresponding k is 4.222301e+00 The corresponding k is 4.247433e+00 The corresponding k is 4.272566e+00 The corresponding k is 4.297699e+00 The corresponding k is 4.322831e+00 The corresponding k is 4.347964e+00 The corresponding k is 4.373097e+00 The corresponding k is 4.398230e+00 The corresponding k is 4.423362e+00 The corresponding k is 4.448495e+00 The corresponding k is 4.473628e+00 The corresponding k is 4.498761e+00 The corresponding k is 4.523893e+00 The corresponding k is 4.549026e+00 The corresponding k is 4.574159e+00 The corresponding k is 4.599292e+00 The corresponding k is 4.624424e+00 The corresponding k is 4.649557e+00 The corresponding k is 4.674690e+00 The corresponding k is 4.699823e+00 The corresponding k is 4.724955e+00 The corresponding k is 4.750088e+00 The corresponding k is 4.775221e+00 The corresponding k is 4.800354e+00 The corresponding k is 4.825486e+00 The corresponding k is 4.850619e+00 The corresponding k is 4.875752e+00 The corresponding k is 4.900885e+00 The corresponding k is 4.926017e+00 The corresponding k is 4.951150e+00 The corresponding k is 4.976283e+00 The corresponding k is 5.001416e+00 The corresponding k is 5.026548e+00 The corresponding k is 5.051681e+00 The corresponding k is 5.076814e+00 The corresponding k is 5.101946e+00 The corresponding k is 5.127079e+00 The corresponding k is 5.152212e+00 The corresponding k is 5.177345e+00 The corresponding k is 5.202477e+00 The corresponding k is 5.227610e+00 The corresponding k is 5.252743e+00 The corresponding k is 5.277876e+00 The corresponding k is 5.303008e+00 The corresponding k is 5.328141e+00 The corresponding k is 5.353274e+00 The corresponding k is 5.378407e+00 The corresponding k is 5.403539e+00 The corresponding k is 5.428672e+00 The corresponding k is 5.453805e+00 The corresponding k is 5.478938e+00 The corresponding k is 5.504070e+00 The corresponding k is 5.529203e+00 The corresponding k is 5.554336e+00 The corresponding k is 5.579469e+00 The corresponding k is 5.604601e+00 The corresponding k is 5.629734e+00 The corresponding k is 5.654867e+00 The corresponding k is 5.680000e+00 The corresponding k is 5.705132e+00 The corresponding k is 5.730265e+00 The corresponding k is 5.755398e+00 The corresponding k is 5.780530e+00 The corresponding k is 5.805663e+00 The corresponding k is 5.830796e+00 The corresponding k is 5.855929e+00 The corresponding k is 5.881061e+00 The corresponding k is 5.906194e+00 The corresponding k is 5.931327e+00 The corresponding k is 5.956460e+00 The corresponding k is 5.981592e+00 The corresponding k is 6.006725e+00 The corresponding k is 6.031858e+00 The corresponding k is 6.056991e+00 The corresponding k is 6.082123e+00 The corresponding k is 6.107256e+00 The corresponding k is 6.132389e+00 The corresponding k is 6.157522e+00 The corresponding k is 6.182654e+00 The corresponding k is 6.207787e+00 The corresponding k is 6.232920e+00 The corresponding k is 6.258053e+00 The corresponding k is 6.283185e+00
%%%Plotting
figure(1)
k = -2*pi : 4*pi/500 : 2*pi
k = 1×501
-6.2832 -6.2581 -6.2329 -6.2078 -6.1827 -6.1575 -6.1324 -6.1073 -6.0821 -6.0570 -6.0319 -6.0067 -5.9816 -5.9565 -5.9313 -5.9062 -5.8811 -5.8559 -5.8308 -5.8057 -5.7805 -5.7554 -5.7303 -5.7051 -5.6800 -5.6549 -5.6297 -5.6046 -5.5795 -5.5543
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
plot(k, Eigenvalues1_vector,k, Eigenvalues2_vector,k,Eigenvalues3_vector)
Warning: Imaginary parts of complex X and/or Y arguments ignored.
axis([-2*pi 2*pi -3 3])
xlabel('k')
ylabel('E_k')

Respuesta aceptada

Christine Tobler
Christine Tobler el 8 de En. de 2025
Sorting the returned eigenvalues by their real part (since you are plotting only the real part of the eigenvalues) should help. You will avoid the discontinuities, but will have a line in the plot that's always the largest real part of the eigenvalues, another for the smallest, and a third for the middle.
If you're looking for smooth tracing of the eigenvalues, you could try https://www.mathworks.com/matlabcentral/fileexchange/22885-eigenshuffle from the File Exchange, which shuffles the eigenvalues in a sequence to form smooth curves.

Más respuestas (1)

Sameer
Sameer el 12 de Dic. de 2024
Editada: Sameer el 12 de Dic. de 2024
It looks like the issue you are encountering is due to the discontinuities in the eigenvalues, which can happen if the eigenvalues are not sorted consistently for each value of "k". This can cause the plot to jump between different eigenvalues, leading to the discontinuities you are seeing.
To fix this, you can sort the eigenvalues at each step of the loop to ensure they are consistently ordered.
Hope this helps!

Categorías

Más información sobre Linear Algebra en Help Center y File Exchange.

Productos


Versión

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by