正方晶格光子晶体能带 下面给出计算正方晶格光子晶体能带的代码。这里划分边界元的方式是从左上到左下,从左下到右下,从右上到右下,从左上到右上。 a=1; r=0.45*a; n=[1,sqrt(13)]; %折射率 Gfunc=@(x,y,k,b,n) -1i/4*besselh(0,k*n*sqrt((x(1)-b).^2+(x(2)-y).^2)); %计算...