理解fisher判别分析

原理

fisher判别法是去找一个过原点的直线,这个直线要达到的效果如图:
在这里插入图片描述
具体实现方法:
在这里插入图片描述
注:这里的Sw为离差阵,也叫做协方差矩阵。
w即为这条直线的方向

直线为Y = wX

临界点y0 = (n0 * u0 + n1 * u1) / (n0 + n1)

判别时计算数据在线上的投影,将此投影与聚类的中心投影相比较得出此数据应属于哪一类。

实战

例:为了解某河段As,Pb污染状况,设在甲,乙两地监测,采样测的这两种元素在水中和底泥中的浓度(如下表)。依据这些数据判别未知样本是从哪个区域采得的。
在这里插入图片描述
注:后两行为需要判别属于哪一类的样本。


clc
clear
X=load("fisher_values.txt");
x1=X(1:5,:);%甲地的已知数据
x2=X(6:10,:);%乙地的已知数据
r1=size(x1,1);%相当于求甲地的样本数
r2=size(x2,1);%相当于求乙地的样本数
sample=X(11:12,:);%待分类的数据
r3=size(sample,1);
%求各类别的均值
ave1=mean(x1);
avel=ave1';
ave2=mean(x2);
ave2=ave2';
%离差矩阵即协方差矩阵。
s1=cov(x1)*(r1-1);
s2=cov(x2)*(r2-1);
sw=s1+s2;%求出协方差矩阵 
w=inv(sw)*(ave1-ave2)*(r1+r2-2);%求出w
y1=mean(w'.*ave1);
y2=mean(w'.*ave2);
y0=(r1*y1+r2*y2)/(r1+r2);
y0=y0';
%利用0,1来区分类别
for i=1:r3
  y(i,:)=w'*sample(1,:)';
   if y(i,:)>y0
      m(i)=0;
   else
      m(i)=1;
   end
end
m %m表示判断结果
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值