一、知识回顾
- 在简单二元假设检验中,我们有假设 { H 0 r(t)=n(t) H 1 r(t)=d+n(t) \begin{cases} H_0&\text{r(t)=n(t)} \\H_1&\text{r(t)=d+n(t)} \end{cases} {H0H1r(t)=n(t)r(t)=d+n(t) ,并且假定 p ( H 1 ) = 1 − ξ , p ( H 0 ) = ξ p(H_1)=1-\xi,p(H_0)=\xi p(H1)=1−ξ,p(H0)=ξ,即已知先验分布。
- 似然比检测器:
- 划分判决域之后,所求得的判决门限值仍然为 λ 0 = ξ 1 − ξ \lambda_0= \frac{\xi}{1-\xi} λ0=1−ξξ,得出最大后验概率准则(MAP) ⇔ \Lrarr ⇔最小错误概率准则
- 当 ξ = 1 2 \xi=\frac12 ξ=21时最大似然准则与最小错误准则等价
- 两类错误概率:
- 虚报概率 P F = P ( D 1 ∣ H 0 ) = ∫ R 1 p 0 ( x ) d x P_F=P(D_1|H_0)=\int_{R_1} p_0(x) dx PF=P(D1∣H0)=∫R1p0(x)dx
- 漏报概率 P M = P ( D 0 ∣ H 1 ) = ∫ R 0 p 1 ( x ) d x P_M=P(D_0|H_1)=\int_{R_0} p_1(x)dx PM=P(D0∣H1)=∫R0p1(x)dx
二、贝叶斯准则
2.1 代价因子的引入
- 在最大后验概率准则和最小错误概率准则中,两类错误都没有进行加权,即并没有对每类错误所造成的后果进行衡量。
- 在实际系统中,这两类错误所造成的后果大多数时候是不一样的,如:雷达系统中漏报的风险代价是高于虚报的风险代价的。
- 虽然代价因子很难确定,但是具有一定的理论意义
2.2 代价因子定义和性质
- 定义
将 H j H_j Hj为真选择 H i H_i Hi的代价定义为: c i j c_{ij} cij,也就是说实际发送的是事件 H j H_j Hj,但是接收机判定为 D i D_i Di。 - 性质
①代价具有非负性
②在任何情况下,判断正确的代价都应该小于判断错误的代价,即有 { c 10 − c 00 > 0 发送0 c 01 − c 11 > 0 发送1 \begin{cases}c_{10} -c_{00}\gt0 &\text{发送0}\\ c_{01} -c_{11}\gt0&\text{发送1}\end{cases} {c10−c00>0c01−c11>0发送0发送1
2.3 代价的计算
在平均错误概率的基础上对错误概率进行加权:
C
‾
=
ξ
[
c
00
(
1
−
P
F
)
+
c
10
P
F
]
+
(
1
−
ξ
)
[
c
01
P
M
+
c
11
(
1
−
P
M
)
]
\overline{C}=\xi[c_{00}(1-P_F)+c_{10}P_F]+(1-\xi)[c_{01}P_M+c_{11}(1-P_M)]
C=ξ[c00(1−PF)+c10PF]+(1−ξ)[c01PM+c11(1−PM)]
其中,
P
F
=
∫
R
1
p
0
(
x
)
d
x
P_F=\int_{R1}p_0(x)dx
PF=∫R1p0(x)dx,
P
M
=
∫
R
0
p
1
(
x
)
d
x
P_M=\int_{R0}p_1(x)dx
PM=∫R0p1(x)dx
从这个式子可以看出,平均代价也是与接收机的判决方式有关的函数,需要做的就是最小化该代价函数
-
最小化代价函数: m i n R 1 C ‾ \underset{R_1}{min}\overline{C} R1minC
C ‾ = ξ [ c 00 ( 1 − ∫ R 1 p 0 ( x ) d x ) + c 10 ∫ R 1 p 0 ( x ) d x ] + ( 1 − ξ ) [ c 01 ∫ R 0 p 1 ( x ) d x + [ c 11 ( 1 − ∫ R 0 p 1 ( x ) d x ) ] \overline{C}=\xi[c_{00}(1-\int_{R1}p_0(x)dx)+c_{10}\int_{R1}p_0(x)dx] \\+(1-\xi)[c_{01}\int_{R0}p_1(x)dx+[c_{11}(1-\int_{R0}p_1(x)dx)] C=ξ[c00(1−∫R1p0(x)dx)+c10∫R1p0(x)dx]+(1−ξ)[c01∫R0p1(x)dx+[c11(1−∫R0p1(x)dx)]
由于 R 0 R_0 R0和 R 1 R_1 R1为互斥完备集,所以有 ∫ R 0 p 1 ( x ) d x = 1 − ∫ R 1 p 1 ( x ) d x \int_{R0}p_1(x)dx=1-\int_{R1}p_1(x)dx ∫R0p1(x)dx=1−∫R1p1(x)dx,将其代入,可得:
C ‾ = ξ [ c 00 ( 1 − ∫ R 1 p 0 ( x ) d x ) + c 10 ∫ R 1 p 0 ( x ) d x ] + ( 1 − ξ ) [ c 01 ( 1 − ∫ R 1 p 1 ( x ) d x ) + c 11 ∫ R 1 p 1 ( x ) d x ] = ξ c 00 + c 01 ( 1 − ξ ) + ∫ R 1 [ ξ ( c 10 − c 00 ) p 0 ( x ) − ( 1 − ξ ) ( c 01 − c 11 ) p 1 ( x ) ] d x \overline{C}=\xi[c_{00}(1-\int_{R1}p_0(x)dx)+c_{10}\int_{R1}p_0(x)dx] \\+(1-\xi)[c_{01}(1-\int_{R1}p_1(x)dx)+c_{11}\int_{R1}p_1(x)dx] \\ = \xi c_{00}+c_{01}(1-\xi)+\int_{R_1}[\xi(c_{10}-c_{00})p_0(x)-(1-\xi)(c_{01}-c_{11})p_1(x)]dx C=ξ[c00(1−∫R1p0(x)dx)+c10∫R1p0(x)dx]+(1−ξ)[c01(1−∫R1p1(x)dx)+c11∫R1p1(x)dx]=ξc00+c01(1−ξ)+∫R1[ξ(c10−c00)p0(x)−(1−ξ)(c01−c11)p1(x)]dx
为了使平均代价最小,只需要在判决域 R 1 R_1 R1上被积函数满足: [ ξ ( c 10 − c 00 ) p 0 ( x ) − ( 1 − ξ ) ( c 01 − c 11 ) p 1 ( x ) ] < 0 [\xi(c_{10}-c_{00})p_0(x)-(1-\xi)(c_{01}-c_{11})p_1(x)]\lt0 [ξ(c10−c00)p0(x)−(1−ξ)(c01−c11)p1(x)]<0利用该思想,可以得出下列关系式: p 1 ( x ) p 0 ( x ) > ξ ( c 10 − c 00 ) ( 1 − ξ ) ( c 01 − c 11 ) \frac{p_1(x)}{p_0(x)}\gt\frac{\xi(c_{10}-c_{00})}{(1-\xi)(c_{01}-c_{11})} p0(x)p1(x)>(1−ξ)(c01−c11)ξ(c10−c00)
即:
分析:根据上式判决方式,如果 λ 0 = = 0 \lambda_0==0 λ0==0,则接收端永远判决为 D 1 D_1 D1,且此时虚报概率 P F = P ( D 1 ∣ H 0 ) = 1 P_F=P(D_1|H_0)=1 PF=P(D1∣H0)=1,漏报概率 P M = P ( D 0 ∣ H 1 ) = 0 P_M=P(D_0|H_1)=0 PM=P(D0∣H1)=0,该情况适合漏报代价 c 01 c_{01} c01非常大的情况;
如果 λ 0 = = ∞ \lambda_0==\infty λ0==∞,则接收机永远判决为 D 0 D_0 D0,此时虚报概率 P F = P ( D 1 ∣ H 0 ) = 0 P_F=P(D_1|H_0)=0 PF=P(D1∣H0)=0,漏报概率 P M = P ( D 0 ∣ H 1 ) = 1 P_M=P(D_0|H_1)=1 PM=P(D0∣H1)=1,该情况适合漏报代价 c 10 c_{10} c10非常大的情况;
2.4 贝叶斯准则的判决门限
-
准则决定门限
贝叶斯准则的思想:使判决的代价最小化
根据2.3节有关代价的计算可以得出,贝叶斯准则的判决门限值为: ξ ( c 10 − c 00 ) ( 1 − ξ ) ( c 01 − c 11 ) \frac{\xi(c_{10}-c_{00})}{(1-\xi)(c_{01}-c_{11})} (1−ξ)(c01−c11)ξ(c10−c00) -
对于数字通信系统而言,正确判决无代价,错误判决代价相同,即 { c 11 = c 00 = 0 c 01 = c 10 = 1 \begin{cases} c_{11}=c_{00}=0 \\ c_{01}=c_{10}=1\end{cases} {c11=c00=0c01=c10=1
此时,贝叶斯准则和最大后验概率准则等价。