Problem in assigning proper eigevalues in a for-loop

5 Ansichten (letzte 30 Tage)
Tsz Tsun
Tsz Tsun am 12 Dez. 2024
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')

Akzeptierte Antwort

Christine Tobler
Christine Tobler am 8 Jan. 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.

Weitere Antworten (1)

Sameer
Sameer am 12 Dez. 2024
Bearbeitet: Sameer am 12 Dez. 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!

Kategorien

Mehr zu Linear Algebra finden Sie in Help Center und File Exchange

Produkte


Version

R2023a

Community Treasure Hunt

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

Start Hunting!

Translated by