设采样点为
y
n
,
n
=
1
,
.
.
.
,
N
y_n,n=1,...,N
yn,n=1,...,N,这些采样点由目标信号
m
=
A
e
i
θ
m=Ae^{i\theta}
m=Aeiθ和I/Q两路功率均为
σ
2
/
2
\sigma^2/2
σ2/2的高斯白噪声
ω
n
\omega_n
ωn组成,对应假设
H
0
H_0
H0和
H
1
H_1
H1下
y
n
y_n
yn的组成分别为:
H
0
:
y
n
=
ω
n
H
1
:
y
n
=
m
+
ω
n
H_0: y_n=\omega_n \\ H_1: y_n=m+\omega_n
H0:yn=ωnH1:yn=m+ωn
设
z
n
=
∣
y
n
∣
z_n=|y_n|
zn=∣yn∣,则对应假设
H
0
H_0
H0和
H
1
H_1
H1下
z
n
z_n
zn分别服从瑞利分布和莱斯分布,其概率密度函数可以表示为:
H
0
:
z
n
=
∣
ω
n
∣
,
f
(
z
n
∣
H
0
)
=
2
z
n
σ
2
e
x
p
(
−
z
n
2
σ
2
)
,
z
n
≥
0
H_0: z_n=|\omega_n|, f(z_n|H_0)=\frac{2z_n}{\sigma^2}{\rm exp}(-\frac{z_n^2}{\sigma^2}),z_n \geq 0
H0:zn=∣ωn∣,f(zn∣H0)=σ22znexp(−σ2zn2),zn≥0
H
1
:
z
n
=
∣
m
+
ω
n
∣
,
f
(
z
n
∣
H
1
)
=
2
z
n
σ
2
e
x
p
(
−
z
n
2
+
A
2
σ
2
)
I
0
(
2
A
z
n
σ
2
)
,
z
n
≥
0
H_1: z_n=|m+\omega_n|, f(z_n|H_1)=\frac{2z_n}{\sigma^2}{\rm exp}(-\frac{z_n^2+A^2}{\sigma^2})I_0(\frac{2Az_n}{\sigma^2}),z_n \geq 0
H1:zn=∣m+ωn∣,f(zn∣H1)=σ22znexp(−σ2zn2+A2)I0(σ22Azn),zn≥0
则
N
N
N个采样点的联合概率密度函数可以表示为:
H
0
:
z
n
=
∣
ω
n
∣
,
f
(
z
∣
H
0
)
=
∏
n
=
0
N
−
1
2
z
n
σ
2
e
x
p
(
−
z
n
2
σ
2
)
,
z
n
≥
0
H_0: z_n=|\omega_n|,f(\boldsymbol{z}|H_0)=\prod_{n=0}^{N-1}{\frac{2z_n}{\sigma^2}{\rm exp}(-\frac{z_n^2}{\sigma^2})},z_n \geq 0
H0:zn=∣ωn∣,f(z∣H0)=n=0∏N−1σ22znexp(−σ2zn2),zn≥0
H 1 : z n = ∣ m + ω n ∣ , f ( z ∣ H 1 ) = ∏ n = 0 N − 1 2 z n σ 2 e x p ( − z n 2 + A 2 σ 2 ) I 0 ( 2 A z n σ 2 ) , z n ≥ 0 H_1: z_n=|m+\omega_n|,f(\boldsymbol{z}|H_1)=\prod_{n=0}^{N-1}{\frac{2z_n}{\sigma^2}{\rm exp}(-\frac{z_n^2+A^2}{\sigma^2})I_0(\frac{2Az_n}{\sigma^2})},z_n \geq 0 H1:zn=∣m+ωn∣,f(z∣H1)=n=0∏N−1σ22znexp(−σ2zn2+A2)I0(σ22Azn),zn≥0
所以,此时似然比检验LRT可以表示为:
Λ
=
f
(
z
∣
H
1
)
f
(
z
∣
H
0
)
=
e
x
p
(
−
A
2
σ
2
)
∏
n
=
0
N
−
1
I
0
(
2
A
z
n
σ
2
)
≷
−
λ
\Lambda=\frac{f(\boldsymbol{z}|H_1)}{f(\boldsymbol{z}|H_0)}={\rm exp}(-\frac{A^2}{\sigma^2})\prod_{n=0}^{N-1}{I_0(\frac{2Az_n}{\sigma^2})} \gtrless -\lambda
Λ=f(z∣H0)f(z∣H1)=exp(−σ2A2)n=0∏N−1I0(σ22Azn)≷−λ
l n ( Λ ) = − A 2 σ 2 + ∑ n = 0 N − 1 l n [ I 0 ( 2 A z n σ 2 ) ] ≷ l n ( − λ ) → ∑ n = 0 N − 1 l n [ I 0 ( 2 A z n σ 2 ) ] ≷ l n ( − λ ) + A 2 σ 2 = T {\rm ln}(\Lambda)=-\frac{A^2}{\sigma^2}+\sum_{n=0}^{N-1}{{\rm ln}[I_0(\frac{2Az_n}{\sigma^2})]} \gtrless {\rm ln}(-\lambda)\rightarrow \sum_{n=0}^{N-1}{{\rm ln}[I_0(\frac{2Az_n}{\sigma^2})]} \gtrless {\rm ln}(-\lambda)+\frac{A^2}{\sigma^2}=T ln(Λ)=−σ2A2+n=0∑N−1ln[I0(σ22Azn)]≷ln(−λ)→n=0∑N−1ln[I0(σ22Azn)]≷ln(−λ)+σ2A2=T
在上面的表达式中
I
0
(
)
I_0()
I0()是一阶修正贝塞尔函数,它满足如下所似情况:
I
0
(
x
)
=
∑
k
=
0
∞
1
(
k
!
)
2
(
x
2
)
2
k
=
1
+
x
2
4
+
x
4
64
+
.
.
.
I_0(x)=\sum_{k=0}^{\infty}{\frac{1}{(k!)^2}(\frac{x}{2})^{2k}}=1+\frac{x^2}{4}+\frac{x^4}{64}+...
I0(x)=k=0∑∞(k!)21(2x)2k=1+4x2+64x4+...
当
x
x
x较小时,
I
0
(
x
)
≈
1
+
x
2
4
I_0(x)\approx 1+\frac{x^2}{4}
I0(x)≈1+4x2,对数函数有如下展开式:
l
n
(
1
+
x
)
=
x
−
x
2
/
2
+
x
3
/
3
−
x
4
/
4
+
.
.
.
{\rm ln}(1+x)=x-x^2/2+x^3/3-x^4/4+...
ln(1+x)=x−x2/2+x3/3−x4/4+...
所以
∑
n
=
0
N
−
1
l
n
[
I
0
(
2
A
z
n
σ
2
)
]
≈
∑
n
=
0
N
−
1
A
2
z
n
2
σ
4
\sum_{n=0}^{N-1}{{\rm ln}[I_0(\frac{2Az_n}{\sigma^2})]} \approx \sum_{n=0}^{N-1}{\frac{A^2 z_n^2}{\sigma^4}}
n=0∑N−1ln[I0(σ22Azn)]≈n=0∑N−1σ4A2zn2
所以,上面的LRT检验可以化简为:
∑
n
=
0
N
−
1
A
2
z
n
2
σ
4
≷
T
→
Υ
=
∑
n
=
0
N
−
1
z
n
2
≷
σ
4
T
A
2
=
T
′
(1)
\sum_{n=0}^{N-1}{\frac{A^2 z_n^2}{\sigma^4}} \gtrless T \rightarrow \Upsilon=\sum_{n=0}^{N-1}{z_n^2} \gtrless \frac{\sigma^4 T}{A^2}=T' \tag{1}
n=0∑N−1σ4A2zn2≷T→Υ=n=0∑N−1zn2≷A2σ4T=T′(1)
其中
Υ
\Upsilon
Υ表示各个采样点包络平方和,是该检验的充分统计量。从
(
1
)
(1)
(1)可知,利用信号包络非相干积累并进行检测,其最佳检测器近似等于对信号包络平方求和进行检测,即近似于平方率检测,一般衡量一个检测器性能优劣的指标是检测概率和虚警概率,两者计算方法如下:
P
F
A
=
∫
T
∞
f
(
Υ
∣
H
0
)
d
Υ
(2)
P_{FA}=\int_{T}^{\infty}{f(\Upsilon|H_0)}d{\Upsilon} \tag{2}
PFA=∫T∞f(Υ∣H0)dΥ(2)
P
D
=
∫
T
∞
f
(
Υ
∣
H
1
)
d
Υ
(3)
P_{D}=\int_{T}^{\infty}{f(\Upsilon|H_1)}d{\Upsilon} \tag{3}
PD=∫T∞f(Υ∣H1)dΥ(3)
上面的表达式可以这么理解:
P
F
A
P_{FA}
PFA表示在没有目标的时候(
H
0
H_0
H0情况下)检测到了目标,即在没有目标的情况下,样本的充分统计
Υ
\Upsilon
Υ超过了检测门限的情况,所以虚警概率是所有在
H
0
H_0
H0情况下
Υ
\Upsilon
Υ超过检测门限的情况之和。
P
D
P_{D}
PD表示在有目标的时候(
H
1
H_1
H1情况下)成功检测出来,所以在有目标的情况样本的充分统计超过门限的所有情况即为检测概率。所以从上面的分析可知,求某个检测器的性能,首先需要确定检测门限,检测门限由似然比检测器求得,如
(
1
)
(1)
(1)所示,该检测器的检测门限为
T
′
T'
T′。然后需要确定样本的充分统计量,如
(
1
)
(1)
(1)所示,该检测器的充分统计为
Υ
=
∑
n
=
0
N
−
1
z
n
2
\Upsilon=\sum_{n=0}^{N-1}{z_n^2}
Υ=∑n=0N−1zn2。最后需要确定样本充分统计的概率密度函数。知道了上面几个条件就可以求得检测器的检测概率和虚警概率。现在不进行具体求解。
从
(
2
)
(2)
(2)可以看出,
P
F
A
P_{FA}
PFA仅由干扰和噪声分布有关(仅与
H
0
H_0
H0情况下的信号分布有关),所以对于起伏目标和非起伏目标来说,
P
F
A
P_{FA}
PFA的最终表达式是一样的,可以推导得到高斯白噪声干扰下的平方率检波器的虚警概率为
P
F
A
=
1
−
I
(
T
N
,
N
−
1
)
I
(
μ
,
M
)
=
∫
0
μ
M
+
1
e
−
τ
τ
M
M
!
d
τ
P_{FA}=1-I(\frac{T}{\sqrt{N}},N-1) \\ I(\mu,M)=\int_{0}^{\mu\sqrt{M+1}}{\frac{e^{-\tau}\tau^M}{M!}}d\tau
PFA=1−I(NT,N−1)I(μ,M)=∫0μM+1M!e−ττMdτ
可以得到虚警概率受检测门限和积累样本数的影响,当积累样本数为1时,对应的虚警概率为:
P
F
A
=
e
−
T
(4)
P_{FA}=e^{-T} \tag{4}
PFA=e−T(4)
需要注意的是,在推导
(
4
)
(4)
(4)结果的过程中我们进行了如下变量代换:
z
n
′
=
z
n
/
σ
z_n^{'}=z_n/\sigma
zn′=zn/σ,其中
z
n
z_n
zn实际采样点的模值,
z
n
′
z_n^{'}
zn′是为了计算方便得到的检波器的输出,对于没进行归一化的采样值
z
n
z_n
zn和平方率检波来说有
P
F
A
=
e
−
T
/
σ
2
P_{FA}=e^{-T/\sigma^2}
PFA=e−T/σ2,其中
σ
2
\sigma^2
σ2表示全部噪声功率(包括I通道和Q通道)。对于恒虚警检测来说,一般手动设置
P
F
A
P_{FA}
PFA大小,根据
P
F
A
=
e
−
T
/
σ
2
→
T
=
−
σ
2
l
n
P
F
A
P_{FA}=e^{-T/\sigma^2}\rightarrow T=-\sigma^2{\rm ln}P_{FA}
PFA=e−T/σ2→T=−σ2lnPFA可以获取检测门限。但是需要注意的是上面的公式在推导过程中用到了很多近似,实际情况下检测门限可以建模为
T
=
α
σ
2
T=\alpha\sigma^2
T=ασ2,
α
\alpha
α是一个需要后续确定的一个参数,它显然与设置的虚警概率有关,当然还可能与其他参数相关。所以准确估计信号中的噪声功率对于设置检测门限是至关重要的。
目前雷达系统通常是对距离多普勒谱上每个单元进行检测的,判断每个单元上是否存在目标。所以距离多普谱上每个单元的数据就构成了采样样本
y
n
y_n
yn,在处理过程中并没有利用多个RD谱上同一个单元的数据进行积累,所以属于单样本检测问题,所以可以不考虑目标起伏模型,把所有目标当成非起伏目标进行处理。在处理过程中同样假设数据只受高斯白噪声干扰,所以其虚警概率模型仍然可以表示为
P
F
A
=
e
−
T
/
σ
2
P_{FA}=e^{-T/\sigma^2}
PFA=e−T/σ2。现在有一个问题就是设置检测门限时需要知道采样样本中的噪声功率,这个功率是很难从单个样本(待检测单元)中获取得到的(很难对待检测单元进行信号分离,获取噪声信号),所以需要利用其它单元对该噪声功率进行估计。此时对其他单元是有一定要求的:
(1)临近单元所含杂波的统计特性与待检测单元的一致;
(2)临近单元内不能包含目标。
高斯噪声背景下,距离多普勒谱上的检测单元服从指数分布,所以可以假设待检测单元
x
i
x_i
xi服从如下概率分布:
f
(
x
i
)
=
1
σ
2
e
x
p
(
−
x
i
/
σ
2
)
f(x_i)=\frac{1}{\sigma^2}{\rm exp}(-x_i/\sigma^2)
f(xi)=σ21exp(−xi/σ2)
利用相邻
N
N
N个单元上的数据进行噪声功率的估计,则这
N
N
N个单元的联合概率密度可以表示为:
f
(
x
)
=
∏
n
=
1
N
1
σ
2
e
x
p
(
−
x
i
/
σ
2
)
=
1
σ
2
e
x
p
[
−
(
∑
i
=
1
N
x
i
)
/
σ
2
]
f(\boldsymbol{x})=\prod_{n=1}^{N}{\frac{1}{\sigma^2}{\rm exp}(-x_i/\sigma^2)}=\frac{1}{\sigma^2}{\rm exp}{[-(\sum_{i=1}^N{x_i})/\sigma^2]}
f(x)=n=1∏Nσ21exp(−xi/σ2)=σ21exp[−(i=1∑Nxi)/σ2]
令
∂
l
n
[
f
(
(
x
)
)
]
∂
σ
2
=
0
\frac{\partial{\rm ln}[f(\boldsymbol(x))]}{\partial\sigma^2}=0
∂σ2∂ln[f((x))]=0,可以求得
σ
^
2
=
1
N
∑
i
=
1
N
x
i
\hat{\sigma}^2=\frac{1}{N}\sum_{i=1}^{N}x_i
σ^2=N1i=1∑Nxi
所以门限可以通过下式进行估计
T
^
=
α
σ
^
2
=
α
N
∑
i
=
1
N
x
i
\hat{T}=\alpha\hat{\sigma}^2=\frac{\alpha}{N}\sum_{i=1}^{N}x_i
T^=ασ^2=Nαi=1∑Nxi
令
z
i
=
α
x
i
/
N
z_i=\alpha x_i/N
zi=αxi/N,则
T
^
=
∑
i
=
1
N
z
i
\hat{T}=\sum_{i=1}^{N}z_i
T^=∑i=1Nzi,其中
z
i
z_i
zi是
x
i
x_i
xi的函数,现在需要利用
x
i
x_i
xi的概率密度函数求
z
i
z_i
zi的概率密度函数,这是求随机变量函数的概率密度函数的问题,可以直接得到
z
i
z_i
zi的概率密度函数为
f
(
z
i
)
=
N
α
σ
2
e
x
p
(
−
N
z
i
α
σ
2
)
f(z_i)=\frac{N}{\alpha\sigma^2}{\rm exp}(\frac{-Nz_i}{\alpha\sigma^2})
f(zi)=ασ2Nexp(ασ2−Nzi)
T
^
=
∑
i
=
1
N
z
i
\hat{T}=\sum_{i=1}^{N}z_i
T^=∑i=1Nzi的概率密度函数直接求取比较不方便,可以利用分布的特征函数来间接求取
T
^
\hat{T}
T^的概率密度函数。这边不进行详细推导,直接给出其概率密度函数:
f
(
T
^
)
=
(
N
α
σ
2
)
N
T
^
N
−
1
(
N
−
1
)
!
e
−
N
T
^
/
α
σ
2
f(\hat{T})=(\frac{N}{\alpha\sigma^2})^N\frac{\hat{T}^{N-1}}{(N-1)!}e^{-N\hat{T}/\alpha\sigma^2}
f(T^)=(ασ2N)N(N−1)!T^N−1e−NT^/ασ2
根据估计的门限值
T
^
\hat{T}
T^算出来的虚警概率可以表示为
P
F
A
=
e
−
T
^
/
σ
2
P_{FA}=e^{-\hat{T}/\sigma^2}
PFA=e−T^/σ2,所以
P
F
A
P_{FA}
PFA的期望
P
ˉ
F
A
\bar{P}_{FA}
PˉFA可以表示为:
P
ˉ
F
A
=
∫
−
∞
∞
e
−
N
T
^
/
α
σ
2
f
(
T
^
)
d
T
^
=
(
1
+
α
/
N
)
−
N
\bar{P}_{FA}=\int_{-\infty}^{\infty}{e^{-N\hat{T}/\alpha\sigma^2}f(\hat{T})d\hat{T}}=(1+\alpha/N)^{-N}
PˉFA=∫−∞∞e−NT^/ασ2f(T^)dT^=(1+α/N)−N
在上面的求解过程中用到了随机变量函数的期望的求解方法。所以可以求解得到:
α
=
N
(
P
ˉ
F
A
−
1
/
N
−
1
)
\alpha=N({\bar{P}_{FA}}^{-1/N}-1)
α=N(PˉFA−1/N−1)
研究表明,对于Swerling1型目标的一个检测单元数据在给定门限
T
^
\hat{T}
T^下,检测概率可以表示为
P
D
=
e
x
p
[
−
T
^
/
(
1
+
χ
)
]
P_D={\rm exp}[-\hat{T}/(1+\chi)]
PD=exp[−T^/(1+χ)],因为在距离多普勒谱上面进行目标检测,每次只包含一个检测样本,所以不存在目标起伏的问题,所以上面的结论也适用Swerling2类型目标,所以检测概率的数学期望可以表示为:
P
ˉ
D
=
∫
−
∞
∞
e
−
T
^
/
(
1
+
χ
)
f
(
T
^
)
d
T
^
=
(
1
+
α
N
(
1
+
χ
)
)
−
N
\bar{P}_D=\int_{-\infty}^{\infty}{e^{-\hat{T}/(1+\chi)}f(\hat{T})d\hat{T}}=(1+\frac{\alpha}{N(1+\chi)})^{-N}
PˉD=∫−∞∞e−T^/(1+χ)f(T^)dT^=(1+N(1+χ)α)−N
需要理解的是上面的结果是基于以下假设进行的,即I/Q通道的噪声均服从高斯分布,平方率检波器,Swerling1 或 Swerling2型目标,并且只有一个待检测单元的数据。
需要理解上面为什么计算检测概率和虚警概率的期望。可以看出这两者都是检测门限
T
^
\hat{T}
T^的函数,如果
T
^
\hat{T}
T^是确定的,那么
P
D
P_D
PD和
P
F
A
P_{FA}
PFA也将是确定的,但是实际上
T
^
\hat{T}
T^是服从一定分布的,但是
P
D
P_D
PD和
P
F
A
P_{FA}
PFA与
T
^
\hat{T}
T^的关系却是确定的,所以这样的问题,显然需要采用统计学的方法来解决。