耦合矩阵计算