平面波展开法计算声子晶体带隙

代码说明：

平面波展开法计算声子晶体带隙clear; clc; tic; epssys=1.0e-6; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%a = 0.02; a1 = a* [1 0 0];a2 = a* [0 1 0];a3 = a* [0 0 1];%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%da = 11600; db = 1300; la1a=1.49e10; la1b=4e4; %la2a=4.20e10; %la2b=6e5; f = 0.3; Au = norm(cross(a1,a2)); Rc = 0.008; Ac = pi*(Rc)^2;