CN102546116A - 基于压缩感知的跳频信号盲检测方法 - Google Patents

基于压缩感知的跳频信号盲检测方法 Download PDF

Info

Publication number
CN102546116A
CN102546116A CN2012100319747A CN201210031974A CN102546116A CN 102546116 A CN102546116 A CN 102546116A CN 2012100319747 A CN2012100319747 A CN 2012100319747A CN 201210031974 A CN201210031974 A CN 201210031974A CN 102546116 A CN102546116 A CN 102546116A
Authority
CN
China
Prior art keywords
frequency hopping
signal
detection
sampled point
under
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012100319747A
Other languages
English (en)
Other versions
CN102546116B (zh
Inventor
刘乃安
吴俊�
付卫红
沈常林
张妍飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201210031974.7A priority Critical patent/CN102546116B/zh
Publication of CN102546116A publication Critical patent/CN102546116A/zh
Application granted granted Critical
Publication of CN102546116B publication Critical patent/CN102546116B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明公开了一种跳频信号盲检测方法,该方法包括的步骤1)通过压缩采样系统获得压缩采样值y;2)分别求出每个跳频采样点与其在H1情况下数学期望的偏差,然后分别平方之后取和,记为H;3)分别求出每个跳频采样点与其在H0情况下数学期望的偏差,然后分别平方之后取和,记为H;4)若H<H,判为H1,即存在跳频信号;否则,判为H0,不存在跳频信号。这种方法可以根据检测要求灵活选择采样点数完成检测。可以不需要使用所有的采样点,因为每一个采样点都有对应的数学期望,都可以作为检测判决依据,而且不需要采用蒙特卡洛方法统计检测阈值。

Description

基于压缩感知的跳频信号盲检测方法
技术领域
本发明属于通信对抗领域,具体涉及一种基于压缩感知的跳频信号盲检测方法。
背景技术
跳频通信因其良好的抗干扰、低截获特性及良好的组网能力,已在军事领域得到了广泛的应用。开展对跳频通信对抗的研究、寻求截获、分选识别和干扰跳频通信信号的方法,已成为了当前通信对抗领域紧迫而困难的重要研究课题之一。
目前跳频信号的检测方法多种多样,主要有信道辐射机检测、压缩接收机检测、空间相关检测、时频分析检测、基于小波变换检测、声光频谱分析、DFT谱分析及其他谱分析方法。大部分的检测方法都是在假设已知部分参数、信道环境较理想的条件下实施的,对于恶劣电磁环境下的信号盲识别效果较差,其中利用相关检测方法实现对跳频信号的检测,其判决门限的确定和信噪比密切相关,在实际实施中存在不便,而时频分析和小波算法在低信噪比条件下虽然可以达到较好的检测效果,但是这两种算法非常复杂,计算量极大,实时性差并且对硬件要求很高。
然而在通信侦察和对抗等非协作通信应用中,通常需要在不能事先获知任何协议和参数的情况下,快速并准确的检测出整个频域段的跳频信号,它是实现跳频信号分离、干扰、监听、欺骗的基础。因此要求检测跳频信号的方法不仅具有简单快速高实时性的基本要求而且能够实现低信噪比复杂电磁环境下的准确识别,这是跳频信号盲检测面临的一个巨大难题。
发明内容
为了解决背景技术中所述的相关检测方法对信噪比敏感、时频分析和小波算法非常复杂、计算量极大、实时性差并且对硬件要求很高的技术问题,本发明提供了一种基于压缩感知的跳频信号盲检测方法。
本发明的技术解决方案是:本发明提供了一种基于压缩感知的跳频信号盲检测方法,
检测的目的是区别两种假设:
H0:x=n
H1:x=s+n
H0表示不存在跳频信号的情况;H1表示存在跳频信号的情况。
式中:s∈RN表示跳频信号,n是加性高斯白噪声,均值为0,方差为σ2IN,x∈RN是长度为奈奎斯特采样点数的接收信号,常规的跳频信号检测方法直接以x为样本进行判决,为了保证高的检测性能,需要A/D采样器以很高的采样频率对接收到的连续信号进行采样,而采用基于压缩感知的跳频信号检测方法解决了这个难题,接收信号x首先经过压缩采样系统从高维映射成低维信号y,也就是y=Φx,其中Φ∈RM×N是一个满足有限等距性质(restricted isometry property,RIP)的M×N观测矩阵,y是压缩感知获得的采样值。
因此,可在不重构原信号的情况下,通过直接处理压缩采样值y来完成检测。具体地,压缩采样值y为:
H0:y=Φn
H1:y=Φ(s+n)
下面分别分析压缩感知采样值在H1和H0两种情况下的数字特征(NumericalCharacteristics,NC):
H1情况下,高斯白噪声的分布满足n~N(0,σ2IN),跳频信号s∈RN,因此s+n~N(s,σ2IN),在每次压缩感知过程中Φ可以认为是一个确定的、已知的M×N的测量矩阵,那么在H1中的过程可以认为是一个N维高斯变量s+n的线性变换过程。由概率论与数理统计中的相关理论可知,N维高斯变量的线性变换仍服从高斯分布,因此可以获得采样值y的数字特征,即
y=Φ(s+n)~N(Φs,Φσ2ΦT)
同理H0情况下y=Φn~N(0,Φσ2ΦT)
由于y是一个M维的变量,因此,可以认为压缩感知的采样点有M个,由y的数字特征可知,每一个采样点yi的数字特征也是可知的,即:
H0:E(yi)=0,
Figure BSA00000669379800031
H1:E(yi)=φis,
Figure BSA00000669379800032
l≤i≤M
式中:φi为矩阵Φ的第i行向量。
由上面分析可知,在H1和H0情况下,压缩感知获得的采样值y方差相同,数学期望不同。H1情况下的数学期望是φS,H0情况下各个采样值的数学期望均为0,两种情况下采样值的数学期望均与噪声强度无关。因此可以利用采样值在不同假设情况下的数学期望不同这一特点,采用将实际采样值与其在两种假设情况下数学期望的偏差作为判决依据的方法,来完成跳频信号检测的目的。
具体步骤如下:
1)接收信号x首先经过压缩采样系统从高维映射成低维信号y,也就是y=Φx,其中Φ∈RM×N是一个M×N的观测矩阵,y是压缩感知获得的M个跳频采样点。
2)在由压缩感知获得的M个跳频采样点中任意选择l(l≤M)个跳频采样点,分别求出每个跳频采样点与其在H1情况下数学期望的偏差,然后分别平方之后取和:
Figure BSA00000669379800033
式中y1,...,yl为从M个跳频采样点中任意选择的l个跳频采样点。
3)使用相同的l(l≤M)个跳频采样点,分别求出每个跳频采样点与其在H0情况下数学期望的偏差,然后分别平方之后取和:
H 0 Δ = Σ i = 1 l ( y i - E ( y i ) | H 0 ) 2 = Σ i = 1 l ( y i - 0 ) 2 ;
4)若H<H,判为H1,即存在跳频信号;否则,判为H0,不存在跳频信号。
本发明与现有方法相比具有如下优点:
1)本发明克服了低信噪比下检测效果不好的问题,因为噪声强度的变化不会改变采样值的数学期望,它只会改变采样值的方差,而采样值在H1和H0情况下的方差又是相同的,而且不需要统计检测门限。
2)降低了数据量及计算复杂度,需要处理的跳频信号由原来长度为N的采样信号变为长度为M的观测信号,而且可以根据检测要求灵活选择采样点数完成检测。
3)与传统的先重构再检测的跳频信号盲检测算法对比,也节省了重构算法本身所需要的计算时间,因此本发明算法具有更快的检测速度,在要求快速且准确的跳频信号检测中具有重要的实际意义。
附图说明
图1是本发明方法的流程框图;
图2是压缩不同情况下跳频信号检测成功率对比分析示意图;
图3不同跳频信号检测方法检测时间对比分析示意图;
具体实施方式
以下结合具体实施例,对本发明进行详细说明。
参照图1,本发明的具体实施方式如下:
1)接收信号x首先经过压缩采样系统从高维映射成低维信号y,也就是y=Φx,其中Φ∈RM×N是一个M×N的观测矩阵,y是压缩感知获得的M个跳频采样点。
2)在由压缩感知获得的M个跳频采样点中任意选择l(l≤M)个跳频采样点,分别求出每个跳频采样点与其在H1情况下数学期望的偏差,然后分别平方之后取和: H 1 Δ = Σ i = 1 l ( y i - E ( y i ) | H 1 ) 2 = Σ i = 1 l ( y i - φ i s ) 2 ; 式中y1,...,yl为从M个跳频采样点中任意选择的l个跳频采样点。
3)使用相同的l(l≤M)个跳频采样点,分别求出每个跳频采样点与其在H0情况下数学期望的偏差,然后分别平方之后取和:
H 0 Δ = Σ i = 1 l ( y i - E ( y i ) | H 0 ) 2 = Σ i = 1 l ( y i - 0 ) 2 ;
4)若H<H,判为H1,即存在跳频信号;否则,判为H0,不存在跳频信号。
由上述流程可以看到,检测决定是通过比较真实采样值与其在两种假设情况下数学期望的偏差程度完成的。由于压缩感知采样值在H1和H0情况下的方差相同,因此,当真实采样值与压缩感知采样值在H1情况下的数学期望的偏差较小时,就表示该采样值与H1情况下的数学期望较接近,由高斯分布的性质可知,就表明真实采样值来自H1情况的概率较高;反之,表明真实采样值来自H0情况的概率较高。
实验中模拟了一个跳速为1000H/s,跳频带宽为240MHz,采样频率为480MHZ,每跳采样N=4.8×105个点,共10个跳变频率点,所以总共采样点为4.8×106个点。考察在不同的压缩比下,NC检测算法针对跳频信号在各个信噪比下的检测效果。令M=N/8、N/16、N/32,信噪比变化范围[-30,10],步进2。NC算法使用全部采样点进行检测,即L=M.其仿真实验结果参照如图2所示。由图2可以看出,本发明检测成功率随着M的减小而下降,尤其在低信噪比时;二是当信噪比高于-6dB时,本发明基本上可以保证跳频信号检测获得成功;三是相比较传统的检测方法,本发明可以保证较高的检测成功率。
参照图3描述了压缩率M/N与检测时间的关系曲线。选择了基于OMP和ROMP的检测算法和本算法进行对比。信号长度N=320,在信噪比SNR=-6dB时,每一个观测值M进行10次检测。从图3中可以看出本检测算法的检测时间远低于基于OMP和ROMP的检测算法。由此可见,与传统的重构检测法相比较,本发明算法具有更快的检测速度,在要求快速且准确的跳频信号检测中具有重要的实际意义。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (1)

1.一种基于压缩感知的跳频信号盲检测方法,其特征在于:该方法包括以下步骤:
1)接收信号x首先经过压缩采样系统从高维映射成低维信号y,也就是y=Φx,其
中Φ∈RM×N是一个M×N的观测矩阵,y是压缩感知获得的M个采样值。
2)在由压缩感知获得的M个跳频采样点中任意选择l(l≤M)个跳频采样点,分别
求出每个跳频采样点与其在H1情况下数学期望的偏差,然后分别平方之后取和:
Figure FSA00000669379700011
式中y1,...,yl为从M个跳频采样点中任意选择的l个跳频采样点。
3)使用相同的l(l≤M)个跳频采样点,分别求出每个跳频采样点与其在H0情况下数学期望的偏差,然后分别平方之后取和:
H 0 Δ = Σ i = 1 l ( y i - E ( y i ) | H 0 ) 2 = Σ i = 1 l ( y i - 0 ) 2 ;
4)若H<H,判为H1,即存在跳频信号;否则,判为H0,不存在跳频信号。
CN201210031974.7A 2012-02-14 2012-02-14 基于压缩感知的跳频信号盲检测方法 Expired - Fee Related CN102546116B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210031974.7A CN102546116B (zh) 2012-02-14 2012-02-14 基于压缩感知的跳频信号盲检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210031974.7A CN102546116B (zh) 2012-02-14 2012-02-14 基于压缩感知的跳频信号盲检测方法

Publications (2)

Publication Number Publication Date
CN102546116A true CN102546116A (zh) 2012-07-04
CN102546116B CN102546116B (zh) 2015-07-29

Family

ID=46352143

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210031974.7A Expired - Fee Related CN102546116B (zh) 2012-02-14 2012-02-14 基于压缩感知的跳频信号盲检测方法

Country Status (1)

Country Link
CN (1) CN102546116B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103051367A (zh) * 2012-11-27 2013-04-17 西安电子科技大学 一种基于聚类的同步正交跳频信号盲源分离方法
CN103532645A (zh) * 2013-10-10 2014-01-22 南京邮电大学 一种观测矩阵优化的压缩频谱感知方法
CN106209703A (zh) * 2016-07-08 2016-12-07 中国人民解放军信息工程大学 一种跳频信号参数盲估计方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101512998A (zh) * 2006-09-11 2009-08-19 艾利森电话股份有限公司 时间-频率跳频图案的检测
CN102253117A (zh) * 2011-03-31 2011-11-23 浙江大学 一种基于压缩感知的新型声学信号采集方法
CN102316061A (zh) * 2010-07-07 2012-01-11 中国科学院微电子研究所 一种跳频正交频分复用系统的时间同步方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101512998A (zh) * 2006-09-11 2009-08-19 艾利森电话股份有限公司 时间-频率跳频图案的检测
CN102316061A (zh) * 2010-07-07 2012-01-11 中国科学院微电子研究所 一种跳频正交频分复用系统的时间同步方法及装置
CN102253117A (zh) * 2011-03-31 2011-11-23 浙江大学 一种基于压缩感知的新型声学信号采集方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王健雄等: "基于循环平稳特性的跳频信号盲检测算法*", 《舰船电子工程》 *
范海宁: "基于循环自相关的跳频信号盲检测算法", 《微计算机信息》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103051367A (zh) * 2012-11-27 2013-04-17 西安电子科技大学 一种基于聚类的同步正交跳频信号盲源分离方法
CN103051367B (zh) * 2012-11-27 2015-08-26 西安电子科技大学 一种基于聚类的同步正交跳频信号盲源分离方法
CN103532645A (zh) * 2013-10-10 2014-01-22 南京邮电大学 一种观测矩阵优化的压缩频谱感知方法
CN103532645B (zh) * 2013-10-10 2015-09-30 南京邮电大学 一种观测矩阵优化的压缩频谱感知方法
CN106209703A (zh) * 2016-07-08 2016-12-07 中国人民解放军信息工程大学 一种跳频信号参数盲估计方法及装置
CN106209703B (zh) * 2016-07-08 2019-06-18 中国人民解放军信息工程大学 一种跳频信号参数盲估计方法及装置

Also Published As

Publication number Publication date
CN102546116B (zh) 2015-07-29

Similar Documents

Publication Publication Date Title
CN106443203B (zh) 一种脉冲信号检测系统及方法
CN105678273B (zh) 射频指纹识别技术瞬态信号的起始点检测算法
CN103399297B (zh) 一种基于机器学习的超宽带非视距鉴别方法
CN108600135A (zh) 一种信号调制方式的识别方法
CN104717106B (zh) 一种基于多变量序贯分析的分布式网络流量异常检测方法
CN106772268B (zh) 一种高斯白噪声下的弱信号盲检测方法
CN106130942A (zh) 一种基于循环谱的无线通信信号调制识别及参数估计方法
CN102510363A (zh) 一种强干扰源环境下的lfm信号检测方法
CN110099019A (zh) 基于深度学习的LoRa调制信号检测方法
CN105785324A (zh) 基于mgcstft的线性调频信号参数估计方法
CN105429719A (zh) 基于功率谱和多尺度小波变换分析强干扰信号检测方法
CN103973383B (zh) 基于Cholesky分解与特征值的协作频谱检测方法
CN106170139B (zh) 一种频谱检测方法及系统
CN103973382B (zh) 基于有限随机矩阵的频谱检测方法
CN109342828A (zh) 一种基于频域恒虚警的雷电脉冲信号检测方法
CN110398722A (zh) 基于随机矩阵有限谱的扩展目标回波检测方法
CN112462355A (zh) 一种基于时频三特征提取的对海目标智能检测方法
CN110764152B (zh) 一种无人机快速检测与识别的装置及方法
CN102546116A (zh) 基于压缩感知的跳频信号盲检测方法
CN102594836A (zh) 一种基于小波能量谱的流量识别方法
CN103561469A (zh) 一种室内定位方法及系统
Martone et al. Fast technique for wideband spectrum sensing
CN108718223B (zh) 一种非合作信号的盲频谱感知方法
CN105680964A (zh) 一种频谱感知方法和频谱感知系统、客户端和服务端
CN102111228A (zh) 一种基于循环对称性的认知无线电频谱感知方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150729

Termination date: 20160214