一种功率谱和统计融合的信号检测方法与流程

本发明属于电子信号侦察领域，具体涉及一种功率谱和统计融合的信号检测方法。

背景技术：

电子信号侦察领域中信号检测是电子侦察信号处理流程中的首要步骤，对于低负信噪比信号或者噪声属于非白类型的检测是信号检测过程中需要面临的一种问题。目前常用的信号检测方法有能量检测法、匹配滤波检测、循环平稳检测法等。

能量检测法应用广泛，但是对低负信噪比信号能量检测接收机往往无法达到令人满意的检测性能，同时能量检测法对于噪声功率变化敏感，当噪声功率变化时检测性能会下降。

匹配滤波检测是一种能够使输出信噪比达到最大的检测器，这种检测器在已知信号的情况下是一种最优检测器，但是匹配滤波检测需要知道信号的先验信息，对于非合作的信号检测这种方法并不适用。

循环平稳检测法利用信号的循环平稳特性，信号的循环谱在循环频率处会有峰值出现，循环平稳检测检测性能好，但是信号循环谱的计算量大是循环检测存在的问题。

因此在非合作的条件下寻到一种对于低负信噪比信号的检测方法成为了现今需要解决的问题。

技术实现要素：

本发明的目的在于提供一种功率谱和统计融合的信号检测方法，适用于对低负信噪比信号进行检测。

实现本发明目的的技术解决方案为：一种功率谱和统计融合的信号检测方法，利用welch功率谱估计方法的优点，通过welch方法提高信号功率谱的信噪比，弥补了低负信噪比的缺陷，通过统计估计方法完成噪声功率估计，在此基础上完成信号检测，具体步骤如下：

步骤1：在盲检测的情况下，对于信号的先验知识知道的很少或者完全不知道，此时令接收到的信号x(n)为

其中s(n)为辐射源发射的信号，高斯白噪声z(n)符合高斯分布n(0,σ²)，h1表示有信号，h0表示无信号，n∈{1,2,3,…n}，n为采样点个数；

步骤2：计算接收到的信号x(n)的welch功率谱；

步骤3：将得到的功率谱φx(k)分成g段，每段q个离散频点：

假设分段后每段的功率谱表示成φx,g(q)，φx,g(q)定义为

φx,g(q)＝φx(q+qg),k1＝1,2,3,...,q；g＝0,1,2,3,...,g-1

其中g表示被划分的第g个频段。

步骤4：选择观测频段g0的功率谱对与其相邻的g0-1和g0+1两个频段的功率谱分别进行统计分析，得到其分布密度曲线，估计噪声功率；

步骤5：设置检测门限，进行信号检测；

设置检测门限，将g0频带内的功率谱与门限进行比较，大于门限的判决为存在信号，小于门限的判决为不存在信号。

本发明与现有技术相比，其显著优点在于：

(1)算法原理简单，计算量小，易于工程实现。

(2)非合作条件下能够完成信号检测。

(3)低信噪比下检测性能较好。

附图说明

图1为本发明的算法流程图。

图2为本发明的功率谱估计对比图，其中横坐标表示频率，单位mhz。

图3为本发明实施例-20db信噪比时临近频带分布密度曲线图，其中横坐标表示观测值大小。

图4为本发明实施例-20db信噪比条件下功率谱检测效果图，其中横坐标表示频率，单位mhz。

图5为本发明与能量检测法的信噪比-检测概率曲线图。

图6为本发明实施例实际信号检测效果图，其中横坐标表示频率，单位mhz。

具体实施方式

下面结合附图对本发明作进一步详细描述。

结合图1，本发明所述的功率谱和统计融合的信号检测方法，步骤如下：

步骤1：在盲检测的情况下，对于信号的先验知识知道的很少或者完全不知道，此时令接收到的信号x(n)为

其中s(n)为辐射源发射的信号，高斯白噪声z(n)符合高斯分布n(0,σ²)，h1表示有信号，h0表示无信号。n∈{1,2,3,…n}，n为采样点个数。因为信号大多都不是确定信号，不只是含有噪声，其中信号的出现与否也含有随机性，所以这类信号的数学模型应该是随机序列。随机信号的每一次实现都是离散时间序列，其能量并不是有限的，因此其dtft并不存在，但是却有有限的平均功率，可以用功率谱密度来表示，因此对信号的讨论也是对随机序列的讨论。

步骤2：计算接收到的信号x(n)的welch功率谱

一般随机信号的功率谱表示如下式，

其中，k表示功率谱的离散频点，k＝1,2,3,...,n，φx(k)表示对信号进行功率谱估计后得到的功率谱。为了提高信噪比，压制噪声干扰，本发明采用welch功率谱估计的方法估计功率谱。welch法计算功率谱的详细步骤如下：

步骤2-1：对信号x(n)分段，xl(m)是对x(n)进行分段后得到的第l段的数据，每段数据和后一段数据有50％的重叠部分，以减少观测时间。

其中m表示第l段数据中第m个样点，每段数据共m个样点，x(n)的整个观测时间数据被分为l段。

步骤2-2：计算每段数据的加窗功率谱

其中k表示功率谱的离散频点，j为虚数，w(m)为窗函数，p是窗函数的功率，表示为

步骤2-3：对各段加窗功率谱求平均得到welch功率谱。

信号x(n)的welch功率谱表示为

welch方法对每段数据加窗后进行fft，取模的平方然后平均作为功率谱，这样减小了功率谱估计的方差。在信噪比为-5db的情况下，经过步骤2处理后得到的信号功率谱与处理前功率谱对比图如图2所示。

步骤3：将得到的功率谱φx(k)分成g段，每段q个离散频点。

假设分段后每段的功率谱表示成φx,g(q)，φx,g(q)定义为

φx,g(k1)＝φx(k1+qg),k1＝1,2,3,...,q；g＝0,1,2,3,...,g-1

其中g表示被划分的第g个频段。

步骤4：选择观测频段g0的功率谱φx,g(q)，对与其相邻的g0-1和g0+1两个频段的功率谱分别进行统计分析，得到其分布密度曲线：

假设观测频段为g0，对g0-1和g0+1两个频段进行统计分析，分别做这两个频段的分布直方图，由分布直方图画分布曲线。在只存在噪声的情况下，根据中心极限定理，的渐进分布趋近于高斯分布，可以证明这种情况下的均值为噪声方差σ²，因此分布峰值对应的横坐标即为对应的噪声功率。因为实际应用中可能会出现噪声非白的情况，我们假设噪声在g0频段及相邻频段内是近似白化的，用相邻频段内的噪声功率表示g0频段内的噪声功率。在由统计直方图得到分布密度曲线的时候本文将原始的直方图进行了平滑，将两临近频带通过直方图估计得到的噪声功率和进行比较，较小的作为噪声的功率σ²。

仿真设置信号为bpsk信号，信噪比为-20db，载频5mhz，采样率12mhz，信息速率为20kbps，此时临近频带的分布密度曲线如图3所示，其中横坐标x表示所统计的功率谱的样本值。

步骤5：设置检测门限，进行信号检测。

将估计得到的噪声功率的提高3db作为检测门限，将感兴趣频带内的功率谱与门限进行比较，大于门限的判决为存在信号，小于门限的判决为不存在信号。-20db信噪比条件下功率谱与检测门限如图4所示。

对本发明进行蒙特卡洛仿真的到的信噪比检测概率曲线如图5所示。

实际接收机对盲采得到的信号进行判决检测的功率谱图如图6所示，图中横线为检测门限。