我正在计算包含符号变量"W"的 8x8 矩阵的特征值。将 8 个特征值绘制为 W 的函数会在图中返回一个奇怪的结果,看起来像有人骑着自行车经过我的图。
对于第二个矩阵,我只是将一些非对角线元素设置为等于 0,一切正常。但我不知道第一个有什么问题。
syms W;
w0=1/780;
wl=1/1064;
h=1; % for now this seems unnecessary, but I want to change this value later on
% This is the 1st matrix which causes some strange plotting results
A=h*[w0+3*wl 2*W 0 0 0 sqrt(3)*W 0 0;
2*W 4*wl 0 0 0 0 0 0;
0 0 2*wl+w0 sqrt(3)*W 0 0 0 sqrt(2)*W;
0 0 sqrt(3)*W 3*wl 0 0 0 0;
0 0 0 0 wl+w0 sqrt(2)*W 0 0;
sqrt(3)*W 0 0 0 sqrt(2)*W 2*wl 0 0;
0 0 0 0 0 0 w0 W;
0 0 sqrt(2)*W 0 0 0 W wl];
% This is the 2nd matrix for which everything is working fine
B=h*[w0+3*wl 2*W 0 0 0 0 0 0;
2*W 4*wl 0 0 0 0 0 0;
0 0 2*wl+w0 sqrt(3)*W 0 0 0 0;
0 0 sqrt(3)*W 3*wl 0 0 0 0;
0 0 0 0 wl+w0 sqrt(2)*W 0 0;
0 0 0 0 sqrt(2)*W 2*wl 0 0;
0 0 0 0 0 0 w0 W;
0 0 0 0 0 0 W wl];
X = eig(A);
X2 = eig(B);
eva22 = X2(1);
eva1 = X(1);
figure(1);
fplot(X2,[-0.002 0.002]);
hold on;
fplot(X,[-0.002 0.002]);
hold off;
xlabel('Rabi frequency [THz]','FontSize',11);
ylabel('dressed states','FontSize',11);
grid on;
box on;
我希望矩阵 A 的图与矩阵 B 的图相似,但不知何故它无法正常工作。我将不胜感激一些如何解决此问题的提示和技巧。
第二个图看起来像这样,因为 B 的特征值是虚数。使用 plot(( 时,它默认绘制复数的实部,但显然 fplot(( 不会。您可以改为fplot(real(X), [-0.002 0.002])
绘制特征值的真实部分(假设这是您想要的(。