遇到了小bug：

R中矩阵相乘错误于A %*% B: 需要数值/复数矩阵/矢量参数

看到网上别人的做法，发现了用class(A)和class(B)之后才发现，是因为读入的时候数据的类型不对，A、B的类型并不是matrix，才导致了这个问题。

用as.matrix来变型一下，就OK了。