【正文】 對(duì)角陣: 且 (324) 且 (325) (326) 其中 為矩陣Z和的非零特征值，和分別是Z和對(duì)應(yīng)于的特征向量，其中,則這兩個(gè)正交矩陣和一個(gè)對(duì)角陣滿足: (327)即可得到協(xié)方差矩陣S的正交歸一特征向量為 (328)其中為的特征值，為其對(duì)應(yīng)的正交歸一特征向量。 然后。根據(jù)閾值選定前K個(gè)特征向量,...構(gòu)成特征子空間，閾值的計(jì)算公式如下: (329) 通常情況下。 (3)計(jì)算訓(xùn)練集的人臉特征向量 人臉特征向量即將人臉圖像向特征空間投影之后得到的系數(shù)向量，即: (330) 將訓(xùn)練集中所有己知人臉圖像的圖像向量、特征向量,和其對(duì)應(yīng)的個(gè)人身份信息保存起來，就構(gòu)成了己知人臉圖像的訓(xùn)練數(shù)據(jù)庫。 (1)建立測(cè)試集 選取人臉庫中的每個(gè)人的部分圖像作為測(cè)試圖像，并進(jìn)行相應(yīng)的預(yù)處理，再將每個(gè)人臉圖像由二維向量轉(zhuǎn)化為一維向量。(2)計(jì)算測(cè)試集的人臉特征向量假設(shè)向量表示一幅測(cè)試圖像，將向特征子空間投影，即可得到該測(cè)試圖相對(duì)應(yīng)的人臉特征向量: (331) (3)選擇距離函數(shù)進(jìn)行分類判別 將人臉圖像投影到特征子空間，得到相應(yīng)的人臉特征向量之后，余下的任務(wù)就是如何判別測(cè)試圖像所屬的類別。首先我們要計(jì)算出圖像之間的相似度，然后選擇合適的分類器進(jìn)行分類判別。那么本章所選擇相似性度量工具是歐氏距離，然后選擇最近鄰分類器進(jìn)行圖像的分類。分類過程如下: 對(duì)所有訓(xùn)練樣本，…，每一個(gè)樣本屬于某個(gè)特定的類別(k=1,2,…,a)，其中a為訓(xùn)練樣本中人的個(gè)數(shù)，它們相對(duì)應(yīng)的特征向量為,… ,。對(duì)于給定的一個(gè)人臉測(cè)試圖像，其對(duì)應(yīng)的特征向量為，如果，并且有，那么分類結(jié)果即為。第四章 實(shí)驗(yàn)過程顯示及分析 引言 ORL人臉數(shù)據(jù)庫是目前使用最廣泛的人臉庫。該數(shù)據(jù)庫是由40個(gè)人，每人10幅，共400幅92X112大小的正面人臉圖像組成，其中包括不同時(shí)期同一背景的人臉面部表情和細(xì)節(jié)的變化圖像。本實(shí)驗(yàn)是從中選取了10個(gè)人，每人選取了表情變化較大的3幅圖像，其中選取每人的一張圖像，組成了有10張圖像的測(cè)試集，另外的20張組成訓(xùn)練集，在運(yùn)用PCA方法進(jìn)行人臉識(shí)。 ，編寫了一個(gè)簡(jiǎn)單的人臉識(shí)別程序，實(shí)現(xiàn)了隨即的單個(gè)樣本的人臉識(shí)別。 實(shí)驗(yàn)過程 本系統(tǒng)在自選的人臉數(shù)據(jù)庫上進(jìn)行簡(jiǎn)單的人臉識(shí)別，通過隨機(jī)輸入想要測(cè)試的圖片，系統(tǒng)會(huì)自動(dòng)讀取，進(jìn)行識(shí)別后，同時(shí)顯示出測(cè)試的圖片和識(shí)別后的圖片。系統(tǒng)界面截圖所示:第一步：，出現(xiàn)瀏覽文件夾對(duì)話框，點(diǎn)擊“確定”。,點(diǎn)擊“確定”。第二步：在出現(xiàn)的輸入對(duì)話框中，輸入要測(cè)試的照片編號(hào)，點(diǎn)擊“ok”，顯示出兩張照片，一個(gè)是測(cè)試的圖片本身，一個(gè)是訓(xùn)練集中對(duì)應(yīng)的照片。顯示出訓(xùn)練集中對(duì)應(yīng)的照片編號(hào)，識(shí)別工作完成。下面是輸入照片編號(hào)“3”和“5”，系統(tǒng)界面截圖：※ ※ ※ ※致 謝 大學(xué)時(shí)光，轉(zhuǎn)瞬及逝。在論文擱筆付梓之際，我的校園生活也即將告一段落，此時(shí)此刻，心中感慨頗多，既有收獲的喜悅，又有離別的難舍，還有對(duì)老師和同學(xué)的滿腔感激之情。 本課題的開展與研究是在我的導(dǎo)師李曉媛教授悉心指導(dǎo)下完成的，從課題的選題到論文的完成，導(dǎo)師都傾注了大量的時(shí)間和精力。導(dǎo)師經(jīng)常為我答疑解惑、開拓思路，提供了許多建設(shè)性的方案意見。她淵博的學(xué)識(shí)、嚴(yán)謹(jǐn)?shù)淖黠L(fēng)、寬廣的胸懷和淡泊明志的精神境界都給我以潛移默化的影響，將使我終生受益?！疤依畈谎裕伦猿甚琛?。對(duì)導(dǎo)師的感激之情非言語所能表達(dá)，導(dǎo)師無私的關(guān)懷和殷切的厚望惟以不斷的努力來報(bào)答。 在課題研究的過程中，得到了黃玲，劉穎等同學(xué)的幫助，在此表示感謝。感謝實(shí)驗(yàn)室諸多同學(xué)的關(guān)心、幫助，經(jīng)常的交流、探討讓我受益匪淺。感謝所有關(guān)心我的朋友。特別要感謝我的父母和家人，他們?nèi)找共賱冢燎诠ぷ?，為我解除了學(xué)習(xí)的后顧之憂。在我成長(zhǎng)的過程中，我的父母給予我許多學(xué)習(xí)上的支持和精神上的鼓勵(lì)，深深感謝他們的無私付出，深深感謝他們的養(yǎng)育之恩!參考文獻(xiàn)[1]徐飛，2002.[2]李剛，[J].計(jì)算機(jī)與現(xiàn)代化，2003, 3： 16.[3]周杰，盧春雨，張長(zhǎng)水，[[J].電子學(xué)報(bào)，2000, 4：2025.[4][D].安徽:中國(guó)科學(xué)技術(shù)大學(xué)，2007.[5]周激流，[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，200711：180184.[6] Yang M H, Ahuja N, Kriegrnan D. A survey on face detection methods[J].IEEE Trans PAMI, 2002,24(l)：3458.[7] Erik H, Boom K. Detecting Faces in Images: A Survey[J].Computer Vision and Image Understanding, 2001, 83：236274.[8]李華勝，張曄，[[J].北京交通大學(xué)學(xué)報(bào)，2001,25(2)：4549.[9]姜軍，張桂林一種基于知識(shí)的快速人臉檢測(cè)方法[[J].中國(guó)圖象學(xué)報(bào)，2002,7(1)：1419.[10]盧春雨，張長(zhǎng)水，[[J].清華大學(xué)學(xué)報(bào)，1999,39(l)：101105.附錄:主程序：exampleclear allclcclose all TrainDatabasePath = uigetdir(39。F:\Applications\work39。, 39。Select training database path39。 )。TestDatabasePath = uigetdir(39。F:\Applications\work39。, 39。Select test database path39。)。 prompt = {39。Enter test image name (a number between 1 to 10):39。}。dlg_title = 39。Input of PCABased Face Recognition System39。num_lines= 1。def = {39。139。}。 TestImage = inputdlg(prompt,dlg_title,num_lines,def)。TestImage = strcat(TestDatabasePath,39。\39。,char(TestImage),39。.jpg39。)。im = imread(TestImage)。 T = CreateDatabase(TrainDatabasePath)。[m, A, Eigenfaces] = EigenfaceCore(T)。OutputName = Recognition(TestImage, m, A, Eigenfaces)。 SelectedImage = strcat(TrainDatabasePath,39。\39。,OutputName)。SelectedImage = imread(SelectedImage)。 imshow(im)title(39。Test Image39。)。figure,imshow(SelectedImage)。title(39。Equivalent Image39。)。 str = strcat(39。Matched image is : 39。,OutputName)。disp(str)%函數(shù)CreateDatabase:function T = CreateDatabase(TrainDatabasePath)%%%%%%%%%%%%%%%%%%%%%%%% File managementTrainFiles = dir(TrainDatabasePath)。Train_Number = 0。 for i = 1:size(TrainFiles,1) if not(strcmp(TrainFiles(i).name,39。.39。)|strcmp(TrainFiles(i).name,39。..39。)|strcmp(TrainFiles(i).name,39。39。)) Train_Number = Train_Number + 1。 endend %%%%%%%%%%%%%%%%%%%%%%%% Construction of 2D matrix from 1D image vectorsT = []。for i = 1 : Train_Number str = int2str(i)。 str = strcat(39。\39。,str,39。.jpg39。)。 str = strcat(TrainDatabasePath,str)。 img = imread(str)。 img = rgb2gray(img)。 [irow icol] = size(img)。 temp = reshape(img39。,irow*icol,1)。 T = [T temp]。 End%函數(shù)EigenfaceCore：function [m, A, Eigenfaces] = EigenfaceCore(T) %%%%%%%%%%%%%%%%%%%%%%%% Calculating the mean image m = mean(T,2)。 % Train_Number = size(T,2)。 %%%%%%%%%%%%%%%%%%%%%%%% Calculating the deviation of each image from mean imageA = []。 for i = 1 : Train_Number temp = double(T(:,i)) m。 A = [A temp]。end L = A39。*A。[V D] = eig(L)。 L_eig_vec = []。for i = 1 : size(V,2) if( D(i,i)1 ) L_eig_vec = [L_eig_vec V(:,i)]。 endend Eigenfaces = A * L_eig_vec。 %函數(shù)Recognition：function OutputName = Recognition(TestImage, m, A, Eigenfaces)ProjectedImages = []。Train_Number = size(Eigenfaces,2)。for i = 1 : Train_Number temp = Eigenfaces39。*A(:,i)。 ProjectedImages = [ProjectedImages temp]。 end %%%%%%%%%%%%%%%%%%%%%%%% Extracting the PCA features from test imageInputImage = imread(TestImage)。temp = InputImage(:,:,1)。 [irow icol] = size(temp)。InImage = reshape(temp39。,irow*icol,1)。Difference = double(InImage)m。ProjectedTestImage = Eigenfaces39。*Difference。 Euc_dist = []。for i = 1 : Train_Number q = ProjectedImages(:,i)。 temp = ( norm( ProjectedTestImage q ) )^2。 Euc_dist = [Euc_dist temp]。end [Euc_dist_min , Recognized_index] = min(Euc_dist)。OutputName = strcat(int2str(Recognized_index),39。.jpg39。)。