CN111541465A - 一种信号捕获方法、装置、电子设备及介质 - Google Patents

一种信号捕获方法、装置、电子设备及介质 Download PDF

Info

Publication number
CN111541465A
CN111541465A CN202010303884.3A CN202010303884A CN111541465A CN 111541465 A CN111541465 A CN 111541465A CN 202010303884 A CN202010303884 A CN 202010303884A CN 111541465 A CN111541465 A CN 111541465A
Authority
CN
China
Prior art keywords
sparse
matrix
value
observation
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
Application number
CN202010303884.3A
Other languages
English (en)
Other versions
CN111541465B (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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN202010303884.3A priority Critical patent/CN111541465B/zh
Publication of CN111541465A publication Critical patent/CN111541465A/zh
Application granted granted Critical
Publication of CN111541465B publication Critical patent/CN111541465B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • H04B1/70751Synchronisation aspects with code phase acquisition using partial detection
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • H04B1/70751Synchronisation aspects with code phase acquisition using partial detection
    • H04B1/70752Partial correlation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • H04B1/70751Synchronisation aspects with code phase acquisition using partial detection
    • H04B1/70753Partial phase search

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明实施例提供了一种信号捕获方法、装置、电子设备及介质,涉及通信定位技术领域,可以减少信号捕获所消耗的资源和时间。本发明包括:利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值,稀疏向量为稀疏矩阵与接收到的信号的自相关结果,稀疏矩阵包括零元素和非零元素,稀疏矩阵包括的非零元素为本地码矩阵包括的元素,稀疏向量包括的每个元素的索引对应本地码矩阵包括的一个本地码。然后依据第一观测值,确定稀疏向量中的稀疏值的索引。进而获得索引对应的目标本地码的部分匹配滤波结果,将部分匹配滤波结果进行快速傅里叶变换运算,得到信号的频谱的峰值,当峰值大于峰值门限时,获得目标本地码的码相位以及峰值对应的残余多普勒。

Description

一种信号捕获方法、装置、电子设备及介质
技术领域
本发明涉及通信定位技术领域,特别是涉及一种信号捕获方法、装置、电子设备及介质。
背景技术
地面通信网络的时分码分正交频分复用(Time and Code Division-OrthogonalFrequency Division Multiplexing,TC-OFDM)信号采用直接扩频序列(Direct SequenceSpread Spectrum,DSSS)作为定位基础。DSSS具有测距性能好(自相关性)和码分多址(互相关性)的优势而被应用于TC-OFDM定位系统。对于定位终端的基带处理部分,首要任务是完成信号的捕获。
现有的基于压缩感知(Compressed Sensing,CS)的信号捕获方法是通过对具有稀疏性的信号进行随机非均匀压缩采样,再对压缩后的信号进行重构获得原始信号。然而这种方法只有对从高维转换得到的低维信号进行信号重构,才能恢复出原始信号,而信号重构需要消耗大量的资源和时间。
发明内容
本发明实施例的目的在于提供一种信号捕获方法、装置、电子设备及介质,以减少信号捕获所消耗的资源和时间。具体技术方案如下:
第一方面,本发明实施例提供了一种信号捕获方法,所述方法包括:
利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值,所述稀疏向量为稀疏矩阵与接收到的信号的自相关结果,所述稀疏矩阵包括零元素和非零元素,所述稀疏矩阵包括的非零元素为本地码矩阵包括的元素,所述本地码矩阵包括多个码相位不同的本地码,所述稀疏向量包括的每个元素的索引对应所述本地码矩阵包括的一个本地码;
依据所述第一观测值,确定所述稀疏向量中的稀疏值的索引,所述稀疏值为非零值;
获得所述索引对应的目标本地码的部分匹配滤波结果,将所述部分匹配滤波结果进行快速傅里叶变换FFT运算,得到所述信号的频谱的峰值,当所述峰值大于峰值门限时,获得所述目标本地码的码相位以及所述峰值对应的残余多普勒。
可选的,所述第一观测矩阵包括的元素为:
Figure BDA0002455040130000021
其中,
Figure BDA0002455040130000022
为所述第一观测矩阵包括的元素,M1为预设压缩比例,HM01为第一哈达玛矩阵,M1=2M01,m∈{1,2,…,M1},l∈{1,2,…,2ML},α=2ML/M1,M为所述信号的分段长度,L为本地码周期。
可选的,所述依据所述第一观测值,确定所述稀疏向量中的稀疏值的索引,包括:
对所述第一观测值进行沃尔什Walsh变换,得到第一沃尔什向量;
确定所述第一沃尔什向量中,数值大于第一沃尔什门限的第一元素向量;
根据所述第一元素向量的索引组成的索引集,确定第二观测矩阵;
利用所述第二观测矩阵,对所述稀疏向量进行压缩感知,获得第二观测值;
对所述第二观测值进行Walsh变换,得到第二沃尔什向量;
确定所述第二沃尔什向量中,数值大于第二沃尔什门限的第二元素向量;
根据第二元素向量的索引,确定所述稀疏向量中的稀疏值的索引。
可选的,所述第二观测矩阵包括的元素为:
Figure BDA0002455040130000023
其中,M2为预设的峰值搜索长度,HM02为第二哈达玛矩阵,A1(i)为所述索引集,M2=2M02
Figure BDA0002455040130000024
第二方面,本发明实施例提供了一种信号捕获装置,所述装置包括:
压缩模块,用于利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值,所述稀疏向量为稀疏矩阵与接收到的信号的自相关结果,所述稀疏矩阵包括零元素和非零元素,所述稀疏矩阵包括的非零元素为本地码矩阵包括的元素,所述本地码矩阵包括多个码相位不同的本地码,所述稀疏向量包括的每个元素的索引对应所述本地码矩阵包括的一个本地码;
确定模块,用于依据所述压缩模块获得的所述第一观测值,确定所述稀疏向量中的稀疏值的索引,所述稀疏值为非零值;
获得模块,用于获得所述索引对应的目标本地码的部分匹配滤波结果,将所述部分匹配滤波结果进行快速傅里叶变换FFT运算,得到所述信号的频谱的峰值,当所述峰值大于峰值门限时,获得所述目标本地码的码相位以及所述峰值对应的残余多普勒。
可选的,所述第一观测矩阵包括的元素为:
Figure BDA0002455040130000031
其中,
Figure BDA0002455040130000032
为所述第一观测矩阵包括的元素,M1为预设压缩比例,HM01为第一哈达玛矩阵,M1=2M01,m∈{1,2,…,M1},l∈{1,2,…,2ML},α=2ML/M1,M为所述信号的分段长度,L为本地码周期。
可选的,所述确定模块,具体用于:
对所述第一观测值进行沃尔什Walsh变换,得到第一沃尔什向量;
确定所述第一沃尔什向量中,数值大于第一沃尔什门限的第一元素向量;
根据所述第一元素向量的索引组成的索引集,确定第二观测矩阵;
利用所述第二观测矩阵,对所述稀疏向量进行压缩感知,获得第二观测值;
对所述第二观测值进行Walsh变换,得到第二沃尔什向量;
确定所述第二沃尔什向量中,数值大于第二沃尔什门限的第二元素向量;
根据第二元素向量的索引,确定所述稀疏向量中的稀疏值的索引。
可选的,所述第二观测矩阵包括的元素为:
Figure BDA0002455040130000041
其中,M2为预设的峰值搜索长度,HM02为第二哈达玛矩阵,A1(i)为所述索引集,M2=2M02
Figure BDA0002455040130000042
第三方面,本发明实施例提供了一种电子设备,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现上述任一信号捕获方法的步骤。
第四方面,本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现上述任一信号捕获方法的步骤。
第五方面,本发明实施例还提供了一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行上述任一信号捕获方法。
本发明实施例的技术方案至少包括以下有益效果:由于本发明实施例在对稀疏向量进行压缩感知之后,依据压缩感知得到的第一观测值,确定稀疏向量中的稀疏值的索引,当信号的频谱的峰值大于峰值门限时,说明该索引对应的目标本地码与信号包括的本地码相同,此时得到目标本地码的码相位以及峰值对应的残余多普勒。本发明实施例在捕获信号的过程中,未对信号进行信号重构,从而减少了信号捕获所消耗的资源和时间。
当然,实施本发明的任一产品或方法并不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种信号捕获方法流程图;
图2为本发明实施例提供的确定稀疏向量中的稀疏值的索引的方法流程图;
图3为本发明实施例提供的一种信号捕获装置结构示意图;
图4为本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了提高信号捕获速度,本发明实施例提供了一种信号捕获方法,应用于电子设备,其中,电子设备可以是手机、计算机、平板电脑等设备,参见图1,该方法包括如下步骤。
步骤101,利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值。
其中,稀疏向量为稀疏矩阵与接收到的信号的自相关结果,稀疏矩阵包括零元素和非零元素,稀疏矩阵包括的非零元素为本地码矩阵包括的元素,本地码矩阵包括多个码相位不同的本地码,稀疏向量包括的每个元素的索引对应本地码矩阵包括的一个本地码。
步骤102,依据第一观测值,确定稀疏向量中的稀疏值的索引。
其中,稀疏值为非零值。
步骤103,获得该索引对应的目标本地码的部分匹配滤波(Partial matchedfiltering,PMF)结果,将获得的部分匹配滤波结果进行快速傅里叶变换(Fast FourierTransform,FFT)运算,得到信号的频谱的峰值,当峰值大于峰值门限时,获得目标本地码的码相位以及峰值对应的残余多普勒。
可以理解的,由于接收到的信号是发送端生成的本地码与载波的叠加信号,经过部分匹配滤波后得到的是信号的载波,步骤103中的峰值为载波的频谱的峰值。
可以理解的,当信号的频谱的峰值大于峰值门限时,说明目标本地码与发送端的本地码相同,目标本地码的码相位是正确的码相位,峰值对应的残余多普勒也是正确的。
信号捕获的目的是获得正确的码相位以及残余多普勒,在获得这两个参数时,信号捕获成功。
本发明实施例的技术方案至少包括以下有益效果:由于本发明实施例在对稀疏向量进行压缩感知之后,依据压缩感知得到的第一观测值,确定稀疏向量中的稀疏值的索引,当信号的频谱的峰值大于峰值门限时,说明该索引对应的目标本地码与信号包括的本地码相同,此时得到目标本地码的码相位以及峰值对应的残余多普勒。本发明实施例在捕获信号的过程中,未对信号进行信号重构,从而减少了信号捕获所消耗的资源和时间。
由压缩感知理论可知,如果信号是稀疏的,那么可以由远低于采样定理要求的采样频率对该信号进行重建恢复。码相位搜索利用的本地码矩阵可以表示为公式(1)。
Figure BDA0002455040130000061
其中,C为本地码矩阵,C中的每一行元素表示一个本地码。
可以理解的,C中的每一行元素表示一个本地码,第一行本地码的码相位为0,第二行本地码的码相位为1,以此类推。示例性的,本地码可以是伪随机码。
为了构建稀疏矩阵,使得稀疏矩阵中非零元素的个数远小于零元素的个数,部分匹配滤波器需要对本地码矩阵C的每一行元素进行分段,得到稀疏矩阵CPMF,分段后每一行分段矩阵可以表示为公式(2)。
Figure BDA0002455040130000071
其中,CiPMF为CPMF的第i行分段矩阵,K和M为预设的稀疏参数。
稀疏矩阵CPMF的矩阵大小为2ML×2L,稀疏向量z=CPMF×s,其中,s为接收到的信号,稀疏向量z中包括连续M个非零值,其余值均为0(将接近0的值也看做值为0)。
可以理解的,当长度为N的信号中非零值的数量为K,其他值都为0(将接近0的值也看做值为0)时,可以认为信号具有稀疏性,其稀疏度为K。通常,信号不具有稀疏性,但信号在进行矩阵变化后,可以具有稀疏性,矩阵变换后的信号可以表示为公式(3)。
x=Ψs (3)
其中,x为矩阵变换后的信号,Ψ为稀疏矩阵,s为矩阵变换前的信号。||s||0≤K,||·||0为零范数。
压缩感知通过对一个矩阵大小为M×N的观测矩阵对信号进行随机非均匀采样,得到压缩采样矩阵y。其中,M<<N,观测矩阵与稀疏矩阵的稀疏基不相关。
y=Φx (4)
其中,y为压缩采样矩阵,Φ为观测矩阵,x为矩阵变换后的信号,Θ=ΦΨ为感知矩阵,Ψ为稀疏矩阵。
由于信号本身不具备稀疏性,但信号具备良好的自相关性能,在码相位搜索域上,当接收到的信号与本地码之间码相位对齐时,接收到的信号与本地码的自相关结果达到峰值,因此自相关结果具备一定的稀疏性。
基于此,在本发明实施例中,在步骤101之前,可以构造一个M1×2ML的第一观测矩阵Φ1,第一观测矩阵的压缩比为α=2ML/M1。第一观测矩阵φ1包括的元素为公式(5)。
Figure BDA0002455040130000072
其中,
Figure BDA0002455040130000081
为第一观测矩阵包括的元素,M1为预设压缩比例,HM01为第一哈达玛(Hadamard)矩阵,M1=2M01,m∈{1,2,…,M1},l∈{1,2,…,2ML},α=2ML/M1,M为信号的分段长度,L为本地码周期。
可以理解的,由哈达玛矩阵的性质可知,哈达玛矩阵为2的整数次幂的方阵,M1=2M01,且HM01的矩阵大小为M1×M1
可选的,可以通过公式(6),利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值。
y1=Φ1z=Φ1CPMF×s=Φ′1s (6)
其中,y1为第一观测值,y1的矩阵大小为M1×1,Φ1为第一观测矩阵,z为稀疏向量,CPMF为稀疏矩阵,Φ′1为第一感知矩阵,s为接收到的信号。
可以理解的,公式(6)将稀疏向量z中连续α个元素的值表征为y1中的一个元素。具体为:将PMF的α段结果累加为一个值,得到大小为2L/M1的自相关结果。
由自相关公式得到当2L/M1大于3时,压缩感知后的自相关结果中包含了自相关结果的大部分峰值。而L为本地码周期,L的大小是已知的,可以依据2L/M1>3,设置M1的大小。
可以理解的,压缩感知使用了M1个并行相关器,与传统的捕获方法相比,压缩感知将本地码与接收到的信号的自相关结果压缩了α倍。
可选的,参见图2,上述步骤102依据第一观测值,确定稀疏向量中的稀疏值的索引,可以包括以下步骤。
步骤201,对第一观测值进行沃尔什(Walsh)变换,得到第一沃尔什向量。
可以理解的,由于哈达玛矩阵具有良好的互相关无关性,对于y1进行Walsh变换,得到Walsh变换后得到的第一沃尔什向量wh1,获得y1中每一列Hadamard对应的系数。
步骤202,确定第一沃尔什向量中,数值大于第一沃尔什门限的第一元素向量。
步骤203,根据第一元素向量的索引组成的索引集,确定第二观测矩阵。
可选的,大于第一沃尔什门限Vth1的沃尔什向量,在wh1中的索引组成的索引集为A1,A1中包含了本地码正确的码相位。A1中的每个索引能够表示α各相位信息,因此可以对A1包括的索引信息展开。
A′i(i)={α(A1(i)-1)+1,α(A1(i)-1)+2,…,min[αA1(i),2L]} (7)
其中,A′i(i)为展开后的索引集,A1(i)为展开前的索引集,α=2ML/M1,M1为预设压缩比例,M为信号的分段长度,L为本地码周期。
A′1中包含以半个码片为精度的码相位位置,第二观测矩阵包括的元素为公式(8)。
Figure BDA0002455040130000091
其中,
Figure BDA0002455040130000092
为第二观测矩阵包括的元素,HM02为第二哈达玛矩阵,M2为第二阶段预设的峰值搜索长度,M2=2M02
Figure BDA0002455040130000093
lb为以2为底的对数,||·||0为零范数,M为信号的分段长度,L为本地码周期,HM02的矩阵大小为M2×M2
步骤204,利用第二观测矩阵,对稀疏向量进行压缩感知,获得第二观测值。
可选的,可以通过公式(9)获得第二观测值。
y2=Φ2z=Φ2CPMF×s=Φ′2s (9)
其中,y2为第二观测值,Φ2为第二观测矩阵,z为稀疏向量,CPMF为稀疏矩阵,Φ′2为第二感知矩阵,Φ′2的矩阵大小为M2×2L,s为接收到的信号。
步骤205,对第二观测值进行Walsh变换,得到第二沃尔什向量。
可以理解的,由于每个部分匹配滤波器对应矩阵y2的一列,且y2的各个列具有正交性,通过对y2进行Walsh变换,可以得到y2的每一列对应的系数。
步骤206,确定第二沃尔什向量中,数值大于第二沃尔什门限的第二元素向量。
由于y2为PMF结果且带有残余多普勒,因此大于第二沃尔什门限Vth2的元素的索引组成的索引集合A2为长度小于M的向量。
步骤207,根据第二元素向量的索引,确定第一观测值中的稀疏值的索引。
可选的,稀疏向量中的稀疏值的索引为
Figure BDA0002455040130000101
即为部分匹配滤波的所有滤波器在CPMF中的索引起始位置,同时对应着码相位结果。
例如,结合公式(1),当
Figure BDA0002455040130000102
时,该索引对应的码相位为1,该索引对应着本地码矩阵的第二行,因此目标本地码为:c2,c3,…,c2L,c1
本发明实施例的技术方案还可以带来以下有益效果:在PMF-FFT捕获算法中应用压缩感知的信号捕获方法,不仅具备PMF-FFT在频域搜索速度快、精度高的优点,还极大的改善了PMF在码域搜索时消耗资源大,耗费时间长的问题。另外,经过理论分析和实验仿真均表明,本发明实施例提供的信号捕获方法可以有效地减少所需相关器的数量至传统PMF-FFT算法的50%以下。
针对上述步骤103、获得该索引对应的目标本地码的部分匹配滤波结果,将获得的部分匹配滤波结果进行FFT运算,得到信号的频谱的峰值,当峰值大于峰值门限时,获得目标本地码的码相位以及频率峰值对应的残余多普勒。
在本发明实施例中,当获得的索引为i时,索引i对应的目标本地码的部分匹配滤波结果为CiPMF×s,其中,CiPMF为第i个本地码对应的分段矩阵,s为接收到的信号。
部分匹配滤波结果进行FFT运算后,得到载波的频谱,频谱纵轴为幅值,横轴为频率,频谱的峰值为频谱的幅值最高值,残余多普勒为频谱的幅值最高值对应的频率。
对应于上述方法实施例,本发明实施例还提供了一种信号捕获装置,参见图3,该装置包括:压缩模块301、确定模块302和获得模块303。
压缩模块301,用于利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值,稀疏向量为稀疏矩阵与接收到的信号的自相关结果,稀疏矩阵包括零元素和非零元素,稀疏矩阵包括的非零元素为本地码矩阵包括的元素,本地码矩阵包括多个码相位不同的本地码,稀疏向量包括的每个元素的索引对应本地码矩阵包括的一个本地码;
确定模块302,用于依据压缩模块301获得的第一观测值,确定稀疏向量中的稀疏值的索引,稀疏值为非零值;
获得模块303,用于获得确定模块302确定的索引对应的目标本地码的部分匹配滤波结果,将部分匹配滤波结果进行快速傅里叶变换FFT运算,得到信号的频谱的峰值,当峰值大于峰值门限时,获得目标本地码的码相位以及峰值对应的残余多普勒。
可选的,第一观测矩阵包括的元素为:
Figure BDA0002455040130000111
其中,
Figure BDA0002455040130000112
为第一观测矩阵包括的元素,M1为预设压缩比例,HM01为第一哈达玛矩阵,M1=2M01,m∈{1,2,…,M1},l∈{1,2,…,2ML},α=2ML/M1,M为信号的分段长度,L为本地码周期。
可选的,确定模块302,具体用于:
对第一观测值进行沃尔什Walsh变换,得到第一沃尔什向量;
确定第一沃尔什向量中,数值大于第一沃尔什门限的第一元素向量;
根据第一元素向量的索引组成的索引集,确定第二观测矩阵;
利用第二观测矩阵,对稀疏向量进行压缩感知,获得第二观测值;
对第二观测值进行Walsh变换,得到第二沃尔什向量;
确定第二沃尔什向量中,数值大于第二沃尔什门限的第二元素向量;
根据第二元素向量的索引,确定稀疏向量中的稀疏值的索引。
可选的,第二观测矩阵包括的元素为:
Figure BDA0002455040130000113
其中,M2为预设的峰值搜索长度,HM02为第二哈达玛矩阵,A1(i)为索引集,M2=2M02
Figure BDA0002455040130000114
本发明实施例还提供了一种电子设备,如图4所示,包括处理器401、通信接口402、存储器403和通信总线404,其中,处理器401,通信接口402,存储器403通过通信总线404完成相互间的通信,
存储器403,用于存放计算机程序;
处理器401,用于执行存储器403上所存放的程序时,实现上述方法实施例中的方法步骤。
上述电子设备提到的通信总线可以是外设部件互连标准(Peripheral ComponentInterconnect,PCI)总线或扩展工业标准结构(Extended Industry StandardArchitecture,EISA)总线等。该通信总线可以分为地址总线、数据总线、控制总线等。为便于表示,图中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。
通信接口用于上述电子设备与其他设备之间的通信。
存储器可以包括随机存取存储器(Random Access Memory,RAM),也可以包括非易失性存储器(Non-Volatile Memory,NVM),例如至少一个磁盘存储器。可选的,存储器还可以是至少一个位于远离前述处理器的存储装置。
上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(Network Processor,NP)等;还可以是数字信号处理器(Digital SignalProcessing,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
在本发明提供的又一实施例中,还提供了一种计算机可读存储介质,该计算机可读存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现上述任一信号捕获方法的步骤。
在本发明提供的又一实施例中,还提供了一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行上述实施例中任一信号捕获方法。
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘Solid State Disk(SSD))等。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的较佳实施例,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。

Claims (10)

1.一种信号捕获方法,其特征在于,所述方法包括:
利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值,所述稀疏向量为稀疏矩阵与接收到的信号的自相关结果,所述稀疏矩阵包括零元素和非零元素,所述稀疏矩阵包括的非零元素为本地码矩阵包括的元素,所述本地码矩阵包括多个码相位不同的本地码,所述稀疏向量包括的每个元素的索引对应所述本地码矩阵包括的一个本地码;
依据所述第一观测值,确定所述稀疏向量中的稀疏值的索引,所述稀疏值为非零值;
获得所述索引对应的目标本地码的部分匹配滤波结果,将所述部分匹配滤波结果进行快速傅里叶变换FFT运算,得到所述信号的频谱的峰值,当所述峰值大于峰值门限时,获得所述目标本地码的码相位以及所述峰值对应的残余多普勒。
2.根据权利要求1所述的方法,其特征在于,所述第一观测矩阵包括的元素为:
Figure FDA0002455040120000011
其中,
Figure FDA0002455040120000012
为所述第一观测矩阵包括的元素,M1为预设压缩比例,HM01为第一哈达玛矩阵,M1=2M01,m∈{1,2,…,M1},l∈{1,2,…,2ML},α=2ML/M1,M为所述信号的分段长度,L为本地码周期。
3.根据权利要求1所述的方法,其特征在于,所述依据所述第一观测值,确定所述稀疏向量中的稀疏值的索引,包括:
对所述第一观测值进行沃尔什Walsh变换,得到第一沃尔什向量;
确定所述第一沃尔什向量中,数值大于第一沃尔什门限的第一元素向量;
根据所述第一元素向量的索引组成的索引集,确定第二观测矩阵;
利用所述第二观测矩阵,对所述稀疏向量进行压缩感知,获得第二观测值;
对所述第二观测值进行Walsh变换,得到第二沃尔什向量;
确定所述第二沃尔什向量中,数值大于第二沃尔什门限的第二元素向量;
根据第二元素向量的索引,确定所述稀疏向量中的稀疏值的索引。
4.根据权利要求3所述的方法,其特征在于,所述第二观测矩阵包括的元素为:
Figure FDA0002455040120000021
其中,M2为预设的峰值搜索长度,HM02为第二哈达玛矩阵,A1(i)为所述索引集,M2=2M02
Figure FDA0002455040120000024
5.一种信号捕获装置,其特征在于,所述装置包括:
压缩模块,用于利用第一观测矩阵,对稀疏向量进行压缩感知,获得第一观测值,所述稀疏向量为稀疏矩阵与接收到的信号的自相关结果,所述稀疏矩阵包括零元素和非零元素,所述稀疏矩阵包括的非零元素为本地码矩阵包括的元素,所述本地码矩阵包括多个码相位不同的本地码,所述稀疏向量包括的每个元素的索引对应所述本地码矩阵包括的一个本地码;
确定模块,用于依据所述压缩模块获得的所述第一观测值,确定所述稀疏向量中的稀疏值的索引,所述稀疏值为非零值;
获得模块,用于获得所述索引对应的目标本地码的部分匹配滤波结果,将所述部分匹配滤波结果进行快速傅里叶变换FFT运算,得到所述信号的频谱的峰值,当所述峰值大于峰值门限时,获得所述目标本地码的码相位以及所述峰值对应的残余多普勒。
6.根据权利要求5所述的装置,其特征在于,所述第一观测矩阵包括的元素为:
Figure FDA0002455040120000022
其中,
Figure FDA0002455040120000023
为所述第一观测矩阵包括的元素,M1为预设压缩比例,HM01为第一哈达玛矩阵,M1=2M01,m∈{1,2,…,M1},l∈{1,2,…,2ML},α=2ML/M1,M为所述信号的分段长度,L为本地码周期。
7.根据权利要求5所述的装置,其特征在于,所述确定模块,具体用于:
对所述第一观测值进行沃尔什Walsh变换,得到第一沃尔什向量;
确定所述第一沃尔什向量中,数值大于第一沃尔什门限的第一元素向量;
根据所述第一元素向量的索引组成的索引集,确定第二观测矩阵;
利用所述第二观测矩阵,对所述稀疏向量进行压缩感知,获得第二观测值;
对所述第二观测值进行Walsh变换,得到第二沃尔什向量;
确定所述第二沃尔什向量中,数值大于第二沃尔什门限的第二元素向量;
根据第二元素向量的索引,确定所述稀疏向量中的稀疏值的索引。
8.根据权利要求7所述的装置,其特征在于,所述第二观测矩阵包括的元素为:
Figure FDA0002455040120000031
其中,M2为预设的峰值搜索长度,HM02为第二哈达玛矩阵,A1(i)为所述索引集,M2=2M02
Figure FDA0002455040120000032
9.一种电子设备,其特征在于,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现权利要求1-4任一项所述的方法步骤。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1-4任一所述的方法步骤。
CN202010303884.3A 2020-04-17 2020-04-17 一种信号捕获方法、装置、电子设备及介质 Active CN111541465B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010303884.3A CN111541465B (zh) 2020-04-17 2020-04-17 一种信号捕获方法、装置、电子设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010303884.3A CN111541465B (zh) 2020-04-17 2020-04-17 一种信号捕获方法、装置、电子设备及介质

Publications (2)

Publication Number Publication Date
CN111541465A true CN111541465A (zh) 2020-08-14
CN111541465B CN111541465B (zh) 2021-12-10

Family

ID=71978763

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010303884.3A Active CN111541465B (zh) 2020-04-17 2020-04-17 一种信号捕获方法、装置、电子设备及介质

Country Status (1)

Country Link
CN (1) CN111541465B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830409A (zh) * 2012-08-30 2012-12-19 西安电子科技大学 一种基于压缩感知的导航信号采集方法
WO2013115481A1 (ko) * 2012-01-30 2013-08-08 한국과학기술원 지엔에스에스 및 대역확산 신호의 신속한 신호 획득을 위한 압축 감지 방법 및 그 장치
CN104199057A (zh) * 2014-09-11 2014-12-10 重庆大学 一种基于压缩感知技术的boc信号无模糊捕获算法
CN104280750A (zh) * 2014-10-23 2015-01-14 西安电子科技大学 基于部分匹配滤波fft算法的长码捕获系统及方法
CN105301608A (zh) * 2015-10-30 2016-02-03 西安烽火电子科技有限责任公司 高动态接收机及其对北斗卫星信号b1频点的捕获方法
CN107315186A (zh) * 2017-06-30 2017-11-03 重庆邮电大学 在残余频偏下基于加窗pmf‑fft的boc信号捕获算法
CN109150235A (zh) * 2017-06-28 2019-01-04 南京理工大学 基于压缩感知的多周期直扩msk信号二维联合捕获方法
CN109765586A (zh) * 2019-01-09 2019-05-17 安徽师范大学 基于压缩感知的glonass信号快速捕获方法
CN110907958A (zh) * 2019-10-23 2020-03-24 深圳华大北斗科技有限公司 信号捕获方法、装置、计算机设备和存储介质

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013115481A1 (ko) * 2012-01-30 2013-08-08 한국과학기술원 지엔에스에스 및 대역확산 신호의 신속한 신호 획득을 위한 압축 감지 방법 및 그 장치
CN102830409A (zh) * 2012-08-30 2012-12-19 西安电子科技大学 一种基于压缩感知的导航信号采集方法
CN104199057A (zh) * 2014-09-11 2014-12-10 重庆大学 一种基于压缩感知技术的boc信号无模糊捕获算法
CN104280750A (zh) * 2014-10-23 2015-01-14 西安电子科技大学 基于部分匹配滤波fft算法的长码捕获系统及方法
CN105301608A (zh) * 2015-10-30 2016-02-03 西安烽火电子科技有限责任公司 高动态接收机及其对北斗卫星信号b1频点的捕获方法
CN109150235A (zh) * 2017-06-28 2019-01-04 南京理工大学 基于压缩感知的多周期直扩msk信号二维联合捕获方法
CN107315186A (zh) * 2017-06-30 2017-11-03 重庆邮电大学 在残余频偏下基于加窗pmf‑fft的boc信号捕获算法
CN109765586A (zh) * 2019-01-09 2019-05-17 安徽师范大学 基于压缩感知的glonass信号快速捕获方法
CN110907958A (zh) * 2019-10-23 2020-03-24 深圳华大北斗科技有限公司 信号捕获方法、装置、计算机设备和存储介质

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
南花妮: "高动态卫星通信系统快速捕获和跟踪算法研究", 《中国优秀硕士学位论文全文数据库》 *
王凯,吴斌,汪勃: "一种利用压缩感知改进的PMF-FFT扩频捕获算法", 《电讯技术》 *
王凯,吴斌,汪勃: "基于压缩感知的扩频信号统一捕获技术", 《无线电工程》 *
程艳合,杨文革: "基于压缩感知的DSSS信号双阶段捕获方法", 《北京航空航天大学学报》 *
郭晓彤: "基于压缩感知的GPS信号捕获算法研究", 《中国优秀硕士学位论文全文数据库》 *

Also Published As

Publication number Publication date
CN111541465B (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
US8750617B2 (en) Signature representation of data having high dimensionality
Bao et al. A compressed sensing approach to blind separation of speech mixture based on a two-layer sparsity model
CN110163061B (zh) 用于提取视频指纹的方法、装置、设备和计算机可读介质
CN110866881A (zh) 图像处理方法及装置、存储介质和电子设备
CN106681716B (zh) 智能终端及其应用程序的自动分类方法
CN111404853A (zh) 一种载波频偏估计方法、装置及计算机存储介质
US20150229499A1 (en) Spectrum sensing method and spectrum sensing device
WO2010080857A2 (en) Scalable media fingerprint extraction
CN111541465B (zh) 一种信号捕获方法、装置、电子设备及介质
CN106411795B (zh) 一种非重构框架下的信号估计方法
US20140198888A1 (en) Discrete signal synchronization based on a known bit pattern
Jabeen et al. Sequency domain signal processing using complex Hadamard transform
CN111931932A (zh) 对抗样本的生成方法及装置、电子设备和可读存储介质
CN115021759B (zh) 一种基于二值最小二乘法的二值稀疏信号恢复方法及系统
Razavi et al. High-resolution cyclic spectrum reconstruction from sub-Nyquist samples
CN111308426B (zh) 一种适用于单天线接收机的低信噪比周期调频信号检测与分离方法
El-Khamy et al. C9. a new wideband spectrum sensing technique based on compressive sensing, wavelet edge detection and parallel processing
US9236056B1 (en) Variable length local sensitivity hash index
CN105811992A (zh) 一种可分离稀疏信号的压缩感知方法及其系统
CN110689133A (zh) 一种训练机器学习引擎的方法、系统及相关装置
CN110336574B (zh) 一种源信号的恢复方法及装置
Khandakar et al. Complexity analysis and accuracy of image recovery based on signal transformation algorithms
Kara Data-Driven Forward–Backward Pursuit for Sparse Signal Reconstruction
Cadavid et al. Simulation and analysis of compressed sensing technique as sampling and data compression and reconstruction of signals using convex programming
CN107908594B (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