第四章作业统计判别.doc

文档介绍：第四章作业
:
设,求该两类模式之间的贝叶斯判别界面的方程式。
绘出判别界面。
解:(a),又因,
故
故判别界面为
(b) 判别界面:
。
function d12=Bayes(W1,W2,pW1,pW2)
m1=mean(W1,1)';
m2=mean(W2,1)';
n=length(W1(1,:));
N1=length(W1(:,1));
N2=length(W2(:,1));
C1=zeros(n,n);
for i=1:N1
C1=C1+W1(i,:)'*W1(i,:);
end;
C1=C1/N1;
C1=C1-m1*m1';
C2=zeros(n,n);
for i=1:N2
C2=C2+W2(i,:)'*W2(i,:);
end
C2=C2/N2;
C2=C2-m2*m2';
syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 real ;%假设每个模式最多12维
X=[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12];
X(n+1:12)=[];
X=X';
d1=log(pW1)-*log(abs(det(C1)))-*((X-m1)'*inv(C1)*(X-m1));
d2=log(pW2)-*log(abs(det(C2)))-*((X-m2)'*inv(C2)*(X-m2));
d12=d1-d2;
d12=simple(d12)