有不同的方法可以做到这一点。
这种方法简单地将它们连接为a的第三个维度
(4 x4xnumel (B))
数组,
d0 = 0.2;%meV δ 0
d1 = 0.18;% #δ1
d2 = 0.05;% #δ2
%#在z和x方向上的电子和空穴有效g因子。
ghx = -0.35;
ghz = -2.2;
gex = -0.65;
gez = -0.8;
α= 0.02;%#meV/T^2抗磁位移系数，
ga = 0.0;
gb = ga;%#与相同极化激子态的相互作用强度。
θ=π/ 3;
wx = 0.01;
wd = wx-d0;
佤邦= wx + 0.5;
世行= wa;
王= wa;
G0a = 0.05;%#是脉冲能量的最大值
G0b = G0a;
tc = 30;%#脉冲中心的时间
tal = 10;%#是高斯均方根宽度
mub = 0.0579;玻尔磁子
B = 0:0.1:10;
V = NaN(4,4,numel(B));% Preallocate
D = v;% Preallocate
为k = 1:数字(B)
英国石油(bp) = mub * B (k)。*罪(θ)* (gez + ghz) * 0.5;%英国石油公司
bm = mub * B (k)。* sin(θ)* (gez-ghz) * 0.5;% bm
= mub * B (k)。* cos(θ)* gex * 0.5;%是
bh = mub * B (k)。* cos(θ)* ghx * 0.5;%黑洞
H1 = [wx d1/2 0,0; d1/2, 0, 0, 0, 0, wd, d2/2; 0, 0, d2/2, wd);
H2 = [bp +α* B (k)。^ 2,0,,黑洞;0,bp +α* B (k)。^ 2,黑洞,;是黑洞,bm +α* B (k)。^ 2 0;黑洞,是0,bm +α* B (k)。^ 2];
HT = H1 H2 +;
[V(:，:，k)，D(:，:，k)] = eig(HT);
结束
另一种选择是将它们存储在单元格数组中。
．
