马尔可夫矩阵法计算