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個文件信息