code

S = np.linalg.inv(Sw)*SB

Sw和SB都是ndarray，那么这里的乘号含义应该是矩阵各个位置相乘。如果我没有搞错的话，应该改为矩阵乘法（.dot）才对吧？

这样改的话，结果跟sklearn里边的也比较像。

reconstruct版里边的实现和这个一样。我不确定仓库里是否有其他类似问题。

（写作业的时候发现对不上，研究了好久。说实话，numpy里边一种乘号还能有两种意义的，实在是令人震惊...）