band.m

用传输矩阵的方法来计算光子晶体带隙

w=[0.01:0.01:2]*2*pi*3e8;
kA=w*3.4/3e8;
kB=w*1/3e8;
dA=0.3;
dB=0.3;
k=acos(cos(kA*dA).*cos(kB*dB)-1/2*(kA./kB+kB./kA).*sin(kA*dA).*sin(kB.*dB));

plot(k,w/2/pi/3e8)