说明：  二维声子晶体能带结构计算程序，采用的是平面波展开法，并且可以用matlab到出能带结构图
(Dimensional phononic crystal band structure calculation program, using the plane wave expansion method, and can be used to matlab energy band structure of FIG.)