optics算法matlab代码