动态因子分析中需要计算不同时间的协差阵的平均数，文章dynamic factor analysis with stata给出一个代码

matrix ST=J(3,3,0)
forvalues i=1(10)40 {
matrix C=A[`i'..(`i'+10-1),1...]
svmat C
matrix accum cov = C1-C3, deviations noconstant
matrix cov=cov/(r(N)-1)
matrix ST=ST+cov
drop C1-C3
}


然后本人用excel验证时，把每个时间的协差阵分别算出来，然后再算这些协差阵的平均数，发现两种方法得出的结果不一致。
因为原文中的stata代码看不懂，也不清楚到底是哪里的问题，有没有朋友能给解释一下？谢谢

附件中是我模拟的数据和计算过程