内容发布更新时间 : 2024/11/8 15:30:39星期一 下面是文章的全部内容请认真阅读。

河 北 大 学 工 商 学 院

Y=sim(net,p);%仿真模拟lvq网 Yc=vec2ind(Y);

test=nnz(Yc-Tc);%测试lvq网

pp=[1.24 1.28 1.40;1.80 1.84 2.04];%等待分类蠓虫的数据 y=sim(net,pp);%仿真模拟lvq网对蠓虫进行分类 index=vec2ind(y)

程序运行后，从命令窗口得到运行结果，

TRAINR, Epoch 0/1000 TRAINR, Epoch 3/1000

TRAINR, Performance goal met. index =

2 2 2

因此，通过LVQ网进行分类，所得三个蠓虫的类属为2，即依然为 Apf 类。

4.42 Bayes判别法

(1) 判别总体的协方差矩阵是否相等 (2) 判断总体是否服从正态分布

首先对每个指标进行一元正态分布的检验，若有一个指标不服从正态分布，则总体不服从正态分布；若每个指标都服从一元正态分布，且各指标不相关则总体服从正态分布。

(3) 利用按比例分配方法估计两个总体的先验概率： p1=6/(6+9)=0.4; p2=9/(6+9)=0.6; (4) 编写程序Bayes.m

apf=[1.14,1.78; 1.18,1.96;1.20,1.86;1.26,2.;1.28,2;1.30,1.96]; af=[1.24,1.72;1.36,1.74;1.38,1.64;1.38,1.82;1.38,1.90; 1.40,1.70;1.48,1.82;1.54,1.82;1.56,2.08]; x= [1.24,1.80;1.28,1.84; 1.4,2.04]; p1=0.4; p2=0.6;

m1=mean(apf);m2=mean(af);s1=cov(apf);s2=cov(af); S=(5*s1+8*s2)/13; for i=1:3,

d(i,:)=log(0.4)-0.5*((x(i,:)-m1)*inv(S)* (x(i,:)-m1)')- (log(0.6)-0.5*((x(i,:)-m2)*inv(S)* (x(i,:)-m2)')); end d

(5) 得到运行结果： d =

第10页

河 北 大 学 工 商 学 院

1.7585 0.9514

1.5747

由于d的结果全部大于零，即p1f1>p2f2，因此，通过Bayes判别法进行分类，三个蠓虫均属于Apf类。

评价：对于以上具体描述的判别方法，每一种方法都存在一个错判概率的问题，至于具体问题哪一种判别方法好，错判概率是一个指标，它应该尽量的小，但是最终的判别结果往往需要综合考虑。

5.总结体会

这次实验是一次较为困难却有意义的实践，我用了将近半个月的时间来完整的解决这个问题。起初并不知道从何下手，后来经过查阅资料了解到这类问题可以有多种解决方案。于是我首先采用了分界线方法进行求解，因为这种方法容易上手，在平常的课堂实验中就有所涉及，但是调试程序的时候还是遇到了些小困难，不过很快便解决了。对于实验中第三个问题该怎样表述，我想了很久，最后只得以例子的形式表达出来，不过自我感觉也还是能自圆其说的。

在查阅了关于这种分类问题的相关资料后，同时了解到人工智能网中的LVQ方法以及Bayes判别法，于是专门去图书馆学习了一下相关知识，发现用比较基础的知识就可以将这一问题解决，于是加以运用，便得到了上文提到的lvq.m及Bayes.m文件，得到的运行结果与使用分界线方法是一致的。

通过这次实验，我更深刻的了解了Matlab软件的强大功能以及数学应用的广泛性，同时也增强了我分析问题的能力。另外，通过这学期的学习与实践，我学 到了一些东西，也体会到了专业知识的运用，在这里要感谢邢老师的认真指导。

第11页