多源传感器检测融合算法及仿真实例(基于Bayes理论)

一、背景介绍
1.二元信号统计检测的基本模型
统计信号处理的理论基础之一是信号的统计检测。信号的统计检测理论主要研究在受噪声干扰的随机信号中信号的有/无？或信号属于哪个状态的最佳判决的概念、方法和性能等问题，其数学基础就是统计判决理论，又称假设检验理论。统计检测理论的基本模型是二元信号检测问题：

• 信源
  检测的信号发生源。

• 概率转移机构
  它是在信源输出其中一个为真的基础上，把噪声干扰背景中的假设$H_j(j=0,1)$为真的信号以一定概率关系映射到观测空间。
  如果没有噪声干扰，信源输出的某一种确知信号将映射到观测空间中的某一确定点，但在噪声干扰的情况下，它将以一定的概率映射到整个观测空间，

• 观测空间
  观测空间R 是在信源输出不同信号状态下，在噪声干扰背景中，由概率转移机构所生成的全部可能观测量的集合

• 判决规则
  观测量落入观测空间后，就可以用来推断哪一个假设成立是合理的，即判决信号属于哪种状态。为此需要建立一个判决规则，以便使观测空间中的每一个点对应着一个相应的假设$H_i(i＝0,1)$。判决结果就是选择假设H0 成立，还是选择假设H1 成立。

统计判决，即统计假设检验的任务，是根据观测量落在观测空间中的位置，按照某种检验规则，作出信号状态是属于哪个假设的判决。因此，统计判决问题实际上是对观测空间$R$的划分问题。

在二元信号检测中，是把整个观测空间$R$划分为$R_0$和$R_1$ 两个子空间(成为判决域)，并满足：
$$R=R_0\cup R_1\Phi =R_0\cap R_1$$
如果观测空间$R$中的某个观测量落入$R_0$域，就判决$H_0$成立，否则就判决假设$H_1$成立。

• 信号检测
  用于判决的准则可以是多种多样的，但信号检测的基本工具却大多利用了数理统计的假设检测方法。

在目标检测中的决策任务实际上属于二元假设检验问题。所谓的假设检验问题是指给定一组假设$H_i，i＝0，1，\dots ，m-1$，通过对已有的数据集$y$ 进行处理，来确定当前哪一个假设$H_j$是成立的。

如何根据已有的观测数据$z$得到最优决策结果，属于最优决策策略问题。
在目标检测中的二元假设检验，例如在监视区域里检测目标的存在与否（传感器网、雷达等），如果假设$H_0$ 表示无目标，$H_1$ 表示有目标，$P(H_0)＝P_0$和$P(H_1)＝P_1$ 分别为$H_0$ 和$H_1$出现的先验概率，且$P_0+P_1＝1$。
单个传感器的决策值u为二元值，一般定义为：
$$u=\{\begin{array}{l}0假设H_0(判定无目标)\\ 现象z(传感器的实际测量值)\\ 1假设H_1(判定有目标)\end{array}$$

二元假设检验问题有两种可能的错误：
$H_0$为真、判决$u=1$为第一类错误，即虚警
$H_1$为真，判决$u=0$为第二类错误，即漏检

$$\begin{array}{rl}& 检测概率：P_d=P(u=1|H_1)=1-P_m\\ & 虚警概率：P_f=P(u=1|H_0)\\ & 漏检概率：P_m=P(u=0|H_1)\end{array}$$

对于二元假设检验问题，数据集$z$空间可划分为两个区域：$R_0、R_1。$如果数据点$z$位于$R_0$区，则认为假设$H_0$成立而作出决策$u=0$；反之若数据点$z$位于$R_1$区，则认为假设$H_1$成立而做出决策$u=1$

信号检测的数学描述
假设$p_0(z)$和$p_1(z)$分别为$z$关于$H_0$和$H_1$的条件概率密度函数，下图中斜线部分红色和天蓝表示了第一类错误概率$P_f$和第二类错误概率$P_m$。

多源检测融合的系统架构

分布式数据融合是每个传感器都先对原始观测数据进行初步分析处理，做出本地判决结论，只把这种本地判决结论及其有关信息，或经初步分析认定可能存在某种结论但又不完全可靠的结论及其有关信息，向数据融合中心呈报；
$$\begin{array}{rl}& u_i=0/1目标有无判定信息观测z_i\\ & u_0=\{\begin{array}{l}0假设H_0(判定无目标)\\ U=(u_1,u_2,u_3,...,u_N,)\\ 1假设H_1(判定有目标)\end{array}\end{array}$$

二.基于似然函数比的融合准则

统计分析中的似然函数($likelihood$):是指在某种非随机参数条件下观测到数据$z$的概率。似然函数就是条件概率密度函数$p(z|\theta )$.
单个传感器检测中，在某种假设下观测信号的条件概率函数为似然函数，即$p(z|H_i)$
融合中心观测数据时各个传感器处理结果，有:多元融合检测中是在某种假设下各种判定输出的概率函数为似然函数，即:$p(u|H_i)$。其中$u=(u_1,u_2,...,u_N)$
注意：这里$u_i$ 是各个传感器各自处理后的数据，直接的测量数据为$z_i$

1.最大后验概率融合检测准则
最大后验概率准则是根据已有的观测数据$z$，分析给定全局决策值$u$下各种假设（可能）的概率，选择最优可能产生该全局决策值的假设.
在二元假设目标检测中就是取两个概率中较大者所对应的假设。即若$P(H_1|u)>P(H_0|u)$,则选择$H_1$;反之选择$H_0$.
上述规则可写为：
$$\frac{P(H_1|u)}{P(H_0|u)}>1?H_1:H_0$$
应用贝叶斯法则得：$$\frac{P(H_1/u)}{P(H_0/u)}=\frac{P(u/H_1)P(H_1)}{P(u/H_0)P(H_0)}>1?H_1:H_0$$
可得似然函数比：
$$\lambda (u)=\frac{P(u/H_1)}{P(u/H_0)}>\frac{P(H_0)}{P(H_1)}?H_1:H_0$$

N个传感器的最大后验概率融合检测
最大后验概率检测准则是选择似然函数比大于$P_0/P_1$ 的数据为$R_1$ 数据空间进行检测。

似然函数比的计算：
假设 u = { u 1 , u 2 , . . . , u N } u=\lbrace u_{1},u_{2},...,u_{N}\rbrace u={u1​,u2​,...,uN​} 中判定为"有目标"的传感器个数为k个，即:
λ ( u ) = P ( u / H 1 ) P ( u / H 0 ) = P ( { u 1 , u 2 , . . . , u N } / H 1 ) P ( { u 1 , u 2 , . . . , u N } / H 0 ) = ∏ ( k ) P ( u i = 1 / H 1 ) ∏ ( N − k ) P ( u i = 0 / H 1 ) ∏ ( k ) P ( u i = 1 / H 0 ) ∏ ( N − k ) P ( u i = 0 / H 0 ) = ∏ ( k ) P d i ∏ ( N − k ) P m i ∏ ( k ) P f i ∏ ( N − k ) ( 1 − P f i ) = ∏ ( k ) ( 1 − P m i ) ∏ ( N − k ) P m i ∏ ( k ) P f i ∏ ( N − k ) ( 1 − P f i ) \begin{aligned}\lambda(u)&=\frac{P(u/H_{1})}{P(u/H_{0})}=\frac{P(\lbrace u_{1},u_{2},...,u_{N} \rbrace/H_{1})}{P(\lbrace u_{1},u_{2},...,u_{N} \rbrace/H_{0})}\\ &=\frac{\prod\limits_{(k)}P(u_{i}=1/H_{1})\prod\limits_{(N-k)}P(u_{i}=0/H_{1})}{\prod\limits_{(k)}P(u_{i}=1/H_{0})\prod\limits_{(N-k)}P(u_{i}=0/H_{0})}\\ &=\frac{\prod\limits_{(k)}P^{i}_{d}\prod\limits_{(N-k)}P^{i}_{m}}{\prod\limits_{(k)}P^{i}_{f}\prod\limits_{(N-k)}(1-P^{i}_{f})}=\frac{\prod\limits_{(k)}(1-P^{i}_{m})\prod\limits_{(N-k)}P^{i}_{m}}{\prod\limits_{(k)}P^{i}_{f}\prod\limits_{(N-k)}(1-P^{i}_{f})} \end{aligned} λ(u)​=P(u/H0​)P(u/H1​)​=P({u1​,u2​,...,uN​}/H0​)P({u1​,u2​,...,uN​}/H1​)​=(k)∏​P(ui​=1/H0​)(N−k)∏​P(ui​=0/H0​)(k)∏​P(ui​=1/H1​)(N−k)∏​P(ui​=0/H1​)​=(k)∏​Pfi​(N−k)∏​(1−Pfi​)(k)∏​Pdi​(N−k)∏​Pmi​​=(k)∏​Pfi​(N−k)∏​(1−Pfi​)(k)∏​(1−Pmi​)(N−k)∏​Pmi​​​
似然函数比求对数：
$$\begin{array}{rl}& \lambda (u)=\frac{\prod _{(k)}(1-P_m^i)\prod _{(N-k)}{P}_m^i}{\prod _{(k)}{P}_f^i\prod _{(N-k)}(1-P_f^i)}\\ & \Rightarrow ln\lambda (u)=\sum _{(k)}ln\frac{1-P_m^i}{P_f^i}+\sum _{(N-k)}ln\frac{P_m^i}{1-P_f^i}\end{array}$$
计算似然函数比，并求对数：
$$\sum _{i=1}^N{w}_i>ln(P_0/P_1)?H_1:H_0$$
其中:$$w_i=\{\begin{array}{l}ln\frac{1-P_m^i}{P_f^i}{u}_i=1\\ ln\frac{P_m^i}{1-P_f^i}{u}_i=0\end{array}$$
一般记为：$u_0=\sum _{i=1}^N{w}_i>ln(P_0/P_{(}1))?1:0$
$u_0=1表示有目标，否则为无目标$

缺点：
在最大后验概率融合检测准则中，实现最佳融合检测准则必须知道假设$H_1或H_0$的先验概率$P_1$或$P_0=1-P_1$以及各个局部检测器的虚警概率$P_f^i$和漏检概率 $P_m^i$。

贝叶斯融合检测基本准则

在最大后验概率融合检测准则中，虚警和漏检两类错误都没有特殊加权，这样相当于假定两类错误是同等危险的。在工程应用中，各类错误的后果并非同等严重，即不同类型的错误所造成的损失或者说要付出的代价是不同的。为了反映这些差别，应对各类错误概率分别规定不同的代价，即所谓代价函数来反映这种损失的不同。

贝叶斯融合检测准则对每一个决策结果分配相应的代价值，基于假设概率得到平均总代价，检测策略是使平均总代价最小。
以单传感器二元假设检验问题为例定义代价：
$C_{ij}$为假设$H_i$成立时作出$u_j$决策的代价，且满足：
$C_{i̸=j,j}>C_{jj}$,平均总代价：
$$\begin{array}{rl}C& =P_0[C_{00}P(u_0|H_0)+C_{10}P(u_1|H_0)]+P_1[C_{01}P(u_0|H_1)+C_{11}P(u_1|H_1)]\\ & =P_0{C}_{00}P(u_0|H_0)+P_0{C}_{10}P(u_1|H_0)+P_1{C}_{01}P(u_0|H_1)+P_1{C}_{11}P(u_1|H_1)\\ & =P_0{C}_{00}(1-P(u_1|H_0))+P_0{C}_{10}P(u_1|H_0)+P_1{C}_{01}(1-P(u_1|H_1))+P_1{C}_{11}P(u_1|H_1)\\ & =P_0{C}_{00}+P_1{C}_{01}+P_0(C_{10}-C_{00})P(u_1|H_0)-P_1(C_{01}-C_{11})P(u_1|H_1))\end{array}$$
结合上面概率密度曲线图：注意到上式中：
$$虚警概率:P(u_1|H_0)={\int }_{R_1}{p}_0(z)dz$$
$$检测概率:P(u_1|H_1)={\int }_{R_1}{p}_1(z)dz$$
其中：$p_0(z)为没有目标的量测分布概率密度函数，p_1(z)为有目标的量测分布概率密度函数$
所以平均总代价：
$$C=P_0{C}_{00}+P_1{C}_{01}+{\int }_{R_1}[P_0(C_{10}-C_{00})p_0(z)-P_1(C_{01}-C_{11})p_1(z)]dz$$
应选择区域$R_1$使平均代价最小.即保证上式积分项小于零。
由此得到二元假设检测的贝叶斯检测准则：
$$\frac{p_1(z)}{p_0(z)}>\frac{P_0(C_{10}-C_{00})}{P_1(C_{01}-C_{11})}?H_1:H_0$$
对比最大后验概率准则：
$$\lambda (u)=\frac{P(u/H_1)}{P(u/H_0)}>\frac{P(H_0)}{P(H_1)}?H_1:H_0$$
显然：按贝叶斯准则与按最大后验概率准则得到的检测系统只是门限不同，【另外，注意表达式中概率密度函数与概率的区别】而当代价的选取满足$C_{10}-C_{00}=C_{01}-C_{11}$时，最大后验概率准则就是贝叶斯判决准则的特例。

最小误差概率准则：
在有一些应用场合，对两类错误没有什么特殊的区别，那么令所有误差的代价函数最小也是一个合理的准则，即令代价$C_{00}=C_{11}=0,C_{01}=C_{10}=1$那么代价函数式为：
$$\begin{array}{rl}C=P_e& =P_{0}P(u_1|H_0)+P_{1}P(u_0|H_1)=P_{0}P(u_1|H_0)+P_1(1-P(u_1|H_1))\\ & =P_{0}P(u_1|H_0)+P_1-P_{1}P(u_1|H_1)\\ & =P_1+{\int }_{R_1}[P_0{p}_0(z)-P_1{p}_1(z)]dz\end{array}$$
这里Pe为误差概率，使C最小的策略就是使误差概率Pe最小。为使Pe最小，则应使上式积分项尽可能小，

于是有$P_0{p}_0(z)-P_1{p}_1(z)<0$等价于

可得：$\frac{p_1(z)}{p_0(z)}>\frac{P_0}{P_1}$

因此最小误差概率准则与最大后验概率准则完全相同。准则为：
$$\frac{p_1(z)}{p_0(z)}>\frac{P_0}{P_1}?H_1:H_0$$

贝叶斯融合检测准则讨论

• $C_{10}-C_{00}$大，即虚警引起的损失大，门限应取大一些使虚警出现的可能性小一点；反之亦然。
• $C_{01}-C_{11}$大，即漏检引起的损失大，门限应取小一些使漏检出现的可能性小一点；反之亦然。

贝叶斯融合检测准则是多传感器系统优化决策的主流技术，是迄今为止理论上最完整的信息融合方法。

在各种先验概率及各种错误决策代价已知的情况下，贝叶斯方法是最优的方法，但是如何获得所需先验概率知识及各种错误决策的代价是应用该方法的一个关键问题。

N个传感器的贝叶斯融合检测***(后面几个公式务必深入研究其推导过程和所代表的物理意义)***

以分布式并行结构为例，对于N个传感器融合检测贝叶斯风险：
$$\begin{array}{rl}{R}_B& =C_{00}P(u_0=0|H_0)P_0+C_{01}P(u_0=0|H_1)+C_{10}P(u_0=1|H_0)P_0+C_{11}P(u_0=1|H_1)P_1\\ & =C_{00}(1-P_F)P_0+C_{01}(1-P_D)P_1+C_{10}{P}_F{P}_0+C_{11}{P}_D{P}_1\\ & =C_{01}{P}_1+C_{00}{P}_0+P_0(C_{10}-C_{00})P_F-P_1(C_{01}-C_{11})P_{11}{P}_D\\ & =C+C_F{P}_F-C_D{P}_D\end{array}\text{\hspace{1em}(1)}$$
其中：$$\begin{array}{rl}& C=C_{01}{P}_1+C_{00}{P}_0\\ & C_F=P_0(C_{10}-C_{00})\\ & C_D=P_1(C_{01}-C_{11})\end{array}\text{\hspace{1em}(2)}$$
注意：上面贝叶斯风险是为整个系统的贝叶斯风险，首先由各个传感器本身作出判决后将判决结果送入融合中心，融合中心的判决结果由$u_0$表示，每个传感器的判决结果统一由 $u=(u_1,u_2,...,u_N)$表示。

融合系统的发现概率和虚警概率的计算公式(注意区别于单个传感器)：
$$P_F=\sum _{u}P(u_0=1|u)P(u|H_0),P_D=\sum _{u}P(u_0=1|u)P(u|H_1)\text{\hspace{1em}(3)}$$

由公式(1)(2)(3)合并可推得融合系统贝叶斯风险为：
$$R_B=C+\sum _{u}P(u_0=1|u)[C_{F}P(u|H_0)-C_{D}P(u|H_1)]$$
由融合系统贝叶斯风险公式可知，为使贝叶斯风险达到最小的最优分布检测必要条件为：
$$P(u_0=1|u)=\{\begin{array}{l}1,C_{F}P(u|H_0)-C_{D}P(u|H_1)<0\\ 0,C_{F}P(u|H_0)-C_{D}P(u|H_1)\ge 0\end{array}$$
由以上必要条件可进一步推得：
$$\lambda (u)=\frac{P(u|H_1)}{P(u|H_0)}>\frac{C_F}{C_D}?H_1:H_0$$
或者：
$$\lambda (u)=\frac{P(u|H_1)}{P(u|H_0)}>\frac{C_F}{C_D}?1:0$$
一般来说，Bayes融合理论推导到这里就算全部介绍完成了，但是如果想要将该理论运用于实际工程上还是很难的，通过上式可以看出，要想运用该理论就必须知道每个传感器的似然函数，而一般在传感器技术文档中，是不会给出似然函数概率的，通常以概率密度的方式给出。下面就推导将上式转化为概率密度表达的方式。

设每个传感器有条件概率密度：$p_i(z_i|H_i),p_i(z_i|H_0)$以及定义：${\tilde{u}}_i=(u_i,...,u_{i-1},u_{i+1},...,u_N)$
$$A(\tilde{u})=P(u_0=1|\tilde{u},u_i=1)-P(u_0=1|{\tilde{u}}_i,u_i=0)$$
定理：则基于贝叶斯融合检测准则可以推出各个传感器的检测准则为：
$${\lambda }_i=\frac{p_i(z_i|H_i)}{p_i(z_i|H_0)}>\frac{\sum _{{\tilde{u}}_i}{C}_{F}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_0)}{\sum _{{\tilde{u}}_i}{C}_{D}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_1)}?H_1:H_0$$
证明：



证毕。
在解决了似然函数用条件概率密度替换的问题后，以上计算还需要先验概率密度知识。在计算上是非常困难的。为此人们重点研究各个传感器满足检测概率大于虚警概率，且彼此统计独立的条件下，贝叶斯融合检测准则。
条件1：决定了最优融合规则的单调性。
\quad 设每个传感器满足：$$P_D^i\ge P_F^{i}i=1,2,...,N$$
$$\Rightarrow P(u_0=1|{\tilde{u}}_i,u_i=1)-P(u_0=1|{\tilde{u}}_i,u_i=0)\ge 0$$
$$\Rightarrow A({\tilde{u}}_i)\ge 0$$
条件2：在各个传感器统计独立的条件下的贝叶斯融合检测准则可以简化为：
$${\lambda }_i=\frac{p_i(z_i|H_i)}{p_i(z_i|H_0)}>\frac{\sum _{{\tilde{u}}_i}{C}_{F}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_0)}{\sum _{{\tilde{u}}_i}{C}_{D}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_1)}?H_1:H_0$$
上式中各个条件概率为：
$$P({\tilde{u}}_i|H_0)=\prod _{{\tilde{u}}_{i,j}=1}{\int }_{R_{1,j}}{p}_j(z_j|H_0)d{z}_j\prod _{{\tilde{u}}_{i,j}=0}[1-{\int }_{R_{1,j}}{p}_j(z_j|H_0)d{z}_j]$$
$$P({\tilde{u}}_i|H_1)=\prod _{{\tilde{u}}_{i,j}=1}{\int }_{R_{1,j}}{p}_j(z_j|H_0)d{z}_j\prod _{{\tilde{u}}_{i,j}=0}[1-{\int }_{R_{1,j}}{p}_j(z_j|H_1)d{z}_j]$$
求解最优贝叶斯融合检测准则，需要在给定$C_F$和$C_D$的条件下求出：
$$P(u|H1),P(u|H_0)$$
由于最优融合检测准则与分布的各个传感器的贝叶斯判决门限是耦合的，所以在求解最优贝叶斯融合检测准则的同时，需要求解贝叶斯融合检测准则最佳判决门限：
$$T_i=\frac{\sum _{{\tilde{u}}_i}{C}_{F}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_0)}{\sum _{{\tilde{u}}_i}{C}_{D}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_1)}$$

求解最优贝叶斯融合检测准则和各个传感器判决门限的迭代算法
第一步：初始化条件。
\quad 设各个传感器的观测分布(已知)：
$$p_i(z_i|H_i),p_i(z_i|H_0)i=1,2,...,N$$
以及假设各个传感器的贝叶斯判决门限:
$$\begin{gathered} T_i^{0}i=1,2,...,N \\ 由p_i(z_i|H_j)算出：P^0(u|H_0),P^0(u|H_1) \end{gathered}$$
第二步：
设置循环变量n，和终止控制量$\xi$，准备循环。估算贝叶斯融合检测准则：
$$P^1(u_0=1|u)=\{\begin{array}{l}1,C_F{P}^0(u|H_0)-C_D{P}^0(u|H_1)<0\\ 0,C_F{P}^0(u|H_0)-C_D{P}^0(u|H_1)\ge 0\end{array}$$
计算对应的贝叶斯融合风险$R_B^{(1)}$:
$R_B^{(1)}=C+\sum _u{P}^{(1)}(u_0=1|u)\times [C_F{P}^0(u|H_0)-C_D{P}^0(u|H_1)]$
第三步：以上一次计算的各个传感器的贝叶斯判决门限为基础计算新的判决门限：
$$A^k(\tilde{u})=P^k(u_0=1|\tilde{u},u_i=1)-P^k(u_0=1|{\tilde{u}}_i,u_i=0)$$
$$T_i^k=\frac{\sum _{{\tilde{u}}_i}{C}_{F}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_0)}{\sum _{{\tilde{u}}_i}{C}_{D}A({\tilde{u}}_i)P({\tilde{u}}_i|z_i,H_1)}i=1,...,N$$
根据各个传感器新的门限计算对应传感器的$A({\tilde{u}}_i),P({\tilde{u}}_i|H_0),P({\tilde{u}}_i|H_1)$
以计算的各个传感器的贝叶斯判决门限为基础：$T_i^{k}i=1,2,...,N$
估算贝叶斯融合检测准则：
$$P^{k+1}(u_0=1|u)=\{\begin{array}{l}1,C_F{P}^k(u|H_0)-C_D{P}^k(u|H_1)<0\\ 0,C_F{P}^k(u|H_0)-C_D{P}^k(u|H_1)\ge 0\end{array}$$
第五步：
计算$k+1$此次迭代的贝叶斯风险$R_B^{(K+1)}:$
$R_B^{(k+1)}=C+\sum _u{P}^{(k+1)}(u_0=1|u)\times [C_F{P}^k(u|H_0)-C_D{P}^k(u|H_1)]$

仿真实例

以上题目只要严格按照如下方式编写matlab代码即可得出仿真结果：

clc;
j = 100;
V1 = 2.2;
V2 = 2.5;
V3 = 2.8;
delta1 = 1;
delta2 = 1.5;
delta3 = 1.2;

%初始化三个传感器的门限值
gama = zeros(j,3);

%代价初始化
C00 = 0;
C11 = 0;
C01 = 1;
C10 = 1;




%传感器1概率初始化
pf1 = zeros(j,1);
pm1 = zeros(j,1);
pd1 = zeros(j,1);
p001 = zeros(j,1);

%传感器2概率初始化
pf2 = zeros(j,1);
pm2 = zeros(j,1);
pd2 = zeros(j,1);
p002 = zeros(j,1);

%传感器3概率初始化
pf3 = zeros(j,1);
pm3 = zeros(j,1);
pd3 = zeros(j,1);
p003 = zeros(j,1);

Rb = zeros(j,1);

%先验概率初始化
P0=1/j:(1/j):1;
P1 = 1 - P0;

C = C01*P1 + C00*P0;
CF = P0*(C10-C00);
CD = P1*(C01-C11);
t = (CF./CD);

%初始判决门限，因为初始情况下，门限公式后面的无法计算

gama(1,1) = t(1);
gama(1,2) = t(1);
gama(1,3) = t(1);




for i = 1:j    
%根据门限计算各个传感器的四类概率
p001(i) = normcdf(gama(i,1),0,delta1);
pm1(i) = normcdf(gama(i,1),V1,delta1);
pd1(i) =1 -  pm1(i);
pf1(i) = 1 - p001(i);


p002(i) = normcdf(gama(i,2),0,delta2);
pm2(i) = normcdf(gama(i,2),V2,delta2);
pd2(i) =1 -  pm2(i);
pf2(i) = 1 - p002(i);

p003(i) = normcdf(gama(i,3),0,delta3);
pm3(i) = normcdf(gama(i,3),V3,delta3);
pd3(i) =1 -  pm3(i);
pf3(i) = 1 - p003(i);

%根据传感器的检测概率和虚假概率，计算融合中心融合准则
Pu000_H1 = pm1(i)*pm2(i)*pm3(i);
Pu000_H0 = p001(i)*p002(i)*p003(i);

Pu001_H1 = pm1(i)*pm2(i)*pd3(i);
Pu001_H0 = p001(i)*p002(i)*pf3(i);

Pu010_H1 = pm1(i)*pd2(i)*pm3(i);
Pu010_H0 = p001(i)*pf2(i)*p003(i);

Pu011_H1 = pm1(i)*pd2(i)*pd3(i);
Pu011_H0 = p001(i)*pf2(i)*pf3(i);

Pu100_H1 = pd1(i)*pm2(i)*pm3(i);
Pu100_H0 = pf1(i)*p002(i)*p003(i);

Pu101_H1 = pd1(i)*pm2(i)*pd3(i);
Pu101_H0 = pf1(i)*p002(i)*pf3(i);

Pu110_H1 = pd1(i)*pd2(i)*pm3(i);
Pu110_H0 = pf1(i)*pf2(i)*p003(i);

Pu111_H1 = pd1(i)*pd2(i)*pd3(i);
Pu111_H0 = pf1(i)*pf2(i)*pf3(i);

Pu0_1_u000 = (Pu000_H1/Pu000_H0) > t(i);
Pu0_1_u001 = (Pu001_H1/Pu001_H0) > t(i);
Pu0_1_u010 = (Pu010_H1/Pu010_H0) > t(i);
Pu0_1_u011 = (Pu011_H1/Pu011_H0) > t(i);
Pu0_1_u100 = (Pu100_H1/Pu100_H0) > t(i);
Pu0_1_u101 = (Pu101_H1/Pu101_H0) > t(i);
Pu0_1_u110 = (Pu110_H1/Pu110_H0) > t(i);
Pu0_1_u111 = (Pu111_H1/Pu111_H0) > t(i);

%计算融合系统的贝叶斯风险

Rb(i) = C(i) + Pu0_1_u000*(CF(i)*Pu000_H0 - CD(i)*Pu000_H1) + ...
    Pu0_1_u001*(CF(i)*Pu001_H0 - CD(i)*Pu001_H1)+...
    Pu0_1_u010*(CF(i)*Pu010_H0 - CD(i)*Pu010_H1)+...
    Pu0_1_u011*(CF(i)*Pu011_H0 - CD(i)*Pu011_H1)+...
    Pu0_1_u100*(CF(i)*Pu100_H0 - CD(i)*Pu100_H1)+...
    Pu0_1_u101*(CF(i)*Pu101_H0 - CD(i)*Pu101_H1)+...
    Pu0_1_u110*(CF(i)*Pu110_H0 - CD(i)*Pu110_H1)+...
    Pu0_1_u111*(CF(i)*Pu111_H0 - CD(i)*Pu111_H1);

%计算A(u)，确定某一个传感器的值后组合两外两个传感器的值
%确定u1
Au1_00 = Pu0_1_u100 - Pu0_1_u000;
Au1_01 = Pu0_1_u101 - Pu0_1_u001;
Au1_10 = Pu0_1_u110 - Pu0_1_u010;
Au1_11 = Pu0_1_u111 - Pu0_1_u011;

%确定u2
Au2_00 = Pu0_1_u010 - Pu0_1_u000;
Au2_01 = Pu0_1_u011 - Pu0_1_u001;
Au2_10 = Pu0_1_u110 - Pu0_1_u100;
Au2_11 = Pu0_1_u111 - Pu0_1_u101;

%确定u3
Au3_00 = Pu0_1_u001 - Pu0_1_u000;
Au3_01 = Pu0_1_u011 - Pu0_1_u010;
Au3_10 = Pu0_1_u101 - Pu0_1_u100;
Au3_11 = Pu0_1_u111 - Pu0_1_u110;

%计算后验概率
Pu1_00_H0 = p002(i)*p003(i);
Pu1_01_H0 = p002(i)*pf3(i);
Pu1_10_H0 = pf2(i)*p003(i);
Pu1_11_H0 = pf2(i)*pf3(i);

Pu2_00_H0 = p001(i)*p003(i);
Pu2_01_H0 = p001(i)*pf3(i);
Pu2_10_H0 = pf1(i)*p003(i);
Pu2_11_H0 = pf1(i)*pf3(i);

Pu3_00_H0 = p001(i)*p002(i);
Pu3_01_H0 = p001(i)*pf2(i);
Pu3_10_H0 = pf1(i)*p002(i);
Pu3_11_H0 = pf1(i)*pf2(i);

Pu1_00_H1 = pm2(i)*pm3(i);
Pu1_01_H1 = pm2(i)*pd3(i);
Pu1_10_H1 = pd2(i)*pm3(i);
Pu1_11_H1 = pm2(i)*pm3(i);

Pu2_00_H1 = pm1(i)*pm3(i);
Pu2_01_H1 = pm1(i)*pd3(i);
Pu2_10_H1 = pd1(i)*pm3(i);
Pu2_11_H1 = pm1(i)*pm3(i);

Pu3_00_H1 = pm1(i)*pm2(i);
Pu3_01_H1 = pm1(i)*pd2(i);
Pu3_10_H1 = pd1(i)*pm2(i);
Pu3_11_H1 = pm1(i)*pm2(i);

%计算三种传感器的融合判决门限

T1 =(CF(i)*(Au1_00*Pu1_00_H0 + Au1_01*Pu1_01_H0 + Au1_10*Pu1_10_H0 + Au1_11*Pu1_11_H0))/...
    (CD(i)*(Au1_00*Pu1_00_H1 + Au1_01*Pu1_01_H1 + Au1_10*Pu1_10_H1 + Au1_11*Pu1_11_H1));

T2 =(CF(i)*(Au2_00*Pu2_00_H0 + Au2_01*Pu2_01_H0 + Au2_10*Pu2_10_H0 + Au2_11*Pu2_11_H0))/...
    (CD(i)*(Au2_00*Pu2_00_H1 + Au2_01*Pu2_01_H1 + Au2_10*Pu2_10_H1 + Au2_11*Pu2_11_H1));

T3 =(CF(i)*(Au3_00*Pu3_00_H0 + Au3_01*Pu3_01_H0 + Au3_10*Pu3_10_H0 + Au3_11*Pu3_11_H0))/...
    (CD(i)*(Au3_00*Pu3_00_H1 + Au3_01*Pu3_01_H1 + Au3_10*Pu3_10_H1 + Au3_11*Pu3_11_H1));
%最后的门限公式和理论公式对应不上
gama(i+1,1) = log(T1)/V1 + V1/2;
gama(i+1,2) = log(T2)/V2 + V2/2;
gama(i+1,3) = log(T3)/V3 + V3/2;

end
%传感器1贝叶斯风险
R1 = (C00*p001+C10*pf1).*P0' + (C01*pm1+C11*pd1).*P1';
%传感器2贝叶斯风险
R2 = (C00*p002+C10*pf2).*P0' + (C01*pm2+C11*pd2).*P1';
%传感器3贝叶斯风险
R3 = (C00*p003+C10*pf3).*P0' + (C01*pm3+C11*pd3).*P1';

plot(P0,R1,'b',P0,R2,'g',P0,R3,'black',P0,Rb,'r'); 
title('贝叶斯风险曲线图');
xlabel('p0');
ylabel('贝叶斯风险值');
legend('传感器1贝叶斯风险曲线','传感器2贝叶斯风险曲线','传感器3贝叶斯风险曲线','融合系统贝叶斯风险曲线');

仿真结果图如下所示。

总结：以上理论和机器学习里面的贝叶斯分类理论相似，都是进行分类。