function?FaceRecognitionByBPneuralnetwork
%%?基于BP神經網絡的人臉識別
%%2014-05-14

clc

global?base?net?????????%定義全局變量，用于?TestofFaceRecognition.m?函數中

%%?第一步，利用主元分析PCA法提取圖像的特征矩陣

allsamples=[];????????????????????????????????%所有訓練圖像
for?i=1:40
????for?j=1:5
??????a=imread（strcat（‘D:\MATLAB\人臉庫\s‘num2str（i）‘\‘num2str（j）‘.bmp‘））;
??????b=a（1:112*92）;??????????????????????????%?b是行矢量?1×N，其中N＝10304，提取順序是先列后行，即從上到下，從左到右
??????b=double（b）;
??????allsamples=[allsamples;?b];?????????????%?allsamples?是一個M?*?N?矩陣，allsamples?中每一行數據代表一張圖片，其中M＝400
????end
end
samplemean=mean（allsamples）;??????????????????%?平均圖片，1?×?N
for?i=1:200
????xmean（i:）=allsamples（i:）-samplemean;????%?xmean是一個M?×?N矩陣，xmean每一行保存的數據是“每個圖片數據-平均圖片”
end

sigma=xmean*xmean‘;????????%?M?*?M?階矩陣
[v?d]=eig（sigma）;
d1=diag（d）;
[d2?index1]=sort（d1）;???????%以升序排序
cols=size（v2）;????????????%?特征向量矩陣的列數
for?i=1:cols
????vsort（:i）?=?v（:?index1（cols-i+1）?）;????%?vsort?是一個M*col（注:col一般等于M）階矩陣，保存的是按降序排列的特征向量每一列構成一個特征向量
????dsort（i）???=?d1（?index1（cols-i+1）?）;?????%?dsort?保存的是按降序排列的特征值，是一維行向量
end????%完成降序排列
%以下選擇85%的能量
dsum?=?sum（dsort）;
????dsum_extract?=?0;
????p=?0;
????while（?dsum_extract/dsum?????????p=p?+?1;
????????dsum_extract?=?sum（dsort（1:p））;
????end
i=1;
%計算特征臉形成的坐標系
while?（i<=p?&&?dsort（i）>0）
????base（:i）?=?dsort（i）^（-1/2）?*?xmean‘?*?vsort（:i）;?????%?base是N×p階矩陣，除以dsort（i）^（1/2）是對人臉圖像的標準化
????i?=?i?+?1;
end

%下面的代碼，將訓練樣本對坐標系上進行投影得到一個?M*p?階矩陣allcoor
allcoor?=?allsamples?*?base;




%%?第二步，創建并訓練BP神經網絡

%生成訓練BP神經網絡的輸入?P?
P=?mapminmax（allcoor）;???????????????%歸一化

%生成目標輸出矢量?T
T=zeros（20040）;
?for?i=1:40
????for?j=1:5
??????T（（i-1）*5+ji）=1;
??end
?end
?
%打亂訓練樣本順序
gx2（:1:47）=P;
gx2（:48:87）=T;
xd=gx2（randperm（numel（gx2）/87）:）;
gx=xd（:1:47）;d=xd（:48:87）;
P=gx‘;
T=d‘;

%創建BP神經網絡
[RQ]=size（P）;
[S2Q]=size（T）;
net=newcf（minmax（P）T[6015]{‘tansig‘‘logsig‘}‘trainscg‘）;

%訓練BP神經網絡
net.trainparam.epochs=5000;
net.trainparam.goal=0.0001;
net.divideFcn?=?‘‘;
net=train（netPT）;

%仿真BP神經網絡
Y=sim（netP）;




%%?第三步，測試BP神經網絡并計算其識別率

%測試BP神經網絡
s=0;
for?i=1:40
????for?j=6:10??????????????????????????????%讀入40x5副測試圖像
?????????a=imread（strcat（‘D:\MATLAB\人臉庫\s‘num2str（i）‘\‘num2str（j）‘.bmp‘））;
?????????b=a（1:10304）;
?????????b=double（b）;
?????????tcoor=?b?*?base;???????????????????%計算坐標，是1×p階矩陣
?????????X?=?mapminmax（tcoor）;??????????????%歸一化
?????????Z=sim（netX‘）;
?????????[ziindex2]=max（Z）;
?????????if?index2==i???
?????????????s=s+1;
?????????else
?????????????i??????????????????????????????%輸出識別出錯的那個人?i
?????????????j??????????????????????????????%輸出識別出錯的那張圖片?j
?????????????index2?????????????????????????%輸出誤識別成的那個人
?????????end
?????end
end

%計算識別率
accuracy=s/Q




end
%%?完成

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????11384??2010-03-24?19:05??人臉庫\s1\1.bmp

?????文件??????11384??2010-03-24?19:43??人臉庫\s1\10.bmp

?????文件??????11384??2010-03-24?19:42??人臉庫\s1\2.bmp

?????文件??????11384??2010-03-24?19:44??人臉庫\s1\3.bmp

?????文件??????11384??2010-03-24?19:44??人臉庫\s1\4.bmp

?????文件??????11384??2010-03-24?19:44??人臉庫\s1\5.bmp

?????文件??????11384??2010-03-24?19:45??人臉庫\s1\6.bmp

?????文件??????11384??2010-03-24?19:45??人臉庫\s1\7.bmp

?????文件??????11384??2010-03-24?19:45??人臉庫\s1\8.bmp

?????文件??????11384??2010-03-24?19:45??人臉庫\s1\9.bmp

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\1.bmp

?????文件??????10318??1994-04-18?14:18??人臉庫\s10\1.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\10.bmp

?????文件??????10318??1994-04-18?14:07??人臉庫\s10\10.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\2.bmp

?????文件??????10318??1994-04-18?14:18??人臉庫\s10\2.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\3.bmp

?????文件??????10318??1994-04-18?14:18??人臉庫\s10\3.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\4.bmp

?????文件??????10318??1994-04-18?14:18??人臉庫\s10\4.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\5.bmp

?????文件??????10318??1994-04-18?14:18??人臉庫\s10\5.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\6.bmp

?????文件??????10318??1994-04-18?14:07??人臉庫\s10\6.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\7.bmp

?????文件??????10318??1994-04-18?14:07??人臉庫\s10\7.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\8.bmp

?????文件??????10318??1994-04-18?14:07??人臉庫\s10\8.pgm

?????文件??????11384??2010-03-25?19:23??人臉庫\s10\9.bmp

?????文件??????10318??1994-04-18?14:07??人臉庫\s10\9.pgm

............此處省略807個文件信息