基于MATLAB的Fisher线性判别代码

clear
clc
close all;
%m1,m2均为10个样本
% m1=[-0.4,0.58,0.089;-0.31,0.27,-0.04;0.38,0.055,-0.035;-0.15,0.53,0.011;-0.35,0.47,0.034;0.17,0.69,0.1;-0.011,0.55,-0.18];
% m2=[0.83,1.6,-0.014;1.1,1.6,0.48;-0.44,-0.41,0.32;0.047,-0.45,1.4;0.28,0.35,3.1;-0.39,-0.48,0.11;0.34,-0.079,0.14];
% a) 对类别m1,m2,利用fisher原则计算投影方向w
%  m1=[0 0 0; 1 0 0; 1 0 2; 1 1 0; 0.5 0 0.5;1 0.5 0;1 0 0.5];
%  m2=[0 0 1; 0 1 0; 0 1 1; 1 2 1;0 0.5 1;0 1 0.5; 0.5 1 2;0,1,2;0,1,1.2;0,1.2,1.3;0,0,2;0,1,2.2;1,2,0.5;1,2,0.5;1,2,0];
% m1=[6.3,3.3,6.0;5.8,2.7,5.1;7.1,3.0,5.9;6.3,2.9,5.6;6.5,3.0,5.8;7.6,3.0,6.6;4.9,2.5,4.5;7.3,2.9,6.3;6.7,2.5,5.8;7.2,3.6,6.1;6.5,3.2,5.1;6.4,2.7,5.3;6.8,3.0,5.5;5.7,2.5,5.0;5.8,2.8,5.1;
%     6.4,3.2,5.3;6.5,3.0,5.5;7.7,3.8,6.7;7.7,2.6,6.9;6.0,2.2,5.0;6.9,3.2,5.7;5.6,2.8,4.9;7.7,2.8,6.7;6.3,2.7,4.9;6.7,3.3,5.7;7.2,3.2,6.0;6.2,2.8,4.8;6.1,3.0,4.9;6.4,2.8,5.6;7.2,3.0,5.8;
%     7.4,2.8,6.1;7.9,3.8,6.4;6.4,2.8,5.6;6.3,2.8,5.1;6.1,2.6,5.6;7.7,3.0,6.1;6.3,3.4,5.6;6.4,3.1,5.5;6.0,3.0,4.8;6.9,3.1,5.4;6.7,3.1,5.6;6.9,3.1,5.1;5.8,2.7,5.1;6.8,3.2,5.9;6.7,3.3,5.7;
%     6.7,3.0,5.2;6.3,2.5,5.0;6.5,3.0,5.2;6.2,3.4,5.4;5.9,3.0,5.1];
% m2=[5.1,3.5,1.4;4.9,3.0,1.4;4.7,3.2,1.3;4.6,3.1,1.5;5.0,3.6,1.4;5.4,3.9,1.7;4.6,3.4,1.4;5.0,3.4,1.5;4.4,2.9,1.4;4.9,3.1,1.5;5.4,3.7,1.5;4.8,3.4,1.6;4.8,3.0,1.4;4.3,3.0,1.1;5.8,4.0,1.2;
%     5.7,4.4,1.5;5.4,3.9,1.3;5.1,3.5,1.4;5.7,3.8,1.7;5.1,3.8,1.5;5.4,3.4,1.7;5.1,3.7,1.5;4.6,3.6,1.0;5.1,3.3,1.7;4.8,3.4,1.9;5.0,3.0,1.6;5.0,3.4,1.6;5.2,3.5,1.5;5.2,3.4,1.4;4.7,3.2,1.6;
%     4.8,3.1,1.6;5.4,3.4,1.5;5.2,4.1,1.5;5.5,4.2,1.4;4.9,3.1,1.5;5.0,3.2,1.2;5.5,3.5,1.3;4.9,3.1,1.5;4.4,3.0,1.3;5.1,3.4,1.5;5.0,3.5,1.3;4.5,2.3,1.3;4.4,3.2,1.3;5.0,3.5,1.6;5.1,3.8,1.9;
%     4.8,3.0,1.4;5.1,3.8,1.6;4.6,3.2,1.4;5.3,3.7,1.5;5.0,3.3,1.4;7.0,3.2,4.7;6.4,3.2,4.5;6.9,3.1,4.9;5.5,2.3,4.0;6.5,2.8,4.6;5.7,2.8,4.5;6.3,3.3,4.7;4.9,2.4,3.3;6.6,2.9,4.6;5.2,2.7,3.9;
%     5.0,2.0,3.5;5.9,3.0,4.2;6.0,2.2,4.0;6.1,2.9,4.7;5.6,2.9,3.6;6.7,3.1,4.4;5.6,3.0,4.5;5.8,2.7,4.1;6.2,2.2,4.5;5.6,2.5,3.9;5.9,3.2,4.8;6.1,2.8,4.0;6.3,2.5,4.9;6.1,2.8,4.7;6.4,2.9,4.3;
%     6.6,3.0,4.4;6.8,2.8,4.8;6.7,3.0,5.0;6.0,2.9,4.5;5.7,2.6,3.5;5.5,2.4,3.8;5.5,2.4,3.7;5.8,2.7,3.9;6.0,2.7,5.1;5.4,3.0,4.5;6.0,3.4,4.5;6.7,3.1,4.7;6.3,2.3,4.4;5.6,3.0,4.1;5.5,2.5,4.0;
%     5.5,2.6,4.4;6.1,3.0,4.6;5.8,2.6,4.0;5.0,2.3,3.3;5.6,2.7,4.2;5.7,3.0,4.2;5.7,2.9,4.2;6.2,2.9,4.3;5.1,2.5,3.0;5.7,2.8,4.1];
% m1=[-300,250,0;-500,150,-100;-700,300,0;-450,120,0;-320,-50,0;-600,100,0;0,-90,0];
% m2=[1200,200,0;1800,-350,0;200,0,0;800,100,0;1600,-100,-200;900,100,0;400,-50,-130;300,-100,0;600,270,300;740,0,0;750,0,-200;480,25,30;550,-120,0];
m1=load('data.txt');
m2=load('data2.txt');
u1=mean(m1);%求均值
u2=mean(m2);
%计算类内散度Si和总类内散度Sw
% S1=(m1(1,:)-u1)'*(m1(1,:)-u1)
% S2=(m2(1,:)-u2)'*(m2(1,:)-u2)
% for i=2:7