pwm法计算二维光子晶体能带的另一个程序.docVIP

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