N-P定理
对于一个给定的
P
F
A
=
α
P_{FA}=\alpha
PFA=α,使
P
D
P_{D}
PD最大的判决为
L
(
x
)
=
p
(
x
;
H
1
)
p
(
x
;
H
0
)
>
γ
L(\mathbf{x})=\frac{p\left(\mathbf{x} ;H_{1}\right)}{p\left(\mathbf{x} ; H_{0}\right)}>\gamma
L(x)=p(x;H0)p(x;H1)>γ
其中,门限
γ
\gamma
γ由
P
F
A
=
∫
{
x
:
L
(
x
)
>
γ
}
p
(
x
;
H
0
)
d
x
=
α
P_{F A}=\int_{\{\mathbf{x}: L(x)>\gamma\}} p\left(\mathbf{x} ; H_{0}\right) d \mathbf{x}=\alpha
PFA=∫{x:L(x)>γ}p(x;H0)dx=α
求出。
推导
NP定理的定义是:对于固定
P
F
A
P_{FA}
PFA使
P
D
P_{D}
PD最大的判决,即
max
P
D
s.t.
P
F
A
=
α
\text { max } P_{D} \text{ s.t. } P_{FA} = \alpha
max PD s.t. PFA=α
采用拉格朗日乘子,构造目标函数
F
=
P
D
+
λ
(
P
F
A
−
α
)
=
∫
R
1
p
(
x
;
H
1
)
d
x
+
λ
(
∫
R
1
p
(
x
;
H
0
)
d
x
−
α
)
=
∫
R
1
(
p
(
x
;
H
1
)
+
λ
p
(
x
;
H
0
)
)
d
x
−
λ
α
\begin{aligned} F&=P_{D}+\lambda\left(P_{F A}-\alpha\right) \\ &=\int_{R_{1}} p\left(\mathbf{x} ; \mathcal{H}_{1}\right) d \mathbf{x}+\lambda\left(\int_{R_{1}} p\left(\mathbf{x} ; \mathcal{H}_{0}\right) d \mathbf{x}-\alpha\right) \\ &=\int_{R_{1}}\left(p\left(\mathbf{x} ; \mathcal{H}_{1}\right)+\lambda p\left(\mathbf{x} ; \mathcal{H}_{0}\right)\right) d \mathbf{x}-\lambda \alpha \end{aligned}
F=PD+λ(PFA−α)=∫R1p(x;H1)dx+λ(∫R1p(x;H0)dx−α)=∫R1(p(x;H1)+λp(x;H0))dx−λα
为了使目标函数F最大,应该在
x
∈
R
1
x \in R1
x∈R1的情况下让
f
(
x
)
f(\mathbf{x})
f(x)大于零
f
(
x
)
=
p
(
x
;
H
1
)
+
λ
p
(
x
;
H
0
)
>
0
,
x
∈
R
1
f(\mathbf{x})=p\left(\mathbf{x} ; \mathcal{H}_{1}\right)+\lambda p\left(\mathbf{x} ; \mathcal{H}_{0}\right)>0 \text{ , } \mathbf{x} \in R1
f(x)=p(x;H1)+λp(x;H0)>0 , x∈R1
此时在
R
1
R1
R1域对
f
(
x
)
f(\mathbf{x})
f(x)的积分最大。
所以在判决时,若
p
(
x
;
H
1
)
p
(
x
;
H
0
)
>
−
λ
\frac{p\left(\mathbf{x} ; \mathcal{H}_{1}\right)}{p\left(\mathbf{x} ; \mathcal{H}_{0}\right)}>-\lambda
p(x;H0)p(x;H1)>−λ
判为
H
1
H_{1}
H1。
令
γ
=
−
λ
\gamma=-\lambda
γ=−λ,得到
p
(
x
;
H
1
)
p
(
x
;
H
0
)
>
γ
\frac{p\left(\mathbf{x} ; \mathcal{H}_{1}\right)}{p\left(\mathbf{x} ; \mathcal{H}_{0}\right)}>\gamma
p(x;H0)p(x;H1)>γ
判
H
1
H_{1}
H1。其中门限
γ
>
0
\gamma>0
γ>0由
P
F
A
=
α
P_{FA}=\alpha
PFA=α求得
P
F
A
=
∫
{
x
:
L
(
x
)
>
γ
}
p
(
x
;
H
0
)
d
x
=
α
P_{F A}=\int_{\{\mathbf{x}: L(x)>\gamma\}} p\left(\mathbf{x} ; H_{0}\right) d \mathbf{x}=\alpha
PFA=∫{x:L(x)>γ}p(x;H0)dx=α
参考《统计信号处理基础》