1702660218
1702660219
1702660220
1702660221
1702660222
图12.129 样本协方差矩阵图
1702660223
1702660224
③估计总体协方差矩阵。操作过程如下：
1702660225
1702660226
1）设置单元格区域。在J3单元格、K2单元格和J4单元格、L2单元格中分别输入“物质1”和“物质2”；合并单元格J1、K1、L1，合并居中后输入“总体协方差矩阵”。单元格区域设置如图12.130所示。
1702660227
1702660228
1702660229
1702660230
1702660231
图12.130 单元格区域设置
1702660232
1702660233
2）在K3单元格中输入公式”=（G8+G13）/（5+5-2）”，按下Enter键即可得到物质1含量的总体方差；同理，在K4单元格中输入公式”=（G9+G14）/（5+5-2）”，按下Enter键即可得到物质1含量与物质2含量的总体协方差；在L3单元格中输入公式”=（G14+G9）/（5+5-2）”，按下Enter键即可得到物质2含量与物质1含量的总体协方差；在L4单元格中输入公式”=（H9+H14）/（5+5-2）”，按下Enter键即可得到物质2含量与物质2含量的总体协方差。至此，可得出总体协方差矩阵图，如图12.131所示。
1702660234
1702660235
1702660236
1702660237
1702660238
图12.131 总体协方差矩阵图
1702660239
1702660240
④求解总体协方差矩阵的逆矩阵。关于求解逆矩阵，Excel 2013提供了一个函数命令，即”MINVERSE”，具体的函数形式如下：
1702660241
1702660242
MINVERSE（array）
1702660243
1702660244
其中，“MINVERSE”函数用来计算某矩阵的逆矩阵；“array”指一个行数等于列数的矩阵区域。
1702660245
1702660246
在本例中运用”MINVERSE”函数来计算总体协方差矩阵的逆矩阵的操作过程如下：
1702660247
1702660248
1）设置单元格区域。在J8、K7单元格和J9、L7单元格中分别输入“物质1”和“物质2”；合并单元格J6、K6、L6，合并居中后输入“总体协方差矩阵的逆矩阵”。单元格区域设置如图12.132所示。
1702660249
1702660250
1702660251
1702660252
1702660253
图12.132 单元格区域设置
1702660254
1702660255
2）在K8单元格中输入函数”=MINVERSE（K3:K3）”，在K9单元格中输入公式”=MINVERSE（K4:L3）”，在L8单元格中输入公式”=MINVERSE（L3:K4）”，在L9单元格中输入公式”=MINVERSE（L4:L4）”，分别按下Enter键即可得到总体协方差矩阵的逆矩阵图，如图12.133所示。
1702660256
1702660257
1702660258
1702660259
1702660260
图12.133 总体协方差矩阵的逆矩阵图
1702660261
1702660262
⑤计算a。根据距离判别法的原理，因为a=总体协方差矩阵的逆矩阵*（肉食动物的平均值矩阵-素食动物的平均值矩阵），所以，此例中将使用Excel 2013中矩阵相乘的函数”MMULT”。该函数的具体形式如下：
1702660263
1702660264
MMULT（array1，array2）
1702660265
1702660266
其中，“MMULT”函数用来计算两个矩阵的乘积；“array1”、“array2”分别指纳入乘法运算中的两个矩阵，第一个矩阵的列数必须等于第二个矩阵的行数，结果矩阵的行数与第一个矩阵相同，列数与第二个矩阵相同。由于是矩阵的运算，所以最后均以按下Ctrl+Shift+Enter键结束。
1702660267