【正文】 empYbuffer + 50*720+i)=0xff。 *(Uint8 *)(tempYbuffer + 200*720+i)=0xff。 } for(j=50。j++) { *(Uint8 *)(tempYbuffer + j*720+260)=0xff。 } 中值濾波處理 中值濾波是一種非線性的信號處理方法。中值濾波在一定的條件下可以克服線性濾波器如最小均方濾波、均直濾波等帶來的圖像細(xì)節(jié)模糊，而且對濾除脈沖干擾及圖像掃描噪聲最為有效。但是對于一些細(xì)節(jié)多，特別是點、線、尖頂細(xì)節(jié)多的圖像不宜采用中值濾波。對于奇數(shù)個元素，中值是指按大 小排序后，中間的數(shù)值；對于偶數(shù)個元素，中值是指排序后中間兩個元素灰度值的平均值。顯然有 IPPT ? ,所以 P 是正交矩陣 .上面的變換過程即??????? 21 ??APP T .由于 ? ?? ? ?????? ???? ????? ????? ????? 22250 225021222221112112212112212222211 s i nac o sas i nac o sas i naa. c o sas i i nas i nac o saAPP T所以只要選擇 ? 滿足 ? ? 02221 121122 ??? ?? c osas inaa 即 22111222 aa atan ??? (6) (當(dāng) 2211 aa ? 時，可選取 4??? ) APPT 就成對角陣，這時 A 的特征值為 ???? 2122222111 s i nas i nac o sa ??? ???? 2122222112 s i nac o sas i na ??? 相應(yīng)的特征向量為 ??????????????? ???? c o ss i nv,s i nc o sv21 雅可比方法的基本思想是通過一系列的由平面旋轉(zhuǎn)矩陣構(gòu)成的正交 變換將實對稱矩陣逐步化為對角陣 ,從而得到 A 的全部特征值及其相應(yīng)的特征向量 .首先引進 nR 中的平面旋轉(zhuǎn)變換 .變換 ???????????j,ikyxco sysinyxsinyco syxkkiijiii???? (7) 記為 yPx ij? ,其中 21 ?ijPjic o ss i ns i nc o s???????????????????????????111?????????? (8) i j ? ?Tn, x...xxx 21? ? ?Tn, y...yyy 21? 則稱 yPx ij? 為 nR 中 ji x,x 平面內(nèi)的一個平面旋轉(zhuǎn)變換 , ijP 稱為 jix,x 平面內(nèi)的平面旋轉(zhuǎn)矩陣 .容易證明 ijP 具有如下簡單性質(zhì)： ① ijP 為正交矩陣 . ② ijP 的主對角線元素中除第 i 個與第 j 個元素為 ?cos 外，其它元素均為 1；非對角線元素中除第 i 行第 j 列元素為 ?sin? ,第 j 行第 i 列 元素為 ?sin 外 ,其它元素均為零 . ③ APT 只改變 A 的第 i 行與第 j 行元素， AP 只改變 A 的第 i 列與第 j 列元素 ,所以 APPT 只改變 A 的第 i 行、第 j 行、第 i 列、第 j 列元素 . 設(shè) ? ? ? ?3?? ? naA nnij 為 n 階實對稱矩陣 , 0?? jiij aa 為一對非對角線元素 .令 ? ?? ? nnijT aAPPA ??? 11 則 1A 為實對稱矩陣 ,且 1A 與 A 有相同的特征值 .通過直接計算知 ??????????????????????????????????j,il,kaaj,ikc o sas i naaaj,iks i nac o saaac o sas i n)aa(aas i nac o sas i naas i nas i nac o saakl)(kljkjk)(jk)(jkjkik)(ki)(ikijiijj)(ji)(ijijjjii)(jjijjjii)(ii1111111221221222122???????????? (9) 當(dāng)取 ? 滿足關(guān)系式 22 jjiiijaa atan ?? 22? （ 10） 時， ? ? ? ? 011 ?? jiij aa ，且 ? ?? ? ? ?? ?? ?? ? ? ?? ?? ?? ??????????????????j,il,kaaaaaaaj,ikaaaaklklijjjiijjiijkikjkik22122221212221212 (11) 由于在正交相似變換下 ,矩陣元素的平方和不變，所以若用 ? ?AD 表示矩陣 A 的對角線元素平方 和，用 ??AS 表示 A 的非對角線元素平方和 ,則由 (11)式得 ?????????212122ijija)A(S)A(Sa)A(D)A(D (12) 這說明用 ijP 對 A 作正交相似變換化為 1A 后 , 1A 的對角線元素平方和比 A 的對角線元素平方和增加了 22ija , 1A 的非對角線元素平方和比 A 的非對角線元素平方和減少了 22ija ,且將事先選定的非對角線元素消去了 (即 ?? 01 ?ija ).因此，只要我們逐次地用這種變換，就可以使得矩陣 A 的非對角線元素平方和趨于零 ,也即使得矩陣 A 逐步化為對角陣 . 這里需要說明一點：并不是對矩陣 A 的每一對非對角線非零元素進行一次這樣的變換就能得到對角陣 .因為在用變換消去 ija 的時 候 ,只有第 i 行、第j 行、第 i 列、第 j 列元素在變化 ,如果 ika 或 kjP 為零 ,經(jīng)變換后又往往不是零了 . 雅可比方法就是逐步對矩陣 A 進行正交相似變換 ,消去非對角線上的非零元素 ,直到將 A 的非對角線元素化為接近于零為止，從而求得 A 的全部特征值 ,把逐次的正交相似變換矩陣乘起來 ,便是所要求的特征向量 . 雅可比方法的計算步驟歸納如下： 第一步 在矩陣 A 的非對角線元素中選取一個非零元素 ija .一般說來，取絕對值最大的非對角線元素 。 第三步 111 APPA T? , 1A 的元素由公式 (9)計算 . 第四步 以 1A 代替 A ，重復(fù)第一、二、三步求出 2A 及 2P ，繼續(xù)重復(fù)這一過程，直到 mA 的非對角線元素全化為充分小 (即小于允許誤差 )時為止 . 第五步 mA 的對角線元素為 A 的全部特征值的近似值 , mP...PPP 21? 的第 j 23 列為對應(yīng)于特征值 j? ( j? 為 mA 的對角線上第 j 個元素 )的特征向量 . 用雅可比方法求得的結(jié)果精度都比較高 ,特別是求得的特征向量正交性很好 ,所以雅可比方法是求實對稱矩陣的全部特征值及其對應(yīng)特征向量的一個較好的方法 .但由于上面介紹的雅可比方法 ,每次迭代都選取絕對值最大的非對角線元素作為消去對象 ,花費很多機器時間 .另外當(dāng)矩陣是稀疏矩陣時 ,進行正 交相似變換后并不能保證其稀疏的性質(zhì) ,所以對階數(shù)較高的矩陣不宜采用這種方法 . 目前常采用一種過關(guān)雅可比方法 .這種方法是選取一個單調(diào)減小而趨于零的數(shù)列 ??na (即 ...aa ?? 21 ,且 0??? nn alim )作為限值，這些限值稱為 ”關(guān) ” ,對矩陣的非對角線元素規(guī)定一個順序 (例如先行后列、自左至右的順序 ).首先對限值 1a 按規(guī)定的順序逐個檢查矩陣的非對角線元素 ,碰到絕對值小于 1a的元素就跳過去，否則就作變換將其化為零 .重復(fù)上述過程 ,直到所有的非對角元素的絕對值都小于 1a 為止 .再取 ,...a,a 32 類似處理，直到所有的非對角線元素的絕對值都小于 ma 時 ,迭代停止 .這時的 ma 應(yīng)小于給定的 誤差限 ? . 實際運算中常用如下的辦法取限值：對于矩陣 A ，計算 0AA? 的非對角線元素平方和 ? ?0AS ,任取 nN? ,取 ? ? ? ??,kNASa kk 211 ?? ? PCA 算法 PCA 算法簡介 主成分分析 (Principal Component Analysis, 簡稱 PCA)是一種常用的基于變量協(xié)方差矩陣對信息進行處理、壓縮和抽提的有效方法。 PCA 方法的基本原理是 :利用 KL 變換抽取人臉的主要成分，構(gòu)成特征臉空間，識別時將測試圖像投影到此空間，得到一組投影系數(shù)， 通過與各個人臉圖像比較進行識別 利用特征臉法進行人臉識別的過程由訓(xùn)練階段和識別階段兩個階段組成 ：訓(xùn)練階段，識別階段。 ? ?Txxxx 20211 , .. .,???????????987654321????? 2 0012 001 i i ix?? ?? ????? ? ?????? ixd iiTiTii AAddC 20212021 2001?? ??? ?2 0 021 , . . . ,dddA ??AATTAA 25 求出 的特征值 及其正交歸一化特征向量 根據(jù)特征值的貢獻率選取前 p 個最大特征向量及其對應(yīng)的特征向量 貢獻率是指選取的特征值的和與占所有特征值的和比，即： . ............................................................ (57) 一般取 即使訓(xùn)練樣本在前 p 個特征向量集上的投影有 99%的能量 求出原協(xié)方差矩陣的特征向量 ........................................... (58) 則 “特征臉 ”空間為： .................................... (59) 第六步 將每一幅人臉與平均臉的差值臉矢量投影到 “特征臉 ”空間，即 ......................... (510) 識別階段 第一步：將待識別的人臉圖像 與平均臉的差值臉投影到特征空間，得到其特征向量表示： ……………………………… (511) 第二步：定義閾值 ……………… . (512) 第三步：采用歐式距離來計算 與每個人臉的距離 ...................... (513) 為了區(qū)分人臉和非人臉，還需要計算原始圖像 與由特征臉空間重建的圖像 AAT i? i?aiiipiii????????20011???%99?a),...,2,1(1 piAvu iii?? ?? ?puuuw ，, . . ., 21?? ?2 0 0, . . . ,2,1??? idw iTi?? ?????? ? Tw? ? 200,...,2,1,max21 , ????? jijiji??? i?? ?200,...,2,122 ????? ? iii?f?? 26 之間的距離 .............................................................