1.一种基于混合误差编码的人脸识别方法，其特征在于：所述人脸识别方法包括以下步骤：步骤1将实数域上的m×n维的人脸图像y作为待识别人脸图像，将来自于K个人的N个不含遮挡的且在正常光照条件下采集到的人脸图像集作为训练样本，这里， 为实数域上的m×n阶矩阵，表示第k个人的第i个训练样本，k＝1，2…K， Nk为第k个人的人脸图像的样本数；

步骤2将y和A分别变换到梯度方向域中：这里，*表示二维卷积，arctan(·)为反正切函数，hd为差分滤波器的系数矩阵，符号T表示矩阵的转置；

步骤3将待识别图像y和训练集A变换到高维特征空间中：其中， 表示待识别图像y在高维特征空间中的表示， 表示训练集A在高维特征空间中的表示，V(·)为拉伸变换，即：将m×n阶矩阵拉伸为D＝m×n维列向量，或者将所有的m×n阶矩阵的集合A中的元素依次拉伸为D＝m×n维列向量，然后再将拉伸后的列向量组装成矩阵，即：其中，

步骤4初始化待识别人脸图像y相对于训练集A的编码系数为 这里，

1N×1表示N维列向量，其全部元素为1；

步骤5初始化参数：t＝1，ζ*＝+∞，这里，t表示迭代次数，ζ*表示历次迭代过程中产生的判别误差的最小质量；

步骤6在像素空间中，按照如下公式，计算待识别图像y的重构图像其中，Mm×n(·)是V(·)的逆运算，即把D＝m×n维列向量还原为m×n阶矩阵，x(t)为实数域上的N维列向量，表示在第t次迭代时所产生的待识别人脸图像y相对于训练集A的编码系数；V(A)为集合A的拉伸变换：步骤7将 变换到梯度方向域中：

(t)

步骤8计算待识别图像y与其重构图像 的梯度方向之差△θ ：(t)

步骤9将△θ 在cos域中做如下变换：步骤10对 进行均值滤波： 这里，hm为均值滤波器的系数；

步骤11将 做如下变换，得到结构化误差： 这里，D＝m×n为 的维数，|·|表示取绝对值，||·||1表示l1范数；

步骤12对结构化误差 进行阈值聚类，得到待识别图像y的有效特征支撑的初步估τ(t)计s ：

τ(t) τ(t)

其中，κτ(·)表示阈值滤波算子，s 和 的下标i∈{1，2，…，m×n}表示s和 中元素的索引，τ为给定的阈值；

(t)

步骤13通过如下公式，进一步估计待识别图像y的有效特征支撑s ：(t) τ(t)

其中，s表示优化式中的与s 相对应的有效特征支撑的变元，s和s 的维数相同，下标i，j∈{1，2，…，m×n}表示s和sτ(t)中元素的索引， 表示i的邻域节点，λs为光滑参数，λu为数据参数；

(t)

步骤14由待识别图像y的有效特征支撑s ，计算y在高维特征空间中的表示 相对(t) (t)于训练集A在高维特征空间中的表示 的稀疏编码系数x ，以及与x 相关联的判别误差其中，s.t.表示约束条件，||·||2表示l2范数， 表示点乘， x表示(t)优化式中的与x 相对应的稀疏编码系数的变元， 表示优化式中的与 相对应的判别误差的变元；

步骤15度量判别误差 的质量ζ(t)： 如果ζ(t)＜ζ*，则令ζ*＝ζ(t)，x*＝x(t)， 这里，ζ*，x*， 分别记录ζ(t)，x(t)， 的最优表示；

步骤16令t＝t+1，迭代上述步骤6～步骤15，直至给定的最大迭代次数；

步骤17在高维特征空间中，按如下公式，计算待识别图像y在训练集中的每一个人的样本空间中的投影残差rk：* *

其中，k∈{1，2，…，K}，δk(x)表示将x中的所有除了对应于第k个人的编码系数外的元素全部置为0；

步骤18按如下公式对待识别图像y进行分类：