第四章微弱信号检测技术

1、第四章 微弱信号检测技术4.1 被动信号检测被动检测是一种常用的检测系统，它已广泛应用于水下引信信号检测及其它工业领域。在被动信号检测中，常用的时域检测方法有以下几种：宽带检测、相干检测、频率随机分布正弦信号的检测技术、时域同步平均检测与波形恢复技术、相关技术等等；而在频域的检测方法主要是基于FFT算法的谱分析技术。4.1.1宽带检测在有些应用场合，干扰噪声和输入信号都是一有限长的限带零均值的高斯分布随机过程，在此情况下一般使用宽带检测技术。4.1.1.1最佳宽带检测器最佳宽带检测器的结构框图如下：Z图4.1 在高斯噪声中检测高斯信号的最佳系统结构图4.1中是信号的功率谱密度，是干扰噪声的功率

2、谱密度。而表示预选滤波的频率响应。当信号和噪声都是限带高斯分布白噪声时，信号和噪声的差别是信号和噪声的功率级不同，为常值，最佳检测器是一个平均功率检测器。从理论上说无论噪声多强，信号多弱，只要他们是平稳的，且他们的方差可准确求出来，那么总可通过比较N和N+S,发现信号。如果过程是各态遍历的,那么方差可通过下式计算出来。 (4.1.1)不难看出，由于截取的样本时间是滑动的，从而图4.1可简化为平方积分系统。由于截断T不是无限长的，所以输出并不等于,而是随t在的均值附近起伏。对于限带白谱：起伏的存在将掩盖信号加噪声（H1）与噪声（H0）的差别。所以系统的信噪比计算公式如下： (4.1.2)在各态遍

3、历条件下，T越长系统的最佳性越好。 当信号和噪声的功率谱不是白谱时，可利用的信息不仅有能量差异，而且还有谱形状的差异。此时的预选滤波器的传输函数的幅度特性如下： (4.1.3)在小输入信噪比情况下： (4.1.4)式(4.1.4)所描述的滤波器称为厄卡特滤波器。若假设信号和噪声有相同的谱形状，则： (4.1.5)上式所描述的是一个白化滤波器，信号和噪声通过后一律变成白噪声。非白谱小信号情况下，其相当于一个白化滤波器和一个匹配滤波器的级联。当信号与噪声有相同形状功率谱时，匹配网络的频率传输函数等于常数，厄卡特滤波器退化为一个白化滤波器,此时虽然不能提高系统输出端的信噪比，但却通过改善噪声谱的形状

4、（白化）提高了系统的等效噪声谱宽。4.1.1.2实用宽带检测器在实际应用中，由于信号和噪声的功率谱很难知道，因此预选滤波器一般没有白化和对信号进行匹配的能力，因此它对系统的输出信噪比影响很小。在实用的宽带检测系统中，主要研究的是宽带能量检测器，对这种接收机一般以系统的输出信噪比的大小或系统处理增益作为衡量系统性能的指标。宽带能量检测器在判决检测前都相应有一个等效积分器，为使讨论具有一般性，可将积分器理解为一个低通滤波器，积分器的传输函数记为H(w)，输入端Y处与输出端Z处的信噪比可按如下公式计算： (4.1.6) (4.1.7)它们和系统参数的关系如下： (4.1.8)其中Wy为积分器等效谱宽

5、；T为积分器的等效积分时间，可表示为： (4.1.9) (4.1.10)式(4.1.9)中为Y(t)的归一化自相关函数，有。可以看出：积分器输出信噪比与积分器的等效积分时间、积分器输入噪声过程的等效谱宽和积分器输入端的输入信噪比有关。对于理想积分器，等效积分时间就等于积分器的积分时间，也等效于被观察信号的作用时间，而对于RC积分器的等效积分时间T=2RC。噪声过程Y(t)的等效噪声谱宽就是其频带内的功率谱密度与的比值。宽带能量检测器（平方律检波器）是很难实现的，所以在实际的接收机中常用易实现的线性检波器代替它，但相应的输出信噪比有所下降。在小输入信噪比条件下，下降只有1.1416倍(0.57d

6、B)，显然是微不足道的，在大输入信噪比条件下，输出信噪比的损失明显增大，但在这种条件下不造成检测的困难。4.1.2相干检测当信号为一弱周期性信号、且伴随着很强的干扰噪声时，相干检测是一种常用的信号检测方法，相干检测的原理框图如图4.2所示：图4.2 相干检测原理图其中SC表示输入信号通道；RC表示参考信号通道；PSD表示相敏检波器；LPF表示低通滤波器。相敏检波器实际上是一个乘法器。假设信号的角频率为，噪声的角频率为，信号与参考信号的相位差为，而噪声与参考信号的相位差为，且，则的输出可表示为： (4.1.11) (4.1.12) (4.1.13)当其通过LPF时，只要LPF的截止频率（或），则

7、噪声分量被滤除。其信号形式如下： (4.1.14)当输入信号的频率有的偏差时，则LPF的输出为： (4.1.15)由上式可看出，当,同时为零时，相干检测的输出取得最大值，因此在相干信号检测过程中，参考信号通道应包括频移或延时环节，可完成从的相移变化。4.1.3频率随机分布正弦信号的检测技术如果信号的频率未知，则可使用扫频检测器确定被检测信号的频率，其原理框图示于图4.3：图4.3 扫频检测器原理图中压控振荡器产生线性调频振荡信号，其在某时刻的输出为：e0=E0cosw0t，当时，差频窄带滤波器的输出为： (4.1.16)如果在噪声背景中存在角频率为的正弦信号，扫频检测器做出存在正弦信号的判决。

8、4.1.4 时域同步平均检测技术 时域同步平均方法是噪声中提取周期性信号的有效方法，是一种积累平均抗干扰检测过程，也称相干检波。对周期为T的信号以KT为间隔截取M断，然后将各段信号对应的点相加后取算术平均。设输入X(t)为信号S(t)和噪声N(t)的合，即，则平均后的信号为，从此式可以看出，平均之后的噪声幅度是平均前信号中噪声幅度的倍，因此信噪比提高倍（有效值），这就是时域同步平均的法则。M次平均后获得的处理增益为：GM=10log10M。显然信号时域同步平均检测把原始输入信号的功率信噪比提高了M倍，即经过多次平均后噪声逐步被抑制，周期信号被显现出来。4.1.5 相关技术相关是测量两个信号的相

9、似性。相关分析分为自相关分析和互相关分析两种，其本质是一种线性滤波，是抑制随机干扰，提高信噪比的一种方法。相关技术在微弱信号检测技术中的主要应用是寻找隐匿于随机噪声中的规律信号。设s(t)表示隐匿在噪声中的正弦信号，N(t)表示干扰噪声，且二者不相关，则其自相关函数如下： 关于相关检测的原理示于图4.4，在实际的检测系统中，常被用来检测周期信号的周期。图4.4 相关检测原理在小型近感探测系统中实现相关器具有实用价值，基本方法包括软件实现技术和硬件实现技术。相关器的软件实现：对于样本长度为N的时间序列信号x(n)，相关器可用下述公式通过编程实现： 相关器的硬件实现：相关器的硬件实现如图4.5所示

10、：图4.5 横向滤波器横向滤波器的输出为：。已出现一些由大规模集成电路构成的相关器，其性能优于横向滤波器。4.1.6 频域检测综述（增加） 谱估计是信号处理的一个重要方面。它在不同领域都有广泛应用。当信号的时域特征和干扰背景有明显差异时，谱分析方法是检测微弱信号的有效方法。随着大规模集成电路的飞速发展，使得小型化的谱分析器实用性越来越好。传统的谱分析方法是以傅立叶分析为基础的，由于FFT算法的出现，使得傅立叶分析的性能大大提高，并迅速进入工业应用领域。传统的谱估计方法主要有自相关法和周期图法，后者是最常用的估计方法。在用谱估计方法检测微弱信号时又细化为：功率谱估计、倒谱估计、zoomFFT分析

11、等。由于短时间信号的FFT分析的频率分辨力很低，而很多实际信号又只能在短时间内视为平稳过程，70年代后导致出现非傅氏方法的现代谱分析，并成为研究热点。这些方法可把观测时间间隔内的数据外推到观测间隔之外，增大了有效观测时间，因此分辨率得到了很大提高。常用的现代谱分析方法主要包括：自回归(AR)谱估计、滑动平均(MA)谱估计、自回归滑动平均(ARMA)谱估计、最大熵谱估计、极小方差谱估计等，但因它们运算量大，在小型检测系统中目前多用基于FFT算法的谱分析技术。由于实际检测的信号不满足高斯、平稳信号的特点、相应的分析系统很难满足线性系统、最小相位系统、因果系统，因此针对非高斯、非平稳信号，以及非线性

12、、非因果、非最小相位系统，人们提出了现代信号处理技术的方法，并成为现代研究热点，主要包括：短时傅立叶变换分析、时频分析、小波分析、高阶统计量分析、循环平稳信号分析、Gabor变换分析、分数阶FT&Radon-Wigner变换分析等技术。4.2 主动信号检测 在被动检测中，被检测的信号是未知的；在主动检测中，信号是已知的或至少知道信号的大部分特性，因此主动检测技术有别于被动检测技术。4.2.1在噪声背景中检测确知信号在主动检测系统中，讨论在噪声背景中检测确知信号需要从最大似然比准则出发。假设信号加噪声的样本矢量为，其复包络矢量为，和分别为信号与噪声的样本矢量和复包络矢量，为零均值复高斯矢

13、量，其协方差矩阵为。 有信号时为： 其概率分布密度函数为： (4.2.1)无信号时为： 其概率分布密度函数为： (4.2.2)则相应的对数似然比为：上式中的第一项与r无关，所以只取第二项作为假设检验统计量，判决方程为： (4.2.3)因为是确知的，为实数，上式可简写为： (4.2.4)4.2.1.1在白噪声背景中检测确知信号：在白噪声背景下，假设噪声样本间相互独立，其方差相等，归一化后取为单位矩阵I。则其检验统计量和判决方程为： (4.2.5)这就是相关接收机的离散型表达式。它计算接收样本矢量和参考样本矢量的内积作为检验统计量。这种离散接收机的框图示于图4.6。图4.6 离散型相关接收机框图信

14、号R(t)先经正交差频低通滤波后得到复包络，然后在信号到达时刻进行采样，在整个信号持续时间内进行存贮，得到整个信号的样本矢量。因假设信号是已知的，中包含的信号必然与相同。将与内积，可得到检验统计量e(tm)，将它和门限相比较即可判断有无目标。在实际中，常用匹配滤波器来代替互相关器，匹配滤波器的脉冲响应为。匹配滤波器的输入高斯噪声加有规信号时，经匹配滤波器这一线性系统，输出仍为高斯噪声加有规信号。如果输入信号的能量为E，输入白噪声的功率谱密度为N0/2，则可给出接收机的工作特性曲线如图4.7所示。图中给出在不同的输出信噪比d=2E/N0时，虚警概率与检测概率的关系。在主动信号检测的相关接收中，为

15、了简化处理运算和克服水声信号的不平稳性，有时采用极性相关器来代替线性相关器。在极性相关器中，输入信号先经强限幅取极性后，再计算相关系数，这相当于1比特量化相关。图4.7 接收机工作特性曲线4.2.1.2 在非白高斯噪声中检测确知信号在非白高斯噪声情况下，相关矩阵不再是单位阵。输入样本矢量与相关阵相乘，使噪声白化后再计算检验统计量。实际分析时，将认为是正定的，常把表示成两个三角阵的乘积（），则其检验统计量和判决方程变为： (4.2.6)由于c是三角阵，它的计算相当于一个物理可实现的滤波器。在的每一个分量的输入过程中，即可计算，特别方便级联实时处理。用三角矩阵解相关的最佳接收机如图4.8所示:图4

16、.8 基于三角阵解相关的最佳接收机框图在连续时间函数情况下，假设输入信号的功率谱为: (4.2.7)则对非白噪声进行白化的传输函数为，白化后的信号频谱变为，因此匹配滤波器的传输函数为：则系统的传输函数为： (4.2.8)在主动信号检测时，s(t)为窄带信号，分布在频带内。弱在内是均匀的，则为矩形，相关函数为型。实际上，海洋噪声或舰艇噪声的功率谱密度有每倍频程约56分贝的下降斜率。对这种情况进行白化处理，其增益不到0.5分贝。一般说来，在主动信号频带内，噪声的频谱如无明显波动，进行修正意义不大。4.2.2在噪声背景中检测参数未知信号在主动信号检测时，事先没有关于目标的先验知识，回波到达时间、初相

17、与多普勒频移均是未知的，而只有回波信号的形式是已知的，并与反射信号相同。这种情况下，主动信号检测常使用正交接收机。4.2.2.1频率已知初相位未知的信号检测：上节讨论确知信号检测时用最佳接收机，因初相位已知，而使相关器的输出为实数。而初相位未知时，相关器的输出为复数。从检测的观点看，可将相关器输出的模做为检验统计量与门限比较，判决有无目标，判决方程为： (4.2.9)对于非白噪声干扰，检测初相位未知信号的最佳接收机可由预白化处理和正交接收机构成，如图4.9所示，图中为将探测距离范围分解的许多距离元，为可能的多普勒频移范围分解的许多频率元。图4.9 初相未知时的最佳接收机是以为中心频率对r(t)

18、取复包络后的样本，而参考信号是以为中心频率对s(t)取复包络后得到的。这样就设置了频率分辨元，距离分辨元是通过对发射信号延时来设置的。这种接收机每次只能完成对一个距离元、一个频率元得判决。4.2.2.2频率和初相位未知的信号检测：由于频率未知要在频率轴上进行搜索，信号到达时间未知要在距离元上搜索，这显然太复杂了。在某些情况下，由于探测目标是移动的，因此回波信号将产生多普勒频移。在此情况下，可计算检验统计量的DFT，输出最大的那个频率即代表目标的多普勒频移。为方便起见，这里讨论白噪声情况，这时检验统计量为： (4.2.10)当没有多普勒频移时，应取为，即等于反射信号得包络。当回波有多普勒频移时，应取为，当多普勒频移小于宽带多普勒容限，即多普勒效应引起的色散效应可忽略时（复包络不变），则有： (4.2.11)式(4.2.10)变为： (4.2.12)上式可视为是计算在频率上的分量。通