meanclass.m

这是人脸识别代码

function [XMean,XMeanClass]= MeanClass(X,XClass)

% 输入样本按列堆积的矩阵X，以及每一列的所属的类别向量XClass
% 计算每一类的均值XMean以及XMean(:,i)的类别XMeanClass(i)。

XCount = size(X,2); % 样本数

% 初始化
ClassCount = 1; % 类别数
XMeanClass(1) = XClass(1); % 第1类的类别
XMean(:,1) = X(:,1); % 第1类的样本总和
XMeanCount(1) = 1; % 第1类的样本数量

for i=2:XCount
    if length(strfind(XMeanClass,XClass(i))) == 0 % 如果第i个样本的类别不在已有的类别中
        ClassCount = ClassCount + 1;
        XMeanClass(ClassCount) = XClass(i); % 第ClassCount类的类别
        XMean(:,ClassCount) = X(:,i); % 第ClassCount类的样本总和
        XMeanCount(ClassCount) = 1; % 第ClassCount类的样本数量
    else
        XMean(:,ClassCount) = XMean(:,ClassCount) + X(:,i);  % 第ClassCount类的样本总和
        XMeanCount(ClassCount) = XMeanCount(ClassCount) + 1;  % 第ClassCount类的样本数量
    end
end;

% 每一类分别除去各自的数量得到每一类的均值
for i=1:ClassCount
    XMean(:,i) = XMean(:,i) / XMeanCount(i);
end;