bayes.m

来自「用Matlab编写的Bayes 分类器」· M 代码 · 共 15 行

function rcount=Bayes(M1,E1,D1,M2,E2,D2)
C=log(det(E1))+det(M1'*E1^(-1)*M1)-log(det(E2))-det(M2'*E2^(-1)*M2);
rcount=0;
for i=1:50
    A=D1(i,:)*(E1^(-1)-E2^(-1))*D1(i,:)';
    B=(E1^(-1)*M1-E2^(-1)*M2)'*D1(i,:)';
    if(-0.5*A+B-0.5*C)>=0
        rcount=rcount+1;
    end
    A=D2(i,:)*(E1^(-1)-E2^(-1))*D2(i,:)';
    B=(E1^(-1)*M1-E2^(-1)*M2)'*D2(i,:)';
    if(-0.5*A+B-0.5*C)<0
        rcount=rcount+1;
    end
end