模式识别实验报告(一二)

信息与通信工程学院模式识别实验报告班级：姓名：学号：日期：2011年12月实验一、Bayes分类器设计一、实验目的：1.对模式识别有一个初步的理解2.能够根据自己的设计对贝叶斯决策理论算法有一个深刻地认识3.理解二类分类器的设计原理二、实验条件：matlab软件三、实验原理：最小风险贝叶斯决策可按下列步骤进行：1)在已知)(iP，)(iXP，i=1,…，c及给出待识别的X的情况下，根据贝叶斯公式计算出后验概率：cjiiiiiPXPPXPXP1)()()()()(j=1,…，x2)利用计算出的后验概率及决策表，按下面的公式计算出采取ia,i=1,…，a的条件风险cjjjiiXPaXaR1)(),()(,i=1,2,…,a3)对(2)中得到的a个条件风险值)(XaRi,i=1,…，a进行比较，找出使其条件风险最小的决策ka，即1,minkiiaRaxRax则ka就是最小风险贝叶斯决策。四、实验内容假定某个局部区域细胞识别中正常（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.97800.79321.18823.0682-1.5799-1.4885-0.7431-0.4221-1.11864.2532已知先验概率是的曲线如下图：-6-4-2024600.10.20.30.40.50.60.70.8)|(1xp)|(2xp类条件概率分布正态分布分别为（-2，0.25）（2,4）试对观察的结果进行分类。五、实验步骤：1.用matlab完成分类器的设计，说明文字程序相应语句，子程序有调用过程。2.根据例子画出后验概率的分布曲线以及分类的结果示意图。3.最小风险贝叶斯决策，决策表如下：状态决策12α106α210重新设计程序，完成基于最小风险的贝叶斯分类器，画出相应的后验概率的分布曲线和分类结果,并比较两个结果。六、实验代码1.最小错误率贝叶斯决策(m1.m)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.48850.7431-0.4221-1.11864.2532]pw1=0.9;pw2=0.1;e1=-2;a1=0.5;e2=2;a2=2;m=numel(x);%得到待测细胞个数pw1_x=zeros(1,m);%存放对w1的后验概率矩阵pw2_x=zeros(1,m);%存放对w2的后验概率矩阵results=zeros(1,m);%存放比较结果矩阵fori=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));endfori=1:mifpw1_x(i)pw2_x(i)%比较两类后验概率result(i)=0;%正常细胞elseresult(i)=1;%异常细胞endenda=[-5:0.05:5];%取样本点以画图n=numel(a);pw1_plot=zeros(1,n);pw2_plot=zeros(1,n);forj=1:npw1_plot(j)=(pw1*normpdf(a(j),e1,a1))/(pw1*normpdf(a(j),e1,a1)+pw2*normpdf(a(j),e2,a2));%计算每个样本点对w1的后验概率以画图pw2_plot(j)=(pw2*normpdf(a(j),e2,a2))/(pw1*normpdf(a(j),e1,a1)+pw2*normpdf(a(j),e2,a2));endfigure(1);holdonplot(a,pw1_plot,'co',a,pw2_plot,'r-.');fork=1:mifresult(k)==0plot(x(k),-0.1,'cp');%正常细胞用五角星表示elseplot(x(k),-0.1,'r*');%异常细胞用*表示end;end;legend('正常细胞后验概率曲线','异常细胞后验概率曲线','正常细胞','异常细胞');xlabel('样本细胞的观察值');ylabel('后验概率');title('后验概率分布曲线');gridonreturn%实验内容仿真：x=[-3.9847,-3.5549,-1.2401,-0.9780,-0.7932,-2.8531,-2.7605,-3.7287,-3.5414,-2.2692,-3.4549,-3.075,-3.9934,2.8792,-0.9780,0.7932,1.1882,3.0682,-1.5799,-1.4885,-0.7431,-0.4221,-1.1186,4.2532]disp(x);pw1=0.9;pw2=0.1;[result]=bayes(x,pw1,pw2);2.最小风险贝叶斯决策(m2.m)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.48850.7431-0.4221-1.11864.2532]pw1=0.9;pw2=0.1;m=numel(x);%得到待测细胞个数R1_x=zeros(1,m);%存放把样本X判为正常细胞所造成的整体损失R2_x=zeros(1,m);%存放把样本X判为异常细胞所造成的整体损失result=zeros(1,m);%存放比较结果e1=-2;a1=0.5;e2=2;a2=2;%类条件概率分布px_w1:（-2，0.25）px_w2（2,4）r11=0;r12=2;r21=4;r22=0;%风险决策表fori=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*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));endfori=1:mifR2_x(i)R1_x(i)%第二类比第一类风险大result(i)=0;%判为正常细胞（损失较小），用0表示elseresult(i)=1;%判为异常细胞，用1表示endenda=[-5:0.05:5];%取样本点以画图n=numel(a);R1_plot=zeros(1,n);R2_plot=zeros(1,n);forj=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(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);holdonplot(a,R1_plot,'co',a,R2_plot,'r-.');fork=1:mifresult(k)==0plot(x(k),-0.1,'cp');%正常细胞用五角星表示elseplot(x(k),-0.1,'r*');%异常细胞用*表示end;end;legend('正常细胞','异常细胞','Location','Best');xlabel('细胞分类结果');ylabel('条件风险');title('风险判决曲线');gridonreturn%实验内容仿真：x=[-3.9847,-3.5549,-1.2401,-0.9780,-0.7932,-2.8531,-2.7605,-3.7287,-3.5414,-2.2692,-3.4549,-3.075,-3.9934,2.8792,-0.9780,0.7932,1.1882,3.0682,-1.5799,-1.4885,-0.7431,-0.4221,-1.1186,4.2532]disp(x);pw1=0.9;pw2=0.1;[result]=bayes(x,pw1,pw2);七、实验结果1.最小错误率贝叶斯决策后验概率曲线与判决显示在上图中后验概率曲线：带红色虚线曲线是判决为异常细胞的后验概率曲线青色实线曲线是为判为正常细胞的后验概率曲线根据最小错误概率准则，判决结果显示在曲线下方：五角星代表判决为正常细胞，*号代表异常细胞各细胞分类结果（0为判成正常细胞，1为判成异常细胞）：0000000000000101110001012.最小风险贝叶斯决策风险判决曲线如上图所示：带红色虚线曲线是异常细胞的条件风险曲线；青色圆圈曲线是正常细胞的条件风险曲线根据贝叶斯最小风险判决准则，判决结果显示在曲线下方：五角星代表判决为正常细胞，*号代表异常细胞各细胞分类结果（0为判成正常细胞，1为判成异常细胞）：100000000000110111000101八、实验分析由最小错误率的贝叶斯判决和基于最小风险的贝叶斯判决得出的图形中的分类结果可看出，样本-3.9934、-3.9847在前者中被分为“正常细胞”，在后者中被分为“异常细胞”，分类结果完全相反。分析可知在最小风险的贝叶斯判决中，影响结果的因素多了一个“损失”。在第一张图中，这两个样本点下两类决策的后验概率相差很小，当结合最小风险贝叶斯决策表进行计算时，“损失”起了主导作用，导致了相反的结果的出现。同时，最小错误率贝叶斯决策就是在0-1损失函数条件下的最小风险贝叶斯决策，即前者是后者的特例。九、实验心得通过本次实验，我对模式识别有了一个初步的理解，开始对模式识别的相关知识从书本上转移到了实践中，并跟据自己的设计对贝叶斯决策理论算法有一个深刻地认识，同时也理解二类分类器的设计原理。同时，之前只学过浅显的Matlab知识，用Matlab实现数值计算的能力又一次得到了训练，对以后的学习和实验都有极大的帮助。实验二、基于Fisher准则线性分类器设计一、实验目的：1.进一步了解分类器的设计概念2.能够根据自己的设计对线性分类器有更深刻地认识3.理解Fisher准则方法确定最佳线性分界面方法的原理及Lagrande乘子求解的原理二、实验条件：matlab软件三、实验原理：线性判别