S=1/(n-1)*(X-X?)?*(X-X?)。其中，n表示样本容量，X表示n行p列的样本数据矩阵，X?表示p维向量，表示样本数据各变量的均值向量。(X-X?)表示样本数据矩阵每个元素减去对应变量的均值，得到的是一个n行p列的矩阵，(X-X?)?表示该矩阵的转置矩阵。最后，将两个矩阵相乘，再除以n-1，得到样本离差阵S。