某地区经勘探证明,A盆地是一个钾盐矿区,B盆地是一个钠盐(不含钾盐)矿区,其他盆地是否含钾盐有待判断。从A和B两个盆地各取5个盐泉样本,从其他盆地抽得8个盐泉样本,具体数据如下表所示,试对后8个待判盐泉进行钾性判别。
思路
距离判别是定义一个样本到某个总体的“距离”的概念,然后根据样本到各个总体的“距离”的远近来判断样本的归属。最常用的是马氏距离,其定义如下:
代码
% 综合实例——盐泉的钾性判别
clear;
clc;
load './res/spring.mat';
miu1 = mean(X1, 2);
miu2 = mean(X2, 2);
miu = (miu1+miu2)/2;
[~, n1] = size(X1);
[m, n2] = size(X2);
for i=1:m
ss1(i,:) = X1(i, :) - miu1(i);
ss2(i,:) = X2(i, :) - miu2(i);
end
s1 = ss1*ss1';
s2 = ss2*ss2';
V = (s1+s2)/(n1+n2-2)