光子矩阵计算