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
1702660268
因此，此例中，分别在K12单元格和K13单元格中输入公式”=MMULT（K8:L8，J12:J13）“和”=MMULT（K9:L9，J12:J13）”，按下Ctrl+Shift+Enter组合键即可得到a的计算结果，如图12.134所示。其中，J12=G4-G3，J13=H4-H3。
1702660269
1702660270
1702660271
1702660272
1702660273
图12.134 指标平均类型之差及a图
1702660274
1702660275
⑥计算判别函数。根据距离判别法的原理，可知判别函数W=a的转置*（待判样本矩阵-样本平均值矩阵）。所以，此例中可以直接计算判别函数。过程如下：
1702660276
1702660277
1）首先计算“待判样本矩阵-样本平均值矩阵”，以计算样本B1为例，在B17单元格中输入公式”=C12-0.5*（G3+G4）”，并在B18单元格中输入公式”=D12-0.5*（H3+H4）”，分别按下Enter键即可得到样本B1的“待判样本矩阵-样本平均值矩阵”。同理，在C17单元格和C18单元格中分别输入公式”=C13-0.5*（G3+G4）“和”=D13-0.5*（H3+H4）”，分别按下Enter键即可得到样本B2的“待判样本矩阵-样本平均值矩阵”。在D17单元格和D18单元格中分别输入公式”=C14-0.5*（G3+G4）“和”=D14-0.5*（H3+H4）”，分别按下Enter键即可得到样本B3的“待判样本矩阵-样本平均值矩阵”。“待判样本矩阵-样本平均值矩阵”计算结果如图12.135所示。
1702660278
1702660279
