%利用類間離散度矩陣進行K-L變換的處理

clear;?close?all;clc
load?20064817924orl_faces_112x92;
row=112;
col=92;

total=40;
g=zeros（row*coltotal）;
?
for?i=1:total
???cur?=?faces（:（i-1）*10+1:（i-1）*10+5）;
???g（:i）=mean（cur2）;????%?計算第i個人前5幅圖像的均值
end

d=mean（g2）;?%?計算所有人前5幅圖像總的平均值

for?i=1:total
????X（:i）=g（:i）-d;
end
R=X‘*X;

[VE]=eig（R）;
E=fliplr（E）;?%?矩陣左右顛倒，將特征值按大小順序排列
V=fliplr（V）;?%?相應地調整特征向量

U=X*V;
U=orth（U）;?%?特征向量正交歸一化

lmE=mean（E）*total;??????%?改變特征值的存儲形式，由矩陣（對角元素）改存為向量
totalE=mean（lmE）*total;?%?特征值之和
for?i=1:total
????curTotalE=mean（lmE（1:i））*i;
????ratio=curTotalE/totalE;
????if?ratio>0.9
????????break;
????end
end

sprintf（‘利用了%d個特征向量組成特征臉空間‘i）

A=U（:1:i）‘;
y=A*g;%每個人的均值圖像向特征空間投影，作為個人識別的依據

%?顯示“特征臉”圖像

%?[rrcc]=size（U）;
%?dd=zeros（rowcol）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1619??2009-06-25?09:32??faceRecogSb.m
?????文件????????1648??2011-06-05?22:34??faceRecogSigma.asv
?????文件????????1648??2011-06-05?22:37??faceRecogSigma.m
?????文件?????????153??2009-06-24?17:53??faceshow.m
?????文件?????????843??2009-06-24?20:48??faceSigma.m
?????文件?????4121792??2006-09-06?08:25??orl_faces_112x92.mat
?????文件?????????392??2009-06-24?17:43??說明.txt