Matlab综合实例——盐泉的钾性判别

最新推荐文章于 2023-02-15 22:29:16 发布

原创

最新推荐文章于 2023-02-15 22:29:16 发布 · 1.6k 阅读

7 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #数据分析

利用Matlab，通过马氏距离对13个盐泉样本进行钾性判别，旨在确定含钾盐泉。代码实现包括距离计算与结果可视化，最终得出的判别公式准确无误，所有样本判别正确。

在这里插入图片描述

某地区经勘探证明，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+