线性判别式分析（Linear Discriminant Analysis LDA）MATLAB实现

线性判别式分析(Linear Discriminant Analysis, LDA)算法，利用MATLAB实现。附实例原始代码以及结果图

clc;
clear;
data = load （‘data_LDA.txt‘）;
cls1_data=[data（1:1:2011）data（1:1:2012）];%第一个类的训练集
cls2_data=[data（202:1:4021）data（202:1:4022）];%第二个类的训练集

%求期望
E_cls1=mean（cls1_data）;%第一类数据的期望矩阵
E_cls2=mean（cls2_data）;%第二类数据的期望矩阵
E_all=mean（[cls1_data;cls2_data]）;%所有训练集的期望矩阵

%%%%%%%%%%%%%%%%%%%%分类前画图%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:201
    plot（cls1_data（i1）cls1_data（i2）‘*r‘）;
    hold on;
end
  plot（E_cls1（1）E_cls1（2）‘^r‘）;
  hold on;
for i=1:201
    plot（cls2_data（i1）cls2_data（i2）‘*b‘）;
    hold on;
end
  plot（E_cls2（1）E_cls2（2）‘^b‘）;
  plot（E_all（1）E_all（2）‘vc‘）;  
  hold on;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 计算类间离散度矩阵：
  x1=E_all-E_cls1;
  x2=E_all-E_cls2;
  Sb=201*x1‘*x1/402+201*x2‘*x2/402;
%% 计算类内离散度矩阵
  y1=0;
  for

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      10827  2018-06-01 16:44  LDA\data_LDA.txt

     文件       2027  2018-10-30 13:07  LDA\LDA.m

     文件      29756  2018-10-29 17:59  LDA\Original_data.jpg

     文件        284  2018-10-29 17:58  LDA\Original_data.m

     文件      28653  2018-10-30 11:47  LDA\投影子空间.jpg

     文件      30130  2018-10-30 11:47  LDA\样本投影到子空间.jpg

     目录          0  2018-11-05 19:24  LDA

----------- ---------  ---------- -----  ----

               101677                    7