平面波方法计算声子晶体能带结构的源程序（复式格子圆柱体按正方形排列,水、四氯化碳为散射体,水银为基体）-Plane wave method phononic crystal band structure of the source (double-cylinder square array grid, water, carbon tetrachloride, the scattering body, mercury as a matrix)