《模式识别实验报告-实验一Bayes分类器设计说明》由会员分享，可在线阅读，更多相关《模式识别实验报告-实验一Bayes分类器设计说明（16页珍藏版）》请在金锄头文库上搜索。
1、实验一 Bayes分类器设计【实验目的】对模式识别有一个初步的理解， 能够根据自己的设计对贝叶斯决策理论算法有一个深刻 地认识，理解二类分类器的设计原理。【实验原理】最小风险贝叶斯决策可按下列步骤进行:在已知P( i)，P(X i)，i=1,，c及给出待识别的X的情况下，根据贝叶斯公式计算出后验概率：P( iX)P(X i)P( i)j=1,，XP(X i)P( i)j 1(2) 利用计算出的后验概率及决策表，按下面的公式计算出采取 aj,i=l,，a的条件风险cR(ai X)(ai , j)P( j X ) ,i=1,2,aj 1(3) 对中得到的a个条件风险值R(aX),i=1,，a进行比
2、较，找出使其条件风险最小的决策ak，即R ak x m|n R ai x则ak就是最小风险贝叶斯决策。【实验容】假定某个局部区域细胞识别中正常(1 )和非正常(2 )两类先验概率分别为正常状态：P (1) =0.9 ;异常状态：P (2 ) =0.1。现有一系列待观察的细胞，其观察值为X:-3.9847 -3.5549-1.2401 -0.9780-0.7932-2.8531-2.7605 -3.7287 -3.5414 -2.2692-3.4549-3.0752-3.99342.8792 -0.9780 0.79321.18823.0682-1.5799 -1.4885 -0.7431 -0
3、.4221-1.11864.2532已知类条件概率是的曲线如下图:P(x| 1)p(x| 2)类条件概率分布正态分布分别为N( -2 , 0.25 )、N( 2,4 )试对观察的结果进行分类。【实验要求】1) 用matlab完成基于最小错误率的贝叶斯分类器的设计，要求程序相应语句有说明文字，要求有子程序的调用过程。2) 根据例子画出后验概率的分布曲线以及分类的结果示意图。3) 如果是最小风险贝叶斯决策，决策表如下：最小风险贝叶斯决策表：、一状态 决策12a 104a 220请重新设计程序，完成基于最小风险的贝叶斯分类器，画出相应的条件风险的分布曲线和分类结果，并比较两个结果。【实验程序】最小错
4、误率贝叶斯决策分类器设计x=-3.9847-3.5549 -1.2401-0.9780-0.7932-2.8531-2.7605-3.7287-3.5414-2.2692-3.4549-3.0752-3.99342.8792-0.97800.79321.18823.0682-1.5799-1.4885-0.7431-0.4221-1.11864.2532 pw1=0.9 ; pw2=0.1e1=-2; a仁 0.5e2=2;a2=2m=numel(x) %得到待测细胞个数pw1_x=zeros(1,m) %存放对w1的后验概率矩阵pw2_x=zeros(1,m) % 存放对w2的后验概率矩阵r
5、esults=zeros(1,m) % 存放比较结果矩阵for i = 1:m%十算在w1下的后验概率pw1_x(i)=(pw1* normpdf(x(i),e1,a1)/(pw1* normpdf(x(i),e1,a1)+pw2* normpdf(x(i),e2,a2)%十算在w2下的后验概率pw2_x(i)=(pw2* normpdf(x(i),e2,a2)/(pw1* normpdf(x(i),e1,a1)+pw2* normpdf(x(i),e2 ,a2)endfor i = 1:mif pw1_x(i)pw2_x(i) %比较两类后验概率result(i)=0 %正常细胞elsere
6、sult(i)=1 %异常细胞endenda=-5:0.05:5 % 取样本点以画图n=nu mel(a)pw1_plot=zeros(1, n)pw2_plot=zeros(1, n)for j=1: npw1_plot(j)=(pw1* normpdf(a(j),e1,a1)/(pw1* normpdf(a(j),e1,a1)+pw2* no rmpdf(a(j),e2,a2)%十算每个样本点对 w1的后验概率以画图pw2_plot(j)=(pw2* no rmpdf(a(j),e2,a2)/(pw1* normpdf(a(j),e1,a1)+pw2* normpdf(a(j) ,e2,a
7、2)end figure(1) hold on plot(a,pw1_plot,k-,a,pw2_plot,r-.)for k=1:mif result(k)=0 plot(x(k),-0.1,b*) % else正常细胞用*表示异常细胞用五角星表示en d;en d;lege nd(正常细胞后验概率曲线xlabel(样本细胞的观察值)ylabel(后验概率)title(后验概率分布曲线)grid onreturn ;plot(x(k),-0.1,rp) %实验容仿真异常细胞后验概率曲线正常细胞,异常细胞)x = -3.9847 , -3.5549 , -1.2401 , -0.9780 ,
8、-0.7932 , -2.8531 ,-2.7605 ,-3.7287 , -3.5414, -2.2692 , -3.4549 , -3.0752, -3.9934 , 2.8792 ,-0.97800.7932 , 1.1882 , 3.0682,-1.5799 ,-1.4885 ,-0.7431 ,-0.4221 , -1.1186 ,4.2532 disp(x)pw1=0.9 pw2=0.1result=bayes(x,pw1,pw2)最小风险贝叶斯决策分类器设计fun ctio n R1_x,R2_x,result=da nger(x,pw1,pw2)m=numel(x) %得到待测
9、细胞个数R1_x=zeros(1,m)%存放把样本X判为正常细胞所造成的整体损失R2_x=zeros(1,m)%存放把样本X判为异常细胞所造成的整体损失result=zeros(1,m)%存放比较结果e1=-2a仁 0.5e2=2a2=2%类条件概率分布 px_w1: ( -2，0.25) px_w2 ( 2,4)r11=0 r12=2 r21=4 r22=0液险决策表for i=1:m%十算两类风险值R1_x(i)=r11*pw1* normpdf(x(i),e1,a1)/(pw1* normpdf(x(i),e1,a1)+pw2* normpdf(x(i),e2,a2)+r21*pw2*
10、normpdf(x(i),e2,a2)/(pw1* normpdf(x(i),e1,a1)+pw2* normpdf(x(i),e2 ,a2)R2_x(i)=r12*pw1* normpdf(x(i),e1,a1)/(pw1* normpdf(x(i),e1,a1)+pw2* normpdf(x(i),e2,a2)+r22*pw2* normpdf(x(i),e2,a2)/(pw1* normpdf(x(i),e1,a1)+pw2* normpdf(x(i),e2 ,a2)endfor i=1:mif R2_x(i)R1_x(i) result(i)=0 elseresult(i)=1end%
11、第二类比第一类风险大%判为正常细胞（损失较小），用0表示%判为异常细胞，用1表示enda=-5:0.05:5%取样本点以画图n=nu mel(a)R1_plot=zeros(1, n)R2_plot=zeros(1, n)for j=1: nR1_plot(j)=r11*pw1* normpdf(a(j),e1,a1)/(pw1* normpdf(a(j),e1,a1)+pw2* normpdf(a(j ),e2,a2)+r21*pw2* normpdf(a(j),e2,a2)/(pw1* normpdf(a(j),e1,a1)+pw2* normpdf(a(j) ,e2,a2)R2_plot
12、(j)=r12*pw1* normpdf(a(j),e1,a1)/(pw1* normpdf(a(j),e1,a1)+pw2* normpdf(a(j ),e2,a2)+r22*pw2* normpdf(a(j),e2,a2)/(pw1* normpdf(a(j),e1,a1)+pw2* normpdf(a(j) ,e2,a2)%十算各样本点的风险以画图endfigure(1) hold onplot(a,R1_plot,b-,a,R2_plot,g*-)for k=1:mif result(k)=0 plot(x(k),-0.1,bA) else%正常细胞用上三角表示end;en d;leg
13、e nd(正常细胞,异常细胞,Location,Best)xlabel(细胞分类结果)ylabel(条件风险)title(风险判决曲线)grid onreturn实验容仿真-0.79323.0682,-1.5799x = -3.9847, -3.5549 ,-1.2401,-0.9780-2.8531 ,-2.7605 , -3.7287 , -3.5414 , -2.2692 , -3.4549 , -3.0752 , -3.9934,2.8792, -0.9780,0.7932,1.1882-1.4885 , -0.7431 , -0.4221 , -1.1186 , 4.2532 dis
14、p(x)pw1=0.9pw2=0.1R1_x,R2_x,result=da nger(x,pw1,pw2)【实验结果和数据】最小错误率贝叶斯决策后验概率曲线与判决结果在一图上：后验概率曲线如图所示，带*的绿色曲线为判决成异常细胞的后验概率曲线；另一条平滑的蓝色曲线为判为正常细胞的后验概率曲线。根据最小错误概率准则，判决结果见曲线下方，其中“上三角”代表判决为正常细胞，“圆圈”代表异常细胞。各细胞分类结果：0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 0 1 0 10为判成正常细胞，1为判成异常细胞率概验后1.21基于最小错误率的贝叶斯判决图最小风险贝叶斯决策风险判决曲线如图2所示，其中带*的绿色曲线代表异常细胞的条件风险曲线；另一条光滑的蓝色曲线为判为正常细胞的条件风险曲线。根据贝叶斯最小风险判决准则，判决
