CN106872778A - 基于压缩感知的微弱信号检测方法 - Google Patents
基于压缩感知的微弱信号检测方法 Download PDFInfo
- Publication number
- CN106872778A CN106872778A CN201710111725.1A CN201710111725A CN106872778A CN 106872778 A CN106872778 A CN 106872778A CN 201710111725 A CN201710111725 A CN 201710111725A CN 106872778 A CN106872778 A CN 106872778A
- Authority
- CN
- China
- Prior art keywords
- signal
- noise
- detected
- compressed sensing
- signals based
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于信号处理领域,尤其涉及基于压缩感知的微弱信号检测方法。本发明的目的是提过一种基于压缩感知的微弱信号检测方法,在低信噪比下对压缩后的接收信号检测待检测信号是否存在。因为待检测信号具有稀疏性,在进行稀疏表示时,x=ΨP,Ψ是参考信号x投影的字典基,P=[0,0,…,1,0,…,1,…]是稀疏向量,x在字典基Ψ下投影向量非零位置是固定的。而高斯白噪声在频谱均匀分布,在字典基中出现非零位置是随机出现并且服从均匀分布,利用噪声投影位置分布的独立性可设计方案,实现检测目标信号的目的。
Description
技术领域
本发明属于信号处理领域,尤其涉及基于压缩感知的微弱信号检测方法。
背景技术
近年来,压缩感知理论在雷达信号处理、无线通信中得到了日趋成熟的应用,比如对信号进行检测,识别和参数估计等。与传统的Nyquist采样定理不同,引入压缩感知之后,在明显减少了对采样数据量的分析的情况下,基本不丢失压缩之前信号的特征。目前雷达系统中采用的压缩感知技术,通常需要利用重构算法对信号进行重构,进而实现目标检测与参数估计。但是,重构算法一般要求较高的信噪比,才能得到满足目标检测与参数估计的重构信号。对于像外辐射源雷达由于其本质是双基地非合作雷达,目标回波功率低,通常不满足信号重构的信噪比要求。因此采用压缩感知技术实现微弱信号检测需要新方法。本发明主要利用了压缩感知理论对信号稀疏表示时稀疏向量位置的统计特性,实现了在低信噪比下对待检测信号的有效检测。
信号检测是通过对接收信号分析判断目标信号是否存在,采用压缩感知技术后,目标检测的数学模型为其中,Φ是传感矩阵,n是高斯白噪声,x是待检测信号。
假设H0是待检测信号不存在,假设H1是待检测信号存在的情况。目前存在的检测算法主要有下面几种方法:(1)2010年刘冰等人提出均值比较算法,其核心思想是噪声是均值为0的白噪声,那么H0情况下E(y)=E(Φn)=0,H1情况下E(y)=E(Φ(x+n))=Φx,而两种情况下方差不改变。所以他采用的方法是将实际采样值与其存在两种假设情况下数学期望的偏差作为判决依据,完成检测;(2)2014年王康等人提出稀疏系数幅度比较算法,其核心思想是对用两个观测测量矩阵同时接收待检测信号,利用OMP算法得到的幅度和位置信息联合设计判决门限,文中提出门限系数的设置是根据经验取值,没有给出具体确定计算方法;(3)2015年Alireza Hariri等人提出最大似然比值算法,其核心思想是没有对待检测信号重构,得到两种情况的似然函数,然后做比值,通过似然比函数概率密度特性得到门限公式,进行判决。但是,上述三种检测算法在信噪比为-5dB以下就不能实现目标检测。这就限制了检测器的应用范围,为了探测更大的范围,就需要解决在更低信噪比下对待检测信号的检测。因此,进一步研究低信噪比下对待检测信号的检测是十分有必要的。
发明内容
本发明的目的是提过一种基于压缩感知的微弱信号检测方法,在低信噪比下对压缩后的接收信号检测待检测信号是否存在。因为待检测信号具有稀疏性,在进行稀疏表示时,x=ΨP,Ψ是参考信号x投影的字典基,P=[0,0,…,1,0,…,1,…]是稀疏向量,x在字典基Ψ下投影向量非零位置是固定的。而高斯白噪声在频谱均匀分布,在字典基中出现非零位置是随机出现并且服从均匀分布,利用噪声投影位置分布的独立性可设计方案,实现检测目标信号的目的。
本发明方法为:已知参考信号x,通过OMP算法得到的稀疏向量P=[0,0,…,1,0,…,1,…],所述稀疏向量P同时也是待检测信号x的稀疏向量,记P中非0向量位置为P1,当只有噪声存在时利用OMP算法得到稀疏向量P,记P中非0位置为P2、P1、P2相互独立,那么P1、P2中存在相同的列数的期望是组合数计算出来的,再由中心极限定理可知独立同分布的统计量服从正态分布且可以通过统计平均代替算术平均,当重复次数越多时其正态分布的均值计算越稳定,从而对采样的数据进行多次蒙特卡洛计算,找出P1、P2相关位置个数的均值,得到判决门限,最终实现检测。结果显示具有很好的检测性能,具体检测方法可分为两部分:(一)判决门限设定;(二)对待检测信号检测。
一种基于压缩感知的微弱信号检测方法,具体步骤为:
S1、确定检测判决门限:
S11、通过测量矩阵Φ对数据压缩处理得到y=Φx,其中,x为已知参考信号;
S12、利用OMP算法计算得到字典基投影稀疏向量P;
S13、得到P中非0位置向量P1,以P1作为参考向量,利用压缩感知内部接收机存在的高斯白噪声,得出判决门限Sμ=1.5*Sβ:
S2、对待检测信号进行检测:通过采样得到待检测信号x和噪声n,通过测量矩阵Φ对数据压缩处理得到y=Φ(x+n),检测步骤具体为:
S21、初始化,令j=1,i=0,采样得到待检测信号x和噪声n,设置重复计算次数;
S22、生成随机测量矩阵Φj对数据压缩处理得到yj=Φj(x+n),利用OMP算法得到输入待检测信号时在字典基投影中的稀疏向量P中非0位置向量Pj;
S23、进行判断,如果P1(t)==Pj(s),i=i+1,直至P1、Pj遍历完,s=1:2*K,t=1:2*K,其中,K为待检测信号的稀疏度;
S24、进行判断,如果j≤N,j=j+1,则返回S22,否则
S25、如果SX>Sμ则判为有信号,否则判为信号不存在。
进一步地,S13所述得出判决门限的具体步骤为:
S131、初始化,令j=1,i=0,采样得到噪声信号n,设置重复计算次数N,其中,N为经验值;
S132、生成随机测量矩阵Φj对数据压缩处理得到yj=Φjn,利用OMP算法得到只存在白噪声时在字典基投影中的稀疏向量P中非0位置向量Pj;
S133、判断如果P1(t)==Pj(s),i=i+1,直至P1、Pj遍历完,s=1:2*K,t=1:2*K,其中,K为待检测信号的稀疏度;
S134、进行判断,如果j≤N,j=j+1,则返回S132,否则取门限值Sμ=1.5*Sβ。
进一步地,S131和S21所述N=1500。
进一步地,S2所述待检测信号x和噪声n的信噪比<-20dB。
本发明的有益效果是:
首先分析白噪声频谱服从均匀分布和稀疏待检测信号在字典基最大投影原子位置固定的特性,为检测微弱信号提供了思路,然后,利用OMP算法分别得到参考信号和噪声信号的位置信息得到判决门限,最后,输入待检测信号得到待检测信号的位置向量与参考位置向量相关得到结果与门限比较判决,实现对微弱信号进行比较检测,对单频正弦信号做待检测信号,可以在信噪比dB时可以成功检测,在待检测探测方向具有很好的参考和实际应用。
附图说明
图1本发明方法门限确定流程图
图2本发明方法对待检测信号检测流程图
图3估计成功概率与信噪比关系。
具体实施方式
下面结合附图和具体实施例进一步阐述本发明。
实施例、
本发明在接收待检测信号信噪比SNR=-30dB到-20dB步进1dB时对信号的存在与否进行检测:
如图1、图2所示:
(一)确定检测判决门限:
1.生成参考位置向量
已知待检测信号是单频正弦信号x=cos(2*pi*f4*Ts*ts),f4=400Hz,n为高斯白噪声,数据长度1000点,信噪比为20dB,通过傅里叶变换进行频谱分析可知,单频信号的稀疏度是2,采样频率fs=1000,ts=1/fs,L=1000,Ts=0:L-1。测量矩阵Φ∈300×1000的高斯随机矩阵。字典基其中Ψ是1000×1000阶的单位阵进行傅里叶变换后的矩阵,经过测量之后得到压缩数据y=Φ(x+n),共300个数据。这里,求解参考位置向量P1的具体方法是:根据2007年IEEE TRANSACTIONS ON INFORMATION THEORY,VOL.53,NO.12,DECEMBER 2007发表的Signal Recovery From Random Measurements Via OrthogonalPursuit文章中给出了OMP算法重构压缩信号的步骤,本发明只用到部分步骤提取位置向量信息,下面门限确定方法和待检测信号检测都用到同样步骤,具体如下:
步骤一:初始化参数,r0=y,t=1,P=[],Aug_t=[];
步骤二:寻找最匹配原子的索引λt为最大内积对应的原子位置;
步骤三:更新位置信息P=[P,λt];
步骤四:更新索引集
步骤五:去除选中的列,把该列置零
步骤六:估计系数α=(Aug_tT*Aug_t)-1*Aug_tT*y;
步骤七:计算残差rt=y-Aug_t*α;
步骤八:判断t≤4,若是则t=t+1,返回步骤二;否则输出P1=P。
这样通过计算之后得到参考的位置信息P1。
2.只存在高斯白噪声得出判决门限
只生成1000个高斯白噪声点记作n,测量矩阵Φ∈300×1000的高斯随机矩阵。其中Ψ是1000×1000阶的单位阵进行傅里叶变换后的矩阵,经过测量之后得到压缩数据y=Φn,共300个数据。
步骤一:只输入噪声n,设置循环次数N=1250,i=0,j=1
步骤二:生成随机测量矩阵Φj,yj=Φjn,初始化参数r0=yj,t=1,P=[],Aug_t=[];
步骤三:寻找最匹配原子的索引λt为最大内积对应的原子位置;
步骤四:更新位置信息P=[P,λt];
步骤五:更新索引集
步骤六:去除选中的列,简单起见,这里把该列置零
步骤七:估计系数α=(Aug_tT*Aug_t)-1*Aug_tT*yj;
步骤八:计算残差rt=yj-Aug_t*α;
步骤九:判断t≤4,若是则t=t+1,返回步骤五;否则输出Pj=P;
步骤十:如果P1(s)==Pj(t),i=i+1,其中s=1:4,t=1:4,直至P1、Pj遍历完;
步骤十一:判断如果j≤N,j=j+1返回步骤二,否则取门限值Sμ=1.5*Sβ。
仿真得到Sβ=0.37与理论值几乎相等,所以取门限值Sμ=0.56。
(二)对待检测信号进行检测:
待检测信号是单频正弦信号x=cos(2*pi*f4*(Ts-td)*ts),f4=400Hz,td为传播时延,本实验td=100,延时100个数据点,采样频率fs=1000,ts=1/fs,L=1000,Ts=0:L-1。n为高斯白噪声数据长度1000点,信噪比SNR=-30dB到-20dB步进1dB,每个信噪比下重复100次实验,测量矩阵Φ∈300×1000的高斯随机矩阵。其中Ψ是1000×1000阶的单位阵进行傅里叶变换后的矩阵,经过测量之后得到压缩数据y=Φ(x+n),共300个数据。
步骤一:
步骤二:forj3=1:100;
步骤三:在不同信噪比下输入信号x和噪声n,设置循环次数N=1250,i=0,j=1;
步骤四:生成随机测量矩阵Φj,yj=Φj(x+n),初始化参数r0=yj,t=1,P=[],Aug_t=[];
步骤五:寻找最匹配原子的索引λt为最大内积对应的原子位置;
步骤六:更新位置信息P=[P,λt];
步骤七:更新索引集
步骤八:去除选中的列,简单起见,这里把该列置零
步骤九:估计系数α=(Aug_tT*Aug_t)-1*Aug_tT*yj;
步骤十:计算残差rt=yj-Aug_t*α;
步骤十一:判断t≤4,若是则t=t+1,返回步骤五,否则输出Pj=P;
步骤十二:如果P1(s)==Pj(t),i=i+1,其中s=1:4,t=1:4,直至P1、Pj遍历完;
步骤十三:判断如果j≤1250,j=j+1,返回步骤四,;否则
步骤十四:判断if(SX>0.56),ii=ii+1,直至100次检测结束;
步骤十五:直至信噪比从-30:-20遍历结束;
图3是对实施例的仿真,结果可以看出信噪比在-26dB时检测概率为95%,很精确的判决出待检测信号的待检测是否存在,本发明具有很好的检测性能。
Claims (4)
1.一种基于压缩感知的微弱信号检测方法,其特征在于,具体步骤为:
S1、确定检测判决门限:
S11、通过测量矩阵Φ对数据压缩处理得到y=Φx,其中,x为已知参考信号;
S12、利用OMP算法计算得到字典基投影稀疏向量P;
S13、得到P中非0位置向量P1,以P1作为参考向量,利用压缩感知内部接收机存在的高斯白噪声,得出判决门限Sμ=1.5*Sβ:
S2、对待检测信号进行检测:通过采样得到待检测信号x和噪声n,通过测量矩阵Φ对数据压缩处理得到y=Φ(x+n),检测步骤具体为:
S21、初始化,令j=1,i=0,采样得到待检测信号x和噪声n,设置重复计算次数;
S22、生成随机测量矩阵Φj对数据压缩处理得到yj=Φj(x+n),利用OMP算法得到输入待检测信号时在字典基投影中的稀疏向量P中非0位置向量Pj;
S23、进行判断,如果P1(t)==Pj(s),i=i+1,直至P1、Pj遍历完,s=1:2*K,t=1:2*K,其中,K为待检测信号的稀疏度;
S24、进行判断,如果j≤N,j=j+1,则返回S22,否则
S25、如果SX>Sμ则判为有信号,否则判为信号不存在。
2.根据权利要求1所述一种基于压缩感知的微弱信号检测方法,其特征在于,S13所述得出判决门限的具体步骤为:
S131、初始化,令j=1,i=0,采样得到噪声信号n,设置重复计算次数N,其中,N为经验值;
S132、生成随机测量矩阵Φj对数据压缩处理得到yj=Φjn,利用OMP算法得到只存在白噪声时在字典基投影中的稀疏向量P中非0位置向量Pj;
S133、判断如果P1(t)==Pj(s),i=i+1,直至P1、Pj遍历完,s=1:2*K,t=1:2*K,其中,K为待检测信号的稀疏度;
S134、进行判断,如果j≤N,j=j+1,则返回S132,否则取门限值Sμ=1.5*Sβ。
3.根据权利要求1所述一种基于压缩感知的微弱信号检测方法,其特征在于,S131和S21所述N=1500。
4.根据权利要求1所述一种基于压缩感知的微弱信号检测方法,其特征在于,S2所述待检测信号x和噪声n的信噪比<-20dB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111725.1A CN106872778B (zh) | 2017-02-28 | 2017-02-28 | 基于压缩感知的微弱信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111725.1A CN106872778B (zh) | 2017-02-28 | 2017-02-28 | 基于压缩感知的微弱信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106872778A true CN106872778A (zh) | 2017-06-20 |
CN106872778B CN106872778B (zh) | 2019-09-13 |
Family
ID=59168003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710111725.1A Active CN106872778B (zh) | 2017-02-28 | 2017-02-28 | 基于压缩感知的微弱信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106872778B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107769864A (zh) * | 2017-09-25 | 2018-03-06 | 成都正扬博创电子技术有限公司 | 一种基于强信号和噪声背景下弱信号提取的方法 |
CN107884752A (zh) * | 2017-11-08 | 2018-04-06 | 电子科技大学 | 一种基于压缩感知的外辐射源雷达对目标检测方法 |
CN107976663A (zh) * | 2018-01-24 | 2018-05-01 | 电子科技大学 | 一种基于子空间投影的外辐射源雷达对目标压缩检测方法 |
CN108037422A (zh) * | 2017-12-07 | 2018-05-15 | 国网冀北电力有限公司秦皇岛供电公司 | 一种光伏系统直流电弧故障特征提取和识别方法及装置 |
CN108872945A (zh) * | 2018-07-16 | 2018-11-23 | 电子科技大学 | 一种基于压缩感知的外辐射源雷达对目标检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103117819A (zh) * | 2013-01-18 | 2013-05-22 | 宁波大学 | 一种基于压缩感知的脉冲超宽带信号检测方法 |
CN103983850A (zh) * | 2014-05-13 | 2014-08-13 | 天津大学 | 基于压缩感知的电力系统谐波压缩信号重构与检测方法 |
CN104122540A (zh) * | 2014-07-31 | 2014-10-29 | 电子科技大学 | 基于匹配字典和压缩感知的雷达一维距离像目标识别方法 |
US9398587B1 (en) * | 2011-04-20 | 2016-07-19 | Hrl Laboratories, Llc | Real-time signal detection over very wide band in low SNR environments |
-
2017
- 2017-02-28 CN CN201710111725.1A patent/CN106872778B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9398587B1 (en) * | 2011-04-20 | 2016-07-19 | Hrl Laboratories, Llc | Real-time signal detection over very wide band in low SNR environments |
CN103117819A (zh) * | 2013-01-18 | 2013-05-22 | 宁波大学 | 一种基于压缩感知的脉冲超宽带信号检测方法 |
CN103983850A (zh) * | 2014-05-13 | 2014-08-13 | 天津大学 | 基于压缩感知的电力系统谐波压缩信号重构与检测方法 |
CN104122540A (zh) * | 2014-07-31 | 2014-10-29 | 电子科技大学 | 基于匹配字典和压缩感知的雷达一维距离像目标识别方法 |
Non-Patent Citations (2)
Title |
---|
CHING-HSIEN CHEN等: "Amplitude-Aided 1-Bit Compressive Sensing Over Noisy Wireless Sensor Networks", 《IEEE WIRELESS COMMUNICATIONS LETTERS》 * |
文方青等: "面向低信噪比的自适应压缩感知方法", 《物理学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107769864A (zh) * | 2017-09-25 | 2018-03-06 | 成都正扬博创电子技术有限公司 | 一种基于强信号和噪声背景下弱信号提取的方法 |
CN107884752A (zh) * | 2017-11-08 | 2018-04-06 | 电子科技大学 | 一种基于压缩感知的外辐射源雷达对目标检测方法 |
CN108037422A (zh) * | 2017-12-07 | 2018-05-15 | 国网冀北电力有限公司秦皇岛供电公司 | 一种光伏系统直流电弧故障特征提取和识别方法及装置 |
CN107976663A (zh) * | 2018-01-24 | 2018-05-01 | 电子科技大学 | 一种基于子空间投影的外辐射源雷达对目标压缩检测方法 |
CN108872945A (zh) * | 2018-07-16 | 2018-11-23 | 电子科技大学 | 一种基于压缩感知的外辐射源雷达对目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106872778B (zh) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106872778A (zh) | 基于压缩感知的微弱信号检测方法 | |
CN109597043B (zh) | 基于量子粒子群卷积神经网络的雷达信号识别方法 | |
CN106646344B (zh) | 一种利用互质阵的波达方向估计方法 | |
CN101764652B (zh) | 基于正交匹配追踪的具有压缩感知过程的信号检测方法 | |
CN110632572A (zh) | 基于无意调相特性的雷达辐射源个体识别方法及装置 | |
CN110673109B (zh) | 一种星载大光斑激光雷达全波形数据分解方法 | |
CN101881826A (zh) | 扫描模式海杂波局部多重分形目标检测器 | |
Chen et al. | Impulse response estimation with binary measurements: A regularized FIR model approach | |
CN110231617A (zh) | 目标障碍物方位探测方法、装置、车载雷达及存储介质 | |
CN110207689A (zh) | 一种基于小波熵的脉冲星信号去噪及辨识方法 | |
CN110441752A (zh) | 海上目标检测方法及装置 | |
CN104320144A (zh) | 稀疏度自适应信号重构方法 | |
CN107884752A (zh) | 一种基于压缩感知的外辐射源雷达对目标检测方法 | |
CN108680832B (zh) | 一种基于sa-apso及降维投影动态收缩圆的局部放电定位方法 | |
Bharti et al. | A general method for calibrating stochastic radio channel models with kernels | |
Mallik et al. | Asymptotics for p-value based threshold estimation in regression settings | |
CN109917330B (zh) | 一种存在相位误差时基于稀疏正交匹配追踪理论的到达角估计方法 | |
CN109117698B (zh) | 一种基于最小均方误差准则的噪声背景估计方法 | |
Mroué et al. | Automatic radar target recognition of objects falling on railway tracks | |
CN106019250A (zh) | 基于角闪烁转发式假目标鉴别方法 | |
CN102075468A (zh) | 一种基于磨光函数与Parzen窗估计的ICA盲信号分离方法及其系统 | |
Kang et al. | A class of transformed hazards models for recurrent gap times | |
Huang et al. | Novel modulation recognizer for frequency-hopping signals based on persistence diagram | |
Qu et al. | WGAN-GP-based synthetic radar spectrogram augmentation in human activity recognition | |
CN105548986A (zh) | 海杂波背景下基于预白化比率中值检测器的目标检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |