本发明涉及通信信号检测技术领域，尤其涉及基于双门限判决的窄带干扰位置检测方法及系统。
背景技术
为了实现无线通信的抗干扰，干扰检测技术是信号检测领域中的重要环节，目的是判断干扰是否存在，为后续的抗干扰措施提供必需的先验信息。
典型的干扰检测技术包括：能量检测算法、匹配滤波检测算法、循环平稳检测方法、协方差矩阵分解算法等。
能量检测算法，首先根据干扰信号存在与否做两种检验假设（H1/H0），然后计算接收信号的检测统计量，与设定的门限值进行比较，判断干扰信号是否存在。能量检测作为一种二元检测方法，不需要干扰信号的先验信息，是一种较为简便有效的盲干扰检测方法。能量检测不能给出干扰信号的频点位置以及属于什么类型的干扰，受背景噪声的影响较大，若噪声不平稳则检测性能会下降。
匹配滤波检测算法，根据主用户信号的某些先验信息，利用线性滤波器与待检测信号完成在时域以及频域上的同步，从而对主用户信号进行解调。如果主用户信号的先验信息未知，则该检测方法无效。
循环平稳检测方法，利用信号的循环平稳特性进行检测，主用户信号和有意干扰信号一般都存在循环频率，属于循环平稳信号，在非零的循环频率处，具有明显的谱相关函数幅值，并且二者的循环频率相异。噪声的自相关函数不具有周期性，不存在循环平稳特性，其谱相关函数在非零循环频率处幅值为零。因此，可以利用干扰信号与主用户信号的循环频率相异特性，实现干扰信号检测和分离。该方法是一种有效的干扰检测算法，但算法复杂度过高。
协方差矩阵分解算法，首先对接收信号进行采样，并求得自相关函数，根据得到的自相关函数组成协方差矩阵，对协方差矩阵进行特征分解，求得最大特征值λmax和最小特征值λmin。将λmax/λmin作为检测统计量，若大于门限值，则存在主用户信号，否则不存在主用户信号。该方法局限于相关性良好的主用户信号，对独立信号的检测效果欠佳。
上述干扰检测方法，对于干扰信号的频点位置识别效果均欠佳，能量检测是一种适用范围较广的信号检测方法，该算法实现简便并且在一定条件下可以较为准确的检测出干扰信号的有、无。
单纯的能量检测存在两个缺点：一方面，无线通信系统存在远近效应，单纯的能量检测不能有效区分干扰信号和非干扰信号，造成干扰误判；另一方面，作为二元检测，单纯的能量检测不能准确地定位窄带干扰的具体频点位置。
技术实现要素：
因此，本发明的目的在于提供一种基于双门限判决的窄带干扰位置检测方法及系统，应对无线通信的远近效应，在接收信号能量大范围波动条件下，对宽带无线通信系统中干扰信道的具体频点位置进行定位，采用连续均值检测量估计和双门限检测算法进行窄带噪声能量估计和干扰位置检测，从而获得准确的信道干扰位置检测结果。
为了实现上述目的，本发明的一种基于双门限判决的窄带干扰位置检测方法，包括以下步骤：
捕获多个窄带信道的并行信号，获得工作频带内各窄带信道的接收信号能量测量值RSSIf和解调信干噪比测量值SINRf；
对各窄带信道的接收信号能量测量值RSSIf进行独立估计，并进行连续均值滤波，得到稳定的接收信号能量水平估计值，作为当前信道的干扰检测量估计值RSSIf_est；
对各窄带信道的解调信干噪比测量值SINRf，进行独立估计，并进行连续均值滤波，得到当前信道的解调信干噪比估计值SINRf_est；
将各个信道的干扰检测量估计值RSSIf_est与预设的第一门限值进行比较，判断干扰信号是否存在；将解调信干噪比估计值SINRf_est与预设的第二门限值进行比较，判断干扰信号是否为有效信号；将存在干扰信号，且干扰信号为有效信号的信道作为干扰位置检测结果输出。
进一步优选的，所述预设的第一门限值为能量检测阈值RSSI_thr，所述预设的第二门限值为信干噪比检测阈值SINR_thr。
进一步优选的，所述各个信道的干扰检测量估计值RSSIf_est，采用如下公式表示：
式中，RSSIf_est为按频点统计的各个信道的干扰检测量估计值，Wdetect为检测窗口长度，N_pulse为每单位检测窗口内的脉冲数量；RSSIif为第f个频点的第i个脉冲的接收信号能量测量值。
进一步优选的，所述能量检测阈值按照如下公式进行计算：
式中，RSSI_thr为能量检测阈值；Pfa为虚警概率，取值为0.001；N为显著统计量的统计次数；RSSI_signal为接收信号额定灵敏度；Q函数表达式为。
进一步优选的，将当前信道的干扰检测量估计值RSSIf_est与预设的第一门限值进行比较，包括：
若干扰检测量估计值RSSIf_est大于能量检测阈值RSSI_thr则判定存在干扰信号，否则判定为不存在干扰信号。
进一步优选的，在判断干扰信号是否为有效信号时，包括将解调信干噪比估计值SINRf_est与预设的信干噪比检测阈值SINR_thr进行对比，超过信干噪比检测阈值SINR_thr的信号作为有效干扰信号，低于信干噪比检测阈值的信号为无效干扰信号。
进一步优选的，还包括对干扰位置检测结果进行统计滤波输出，所述统计滤波输出包括以下步骤：
在检测周期内设定多个检测窗口，对多次干扰位置检测结果进行统计滤波，采用相对统计量进行异常值过滤。
本发明还提供一种基于双门限判决的窄带干扰位置检测系统，包括AD变换模块、RSSI测量模块、RSSI估计模块、SINR测量模块、SINR估计模块、双门限干扰位置检测模块；
所述AD变换模块用于接收射频前端输入的宽带无线信号，并输出数字域的宽带无线信号；
所述RSSI测量模块，用于对数字域的宽带无线信号进行多信道并行能量检测，获得工作频带内多个窄带信道的接收信号能量测量值RSSIf。
RSSI估计模块，用于对各窄带信道的接收信号能量测量值RSSIf进行独立估计，并进行连续均值滤波，得到稳定的接收信号能量水平估计值，作为当前信道的干扰检测量估计值RSSIf_est；
SINR测量模块，用于对各信道输入的低速脉冲信号进行信干噪比测量，得到解调信干噪比测量值SINRf；
SINR估计模块，用于对各窄带信道的解调信干噪比测量值SINRf，进行独立估计，并进行连续均值滤波，得到当前信道的解调信干噪比估计值SINRf_est；
双门限干扰位置检测模块，将各窄带信道的干扰检测量估计值RSSIf_est与预设的第一门限值进行比较，判断干扰信号是否存在；所述第一门限值为能量检测阈值RSSI_thr；将解调信干噪比估计值SINRf_est与预设的第二门限值进行比较；判断干扰信号是否为有效信号；将存在干扰信号，且干扰信号为有效信号的信道作为干扰位置检测结果输出；所述第二门限值为信干噪比检测阈值SINR_thr。
进一步优选的，还包括同步捕获模块和多路信号抽取模块；
所述同步捕获模块用于对输入的宽带无线信号进行滑动自相关，获得各信道的最佳采样时刻，用于各信道的脉冲信号抽取；
所述多路信号抽取模块用于获得最佳采样时刻后，以信息速率抽取输入的高速数字信号，获得各窄带信道的低速数字信号，进行多路窄带信道的接收信号能量测量和估计。
进一步优选的，还包括统计滤波模块，所述统计滤波模块用于在检测周期内设定多个检测窗口，对多次干扰位置检测结果进行统计滤波，采用相对统计量进行异常值过滤。
本申请公开的基于双门限判决的窄带干扰位置检测方法及系统，相比于现有技术至少包括以下优点：
1、本发明提供的基于双门限判决的窄带干扰位置检测方法及系统，采用信道化接收机捕获多个窄带信道的并行信号，可以同时获得工作频带内多个窄带信道的接收信号能量测量值和解调信干噪比测量值，并检测频域干扰点，相比于扫描式信道检测，为干扰位置检测提供有效的频域指示，并行信道检测的效率更高。
2、本发明提供的基于双门限判决的窄带干扰位置检测方法及系统，相较于单纯的能量检测，本发明所述动态阈值检测算法具备优良的环境适应性，能够检测无信号、小信号和干扰信号条件下的窄带信道干扰，同时给出干扰检测和干扰位置检测结果。
3、本发明提供的基于双门限判决的窄带干扰位置检测方法及系统，采用双门限干扰位置检测方法，将各窄带信道的干扰检测量估计值与设定的能量检测阈值进行比较，判断干扰信号是否存在；同时设定解调信噪比加干扰容限作为信干噪比检测阈值，排除干扰信号为有效信号的可能。在信号能量波动范围较大的环境中，双门限干扰位置检测比单纯的能量检测更为准确。
4、本发明提供的基于双门限判决的窄带干扰位置检测方法及系统，对各窄带信道的接收信号能量进行独立估计，并进行连续均值滤波，得到相对稳定的接收信号能量水平估计。干扰检测量的估计窗口长度可以根据检测的实时性和准确性需求进行自适应调整。能够检测无信号、小信号和大信号条件下的窄带信道干扰，同时给出干扰检测和干扰位置检测结果。
附图说明
图1为本发明的基于双门限判决的窄带干扰位置检测方法的流程图；
图2为本发明提供的基于双门限判决的窄带干扰位置检测系统的结构示意图；
图3为本发明实施例的基于双门限判决的窄带干扰位置检测方法的流程图；
图4为本发明提供的基于双门限判决的窄带干扰位置检测系统的窄带干扰位置检测示意图。
具体实施方式
以下通过附图和具体实施方式对本发明作进一步的详细说明。
如1图所示，本发明一方面实施例提供的本发明的一种基于双门限判决的窄带干扰位置检测方法，包括以下步骤：
S1、捕获多个窄带信道的并行信号，获得工作频带内各窄带信道的接收信号能量测量值RSSIf和解调信干噪比测量值SINRf；
S2、对各窄带信道的接收信号能量测量值RSSIf进行独立估计，并进行连续均值滤波，得到稳定的接收信号能量水平估计值，作为当前信道的干扰检测量估计值RSSIf_est；
对各窄带信道的解调信干噪比测量值SINRf，进行独立估计，并进行连续均值滤波，得到当前信道的解调信干噪比估计值SINRf_est。
S3、根据当前信道的干扰检测量估计值RSSIf_est与预设的第一门限值进行比较，判断干扰信号是否存在；将解调信干噪比估计值SINRf_est与预设的第二门限值进行比较，判断干扰信号是否为有效信号；将存在干扰信号，且干扰信号为有效信号的信道作为干扰位置检测结果输出。
进一步，所述第一门限值为能量检测阈值RSSI_thr；所述第二门限值为信干噪比检测阈值SINR_thr。
在S1中，捕获多个窄带信道的并行信号，获得工作频带内多个窄带信道的接收信号能量测量值RSSIf和解调信干噪比测量值SINRf；
并行信道测量，是基于宽带信道化接收机进行多个窄带信道的并行信号捕获和信道测量，获得工作频带内多个窄带信道的接收信号能量测量值RSSIf和解调信干噪比测量值SINRf。其中，
a. 接收信号能量测量值
窄带信道CHf的接收信号能量测量值RSSIf，是对经过射频放大、AD输出的接收信号Yf(t)进行能量测量。可选地，RSSIf可以采用时域积分或频域变换的方式进行能量测量得出，时域积分是对接收信号的多个样点值进行积分和均值量化；频域变换的方式是将接收信号进行FFT变换，通过频域的窄带信号能量均值表示该信道的接收信号能量。
接收信号Yf(t)的能量检测过程，用公式表示如下：
Yf(t) =s(t) n(t) j(t) （公式1）
（公式2）
（公式3）
式中，s(t)、n(t)、j(t)分别为有用信号，背景噪声和干扰信号；L为采样点数量，i为采样点序号，即任意一个采样点；为任意采样点i处捕获的接收信号；RSSIf，为接收信号能量测量值；在干扰信号j(t)足够大时，可以认为接收信号的总功率近似为干扰信号强度。
b. 解调信干噪比测量值
对于数字化通信系统，解调信干噪比测量值可以基于已知的脉冲信号结构信息，分别进行信号能量和噪声能量的统计，二者的比值即为信干噪比。
可选地，跳频通信系统中的解调信干噪比测量值SINRf的测量方式如下：
(公式4)
式中，SINRf为按频点统计的解调信干噪比测量值，Efduty为脉冲占用期能量均值，Efempty为脉冲空闲期能量均值，K为频谱检测周期内频率脉冲f出现的次数。
S2、对各窄带信道的接收信号能量测量值RSSIf进行独立估计，并进行连续均值滤波，得到相对稳定的接收信号能量水平估计，作为当前信道的干扰检测量估计值RSSIf_est。
干扰检测量估计，是对各窄带信道的接收信号能量测量值RSSIf进行独立估计，并进行连续均值滤波，得到相对稳定的接收信号能量水平估计，作为当前信道的干扰检测量估计值RSSIf_est。
信道干扰检测量的估计过程，用公式表示如下：
（公式5）
式中，RSSIf_est为按频点统计的干扰检测量估计值，Wdetect为检测窗口长度，N_pulse为每单位检测窗口内的脉冲数量。RSSIif为第f个频点的第i个脉冲的接收信号能量测量值。
S3、根据当前信道的干扰检测量估计值RSSIf_est与预设的第一门限值进行比较，判断干扰信号是否存在；将解调信干噪比估计值SINRf_est与预设的第二门限值进行比较，判断干扰信号是否为有效信号；将存在干扰信号，且干扰信号为有效信号的信道作为干扰位置检测结果输出。
进一步，所述第一门限值为能量检测阈值RSSI_thr；所述第二门限为信干噪比检测阈值SINR_thr。
优选的，所述能量检测阈值RSSI_thr为背景噪声固定增益的干扰检测量估计值加虚警容限设定的。
若干扰检测量估计值RSSIf大于能量检测阈值RSSI_thr，且信干噪比SINRf小于门限SINR_thr，则判决该信道有干扰H1；其他情况下，判决该信道无干扰H0。
在实际通信环境中，无线信号和背景噪声的能量水平会存在波动，干扰检测门限值若完全取决于当前信号能量的估计值会产生一定概率的误检问题，即信号和噪声能量值大于干扰检测门限。为了减少虚警概率Pfa，能量检测阈值可以用由虚警概率的容限值给出，如下式：
（公式6）
式中，RSSI_thr为能量检测阈值；Pfa为虚警概率，通常取为0.001；N为显著统计量的统计次数；RSSI_signal为接收信号额定灵敏度；Q函数表达式为。
窄带信道干扰判决的先验假设H1/H0如下式：
（公式7）
式中，CQIf为干扰位置检测结果，表示第f个信道是否存在干扰信号。
进一步优选的，还包括对干扰位置检测结果进行统计滤波输出，所述统计滤波输出包括以下步骤：在检测周期内设定多个检测窗口，对多次干扰位置检测结果进行统计滤波，滤除异常检测值。
为了防止误检，在检测周期内设定多个检测窗口，对M次干扰位置检测结果进行统计滤波，滤除异常检测值，进一步保证干扰位置检测的可靠性。
统计滤波采用相对统计量rel_Stat进行异常值过滤，如下式：
（公式8）
（公式9）
式中，rel_Stat为相对统计量，rel_thr为相对统计量阈值，通常取为0.9~1的数量值；M为检测周期内的干扰检测次数，为第f个干扰检测点第i次干扰检测位置的干扰位置检测结果。
如图2所示，本发明还提供一种基于双门限判决的窄带干扰位置检测系统，用于实施上述实施例中的基于双门限判决的窄带干扰位置检测方法，包括AD变换模块、RSSI测量模块、RSSI估计模块、SINR测量模块、SINR估计模块、双门限干扰位置检测模块。
所述AD变换模块用于接收射频前端输入的宽带无线信号，并输出数字域的宽带无线信号；AD变换接收射频前端输入的宽带无线信号，经过低通滤波和模数变换后，输出数字域的宽带无线信号。
所述RSSI测量模块，用于对数字域的宽带无线信号进行多信道并行能量检测，获得工作频带内多个窄带信道的接收信号能量测量值RSSIf。
RSSI测量模块对数字域的宽带无线信号进行多信道并行能量检测，输出多信道能量测量结果RSSIf， RSSI测量可以采用时域积分和频域变换方法。
RSSI估计模块，用于对各窄带信道的接收信号能量测量值RSSIf进行独立估计，并进行连续均值滤波，得到相对稳定的接收信号能量水平估计，作为当前信道的干扰检测量估计值RSSIf_est。
还包括，同步捕获模块，在已知通信系统的先验信息（扩/跳频同步序列）条件下，同步捕获模块对输入的宽带无线信号进行滑动自相关，获得各信道的最佳采样时刻，用于各信道的脉冲信号抽取。
多路信号抽取模块，在获得最佳采样时刻后，多路信号抽取模块以信息速率抽取输入的高速数字信号，获得各窄带信道的低速数字信号，用于多路信号测量和估计。
SINR测量模块，用于对各信道输入的低速脉冲信号进行信干噪比测量，得到解调信干噪比测量值SINRf。
SINR估计模块，用于对各窄带信道的解调信干噪比测量值SINRf，进行独立估计，并进行连续均值滤波，得到当前信道的解调信干噪比估计值SINRf_est；
双门限干扰位置检测模块，用于对各个信道的干扰检测量估计值RSSIf_est进行统计，将各窄带信道的干扰检测量估计值RSSIf_est与设定的能量检测阈值进行比较，判断干扰信号是否存在；并将解调信干噪比估计值SINRf_est与预设的信干噪比检测阈值SINR_thr进行比较；判断干扰信号是否为有效信号。
进一步优选的，还包括统计滤波模块，所述统计滤波模块用于在检测周期内设定多个检测窗口，对多次干扰位置检测结果进行统计滤波，滤除异常检测值。对M次检测结果进行相对统计量估计，若相对统计量大于门限值才输出确定的检测结果，去除异常检测值；经过上述测量模块和滤波模块的处理，干扰位置检测系统能够输出相对灵敏、可靠的干扰位置检测结果。
如图3为在本发明提供的检测方法及系统的典型实施例的实施流程图，包括以下步骤：开始后先检测统计滤波窗口是否结束，若统计滤波窗口为Y时，则多次检测，统计滤波；检查信道检测窗口是否结束；若信道检测窗口为Y，则对干扰检测量估计值，进行动态阈值估计，实施干扰位置检测，若都未结束，则开始进行以下步骤：包括以下步骤：
1）接收信号能量测量
对数字域的宽带无线信号进行多信道并行能量检测，输出各窄带信道的信号能量测量值RSSIf， 接收信号能量测量可以采用时域积分和频域变换方法，时域积分和频域变换方法如前述方法实施例所述。
2）信道干扰检测量估计
在检测窗口内信道干扰检测量估计是对各信道输入的能量测量进行连续均值滤波，得到相对稳定的接收信号能量水平估计，作为当前信道的干扰检测量估计值RSSIf_est。估计方法如前述方法实施例所述。
3）信号同步捕获
在已知通信系统的先验信息（扩/跳频同步序列）条件下，同步捕获对输入的宽带无线信号进行滑动自相关，获得各信道的最佳采样时刻，用于各信道的脉冲信号抽取。
4）多路信号抽取
在获得最佳采样时刻后，多路信号抽取以信息速率抽取输入的高速数字信号，获得各窄带信道的低速数字信号，用于后续的多路信号测量和估计。
5）解调信干噪比测量
解调信干噪比测量是对各信道输入的低速脉冲信号进行信干噪比测量，基于已知的脉冲信号结构信息，分别进行信号能量和噪声能量的统计，二者的比值即为解调信干噪比测量值SINRf。可选地，解调信干噪比测量方法如前述方法实施例所述。
6）SINR检测量估计
在检测窗口内，SINR估计模块对各信道输入的信干噪比测量进行连续均值滤波，作为解调信干噪比检测量估计值SINRf_est。
7）双门限干扰位置检测
在检测窗口内，将当前信道的干扰检测量估计值RSSIf_est与预设的第一门限值进行比较，判断干扰信号是否存在；将解调信干噪比估计值SINRf_est与预设的第二门限值进行比较，判断干扰信号是否为有效信号。干扰位置检测方法，如前述方法实施例中所述。
8）多次检测、统计滤波
为了最大程度地消除干扰误检和虚警问题，本发明还提供了统计滤波技术，对M次检测结果进行相对统计量估计，若相对统计量大于门限值才输出确定的检测结果，去除异常检测值后结束。统计滤波方法，如前述方法实施例所述。
经过上述测量模块和滤波模块的处理，干扰位置检测系统能够输出相对灵敏、可靠的干扰位置检测结果。
图4为该实施例实施之后的窄带干扰位置检测波形示意图。图中深灰色带3个波峰的信号为检测的干扰信号，虚线下方均匀的浅灰色波形为有效信号，干扰信号与有效信号重叠处，为干扰位置。
显然，上述实施例仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
本文用于企业家、创业者技术爱好者查询，结果仅供参考。
