和前面体光子晶体一样，代码也是从Photonic Crystals Physics and Practical Modeling (Igor A. Sukhoivanov, Igor V. Guryev 这本书上扒下来的(page 197) 这个代码同样用的数值法算的相对介电函数的傅里叶系数 书上的原始代码似乎有些问题，它的结果和FEM的结果符合得不是很好 原始代码为： function