如何修复第一个矩阵(matlab)的符号特征值图?

问题描述 投票:0回答:1

我正在计算包​​含符号变量“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的图表类似,但不知何故它不能正常工作。我很感激一些提示和技巧如何解决这个问题。

matlab plot linear-algebra symbolic-math
1个回答
0
投票

第二个图看起来像是因为B的特征值是虚构的。当使用plot()时,默认情况下会绘制复数的实部,但显然fplot()不会。你可以用fplot(real(X), [-0.002 0.002])来绘制特征值的真实部分(假设你想要的那样)。

© www.soinside.com 2019 - 2024. All rights reserved.