- 基于MATLAB的人工智能模式识别
- 周润景 武立群 蔺雨露编著
- 2240字
- 2021-06-10 14:48:36
2.3 最小风险贝叶斯决策
2.3.1 最小风险贝叶斯决策理论
决策理论就是为了实现特定的目标,根据客观可能性,在一定信息和经验的基础上,借助一定的工具、技巧和方法,对影响未来目标实现的诸多因素进行准确的计算和判断优选后,对未来行动做出决定。在某些情况下,引入风险的概念,以求风险最小的决策更为合理,如对癌细胞的识别,因为识别的正确与否直接关系到患者的身体健康甚至生命。风险的概念常与损失相联系。损失函数用于计算当参数的真值和决策结果不一致时带来的损失,这种损失作为参数的真值和决策结果的函数,称为损失函数。而损失函数的期望值便称为风险函数。为了分析,引入损失函数λ(αi,ωj),i=1,2,…,a,j=1,2,…,m。这个函数表示当处于状态ωj时采取决策为αi所带来的损失。在决策理论中,常用决策表一目了然地表示各种情况下的决策损失,如表2-2所示。这是在已知先验概率P(ωj)及类条件概率密度P(X|ωj),j=1,2,…,m的条件下进行讨论的。
表2-2 贝叶斯决策表
![img](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_64.jpg?sign=1734454484-O3ejdpDfdPjTlP5ppgAKgkgwr91WtG6T-0-08fa9ce7703486c2a5ba644e60af3ae9)
根据贝叶斯公式,后验概率为
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_65.jpg?sign=1734454484-j3MyYWT1UDUEOx61qYklXmwY35B3gw78-0-318307d6877380569ec832d89934b51c)
(2-25)
引入“损失”的概念后,当考虑错判所造成的损失时,就不能只根据后验概率的大小做决策了,还必须考虑所采用的决策是否损失最小。对于给定的X,如果采取决策αi,λ可以在m个λ(αi,ωj)(i=1,2,…,m)当中任取一个,其相应概率为P(ωj|X)。因此采取决策αi时的条件期望损失R(αi/X)为
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_66.jpg?sign=1734454484-bN40xVHaa7DPqqTHtXd9oDXGcZGCje8A-0-03e93ab6f21503441d5cb9cabc2f7447)
(2-26)
在决策论中又把采取决策的αi的条件期望损失R(αi|X)称为条件风险。由于X是随机向量的观察值,对X的不同观察值,采取αi决策时,其条件风险的大小是不同的。所以究竟采取哪一种决策将随X的取值而定。决策α可以看成随机向量X的函数,记为α(X),这里定义期望风险R为
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_67.jpg?sign=1734454484-tzrXRL4vXxjyT34TAHaOP3L6LEVfkVHY-0-084fec6312f06ac447debf9e7e9332f3)
(2-27)
式中,dx是特征空间的体积元,积分在整个特征空间进行。期望风险R反映对整个特征空间所有X的取值都采取相应的决策α(X)所带来的平均风险;而条件风险R(αi|X)只是反映了对某一X的取值采取决策αi时所带来的风险。显然,需要采取一系列决策α(X)使期望风险R最小。在考虑错判带来的损失时,我们希望损失最小。如果在采取每一个决策或行动时,都使其风险最小,则对所有的X做出决策时,其期望风险也必然最小,这样的决策就是最小风险贝叶斯决策。
最小风险贝叶斯决策规则为:如果
R(αk|X)=minR(αi|X) i=1,2,…,a
(2-28)
则有α=α(k)(即采取决策αk)。对于实际问题,最小风险贝叶斯决策可按下列步骤实现。
(1)在已知P(ωj),P(X|ωj),j=1,2,…,m,并给出待识别的X的情况下,根据贝叶斯公式可以计算出后验概率
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_68.jpg?sign=1734454484-VUV1D1Xq4hibgBp3Dr1q7MVix7hh3NFw-0-d91edefd24b6886c5dff2f24e1646fdf)
(2-29)
(2)利用计算出的后验概率及决策表,按式(2-29)有
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_69.jpg?sign=1734454484-XBZPgxWZYnHJwonqElajCoJwBPEw54a8-0-1a53b7d98d15d288a9c83e7b2fd00607)
(2-30)
计算出αi(i=1,2,…,a)的条件风险R(αi|X)。
(3)对步骤(2)中得到的a个条件风险值R(αi|X),i=1,2,…,a,进行比较,找出使条件风险最小的决策αk,即
R(αk|X)=minP(αi|X),i=1,2,…,a
(2-31)
则αk就是最小风险贝叶斯决策。
说明
最小风险贝叶斯决策除了要满足符合实际情况的先验概率P(ωj)及类条件概率密度P(X|ωj),j=1,2,…,m的条件,还必须有适合的损失函数λ(αi,ωj),i=1,2,…,a,j=1,2,…,m。实际工作中要列出合适的决策表很不容易,往往要根据研究的具体问题,分析错误决策所造成损失的严重程度,并与有关专家共同商讨来确定。
2.3.2 最小错误率与最小风险贝叶斯决策的比较
最小错误率贝叶斯决策规则与最小风险贝叶斯决策规则有着某种联系。这里再讨论一下两者的关系。设损失函数为
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_70.jpg?sign=1734454484-gtSTanjkNRXoOMwBnJ5oppzwPrg1rvyU-0-23861a9851a31d9d49f2f760dbe374a3)
(2-32)
式(2-32)中假设对m类只有m个决策,即不考虑“拒绝”的情况,对正确决策(即i=j)来说,λ(αi,ωj)=0,就是没有损失;对任何错误决策来说,其损失为1。这样定义的损失函数称为0-1损失函数。此时条件风险为
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_71.jpg?sign=1734454484-c69F2GHEd94YNmjiGrhelNhhOprjhQ7Y-0-3bd5496a77210a911a553fe0ca52c704)
(2-33)
式中,表示对X采取决策ωj(此时ωj就相当于αi)的条件错误概率。所以在采取0-1损失函数时,使
R(αk|X)=minP(αi|X) i=1,2,…,a
(2-34)
最小风险贝叶斯决策这时就等价于式(2-35)的最小错误率贝叶斯决策。
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_73.jpg?sign=1734454484-0W2J8c7IXL3V3wrWWgUjzd0cgKvKX7Nf-0-63fd475df0b056824d32f83602dfd934)
(2-35)
由此可见,最小错误率贝叶斯决策就是采用0-1损失函数条件下的最小风险贝叶斯决策,即前者是后者的特例。
2.3.3 贝叶斯算法的计算过程
(1)输入类数M、特征数n、待分样本数m。
(2)输入训练样本数N和训练集矩阵X(N×n),并计算有关参数。
(3)计算待分析样本的后验概率。
(4)若按最小风险原则分类,则输入各值,并计算各样本属于各类时的风险并判定各样本类别。
2.3.4 最小风险贝叶斯分类的MATLAB实现
1.初始化
初始化程序如下:
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_74.jpg?sign=1734454484-xInpwhOS2PsuRdMPUrVfI9RrgIajU2XZ-0-cbea894463e6d14f6064e7b08383186f)
2.参数计算
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_75.jpg?sign=1734454484-iDmKVneJBfC9gxFBfwTTSuIaIqIelHJ9-0-c9340c30e7f11ef5e87c08d5e72222aa)
3.完整程序及仿真结果
程序代码如下:
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_76.jpg?sign=1734454484-wSvaZL4XIz5qwaOwW7ufB5uixx70nm96-0-b719a447c49d24bce78e1d2a4ad1c1ca)
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_77.jpg?sign=1734454484-g8E721x9RESxO75Ug66AogZHSUKciiXz-0-8527979b0cb998705d6bcdf0d21ee2ed)
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_78.jpg?sign=1734454484-DibKtRfP2OgZG58Xbjx1VY9hPuux9MdN-0-81ec55bc3a3ca19f61116e6c5876d490)
运行程序得到的损失函数的矩阵为:
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_79.jpg?sign=1734454484-ZsuclAHBaXm7bfV1mk898QnQl8WChxGp-0-c809e9af5f34ff1910e8f50928d00059)
得到贝叶斯决策表,如表2-3所示。
表2-3 贝叶斯决策表
![img](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_80.jpg?sign=1734454484-uVjrPvTQvDOIWhnuE4UbQT2jEGmHiJXs-0-0e69cbb89475b107166365900dc5261f)
损失函数图如图2-7所示。
![img](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_81.jpg?sign=1734454484-yN6a4BO27ztAR7ZvDYQcPnGDI9zGFCsf-0-3b204298551eac08898da12d0c98d3f0)
图2-7 损失函数图
继续运行程序,MATLAB命令窗口显示的结果如下:
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_82.jpg?sign=1734454484-xprRGX3b6KvKGy2luXXpwqdiy1252jwK-0-77fa853a1fdda60eb3b399c7e4fc05ca)
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_83.jpg?sign=1734454484-Usf5CifGhF7AtYrZHvOWGkr87zoupe7a-0-10fb02acca46e60a1c7344c2373af9db)
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_84.jpg?sign=1734454484-JlDWrhKGeTUppOl6Zj4HKso1b9krsLPO-0-2ac86ffc02c14b64621ed2f2a6585bbb)
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_85.jpg?sign=1734454484-LblEAUY9DoPDNfqZe8bJnDwTwwzI6yvh-0-f6d269a609599869d817d987fda97cbe)
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_86.jpg?sign=1734454484-NaPIQvRGcL4JT4w0qEHDRa621ewdaYDs-0-4b8c0ba105a56f45f5df7c7bfdee2cb7)
运行m个文件程序,将得到的结果整理成表,如表2-4所示。
表2-4 待测样本分类表
![img](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_87.jpg?sign=1734454484-1THWFI57Mv6uiEX0iUyGKBqoNNj23cx8-0-288956494a670695bb6420206172402b)
续表
![img](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_88.jpg?sign=1734454484-oBlhhhb9Wfi7wFqfZjOSe5tVmodZLiSV-0-29a258aa5ee12a7d79f5de1d5d0f9153)
该分类结果与标准分类存在一个数据类别不同的差异:即数据“1494.63,2072.59,2550.51”。用上述方法得到的分类结果中该数据属于类别3,而在标准分类中该数据属于类别1。出现这种情况的原因可能是基于最小风险的贝叶斯分类的分类方法存在误差。
反过来验证分类结果的正确性。
修改MATLAB循环语句的循环次数与后验概率的输入向量:
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_89.jpg?sign=1734454484-E9x5OOumWQVWIBjKtpOuB5VyXCvE40UR-0-85270fb8178ce06fa6d6a6da47e3863d)
得到MATLAB命令窗显示的结果,部分结果如下:
![](https://epubservercos.yuewen.com/5B96F7/20205398108552606/epubprivate/OEBPS/Images/txt002_90.jpg?sign=1734454484-TS7xnEbv3HqifFEombP2Hoc8lFShwMq7-0-fd631301081834dbfe942faac82813e0)
可以看出,其与训练样本的分类结果是完全吻合的。
总结
以贝叶斯决策理论为核心内容的统计决策理论是统计模式识别的重要基础,理论上该分类理论有最优性能:即分类错误或风险在所有分类器中是最小的,常可以作为衡量其他分类器设计方法优劣的标准。
但是该方法明显的局限性体现为:该方法需要已知类别数,以及各类别的先验概率和类条件概率密度,也就是说要分两步来解决模式识别问题——先根据训练样本设计分类器,再对测试样本进行分类。因此有必要研究直接从测试样本出发设计分类器的其他方法。
思考与练习
(1)什么是最小错误率贝叶斯决策?
(2)什么是最小风险贝叶斯决策?
(3)最小错误率贝叶斯决策与最小风险贝叶斯决策的区别是什么?