基于小波变换的人脸识别
近年来，小波变换在科技界备受重视，不仅形成了一个新的数学分支，而且被广泛地应用于模式识别、信号处理、语音识别与合成、图像处理、计算机视觉等工程技术领域。小波变换具有良好的时频域局部化特性，且其可通过对高频成分采取逐步精细的时域取样步长，从而达到聚焦对象任意细节的目的，这一特性被称为小波变换的“变聚焦”特性，小波变换也因此被人们冠以“数学显微镜”的美誉。
具体到人脸识别方面，小波变换能够将人脸图像分解成具有不同分辨率、频率特征以及不同方向特性的一系列子带信号，从而更好地实现不同分辨率的人脸图像特征提取。
4.1 小波变换的研究背景
法国数学家傅立叶于1807年提出了著名的傅立叶变换，第一次引入“频率”的概念。傅立叶变换用信号的频谱特性来研究和表示信号的时频特性，通过将复杂的时间信号转换到频率域中，使很多在时域中模糊不清的问题，在频域中一目了然。在早期的信号处理领域，傅立叶变换具有重要的影响和地位。定义信号f(t)为在(-∞，+∞)内绝对可积的一个连续函数，则f(t)的傅立叶变换定义如下：
F???????f?t?e?j?tdt ??(4-1)
傅立叶变换的逆变换为：
1f?t??2??????F???ej?td? (4-2)
从上面两个式子可以看出，式（4-1）通过无限的时间量来实现对单个频率的频谱计算，该式表明F(?)这一频域过程的任一频率的值都是由整个时间域上的量所决定的。可见，式（4-1）和（4-2）只是同一能量信号的两种不同表现形式。
尽管傅立叶变换可以关联信号的时频特征，从而分别从时域和频域对信号进
1
行分析，但却无法将两者有效地结合起来，因此傅立叶变换在信号的局部化分析方面存在严重不足。但在许多实际应用中，如地震信号分析、核医学图像信号分析等，研究者们往往需要了解某个局部时段上出现了哪个频率，或是某个频率出现在哪个时段上，即信号的时频局部化特征，傅立叶变换对于此类分析无能为力。
因此需要一种如下的数学工具：可以将信号的时域和频域结合起来构成信号的时频谱，描述和分析其时频联合特征，这就是所谓的时频局部化分析方法，即时频分析法。1964年，Gabor等人在傅立叶变换的基础上引入了一个时间局部化“窗函数”g(t)，改进了傅立叶变换的不足，形成窗口化傅立叶变换，又称“Gabor变换”。
g(t)在有限的区间外恒等于零或很快地趋于零，定义“窗函数”用函数g(t??)乘以f(t)，其效果等同于在t??附近打开一个窗口，即：
Gf??,????f?t?g?t???e?j?tdt (4-3)
????式（4-3）即为函数f(t)关于g(t)的Gabor变换。由定义可知，信号f(t)的Gabor变换可以反映该信号在t??附近的频谱特性。其逆变换公式为：
f?t???12??????d??ej?tg?t???Gf??,??d? (4-4)
可见Gf??,??的确包含了信号f(t)的全部信息，且Gabor窗口位置可以随着?的变化而平移，符合信号时频局部化分析的要求。
虽然Gabor变换一定程度上克服了傅立叶变换缺乏时频局部分析能力的不足，且兼顾信号的时频分辨率，但其本身仍存在不可克服的局限性，即Gabor窗口不具有自适应性，其大小是固定不变的，因此Gabor变换只能进行单一分辨率的分析。而在实际研究中，我们常常希望窗口的大小会随着频率的高低而改变，比如在研究高频信号时，希望窗口开得小一点；反之，在研究低频信号时，则希望窗口开得大一点，这样才更符合实际研究中低频信号分辨率比高频信号分辨率低的特点，因此需要研究更好的解决办法来改善Gabor变换的不足。
为了克服前面所描述的傅立叶变换和Gabor变换存在的不足，学者们提出了小波变换(Wavelet Transform)。
2
4.2 小波变换与逆变换
4.2.1 连续小波变换和离散小波变换
小波变换是一种窗El面积(即窗口大小)固定但窗口形状可变的时频局部化分析方法，其高频部分的时间分辨率较高而频率分辨率较低，而低频部分的频率分辨率较高而时间分辨率较低，因此对信号具有良好的自适应性，被冠以“数学显微镜”的美誉。小波变换又可分为连续小波变换和离散小波变换两种。
连续小波变换的概念是由Morlet等人提出[1]。设??x?为小波变换的核函数，若核函数??x??L2?R?若满足容许性条件 ：
C????F?????2R?d???? (4-5)
则称该函数?(x)为基小波。一维信号f(t)?L2(R)的连续小波变换可定义为：
?Wf??a,b???1a?R?t?b?f?x????dx (4-6)
?a?根据上述对基小波定义，可知：?R??x?dx??，且??x?在无穷远处趋近于零。一般的，记：
?a,b?x?????x?b?? a,b?R,a<0 (4-7) ?a?函数?a,b?x?是由基小波函数??x?经过尺度a的伸缩和b的平移之后所得。常用的连续小波包括：Morlet小波、Daubechies小波、三次样条小波、Meyer小波和Simoncelli小波等。
与连续小波变换相对应的是离散小波变换，其一般形式为：
?W?f??m,n???f,?m?x?nb???0其中??m,n?x??a02??m?a??0?m,n??a0?m2?R?x?nb0f?x????am?0??dx (4-8) ?????m,n?Z??为小波基，a0、b0为两个常量且????a0>0。离散小波最具代表性的为二进小波，即a0为2的幂次2J，b0取整。
3
由于小波变换在时域和频域同时兼有局部化能力，且能逐步聚焦到对象的任何细节进行分析，因此在人脸识别方面得到众多研究者的关注。 4.2.2 几种常见的小波
同傅立叶分析不同，小波分析的基（小波函数）不是唯一存在的，所有满足小波条件的函数都可以作为小波函数，那么小波函数的选取就成了十分重要的问题[8]。
(1) Haar小波
A.Haar于1990年提出一种正交函数系，定义如下：
?10?x?1/2??H(x)???11/2?x?1?0? 其它 （4-9）
这是一种最简单的正交小波，即
?????(t)?(x?n)dx?0 n??1,?2, （4-10）
(2) Daubechies(dbN)小波系
该小波是Daubechies从两尺度方程系数?hk?出发设计出来的离散正交小波。一般简写为dbN，N是小波的阶数。小波?和尺度函数阈中的支撑区为2N-1。?的消失矩为N。除N＝1外(Haar小波)，dbN不具对称性〔即非线性相位〕；dbN没有显式表达式(除N＝1外)。但
P(y)??CkN?1?kykk?0N?1?hk?的传递函数的模的平方有显式表达式。假
设
N?1?kCk，其中，为二项式的系数，则有
12N?1?ik?m0(?)?hke?2k?0其中
m0(?)?(cos22?2)NP(sin2?)2 （4-11）
(3) Biorthogonal(biorNr.Nd)小波系
Biorthogonal函数系的主要特征体现在具有线性相位性，它主要应用在信号与图像的重构中。通常的用法是采用一个函数进行分解，用另外一个小波函数
4
进行重构。Biorthogonal函数系通常表示为biorNr.Nd的形式：
Nr=1 Nd=1，3，5 Nr=2 Nd=2，4，6，8 Nr=3 Nd=1，3，5，7，9 Nr=4 Nd=4 Nr=5 Nd=5 Nr=6 Nd=8
其中，r表示重构，d表示分解。 (4) Coiflet(coifN)小波系
coiflet也是函数由Daubechies构造的一个小波函数，它具有coifN(N=1，2，3，4，5)这一系列，coiflet具有比dbN更好的对称性。从支撑长度的角度看，coifN具有和db3N及sym3N相同的支撑长度；从消失矩的数目来看，coifN具有和db2N及sym2N相同的消失矩数目。
(5) SymletsA(symN)小波系
Symlets函数系是由Daubechies提出的近似对称的小波函数，它是对db函数的一种改进。Symlets函数系通常表示为symN(N=2，3，?，8)的形式。
(6) Morlet(morl)小波
?x?(x)?CeMorlet函数定义为
2/2cos5x，它的尺度函数不存在，且不具有正
交性。
(7) Mexican Hat(mexh)小波 Mexican Hat函数为
?(x)?22?1/4?(1?x2)e?x/23 （4-12）
它是Gauss函数的二阶导数，因为它像墨西哥帽的截面，所以有时称这个函数为墨西哥帽函数。墨西哥帽函数在时间域与频率域都有很好的局部化，并且满足
?????(x)dx?0 （4-13)
由于它的尺度函数不存在，所以不具有正交性。 (8) Meyer函数
5