雷达信号检测问题与奈曼-皮尔逊准则
雷达信号检测问题属于一个二元检测问题,即或者为目标加噪声情况,或者为仅含有噪声情况。对观察信号空间DDD的划分问题,即将DDD空间划分为D1D_1D1 (有信号)和D0D_0D0 (无信号)两个子空间,并满足D=D0∪D1D=D_0\cup D_1D=D0∪D1,D0∩D1=∅D_0\cap D_1=\emptysetD0∩D1=∅。子空间D0D_0D0和D1D_1D1称为判决域。若某个观测量(x∣Hi)(i=0,1)(x|{{H}_{i}})(i=0,1)(x∣Hi)(i=0,1)落入D1D_1D1域,就判决假设H1H_1H1成立,否则就判决假设H0H_0H0成立。
在不考虑杂波与干扰的问题下,认为接收回波中仅可能存在信号加噪声或纯噪声这两种情况。将仅含有噪声的情况设为H0H_0H0假设,含有信号加噪声的情况为H1H_1H1假设,有
{H0:x(t)=n(t)H1:x(t)=s(t)+n(t)
\left\{ \begin{aligned}
& {{H}_{0}}:x(t)=n(t) \\
& {{H}_{1}}:x(t)=s(t)+n(t) \\
\end{aligned} \right.
{H0:x(t)=n(t)H1:x(t)=s(t)+n(t)
其中n(t)n(t)n(t)为均值为0,方差为σn2\sigma_n^2σn2的高斯白噪声,s(t)s(t)s(t)为确知信号。对二元检测问题,存在两种正确判决和两种错误判决:
Pd=P(H1∣H1)Pn=P(H0∣H0)Pm=P(H0∣H1)Pfa=P(H0∣H0)
\begin{aligned}
& {{P}_{d}}=P({{H}_{1}}|{{H}_{1}}) \\
& {{P}_{n}}=P({{H}_{0}}|{{H}_{0}}) \\
& {{P}_{m}}=P({{H}_{0}}|{{H}_{1}}) \\
& {{P}_{fa}}=P({{H}_{0}}|{{H}_{0}}) \\
\end{aligned}
Pd=P(H1∣H1)Pn=P(H0∣H0)Pm=P(H0∣H1)Pfa=P(H0∣H0)
其中,Pd{{P}_{d}}Pd为检测概率,Pn{{P}_{n}}Pn为正确不发现概率,Pm{{P}_{m}}Pm为漏警概率,Pfa{{P}_{fa}}Pfa为虚警概率。
检测概率与漏警概率的关系为:
Pd+Pm=P(H1∣H1)+P(H0∣H1)=1
{{P}_{d}}+{{P}_{m}}=P({{H}_{1}}|{{H}_{1}})+P({{H}_{0}}|{{H}_{1}})=1
Pd+Pm=P(H1∣H1)+P(H0∣H1)=1
正确不发现概率与虚警概率的关系为:
Pn+Pfa=P(H0∣H0)+P(H1∣H0)=1
{{P}_{n}}+{{P}_{fa}}=P({{H}_{0}}|{{H}_{0}})+P({{H}_{1}}|{{H}_{0}})=1
Pn+Pfa=P(H0∣H0)+P(H1∣H0)=1
假设H0{{H}_{0}}H0出现的先验概率为P(H0)P({{H}_{0}})P(H0),H1{{H}_{1}}H1出现的先验概率为P(H1)P({{H}_{1}})P(H1),观测信号x(t)x(t)x(t)在H1{{H}_{1}}H1假设下的条件概率密度函数为:
p(x∣H0)=1(2π)σne−x22σn2
p(x|{{H}_{0}})=\frac{1}{\sqrt{(2\pi )}{{\sigma }_{n}}}{{e}^{-\frac{{{x}^{2}}}{2\sigma _{n}^{2}}}}
p(x∣H0)=(2π)σn1e−2σn2x2
x(t)x(t)x(t)在H0{{H}_{0}}H0假设下的条件概率密度函数为:
p(x∣H1)=1(2π)σne−(x−s)22σn2
p(x|{{H}_{1}})=\frac{1}{\sqrt{(2\pi )}{{\sigma }_{n}}}{{e}^{-\frac{{{(x-s)}^{2}}}{2\sigma _{n}^{2}}}}
p(x∣H1)=(2π)σn1e−2σn2(x−s)2
设定判决门限VT{{V}_{T}}VT,若输入信号的包络超过了该门限,则认为有目标,否则认为没有目标。故根据H1{{H}_{1}}H1假设下的条件概率密度函数可得检测概率为:
Pd=∫VT∞p(x∣H1)dx
{{P}_{d}}=\int_{{{V}_{T}}}^{\infty }{p(x|{{H}_{1}})dx}
Pd=∫VT∞p(x∣H1)dx
根据H0{{H}_{0}}H0假设下的条件概率密度函数可得虚警概率为:
Pfa=∫VT∞p(x∣H0)dx
{{P}_{fa}}=\int_{{{V}_{T}}}^{\infty }{p(x|{{H}_{0}})dx}
Pfa=∫VT∞p(x∣H0)dx
判决门限VT{{V}_{T}}VT的选择与所采用的最佳检测检测准则有关。常用的最佳检测准则有最大后验概率准则、贝叶斯准则、极小极大化准则、最小错误概率准则、奈曼-皮尔逊准则等。在雷达信号检测领域,由于预先未知目标出现的概率,很难确定一次漏警所造成的损失。在一定的虚警概率下,使漏警概率最小或使检测概率最大,为奈曼-皮尔逊准则。
数学上,用如下表达式来描述奈曼-皮尔逊准则:
Pe=Pm+λ0Pfa=P(H0∣H1)+λ0P(H1∣H0)
{{P}_{e}}={{P}_{m}}+{{\lambda }_{0}}{{P}_{fa}}=P({{H}_{0}}|{{H}_{1}})+{{\lambda }_{0}}P({{H}_{1}}|{{H}_{0}})
Pe=Pm+λ0Pfa=P(H0∣H1)+λ0P(H1∣H0)
由于检测概率与漏警概率之间满足Pd+Pm=1{{P}_{d}}+{{P}_{m}}=1Pd+Pm=1,有
Pe=1−Pd+λ0Pfa=1−P(H1∣H1)+λ0P(H1∣H0)
{{P}_{e}}=1-{{P}_{d}}+{{\lambda }_{0}}{{P}_{fa}}=1-P({{H}_{1}}|{{H}_{1}})+{{\lambda }_{0}}P({{H}_{1}}|{{H}_{0}})
Pe=1−Pd+λ0Pfa=1−P(H1∣H1)+λ0P(H1∣H0)
为提高判决质量,减小噪声干扰随机性的影响,一般需要对接收信号进行多次观测采样,对于NNN次独立取样,输入信号为NNN维空间,接收样本矢量表示为
x=[x1,x2,...,xN]T
\mathbf{x}={{[{{x}_{1}},{{x}_{2}},...,{{x}_{N}}]}^{\text{T}}}
x=[x1,x2,...,xN]T
当输入为x(t)=s(t)+n(t)x(t)=s(t)+n(t)x(t)=s(t)+n(t)时,NNN个取样点的联合概率分布密度函数为p(x1,x2,⋯ ,xN∣H1)p({{x}_{1}},{{x}_{2}},\cdots ,{{x}_{N}}|{{H}_{1}})p(x1,x2,⋯,xN∣H1),当输入为x(t)=n(t)x(t)=n(t)x(t)=n(t)时,其联合概率分布密度函数为p(x1,x2,⋯ ,xN∣H0)p({{x}_{1}},{{x}_{2}},\cdots ,{{x}_{N}}|{{H}_{0}})p(x1,x2,⋯,xN∣H0),检测概率和虚警概率分别为
Pd=∬D1⋅⋅⋅∫p(x1,x2,⋅⋅⋅⋅xN∣H1)dx1dx2⋅⋅⋅dxN
{{P}_{d}}=\iint_{{{D}_{1}}}{\cdot }\cdot \cdot \int{p}({{x}_{1}}\text{,}{{x}_{2}}\text{,}\cdot \cdot \cdot \cdot {{x}_{N}}\mid {{H}_{1}})\text{d}{{x}_{1}}\text{d}{{x}_{2}}\cdot \cdot \cdot \text{d}{{x}_{N}}
Pd=∬D1⋅⋅⋅∫p(x1,x2,⋅⋅⋅⋅xN∣H1)dx1dx2⋅⋅⋅dxN
Pfa=∬D1⋅⋅⋅∫p(x1,x2,⋅⋅⋅xN∣H0)dx1dx2⋅⋅⋅dxN {{P}_{fa}}=\iint_{{{D}_{1}}}{\cdot }\cdot \cdot \int{p}({{x}_{1}}\text{,}{{x}_{2}}\text{,}\cdot \cdot \cdot {{x}_{N}}\mid H_0)\text{d}{{x}_{1}}\text{d}{{x}_{2}}\cdot \cdot \cdot \text{d}{{x}_{N}} Pfa=∬D1⋅⋅⋅∫p(x1,x2,⋅⋅⋅xN∣H0)dx1dx2⋅⋅⋅dxN
则总错误概率与联合概率分布密度函数的关系为
Pe=1−∬D1⋅⋅⋅∫[p(x1,x2,⋅⋅⋅,xN∣H1)−λ0p(x1,x2,⋅⋅⋅xN∣H0)]dx1dx2⋅⋅⋅dxN
{{P}_{\text{e}}}=1-\iint_{{{D}_{1}}}{\cdot }\cdot \cdot \int{\left[ p\left( {{x}_{1}},{{x}_{2}},\cdot \cdot \cdot ,{{x}_{N}}\mid {{H}_{1}} \right)-{{\lambda }_{0}}p\left( {{x}_{1}},{{x}_{2}},\cdot \cdot \cdot {{x}_{N}}\mid {{H}_{0}} \right) \right]}\text{d}{{x}_{1}}\text{d}{{x}_{2}}\cdot \cdot \cdot \text{d}{{x}_{N}}
Pe=1−∬D1⋅⋅⋅∫[p(x1,x2,⋅⋅⋅,xN∣H1)−λ0p(x1,x2,⋅⋅⋅xN∣H0)]dx1dx2⋅⋅⋅dxN
观察空间的划分应保证总错误概率Pe{{P}_{\text{e}}}Pe最小,即后面的积分值最大。因此有
p(x1,x2,⋅⋅⋅,xN∣H1)−λ0p(x1,x2,⋅⋅⋅xN∣H0)≥0
p\left( {{x}_{1}},{{x}_{2}},\cdot \cdot \cdot ,{{x}_{N}}\mid {{H}_{1}} \right)-{{\lambda }_{0}}p\left( {{x}_{1}},{{x}_{2}},\cdot \cdot \cdot {{x}_{N}}\mid {{H}_{0}} \right)\ge 0
p(x1,x2,⋅⋅⋅,xN∣H1)−λ0p(x1,x2,⋅⋅⋅xN∣H0)≥0
的点判定为有信号,否则判定为无信号。即
p(x1,x2…xN∣H1)p(x1,x2…xN∣H0){≥λ0<λ0
\frac{p({{x}_{1}}\text{,}{{x}_{2}}\ldots {{x}_{N}}\mid {{H}_{1}})}{p({{x}_{1}}\text{,}{{x}_{2}}\ldots {{x}_{N}}\mid {{H}_{0}})}\left\{ \begin{aligned}
& \ge {{\lambda }_{0}}\quad \\
& <{{\lambda }_{0}}\quad \\
\end{aligned} \right.
p(x1,x2…xN∣H0)p(x1,x2…xN∣H1){≥λ0<λ0
参考文献:
[1] 陈伯孝等. 现代雷达系统分析与设计[M]. 陕西:西安电子科技大学出版社. 2012.
1万+

被折叠的 条评论
为什么被折叠?



