本文主要内容来自文献:
H.Urkowitz (1967). “Energy Detection of Unknown Deterministic Signals.” Proceedings of the IEEE 55(4)
1、白噪声中的检测:低通过程
1.1 系统模型
信号处理过程如图1所示。
1) H 0 \mathcal H_0 H0:
y ( t ) = n ( t ) E [ n ( t ) ] = 0 n o i s e P S D = N 0 ( 单 边 功 率 谱 ) n o i s e B W = W \begin{aligned} &y(t)=n(t)\\ &{\rm E}[n(t)]=0\\ &{\rm noise\ PSD}=N_{0}{\rm (单边功率谱)}\\ &{\rm noise\ BW} = W \end{aligned} y(t)=n(t)E[n(t)]=0noise PSD=N0(单边功率谱)noise BW=W
2) H 1 \mathcal H_1 H1:
y ( t ) = s ( t ) + n ( t ) E [ s ( t ) + n ( t ) ] = s ( t ) \begin{aligned} &y(t)=s(t)+n(t)\\ &{\rm E}[s(t)+n(t)]=s(t)\\ \end{aligned} y(t)=s(t)+n(t)E[s(t)+n(t)]=s(t)
为了方便计算检测概率和虚警概率,我们将
V ′ = 2 N 0 ∫ 0 T [ y ( t ) ] 2 d t (1) \tag{1} V'=\frac{2}{N_{0}}\int_0^T[y(t)]^2dt V′=N02∫0T[y(t)]2dt(1)作为检测变量。
进一步,我们知道,对于带宽为 W W W的信号,根据低通抽样定理,抽样速率为 2 W 2W 2W,那么对于持续时间为 T T T的信号,就一共需要 2 T W 2TW 2TW个样值来表示它。如果是窄带信号,则对同相和正交分量,分别用 W W W的速率抽样。
1.2 判决变量推导
- H 0 \mathcal H_0 H0
由于通过噪声预滤波,因此 n ( t ) n(t) n(t)为低通信号。我们用 2 W 2W 2W速率对 n ( t ) n(t) n(t)进行抽样,显然可以得到抽样后的样值 { n i = n ( i 2 W ) , i = 1 , 2 , … N } \{ n_i=n(\frac{i}{2W}),i=1,2,\ldots N\} {
ni=n(2Wi),i=1,2,…N}。如果利用低通滤波器 H ( f ) = 1 2 W R e c t [ f 2 W ] H(f)=\frac{1}{2W}{\rm Rect}\left[ \frac{f}{2W}\right] H(f)=2W1Rect[2Wf]对抽样后的样值进行重建,由于 h ( t ) = s i n c ( 2 W t ) h(t)={\rm sinc}(2Wt) h(t)=sinc(2Wt),因此重建后的信号为
n ( t ) = ∑ i = − ∞ ∞ n i s i n c ( 2 W t − i ) , (2) \tag{2} n(t)=\sum_{i=-\infty}^\infty n_i{\rm sinc}(2Wt-i), n(t)=i=−∞∑∞nisinc(2Wt−i),(2)显然,这里 n i ∼ N ( 0 , σ n 2 ) n_i\sim {\mathcal N}(0,\sigma_n^2) ni∼N(0,σn2),且 σ n 2 = N 0 W \sigma_n^2=N_0W σn2=N0W。
进一步,根据性质
∫ − ∞ ∞ s i n c ( 2 W t − i ) s i n c ( 2 W t − k ) d t = { 1 2 W , i = k 0 , i ≠ k \int_{-\infty}^{\infty}{\rm sinc}(2Wt-i){\rm sinc}(2Wt-k)dt=\left\{\begin{aligned} \frac{1}{2W}, &\ i=k\\ 0,&\ i \ne k\\ \end{aligned}\right. ∫−∞∞sinc(2Wt−i)sinc(2Wt−k)dt=⎩⎨⎧2W1,0, i=k i=k可以得到
∫ − ∞ ∞ n 2 ( t ) d t = 1 2 W ∑ i = − ∞ ∞ n i 2 . (3) \tag{3} \int_{-\infty}^{\infty}n^2(t)dt=\frac{1}{2W}\sum_{i=-\infty}^{\infty}n_i^2. ∫−∞∞n2(t)dt=2W1