基于最小错误率的贝叶斯决策

a=input('请输入样本值:'); [x,xclass]=ReadIRIS; u=zeros(4,3); for i=1:3; m=(i-1)*50+1; m1=m+49; x1=0; x2=0; x3=0; x4=0; for k=m:m1; a1=x(1,k);x1=x1+a1; a2=x(2,k);x2=x2+a2; a3=x(3,k);x3=x3+a3; a4=x(4,k);x4=x4+a4; end m(1:i)=x1/50; m(2:i)=x2/50; m(3:i)=x3/50; m(4:i)=x4/50; end for i=1:3; m=(i-1)*50+1; m1=m+49; x25=zeros(4,4); for j=m:m1; x22=x(:,j)-u(:,i); x23=x22'; x24=x22*x23; x25=x24+x25; end x25=x25./50; if i==1; c1=x25; elseif i==2; c2=x25; else c3=x25; end end c11=sqrt(det(c1)); c22=sqrt(det(c2)); c33=sqrt(det(c3)); p=(2*pi).^2; a=a'; b1=a-u(:,1); b2=a-u(:,2); b3=a-u(:,3); e1=-0.5.*b1'*(inv(c1))*b1; e2=-0.5.*b2'*(inv(c2))*b2; e3=-0.5.*b3'*(inv(c3))*b3; p1=exp(e1)./(p.*c11); p2=exp(e2)./(p.*c22); p3=exp(e3)./(p.*c33); p=p1+p2+p3; h=zeros(1,3); h(1)=p1./p; h(2)=p2./p; h(3)=p3./p; g=h(1); m=1; for i=1:3; if h(i)>g g=h(i); m=i; end end if m==1; lei='Iris-setosa'; elseif m==2; lei='Iris-versicolor'; else lei='Iris-virginica'; end disp(['所输入属第',num2str(m),'即'，lei]); disp(['属于该类概率为：'，num2str(g)]);