clear?all
clc



prompt?=?{‘Enter?test?image?name?（a?number?between?1?to?10）:‘};
dlg_title?=?‘Input?of?ICA-based?Face?Recognition?System‘;
num_lines=?1;
def?=?{‘1‘};

TestImage??=?inputdlg（promptdlg_titlenum_linesdef）;
TestImage?=?strcat（‘F:\研一課程\圖像分析與計(jì)算機(jī)視覺\ICA?face?recognition\TestDatabase‘‘\‘char（TestImage）‘.BMP‘）;
im?=?imread（TestImage）;


allsamples=[];??%所有訓(xùn)練圖像?
for?i=1:20
????%a=imread（?strcat（?‘D:\Tools\Matlab\matlab?new?project\PCA_based?Face?Recognition?System\TrainDatabase‘?num2str（i）‘.bmp‘）;????%?讀取文件夾里所有圖片
????%?figure（i）imshow（a）;
????str?=?int2str（i）;%將整形轉(zhuǎn)換成字符型
????str?=?strcat（‘\‘str‘.bmp‘）;%str=‘\i.bmp‘
????str?=?strcat（‘F:\研一課程\圖像分析與計(jì)算機(jī)視覺\ICA?face?recognition\TrainDatabase‘str）;%讀入圖片i.bmp
????a?=?imread（str）;??????????????????%讀入圖片i.bmp
????b=a（?1:100*100?）;????????????????????????%?b是行矢量?1×N，其中N＝10000，提取順序是先列后行，即從上到下，從左到右
????b=double（b）;
????allsamples=[?allsamples;?b?];????????????%?allsamples是一個(gè)M*N?矩陣，allsamples中每一行數(shù)據(jù)代表一張圖片，其中M＝10
end

[MT]=size（allsamples）;?%獲取輸入矩陣的行/列數(shù)，行數(shù)為觀測數(shù)據(jù)的數(shù)目，列數(shù)為采樣點(diǎn)數(shù)
average=mean（allsamples‘）‘;%均值
for?i=1:M
????X（i:）=allsamples（i:）-average（i）*ones（1T）;
end

%白化/球化
Cx=cov（X‘1）;%計(jì)算協(xié)方差矩陣Cx
[eigvectoreigvalue]=eig（Cx）;%計(jì)算Cx的特征值和特征向量
W=eigvalue^（-1/2）*eigvector‘;%白化矩陣
Z=W*X;%新的觀測信號

%迭代
Maxcount=10000;%最大迭代次數(shù)
Critical=0.00001;%判斷是否收斂
m=M;%需要估計(jì)的分量的個(gè)數(shù)
W=rand（m）;%產(chǎn)生區(qū)間（0，1）均勻分布的隨機(jī)數(shù)，是m*m的矩陣
for?n=1:m
????WP=W（:n）;%初始權(quán)矢量（任意），只產(chǎn)生了一個(gè)列向量
????%Y=WP‘*Z
????%G=Y.^3;%G為非線性函數(shù)，可取y^3等
????%GG=3Y.^2;G的導(dǎo)數(shù)
????count=0;
????LastWP=zeros（m1）;%全0列向量
????W（:n）=W（:n）/norm（W（:n））;
????while?norm（WP-LastWP）>Critical?&&?norm（WP+LastWP）>Critical?
????????%abs（WP-LastWP）&abs（WP+LastWP）>Critical?%判斷收斂的句子
????????count=count+1;%迭代次數(shù)
????????LastWP=WP;%上次迭代的值
????????%WP=1/T*Z*（（LastWP‘*Z）.^3）‘-3*LastWP;
????????for?i=i:m
????????????WP（i）=mean（Z（i:）.*（tanh（（LastWP）‘*Z）））-（mean（1-（tanh（（LastWP））‘*Z）.^2））.*LastWP（i）;
????????end
????????WPP=zeros（m1）;
????????for?j=1:n-1
????????????WPP=WPP+（WP‘*W（:j））*W（:j）;
????????end
????????WP=WP-WPP;
????????WP=WP/（norm（WP））;
????????
????????if?count==Maxcount
????????????fprint（‘未找到相應(yīng)的信號‘）;
????????????return?;
????????end
????end
????W（:n）=WP;
end
Z=W‘*Z;


allcoor?=allsamples*Z‘?;??
????????????????????
??
????b=im（1:10000）;
????b=double（b）;
???
????tcoor=?b*Z‘;?????????
????for?k=1:20
???????
????????temp（k）=norm（tcoor-allcoor（k:））;
??????????
????end
???[distRecognized_index]=sort（temp）;
?
??
??
??
imshow（im）
title（‘Test?Image‘）;
a=Recognized_index（1）;
matchp=?strcat（‘F:\研一課程\圖像分析與計(jì)算機(jī)視覺\ICA?face?recognition\TrainDatabase‘‘\‘int2str（Recognized_index（1））‘.BMP‘）;?
figureimshow（matchp）;
title（‘Matched?image‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????3038??2013-06-16?17:01??ICA?face?recognition\recognition.m

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\1.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\10.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\2.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\3.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\4.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\5.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\6.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\7.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\8.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TestDatabase\9.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\1.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\10.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\11.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\12.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\13.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\14.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\15.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\16.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\17.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\18.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\19.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\2.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\20.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\3.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\4.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\5.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\6.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\7.bmp

?????文件??????30966??2013-06-16?14:17??ICA?face?recognition\TrainDatabase\8.bmp

............此處省略8個(gè)文件信息