LDA人脸识别自编程序

LDA人脸识别（自编程序）

clc
clear
NumClass=2;
NumEachClass=10;
Trainsamples=[];%所有训练图像 
for i=1:2
   for j=1:5
    a=imread（strcat（‘ORL\s‘num2str（i）‘_‘num2str（j）‘.bmp‘））; %读图像
    c=imresize（a0.5）;
    b=c（1:56*46）; % b是行矢量 1×N，其中N＝10304，提取顺序是先列后行，即从上到下，从左到右
    b=double（b‘）; %
    Trainsamples=[Trainsamples b];
end
end
mean1=sum（Trainsamples（:1:5）2）/5;
mean2=sum（Trainsamples（:6:10）2）/5;
S1=0;
S2=0;
for i=1:5
    S1=S1+（Trainsamples（:i）-mean1）*（Trainsamples（:i）-mean1）‘;
    S2=S2+（Trainsamples（:i+5）-mean2）*（Trainsamples（:i+5）-mean2）‘;
end
Sw=S1+S2;
Wu=inv（Sw）*（mean1-mean2）;
Y1=Wu‘*Trainsamples（:1:5）;
Y2=Wu‘*Trainsamples（:6:10）;
Y1mean=mean（mean（Y1））;
Y2mean=mean（mean（Y2））;
W0=（Y1mean+Y2mean）/2;
if Y1mean    labe1=0;
else 
    label=1;
end
%%%%%%%%%%%%%%%%%样本测试
Testsamples=[];%所有训练图像
for i=1:2
   for j=6:10
    a=imread（strcat（‘ORL\s‘num2str（i）‘_‘num2str（j）‘.bmp‘））; %读图像
    c=imresize（a0.5）;
    b=c（1:56*46）; % b是行矢量 1×N，其中N＝10304，提取顺序是先列后行，即从上到下，从左到右
    b=double（b‘）; %
    Testsamples=[Testsamples b];
end
end
TestLabel=Wu‘*Testsamples;
TestLabel-W0

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

     文件       1271  2018-05-01 17:37  LDA\LDA.asv

     文件       1200  2018-05-01 17:40  LDA\LDA.m

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_1.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_10.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_2.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_3.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_4.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_5.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_6.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_7.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_8.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s10_9.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_1.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_10.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_2.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_3.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_4.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_5.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_6.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_7.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_8.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s11_9.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_1.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_10.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_2.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_3.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_4.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_5.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_6.bmp

     文件      11382  2012-05-12 10:11  LDA\ORL\s12_7.bmp

............此处省略377个文件信息