貝葉斯分類作業(yè)題

1、精選優(yōu)質文檔-傾情為你奉上作業(yè)：在下列條件下，求待定樣本x=(2,0)T的類別，畫出分界線，編程上機。1、二類協(xié)方差相等，2、二類協(xié)方差不等。訓練樣本號k1 2 31 2 3特征x11 1 2-1 -1 -2特征x21 0 -11 0 -1類別 1 21、 二類協(xié)方差不等Matlab程序如下：>> x1=mean(1,1,2),mean(1,0,-1)',x2=mean(-1,-1,-2),mean(1,0,-1)'x1 = 1.3333 0x2 = -1.3333 0>> m=cov(1,1;1,0;2,-1),n=cov(-1,1;-1,0;-2,-

2、1)m = 0.3333 -0.5000 -0.5000 1.0000n = 0.3333 0.5000 0.5000 1.0000>> m1=inv(m),n1=inv(n)m1 = 12.0000 6.0000 6.0000 4.0000n1 = 12.0000 -6.0000 -6.0000 4.0000>> p=log(det(m)/(det(n)p = 0>> q=log(1)q = 0>> x=2,0'x = 2 0>> g=0.5*(x-x1)'*m1*(x-x1)-0.5*(x-x2)'*n1*

3、(x-x2)+0.5*p-qg = -64（說明：g<0,則判定x=2,0T屬于1類）（化簡矩陣多項式0.5*(x-x1)'*m1*(x-x1)-0.5*(x-x2)'*n1*(x-x2)+0.5*p-q，其中x1,x2已知，x設為x= x1,x2T,化簡到(12x1-16+6x2）(x1-4/3)+(6x1-8+4x2)-(12x1+16-6x2)(x1+4/3)-(-6x1-8+4x2)x2，下面用matlab化簡，程序如下）>> syms x2;>> syms x1;>> w=(12*x1-16+6*x2)*(x1-4/3)+(6

4、*x1-8+4*x2)*x2-(12*x1+16-6*x2)*(x1+4/3)-(-6*x1-8+4*x2)*x2,simplify(w) w = (12*x1-16+6*x2)*(x1-4/3)+(6*x1-8+4*x2)*x2-(12*x1+16-6*x2)*(x1+4/3)-(-6*x1-8+4*x2)*x2 ans = -64*x1+24*x2*x1(說明：則-64×x1+24×x2×x1=0,即x1=0，或者x2=8/3,很顯然分界線方程為x1=0，因為x2=8/3連1類與2都分不開)2、二類協(xié)方差相等Matlab程序如下：>> l=m+nl

5、 = 0.6667 0 0 2.0000>> l1= inv(l)l1 = 1.5000 0 0 0.5000>> g1=(x2-x1)'*m1*x+0.5*(x1'*l1*x1-x2'*l1*x2)-qg1 = -64.0000（說明：g1<0,則判定x=2,0T屬于1類）>> (x2-x1)'*m1ans = -32.0000 -16.0000>> syms x11;>> syms x22;>> w1=-32*x11+(-16)*x22+0.5*(x1'*l1*x1-x2&

6、#39;*l1*x2)-q,simplify(w1) w1 = -32*x11-16*x22 ans = -32*x11-16*x22（說明：分界線方程為-32×x1-16×x2=0,即2×x1+x2=0）以下是matlab繪圖程序：>> x1=1;1;2; x2=1;0;-1;plot(x1,x2,'mx','markersize',15);axis(-5,5,-5,5);grid on;hold on>> x1=-1;-1;-2; x2=1;0;-1;plot(x1,x2,'m*','markersize',15);axis(-5,5,-5,5);hold on>> x1=2; x2=0;plot(x1,x2,'mp','markersize',15);axis(-5,5,-5,5);hold on>> x2=-5:0.02:5;x1=0;plot(x1,x2,'b');axis(-5,5,-5,5);>> x1=-5:0.02:5;x2=-2*x1;plot(x1,x2,'-.k&