因为
（JY12）
是400x400和
欧米茄
是200x200
矩阵维度必须同意。
1视图（最近30天）
显示较旧的评论
我有以下代码的Contour图，我面对的“矩阵维度必须同意”。
我该如何解决？
a = 1.855;
EB = -3.276;
EP = -1.979;
T0 = -1.844;
theta = 0;
EE = linspace（-0.3,0.3,200）;
EE = linspace（0,10,200）;
[E，Omega] = meshgrid（ee，ee）;
Sigma = 0.265;
exx = e。*（（cos（theta））。^2-sigma。*（sin（theta））。
exy = e。*（（1+sigma）。*cos（theta）。*sin（theta））;
eyx = e。*（（1+sigma）。*cos（theta）。*sin（theta））;
eyy = e。*（（sin（theta））^2-sigma。*（cos（theta））。
delta1x = a/2;
delta1y = a*sqrt（3）/2;
delta2x = a/2;
delta2y = -a*sqrt（3）/2;
delta3x = -a;
delta3y = 0;
delta11 = a*（1+3.*eyy/4+sqrt（3）。*exy/2+exx/4）;
delta22 = a。*（1+3.*eyy/4-sqrt（3）。*exy/2+exx/4）;
delta33 = a。*（1+exx）;
G1 = -2.6275;
g2 = 2*g1;
del1 = g1。*（exx+eyy）;
del2 = g2。*（exx+eyy）;
on1 = ep+del1;
on2 = eb+del2;
t1 = t0*exp（-3.37。*（delta11/a-1））;
t2 = t0*exp（-3.37。*（delta22/a-1））;
t3 = t0*exp（-3.37。*（delta33/a-1））;
delta2 = 0.2;
ETA = 0.05;
iform =复合物（0.0,1.0）;
x = linspace（-1.5,1.5,200）;
ky = linspace（-2,2,200）;
[x，y] = meshgrid（x，ky）;
U = 0.0;
t = 10;
phi = t1。*exp（iform。*x。*delta1x）。*exp（iform。*y。*delta1y）+t2。*delta2y）+t3。*exp（iform。*x。*delta3x）;
h2 = [on1。*一个（size（phi））phi;conj（phi）on2。*一个（size（phi））];
[v，e] = eig（h2）;
E1 = E（1）;
e2 = e（4）;
v1 = v（：，1）;
v2 = v（：，2）;
phiy = t1。*iform。*delta1y。*exp（iform。*x。*delta1x）。*exp（iform。*y。*delta1y）+t2。*iform。*delta2x）。*exp（iform。*y。*delta2y）;
h2y = [zeros（size（phiy））phiy;
conj（phiy）零（大小（phiy））];
jy12 =（（v1'。*h2y。*v2）。*conj（v1'。*h2y。*v2））;
deltae12 = e2-e1;
fermi1 = 1/（1+exp（（e1-u）/t））;
fermi2 = 1/（1+exp（（e2-u）/t））;
sigma_xx =（iform。*omega。^ - 1）。*（jy12）。*（fermi2-fermi1）。delta2）。^-1）;
轮廓（E，EE，Sigma_xx，500）
