CN101764652A - 基于正交匹配追踪的具有压缩感知过程的信号检测方法 - Google Patents
基于正交匹配追踪的具有压缩感知过程的信号检测方法 Download PDFInfo
- Publication number
- CN101764652A CN101764652A CN201010032485A CN201010032485A CN101764652A CN 101764652 A CN101764652 A CN 101764652A CN 201010032485 A CN201010032485 A CN 201010032485A CN 201010032485 A CN201010032485 A CN 201010032485A CN 101764652 A CN101764652 A CN 101764652A
- Authority
- CN
- China
- Prior art keywords
- signal
- matrix
- interest
- surplus
- column vector
- 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
Images
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
基于正交匹配追踪的具有压缩感知过程的信号检测方法,它涉及信号处理及信号检测领域。它解决了利用现有的压缩感知实现信号检测时需要重构信号所带来的资源浪费的问题,本发明包括步骤一:设定预设迭代次数T;并设定余量rt的初始值r0和空矩阵V0,使r0=y,y为采样值,再计算传感矩阵V;步骤二:在传感矩阵V中选出相关最大列向量vnt;步骤三:更新矩阵Vt-1为Vt;步骤四:获得稀疏系数估计值;步骤五:更新余量rt-1为rt;步骤六:使t=t+1,当t<T时,返回执行步骤二,否则执行步骤七;步骤七:求解稀疏系数估计值的无穷范数,并比较无穷范数与γ进而选择H0假设或选择H1假设,完成对感兴趣信号s的检测。本发明为通信技术的发展奠定了基础。
Description
技术领域
本发明涉及信号处理及信号检测领域,具体涉及基于正交匹配追踪的具有压缩感知过程的信号检测方法。
背景技术
压缩感知(Compressive sensing),也称为“压缩传感”(以下均称为“压缩感知”),是一种非传统的采样方式,每一步观测是通过信号在观测向量上的投影获得的。该理论指出,如果信号是稀疏的或者在某个基下可压缩,那么用少量的采样值就可以保持信号的结构和相关信息。基于该理论,用于精确重构信号的采样需求数量可以远低于观测的维度,这极大缓解了宽带信号处理的压力。
目前,压缩感知的研究多数集中在信号或图像的重构及相关问题。但是,很多信号处理问题并不需要精确重构信号,例如信号检测这样的任务,信号获取的目的并不是为了重构信号,而是为了从采样数据中提取检测目标的信息,完成一个检测决定。因此利用现有的压缩感知重构信号后再进行检测,带来了大量的资源浪费。
发明内容
为了解决利用现有的压缩感知实现信号检测时需要重构信号所带来的资源浪费的问题,本发明提供基于正交匹配追踪的具有压缩感知过程的信号检测方法。
本发明的基于正交匹配追踪的具有压缩感知过程的信号检测方法,它的检测过程为:
步骤一:设定预设迭代次数T,并设定余量rt的初始值r0和空矩阵V0,使r0=y,再计算传感矩阵V=ΦΨ,其中,y为采样值,Φ表示压缩感知采样过程,Ψ表示感兴趣信号s的变换域;
步骤二:根据公式nt=arg在传感矩阵V中选出与余量rt-1相关性最大的相关最大列向量vnt;
步骤三:根据公式Vt=[Vt-1Vnt]更新矩阵Vt-1为Vt,其中,Vt表示第t次迭代后的传感矩阵;
步骤六:使t=t+1,判断t<T,如果是,返回执行步骤二,否则执行步骤七;
步骤七:求解稀疏系数估计值的无穷范数当时,选择H1假设,H1:y=Φ(s+n)=Φs+Φn=ΦΨθ+Φn=Vθ+Φn,所述H1假设说明感兴趣信号s存在;当时,选择H0假设,H0:y=Φn,其中,γ>0,所述γ的具体数值依据感兴趣信号s相对于加性高斯白噪声n的强度设定,H0假设说明感信号s不存在,完成对感兴趣信号s的检测。
说明书附图
图1是本发明的基于正交匹配追踪的具有压缩感知过程的信号检测方法的流程图;图2是具体实施方式中感兴趣信号s的时域波形示意图;图3是应用具体实施方式对感兴趣信号s进行检测的仿真结果示意图。
具体实施方式
具体实施方式一:根据说明书附图1具体说明本实施方式,本实施方式所述的基于正交匹配追踪的具有压缩感知过程的信号检测方法,它的检测过程为:
步骤一:设定预设迭代次数T,并设定余量rt的初始值r0和空矩阵V0,使r0=y,再计算传感矩阵V=ΦΨ,其中,y为采样值,Φ表示压缩感知采样过程,Ψ表示感兴趣信号s的变换域;
步骤三:根据公式Vt=[Vt-1vnt]更新矩阵Vt-1为Vt,其中,Vt表示第t次迭代后的传感矩阵;
步骤六:使t=t+1,判断t<T,如果是,返回执行步骤二,否则执行步骤七;
步骤七:求解稀疏系数估计值的无穷范数当时,选择H1假设,H1:y=Φ(s+n)=ΦsΦn=ΦΨθ+Φn=Vθ+Φn,所述H1假设说明感兴趣信号s存在;当≤γ时,选择H0假设,H0:y=Φn,其中,γ>0,所述γ的具体数值依据感兴趣信号s相对于加性高斯白噪声n的强度设定,H0假设说明感兴趣信号s不存在,完成对感兴趣信号s的检测。
具体实施方式二:本具体实施方式是对具体实施方式一所述的基于正交匹配追踪的具有压缩感知过程的信号检测方法的进一步说明,具体实施方式一中步骤二中所述的根据公式nt=arg<rt-1,vi>在传感矩阵V中选出与余量rt-1相关性最大的相关最大列向量vnt的具体过程为:依次求出余量rt-1与传感矩阵V中各列vi(i=1,…,N)的内积<rt-1,vi>,最大内积所对应的列向量序号为nt,1≤nt≤N,将列向量序号nt所对应的传感矩阵V中的列向量vi选作与第t次迭代中的余量rt-1相关性最大的相关最大列向量vnt,其中,迭代次数t为自然数且小于预设迭代次数T。
本具体实施方式的检测原理为:
检测的目的是区别两种假设:
H0:x=n vs. H1:x=s+n (1)
s表示我们感兴趣的信号,n表示加性高斯白噪声。感兴趣信号s在变换域Ψ内的稀疏度为K,即:
s=Ψsθs||θs||0=k (2)
其中Ψs和θs表示与感兴趣信号s相应的变换基与稀疏系数。由于加性高斯白噪声n在变换域内不稀疏,因此检测问题可以重新描述为区别(3)式中的两种假设:
H0:x=n vs.H1:x=Ψsθs+n (3)
根据(3)式,又可以通过判断θs是否存在来区别两种假设。即:
H0:θs=0vs.H1:θs≠0(4)
基于压缩感知理论,对感兴趣的x通过y=Фx获取,其中Φ为压缩感知采样过程。为了完成检测任务,需要利用一种检测方法从采样值y中提取特征量作为判决依据,用来区别两种假设。如果提取的特征量可以完成(4)式判断,那么就可以区别(4)中的两种假设,达到检测目的。
本具体实施方式,对采样值y进行压缩感知获得传感矩阵,对传感矩阵进行迭代分析,获得稀疏系数估计值对于检测任务,中已经包含了足够的信息,可以从中提取特征量作为判决依据。在感兴趣信号s存在时(H1情况下),该特征量为信号的一个重要分量在变换域内的稀疏系数估计值的绝对值;在感兴趣信号s不存在时(H0情况下),该特征值仅表示加性高斯白噪声在变换域内的一个投影系数。此外,由于加性高斯白噪声对的影响,本具体实施方式利用一个非零阈值γ作为判决条件,该检测用非零阈值γ可以根据先验条件及检测需求通过蒙特卡罗模拟方法选择。
结合本实施方式,对本实施方式的有效性进行仿真验证如下:
设感兴趣信号s长度N=256,由三个正弦分量叠加而成,信号时域波形如图2所示。该信号在频域是稀疏的,满足基于压缩感知信号检测问题对待检测信号的要求。噪声为加性白高斯噪声n。压缩感知过程通过一个M×N的观测矩阵Φ实现,在每次检测实验中,Φ是一个随机产生的元素满足高斯分布的矩阵,M为压缩感知采样点数。仿真实验中,先验概率Pr(H0)=Pr(H1)=1/2,检测成功率为2000次检测实验的统计结果。
令采样点数M=50,信噪比SNR=20,方法迭代次数T=6,依据先验条件检测阈值选择范围为[2.5,3.5],阈值步进为0.1。图3给出了应用本具体实施方式对感兴趣信号s进行检测的仿真结果;由图3的仿真结果可以看出,基于正交匹配追踪的具有压缩感知过程的信号检测方法可以较好的完成检测任务,同时计算量小。
Claims (2)
1.基于正交匹配追踪的具有压缩感知过程的信号检测方法,其特征在于它的检测过程为:
步骤一:设定预设迭代次数T,并设定余量rt的初始值r0和空矩阵V0,使r0=y,再计算传感矩阵V=ΦΨ,其中,y为采样值,Φ表示压缩感知采样过程,Ψ表示感兴趣信号s的变换域;
步骤二:根据公式在传感矩阵V中选出与余量rt-1相关性最大的相关最大列向量vnt;
步骤三:根据公式Vt=[Vt-1vnt]更新矩阵Vt-1为Vt,其中,Vt表示第t次迭代后的传感矩阵;
步骤五:根据公式更新余量rt-1为rt;
步骤六:使t=t+1,判断t<T,如果是,返回执行步骤二,否则执行步骤七;
2.根据权利要求1所述的基于正交匹配追踪的具有压缩感知过程的信号检测方法,其特征在于步骤二所述的根据公式在传感矩阵V中选出与余量rt-1相关性最大的相关最大列向量vnt的具体过程为:依次求出余量rt-1与传感矩阵V中各列vi(i=1,…,N)的内积<rt-1,vi>,最大内积所对应的列向量序号为nt,1≤nt≤N,将列向量序号nt所对应的传感矩阵V中的列向量vi选作与第t次迭代中的余量rt-1相关性最大的相关最大列向量vnt,其中,迭代次数t为自然数且小于预设迭代次数T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100324854A CN101764652B (zh) | 2010-01-18 | 2010-01-18 | 基于正交匹配追踪的具有压缩感知过程的信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100324854A CN101764652B (zh) | 2010-01-18 | 2010-01-18 | 基于正交匹配追踪的具有压缩感知过程的信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101764652A true CN101764652A (zh) | 2010-06-30 |
CN101764652B CN101764652B (zh) | 2012-12-19 |
Family
ID=42495650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010100324854A Expired - Fee Related CN101764652B (zh) | 2010-01-18 | 2010-01-18 | 基于正交匹配追踪的具有压缩感知过程的信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101764652B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908890A (zh) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | 分块大小未知的块稀疏信号的盲重构方法 |
CN101908889A (zh) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | 一种块稀疏度未知的稀疏信号压缩感知重构方法 |
CN102034250A (zh) * | 2010-11-26 | 2011-04-27 | 西安电子科技大学 | 基于边缘结构信息的分块压缩感知重构方法 |
CN102045118A (zh) * | 2010-10-22 | 2011-05-04 | 清华大学 | 基于伪逆相乘的压缩感知重建算法 |
CN102291152A (zh) * | 2011-04-07 | 2011-12-21 | 湖南大学 | 基于奇异值分解的压缩感知含噪信号重构系统 |
CN102565737A (zh) * | 2011-12-12 | 2012-07-11 | 中国科学院深圳先进技术研究院 | 快速磁共振成像方法及系统 |
CN102624399A (zh) * | 2012-03-30 | 2012-08-01 | 北京邮电大学 | 一种压缩感知信号的重构方法 |
CN103036574A (zh) * | 2012-12-13 | 2013-04-10 | 南开大学 | 基于压缩传感的一种自校验稀疏度自适应匹配追踪算法 |
CN103051403A (zh) * | 2013-01-29 | 2013-04-17 | 哈尔滨工业大学 | 基于多个mwc分布式亚奈奎斯特采样联合重构的频谱感知方法 |
CN103532567A (zh) * | 2013-11-01 | 2014-01-22 | 哈尔滨工业大学 | 分布式压缩感知框架下基于快速计算内积的正交匹配追踪算法的信号重构方法 |
CN103679762A (zh) * | 2013-12-13 | 2014-03-26 | 江苏大学 | 一种基于稀疏数据的超声信号重建方法 |
CN103780264A (zh) * | 2014-01-24 | 2014-05-07 | 浙江大学 | 一种基于fpga的压缩感知方法 |
CN103873170A (zh) * | 2014-03-26 | 2014-06-18 | 哈尔滨工业大学 | 一种在盲稀疏条件下压缩感知的频谱检测方法 |
CN104104389A (zh) * | 2013-04-09 | 2014-10-15 | 华为技术有限公司 | 一种信号重建方法及设备 |
CN104215456A (zh) * | 2014-08-20 | 2014-12-17 | 昆明理工大学 | 一种基于平面聚类和频域压缩感知重构的机械故障诊断方法 |
CN104904122A (zh) * | 2013-12-31 | 2015-09-09 | 华为技术有限公司 | 基于压缩感知的信号重建方法及装置 |
CN105099973A (zh) * | 2015-08-25 | 2015-11-25 | 西安电子科技大学 | 用于宽带调制变换压缩采集架构的测量矩阵校准方法 |
CN105471529A (zh) * | 2015-12-30 | 2016-04-06 | 北京邮电大学 | 频谱信号感知方法及装置 |
CN108649961A (zh) * | 2018-05-08 | 2018-10-12 | 北京理工大学 | 一种基于边信息先验的多导联心电信号重构估计方法 |
CN109683051A (zh) * | 2018-12-27 | 2019-04-26 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于稀疏信号重构的小电流接地故障选线方法 |
CN109787924A (zh) * | 2019-03-13 | 2019-05-21 | 重庆邮电大学 | 一种基于压缩感知的LoRa信道估计方法 |
CN109884581A (zh) * | 2019-02-22 | 2019-06-14 | 南昌大学 | 一种基于任意阵列的太赫兹宽空域无源测向方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090110323A (ko) * | 2007-01-04 | 2009-10-21 | 브리티쉬 텔리커뮤니케이션즈 파블릭 리미티드 캄퍼니 | 비디오 신호를 인코딩하는 방법 및 시스템 |
CN101494627B (zh) * | 2009-03-11 | 2013-06-05 | 北京邮电大学 | 宽带移动通信中利用压缩感知减少导频数的信道估计方法 |
CN101572897A (zh) * | 2009-03-24 | 2009-11-04 | 复旦大学 | 分布式认知无线电信道探测方法和系统 |
-
2010
- 2010-01-18 CN CN2010100324854A patent/CN101764652B/zh not_active Expired - Fee Related
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908890B (zh) * | 2010-07-30 | 2012-12-26 | 哈尔滨工业大学 | 分块大小未知的块稀疏信号的盲重构方法 |
CN101908889A (zh) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | 一种块稀疏度未知的稀疏信号压缩感知重构方法 |
CN101908890A (zh) * | 2010-07-30 | 2010-12-08 | 哈尔滨工业大学 | 分块大小未知的块稀疏信号的盲重构方法 |
CN101908889B (zh) * | 2010-07-30 | 2013-03-20 | 哈尔滨工业大学 | 一种块稀疏度未知的稀疏信号压缩感知重构方法 |
CN102045118A (zh) * | 2010-10-22 | 2011-05-04 | 清华大学 | 基于伪逆相乘的压缩感知重建算法 |
CN102045118B (zh) * | 2010-10-22 | 2014-06-18 | 清华大学 | 基于伪逆相乘的压缩感知重建算法 |
CN102034250B (zh) * | 2010-11-26 | 2012-05-23 | 西安电子科技大学 | 基于边缘结构信息的分块压缩感知重构方法 |
CN102034250A (zh) * | 2010-11-26 | 2011-04-27 | 西安电子科技大学 | 基于边缘结构信息的分块压缩感知重构方法 |
CN102291152A (zh) * | 2011-04-07 | 2011-12-21 | 湖南大学 | 基于奇异值分解的压缩感知含噪信号重构系统 |
CN102291152B (zh) * | 2011-04-07 | 2014-06-04 | 湖南大学 | 基于奇异值分解的压缩感知含噪信号重构系统 |
CN102565737A (zh) * | 2011-12-12 | 2012-07-11 | 中国科学院深圳先进技术研究院 | 快速磁共振成像方法及系统 |
CN102624399A (zh) * | 2012-03-30 | 2012-08-01 | 北京邮电大学 | 一种压缩感知信号的重构方法 |
CN102624399B (zh) * | 2012-03-30 | 2014-10-22 | 北京邮电大学 | 一种压缩感知信号的重构方法 |
CN103036574B (zh) * | 2012-12-13 | 2017-03-08 | 南开大学 | 基于压缩传感的一种自校验稀疏度自适应匹配追踪算法 |
CN103036574A (zh) * | 2012-12-13 | 2013-04-10 | 南开大学 | 基于压缩传感的一种自校验稀疏度自适应匹配追踪算法 |
CN103051403A (zh) * | 2013-01-29 | 2013-04-17 | 哈尔滨工业大学 | 基于多个mwc分布式亚奈奎斯特采样联合重构的频谱感知方法 |
CN104104389B (zh) * | 2013-04-09 | 2017-08-18 | 华为技术有限公司 | 一种信号重建方法及设备 |
CN104104389A (zh) * | 2013-04-09 | 2014-10-15 | 华为技术有限公司 | 一种信号重建方法及设备 |
CN103532567A (zh) * | 2013-11-01 | 2014-01-22 | 哈尔滨工业大学 | 分布式压缩感知框架下基于快速计算内积的正交匹配追踪算法的信号重构方法 |
CN103532567B (zh) * | 2013-11-01 | 2016-12-07 | 哈尔滨工业大学 | 分布式压缩感知框架下基于快速计算内积的正交匹配追踪算法的信号重构方法 |
CN103679762A (zh) * | 2013-12-13 | 2014-03-26 | 江苏大学 | 一种基于稀疏数据的超声信号重建方法 |
CN104904122A (zh) * | 2013-12-31 | 2015-09-09 | 华为技术有限公司 | 基于压缩感知的信号重建方法及装置 |
CN103780264A (zh) * | 2014-01-24 | 2014-05-07 | 浙江大学 | 一种基于fpga的压缩感知方法 |
CN103873170B (zh) * | 2014-03-26 | 2016-02-03 | 哈尔滨工业大学 | 一种在盲稀疏条件下压缩感知的频谱检测方法 |
CN103873170A (zh) * | 2014-03-26 | 2014-06-18 | 哈尔滨工业大学 | 一种在盲稀疏条件下压缩感知的频谱检测方法 |
CN104215456B (zh) * | 2014-08-20 | 2017-02-15 | 昆明理工大学 | 一种基于平面聚类和频域压缩感知重构的机械故障诊断方法 |
CN104215456A (zh) * | 2014-08-20 | 2014-12-17 | 昆明理工大学 | 一种基于平面聚类和频域压缩感知重构的机械故障诊断方法 |
CN105099973B (zh) * | 2015-08-25 | 2018-03-20 | 西安电子科技大学 | 用于宽带调制变换压缩采集架构的测量矩阵校准方法 |
CN105099973A (zh) * | 2015-08-25 | 2015-11-25 | 西安电子科技大学 | 用于宽带调制变换压缩采集架构的测量矩阵校准方法 |
CN105471529A (zh) * | 2015-12-30 | 2016-04-06 | 北京邮电大学 | 频谱信号感知方法及装置 |
CN105471529B (zh) * | 2015-12-30 | 2018-02-02 | 北京邮电大学 | 频谱信号感知方法及装置 |
CN108649961A (zh) * | 2018-05-08 | 2018-10-12 | 北京理工大学 | 一种基于边信息先验的多导联心电信号重构估计方法 |
CN108649961B (zh) * | 2018-05-08 | 2019-04-02 | 北京理工大学 | 一种基于边信息先验的多导联心电信号重构估计方法 |
CN109683051A (zh) * | 2018-12-27 | 2019-04-26 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于稀疏信号重构的小电流接地故障选线方法 |
CN109884581A (zh) * | 2019-02-22 | 2019-06-14 | 南昌大学 | 一种基于任意阵列的太赫兹宽空域无源测向方法 |
CN109787924A (zh) * | 2019-03-13 | 2019-05-21 | 重庆邮电大学 | 一种基于压缩感知的LoRa信道估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101764652B (zh) | 2012-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101764652A (zh) | 基于正交匹配追踪的具有压缩感知过程的信号检测方法 | |
CN110632572B (zh) | 基于无意调相特性的雷达辐射源个体识别方法及装置 | |
Tewes et al. | Cosmograil: the cosmological monitoring of gravitational lenses-xi. techniques for time delay measurement in presence of microlensing | |
CN110619296B (zh) | 一种基于奇异分解的信号降噪方法 | |
CN104375976B (zh) | 基于张量正则分解的欠定盲源分离中的混合矩阵识别方法 | |
Li et al. | A mixing matrix estimation algorithm for underdetermined blind source separation | |
White et al. | Methodological tools | |
Xu et al. | Asymptotic properties of Pearson׳ s rank-variate correlation coefficient in bivariate normal model | |
CN107133476B (zh) | 一种基于响应混叠性度量的测试激励与测点的协同优化方法 | |
Kwon et al. | Moderately clipped lasso | |
Meng et al. | Measurement uncertainty of complex-valued microwave quantities | |
GB2476911A (en) | Methods for determining orientation of electromagnetic receivers | |
CN107884752A (zh) | 一种基于压缩感知的外辐射源雷达对目标检测方法 | |
CN101944236A (zh) | 基于采样值分布特性的具有压缩感知过程的信号检测方法 | |
CN105446941A (zh) | 基于联合零对角化的时频域盲信号分离方法 | |
CN104104389B (zh) | 一种信号重建方法及设备 | |
CN107346300B (zh) | 一种基于绝对传递率函数的传递路径分析方法 | |
Jiang | Estimation of effective concentrations from in vitro dose-response data using the log-logistic model | |
Spagnolo et al. | Forensic metrology: uncertainty of measurements in forensic analysis | |
Aouada et al. | Application of the bootstrap to source detection in nonuniform noise | |
Arias-Giraldo et al. | enmpa: An R package for ecological niche modeling using presence-absence data and generalized linear models | |
CN107888298B (zh) | 一种射频隐身信号不确定性的测试方法 | |
CN112285791A (zh) | 一种探地雷达电磁场数据处理方法、装置以及介质 | |
CN112649865B (zh) | 基于特征子波的火山岩储层识别方法及系统 | |
Armijo et al. | Cosmological constraints using Minkowski functionals from the first year data of the Hyper Suprime-Cam |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121219 Termination date: 20140118 |