matlab计算光子晶体光纤带隙,使用FDTD方法,方便好用