本算法进一步深入熟悉平面波展开法
结果如下:与文献molding the flow of light 相一致
ps:当有多余G存在会产生多余线,需要进一步研究
主程序如下:
%This is a simple demo for 1D Photonic Crystals simulation %10 points is considered. %by Gao Haikuo %date:20170411 clear; clc; global NG G f Nkpoints eigenValue modeset kCorner global epsa epsb epssys a b1 b2 epssys=1.0e-6; %设定一个最小量,避免系统截断误差或除0错误 %this is the lattice vector and the reciprocal lattice vector a=