www.pudn.com > facerecognitionmatlab.rar > E.M

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