e.m

来自「基于小波的人脸识别算法的matlab实现」· M 代码 · 共 27 行

function Enum=E(WI1,WI2,width,height);
%计算能量E；
Enum=0;
for x1=1:width
    for x2=1:height
        %计算JI，JM；
        for v=1:5
            for u=1:8
                JI((v-1)*8+u)=WI1(v,u,x1,x2);
                JM((v-1)*8+u)=WI2(v,u,x1,x2);
            end
        end
        %Enum=Enum-JI*rot90(JM)/sqrt(JI*rot90(JI))/sqrt(JM*rot90(JM));
        JIM=0;JMM=0;
        for k=1:40
            JIM=JIM+JI(k)*JI(k);
        end
        for k=1:40
            JMM=JMM+JM(k)*JM(k);
        end
        JIM=sqrt(JIM);JMM=sqrt(JMM);
        for k=1:40
            Enum=Enum-JI(k)*JM(k)/JIM/JMM;%12个中错了2个；如果用平方差，错4个（1/3）；
        end
    end
end