二维光子晶体能带的程序采用平面波传输法.docxVIP

pwm法计算二维光了晶体能带的另一个程序 已经分享了一个，这是另一个編写方法不同滞度更高。程序完整含冇注释，能同时显示TE和TM波的能带 结构。略作修改也可分开显示。 function [ output_args ] = Untitlcd2( input_args ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here %二维光子晶体带结构计算 %计算二维正方格子或三角格子， %介质圆柱(介电常数a)立于背景(介电常数b)之屮 clear; clc; tic; epssys=1.0e-6; %设定一个最小量，避免系统截断误差或除()错误 %定义实际的止空间格了基欠 latticetype=O; %定义格子类型，1计算正方格子，0计算三角格子 switch latticetype case 1 a = 1; al =a* [1 0 0]; a2 = a* [0 1 0]; a3 = a* LOO 1J; case 0 a = 1; al = a* [ sqrt(3)/2 1/2 0]; a2 = a* [-sqrt(3)/2 1/2 0]; a3 = a* [00 1J; end; %定义晶格的参数 epsa = 1; %介质柱的介电常数 epsb = 10.5; %背景的介电