CN106680776A - 对多普勒信息不敏感的低旁瓣波形设计方法 - Google Patents
对多普勒信息不敏感的低旁瓣波形设计方法 Download PDFInfo
- Publication number
- CN106680776A CN106680776A CN201611146183.3A CN201611146183A CN106680776A CN 106680776 A CN106680776 A CN 106680776A CN 201611146183 A CN201611146183 A CN 201611146183A CN 106680776 A CN106680776 A CN 106680776A
- Authority
- CN
- China
- Prior art keywords
- overbar
- formula
- rho
- sigma
- epsiv
- 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了对多普勒信息不敏感的低旁瓣波形设计方法,具体按照以下步骤实施:步骤1:构建“时延‑多普勒”范围内低自相关旁瓣波形设计的复数数学模型,对该数学模型的系数矩阵进行分割;步骤2:使用ADMM方法对目标函数进行变换求解引入变量;步骤3:用分段最小值法求解不等式约束问题;步骤4:使用MM方法得到满足约束条件的探测波形x;步骤5:重复步骤3、4直至x收敛。本发明对多普勒信息不敏感的低旁瓣波形设计方法,所设计的波形在一定“时延‑多普勒”区域内的每个自相关旁瓣很小,均小于‑40dB,满足特定距离和速度范围内运动目标的高分辨率探测需求,同时,本发明方法设计的波形具有多样性,因此还具备反侦察特性。
Description
技术领域
本发明属于波形设计技术领域,具体涉及一种对多普勒信息不敏感的低旁瓣波形设计方法。
背景技术
雷达的工作原理是:发射电磁波对探测区域进行照射并接收回波,通过对比发射信号与回波信号的频率、时延等推测计算探测目标的距离、位置、高度、角度、速度等信息。如图1所示的雷达探测工作原理的示意图。因此,探测波形对于雷达探测性能至关重要。自相关函数旁瓣越低的探测信号,雷达探测的精度越高。
目前,对运动目标的检测普遍采用的是线性或非线性调频信号。线性或非线性调频信号的设计简单。线性调频信号的“时延-多普勒”自相关见图2。图2中可见,其时延因子与多普勒因子在特定范围内的旁瓣会出现尖峰;线性或非线性调频信号的起始频率、终止频率与调频函数一旦被对方确认后,发射信号对于对方将无任何秘密而言。
考虑到实际应用中,只需要探测感兴趣区域的目标信息,对此,Hao He及国际雷达信号处理专家Jian Li教授、PetreStoica教授的著作《Waveform Design for ActiveSensing Systems-A computational approach》(Cambridge University Press,July,2011,pp.88-106)对感兴趣的“时延-多普勒”区域提出了低自相关旁瓣设计方法(本发明称其为He-Li-Stoica方法)。He-Li-Stoica方法可以看做是在“时延-多普勒”区域上开了个窗,使得探测波形在该窗口内的整体旁瓣电平(Integrated Sidelobe Level,ISL)取得极小值。仿真结果见图3。该方法在窗口内的整体旁瓣电平较低,但个别“时延-多普勒”上(尤其是窗口边缘)的自相关旁瓣较大。在实际探测工作中,目标可能出现在该窗口内的任意“时延-多普勒”上,因此,较大的个别“时延-多普勒”上的自相关旁瓣影响实际探测的性能。
因此,上述波形设计方法都存在一定的局限性和不可靠性。
发明内容
本发明的目的是提供一种对多普勒信息不敏感的低旁瓣波形设计方法,解决了现有的目标检测波形设计方法存在的个别“时延-多普勒”上的自相关旁瓣较大的问题。
本发明所采用的技术方案是,对多普勒信息不敏感的低旁瓣波形设计方法,具体按照以下步骤实施:
步骤1:构建“时延-多普勒”范围内低自相关旁瓣波形设计的复数数学模型,对该数学模型的系数矩阵进行分割得到矩阵的实系数矩阵与虚系数矩阵并将非正定的实系数矩阵与虚系数矩阵转换为正定矩阵;
步骤2:使用ADMM方法对目标函数进行变换求解引入变量;
步骤3:用分段最小值法求解不等式约束问题;
步骤4:使用MM方法得到满足约束条件的探测波形x;
步骤5:重复步骤3、4直至x收敛。
本发明的有益效果是:本发明对多普勒信息不敏感的低旁瓣波形设计方法,所设计的波形在一定“时延-多普勒”区域内的每个自相关旁瓣很小,均小于-40dB,满足特定距离和速度范围内运动目标的高分辨率探测需求,同时,本发明方法设计的波形具有多样性,因此还具备反侦察特性。
附图说明
图1是雷达探测工作原理示意图;
图2是线性调频信号的自相关三维图;
图3是He-Li-Stoica方法设计的波形的自相关三维图;
图4是本发明方法的自相关旁瓣迭代结果图;
图5是本发明方法设计的波形的自相关三维图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明对多普勒信息不敏感的低旁瓣波形设计方法,首先,构建“时延-多普勒”范围内低自相关旁瓣波形设计的复数数学模型,该复数模型的目标函数是一元一次函数,约束条件中包含有多元四次方不等式,对该数学模型的系数矩阵进行分割得到矩阵的实系数矩阵与虚系数矩阵并将非正定的实系数矩阵与虚系数矩阵转换为正定矩阵;其次,使用ADMM框架将上述模型分解为若干模型,并分别用分段最小值法与MM方法解决不等式约束与多元四次方优化问题,解出自相关旁瓣约束变量ε、替代变量Flk、Glk与信号变量从而得到满足约束条件的探测波形x。
本发明对多普勒信息不敏感的低旁瓣波形设计方法,就是设计一个长度为N的模为1的信号序列x={x1,...,xN},使得信号序列在指定范围内多普勒与时延的自相关r的模满足|r|≤ε。
具体实施过程为:
步骤1:数据与模型预处理
将“时延-多普勒”范围内的自相关函数的实部与虚部进行分离。根据给定的时延范围与多普勒范围,计算范围内的不同时延与多普勒对应的实部系数矩阵与虚部系数矩阵的正交分解向量。此时,自相关由实部系数矩阵与虚部系数矩阵表示。对得到的各正交分解向量进行实部与虚部的分离,构成新的系数矩阵。
仅考虑时延Kk∈{Kd,...,Ku}(k=1,...,K,0<Kd<Ku<N)时,非周期自相关函数
同时考虑时延Kk与多普勒Ll∈{Ld,...,Lu}(l=1,...,L,Ld<Lu)因素时,非周期自相关函数写为:
rlk=xHAlkx (1)
其中,
x为长度为N的模为1的信号序列,x=[x1,...,xN]T,lk=1,...,LK,LK=L×K
算子diag(·)是将向量扩展为对角矩阵,使得Dl(n,n)=dl(n)。
那么,满足要求的波形设计的数学描述为:
对系数矩阵Alk进行分割得到:
其中,
显然,RElk与IMlk满足:
1)
2)xHRElkx与xHIMlkx均为实数。
那么,自相关函数rlk与RElk、IMlk满足关系:
4|rlk|2=(xHRElkx)2+(xHIMlkx)2 (4)
为了后续使用MM方法,需要将非正定的实系数矩阵与虚系数矩阵转换为正定矩阵,即系数矩阵的特征值均大于零。详细过程见式(5)与式(6)。
对系数矩阵RElk进行特征分解,得到特征根{λ1,...,λN},记MRlk=min{λ1,...,λN},其中min{·}表示求·中的最小元素。
若MRlk≤0,则:
MRlk=MRlk-0.1
若MRlk>0,则:
其中EN为N阶单位矩阵,即:
显然,是正定矩阵。正定矩阵均可分解为一对对称矩阵的乘积,因此正定矩阵分解为:
其中,
同理,系数矩阵IMlk对应的正定矩阵分解为:
其中,
系数矩阵IMlk的特征根的最小值,记为MIlk。
那么式(1)可以表示为:
根据信号序列x的模为1即xHx=1可得
式(8)可以写成另外一种实数形式:
其中,Re(·)与Im(·)分别表示·的实部与虚部。
则是实对称矩阵,即:
证明:由式(6.1)有
则Re(Wlk)T=Re(Wlk),Im(Wlk)T=-Im(Wlk)
同理可证
由式(9),式(4)可以表示为:
令替代变量Flk与Glk分别为则式(10)可以表示为:
则式(2)表示为:
步骤2:基于步骤1得到的系数矩阵,使用ADMM方法对目标函数进行变换求解引入变量,即使用ADMM方法来求解规划式(12)。
引入拉格朗日乘子αlk,βlk与ADMM乘子ρ>0,式(12)可表示为
步骤3:使用分段最小值法求解不等式约束问题,从而更新ε(t+1),Flk(t+1),Glk(t+1),αlk(t+1),βlk(t+1)
①更新Flk(t+1),Glk(t+1)
约束与变量无关,它们仅在迭代求解Flk,Glk的时候起作用,而约束与变量Flk与Glk无关,它们仅在迭代求解的时候起作用。所以对于给定t次的求解ε(t+1),Flk(t+1),Glk(t+1)的问题,式(13)可表示为:
记
那么式(14)可以简化为:
由式(14)的约束条件可知,变量Flk,Glk的取值是受变量ε约束的。而对于给定的ε,式(14)的Flk(t+1),Glk(t+1)为:
其中,
显然,记为Slk,则:
②更新ε(t+1)
给定t次的求解ε(t+1)。将式(17)、(18)代入式(16)得:
忽略掉常数项,式(19)简化为:
其中,
显然,因为式(18)的分段情况,一般情况下式(20)是不连续的,需要分段讨论其极小值点。分段讨论式(20)时,分段的分割点取值范围为其中,εmax是代入原问题得到的ε,ε0是大于0且小于εmax与的实数。
在区间[εd,εu]上,ρ>0且Slk≥0,所以A,B>0,故式(20)取极小值时的ε为:
将式(21)代入式(20)得到L在各段的极小值{L1,...,Lm},选取{L1,...,Lm}中最小值对应的ε为ε(t+1)。
③更新拉格朗日乘子αlk(t+1),βlk(t+1)
步骤4:推导目标函数及目标函数的最优函数,从而将四次型目标函数降幂为一次型函数,然后在一次型最优函数的基础上,使用MM方法求解满足约束条件的信号序列。即已知αlk(t+1),ε(t+1),βlk(t+1),Flk(t+1),Glk(t+1),求解
步骤4.1:推导目标函数
同理于式(14),求解时,式(13)可写为:
将式(23)中的常数项去掉后,得到关于的目标函数:
步骤4.2:求解目标函数(24)的优化函数
显然,目标函数(24)为的四次方函数,为了求解找到式(24)的优化函数,使得该优化函数是的一次函数,再使用MM方法求解
式(24)右端第一项有:
其中,因为是实对称矩阵,故而是实对称矩阵,矩阵
也是对称矩阵。
因为,二次型函数f(x)=xHAx,若A,B是n×n的自共轭矩阵且A≤B,那么
根据式(27)
其中,是的最大特征根,
因为所以式(28)的第二项是常数将式(28)的常数项用CW1表示,即:
则式(28)可写为:
又
其中,
与均为实对称矩阵,故M1是实对称矩阵。
那么根据式(27),若是M1的最大特征根,则:
令
则
同理,式(24)右端第二项
其中, 是M1的最大特征根, 是的最大特征根。
式(24)右端第三项的优化函数的推导过程类似于式(31),优化函数为:
其中,
所以,式(24)的优化函数为:
步骤4.3:更新
式(36)的拉格朗日函数为:
对拉格朗日函数(37)关于求偏导
因为是关于的一次函数,所以其对的偏导是常数,记为dn,令则:
将式(39)代入约束条件得
将式(40)代入式(36)目标函数可知,μn取值为因此
将式(41)计算得到的代入
中,得到本轮迭代的x。
步骤5:重复步骤3-4,直至x收敛。
实施例
为了证明本发明方法的有效性,把He-Li-Stoica方法产生的波形、线性调频信号波形与本发明方法产生的波形做比较,下面结合具体实例以及图表进行详细说明。
我们仿真设计长度为N=128的波形,要求其在“时延-多普勒”区域[Lmin,Lmax]×[Kmin,Kmax]为[-2,2]×[-8,8]内的探测波形的自相关旁瓣低于-40dB,ρ=0.1。
步骤1:构建模型基本框架并初始化数据。根据[Lmin,Lmax]×[Kmin,Kmax]的值设定LK=(Lmax-Lmin)×(Kmax-Kmin),那么非周期的“时延-多普勒”的低自相关旁瓣的探测波形x={x1,...,xN}的设计问题定义为:
其中x=[x1 ... xN]T,Alk是已知的矩阵,其具体形式为
随机初始化x(0),得
根据式(2)得到初始自相关函数的模|rlk(0)|以及初始ε(0)=max{|rlk(0)|,lk=1,...,LK}。其中,MRlk,MIlk,的计算步骤如下:
对系数矩阵Alk进行分割得到为了后续使用MM方法,需要将非正定的实系数矩阵与虚系数矩阵转换为正定矩阵,即系数矩阵的特征值均大于零。对矩阵RElk进行特征分解,得到特征根{λ1,...,λN},记MRlk=min{λ1,...,λN},其中min{·}表示求·中的最小元素。若MRlk≤0,则MRlk=MRlk-0.1,若MRlk>0,则MRlk=0,其中N阶单位矩阵 分解为其中同理,系数矩阵IMlk对应的正定矩阵分解为其中系数矩阵IMlk的特征根的最小值,记为MIlk。
步骤2:构建低自相关旁瓣的“时延-多普勒”波形设计问题的ADMM框架,并初始化ADMM相应的参数。
引入拉格朗日乘子αlk,βlk、替代变量 与ADMM乘子ρ>0,波形设计问题(1)可表示为
其中,
初始化Flk(0),Glk(0):
随机初始化αlk(0),βlk(0)。
步骤3、使用分段最小值法求解不等式约束问题,进而更新ε(t+1),Flk(t+1),Glk(t+1),αlk(t+1),βlk(t+1),lk=1,...,LK
对于给定的ε(t),更新替代变量Flk(t+1),Glk(t+1)为
其中
给定t次迭代的求解ε(t+1),具体步骤如下。
在ε的取值区间[εd,εu]上
其中
因为ρ>0且Slk≥0,所以A,B>0,故取极小值时的ε为
根据上述步骤,计算Lε在ε的所有可能取值区间[εd,εu]上的极小值{L1,...,Lm},选取{L1,...,Lm}中最小值对应的ε为ε(t+1)。
更新拉格朗日乘子αlk(t+1),βlk(t+1),lk=1,...,LK为
步骤4、使用MM方法求解四次方优化问题。已知ε(t+1),αlk(t+1),βlk(t+1),Flk(t+1),Glk(t+1),lk=1,...,LK,更新为
其中,dn是拉格朗日函数
对的偏导。
上述拉格朗日函数中,
和分别是M1和M2的最大特征根,和是和的最大特征根。
步骤5、重复步骤3至步骤4,直至x收敛。
图4是本发明方法的自相关旁瓣迭代结果图。设计得到的波形的“时延-多普勒”[Lmin,Lmax]×[Kmin,Kmax]为[-2,2]×[-8,8]的自相关旁瓣三维图见图5。
Claims (5)
1.对多普勒信息不敏感的低旁瓣波形设计方法,其特征在于,具体按照以下步骤实施:
步骤1:构建“时延-多普勒”范围内低自相关旁瓣波形设计的复数数学模型,对该数学模型的系数矩阵进行分割得到矩阵的实系数矩阵与虚系数矩阵并将非正定的实系数矩阵与虚系数矩阵转换为正定矩阵;
步骤2:使用ADMM方法对目标函数进行变换求解引入变量;
步骤3:用分段最小值法求解不等式约束问题;
步骤4:使用MM方法得到满足约束条件的探测波形x;
步骤5:重复步骤3、4直至x收敛。
2.根据权利要求1所述的对多普勒信息不敏感的低旁瓣波形设计方法,其特征在于,所述步骤1具体为:
同时考虑时延Kk与多普勒Ll∈{Ld,...,Lu}因素时,自相关函数写为:
rlk=xHAlkx (1)
其中,k=1,...,K,0<Kd<Ku<N,l=1,...,L,Ld<Lu,lk=1,...,LK,LK=L×K,
x为长度为N、模为1的信号序列,
算子diag(·)是将向量扩展为对角矩阵,使得Dl(n,n)=dl(n);
那么,满足要求的波形设计的数学描述为:
对系数矩阵Alk进行分割得到:
其中,
1)
2)xHRElkx与xHIMlkx均为实数;
那么,自相关函数rlk与RElk、IMlk满足关系:
4|rlk|2=(xHRElkx)2+(xHIMlkx)2 (4)
对系数矩阵RElk进行特征分解,得到特征根{λ1,...,λN},记MRlk=min{λ1,...,λN},其中min{·}表示求·中的最小元素;
若MRlk≤0,则:
若MRlk>0,则:
其中EN为N阶单位矩阵,即:
正定矩阵分解为:
其中,
同理,系数矩阵IMlk对应的正定矩阵分解为:
其中,
系数矩阵IMlk的特征根的最小值,记为MIlk;
那么式(1)表示为:
根据信号序列x的模为1即xHx=1得:
式(8)可以写成另外一种实数形式:
其中,
Re(·)与Im(·)分别表示·的实部与虚部;
则是实对称矩阵,即:
由式(9),式(4)表示为:
令替代变量Flk与Glk分别为则式(10)表示为:
则式(2)表示为:
3.根据权利要求2所述的对多普勒信息不敏感的低旁瓣波形设计方法,其特征在于,所述步骤2具体为:
引入拉格朗日乘子αlk,βlk与ADMM乘子ρ>0,式(12)表示为:
4.根据权利要求3所述的对多普勒信息不敏感的低旁瓣波形设计方法,其特征在于,所述步骤3使用分段最小值法求解不等式约束问题,从而更新ε(t+1),Flk(t+1),Glk(t+1),αlk(t+1),βlk(t+1),具体为:
①更新Flk(t+1),Glk(t+1)
对于给定t次的αlk(t),βlk(t)求解ε(t+1),Flk(t+1),Glk(t+1)的问题,式(13)表示为:
记
式(14)简化为:
由式(14)的约束条件可知,变量Flk,Glk的取值是受变量ε约束的,而对于给定的ε,式(14)的Flk(t+1),Glk(t+1)为:
其中,
记为Slk,则:
②更新ε(t+1)
给定t次的αlk(t),βlk(t)求解ε(t+1):
将式(17)、(18)代入式(16)得:
忽略掉常数项,式(19)简化为:
其中,
由于式(18)的分段情况,式(20)是不连续的,需要分段讨论其极小值点,分段讨论式(20)时,分段的分割点取值范围为其中,εmax是代入原问题得到的ε,ε0是大于0且小于εmax与的实数;
在区间[εd,εu]上,ρ>0且Slk≥0,所以A,B>0,故式(20)取极小值时的ε为:
将式(21)代入式(20)得到L在各段的极小值{L1,...,Lm},选取{L1,...,Lm}中最小值对应的ε为ε(t+1);
③更新拉格朗日乘子αlk(t+1),βlk(t+1)
5.根据权利要求4所述的对多普勒信息不敏感的低旁瓣波形设计方法,其特征在于,所述步骤4具体为:
步骤4.1:推导目标函数
同理于式(14),求解时,式(13)写为:
将式(23)中的常数项去掉后,得到关于的目标函数:
步骤4.2:求解目标函数(24)的优化函数
式(24)右端第一项有:
其中,
因为,二次型函数f(x)=xHAx,若A,B是n×n的自共轭矩阵且A≤B,那么
根据式(27)
其中,是的最大特征根,
因为所以式(28)的第二项是常数将式(28)的常数项用CW1表示,即:
则式(28)写为:
又
其中,
根据式(27),若是M1的最大特征根,则:
令
则
同理,式(24)右端第二项
其中, 是M1的最大特征根, 是的最大特征根;
式(24)右端第三项的优化函数的推导过程同理于式(31),优化函数为:
其中,所以,式(24)的优化函数为:
步骤4.3:更新
式(36)的拉格朗日函数为:
对拉格朗日函数(37)关于求偏导
因为是关于的一次函数,所以式(37)对的偏导是常数,记为dn,令则:
将式(39)代入约束条件得
将式(40)代入式(36)目标函数,μn取值为因此
将式(41)计算得到的代入
中,得到本轮迭代的x。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611146183.3A CN106680776B (zh) | 2016-12-13 | 2016-12-13 | 对多普勒信息不敏感的低旁瓣波形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611146183.3A CN106680776B (zh) | 2016-12-13 | 2016-12-13 | 对多普勒信息不敏感的低旁瓣波形设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106680776A true CN106680776A (zh) | 2017-05-17 |
CN106680776B CN106680776B (zh) | 2019-06-21 |
Family
ID=58869441
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611146183.3A Expired - Fee Related CN106680776B (zh) | 2016-12-13 | 2016-12-13 | 对多普勒信息不敏感的低旁瓣波形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106680776B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107677999A (zh) * | 2017-09-25 | 2018-02-09 | 西北工业大学 | 一种精确控制相关性旁瓣的序列集设计方法 |
CN107703489A (zh) * | 2017-09-26 | 2018-02-16 | 西北工业大学 | 一种mimo雷达恒模波形和接收机的联合设计方法 |
CN107728118A (zh) * | 2017-09-25 | 2018-02-23 | 西北工业大学 | 无需拟合协方差矩阵的低旁瓣发射波束图设计方法 |
CN109343006A (zh) * | 2018-10-31 | 2019-02-15 | 中国科学院电子学研究所 | 基于增广拉格朗日遗传算法的nflm信号优化方法及装置 |
CN109446665A (zh) * | 2018-10-31 | 2019-03-08 | 中国科学院电子学研究所 | 一种非线性调频信号优化方法及装置、存储介质 |
CN112630732A (zh) * | 2020-12-14 | 2021-04-09 | 中国电子科技集团公司第五十四研究所 | 一种基于isl约束的抗射频干扰设计方法 |
CN113311413A (zh) * | 2021-03-25 | 2021-08-27 | 西北工业大学 | 可控制模糊度函数的声呐波形设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105959015A (zh) * | 2016-04-22 | 2016-09-21 | 西安电子科技大学 | 基于最小多面体模型的ldpc码线性规划译码方法 |
US20160335224A1 (en) * | 2014-03-31 | 2016-11-17 | Los Alamos National Security, Llc | Efficient convolutional sparse coding |
-
2016
- 2016-12-13 CN CN201611146183.3A patent/CN106680776B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160335224A1 (en) * | 2014-03-31 | 2016-11-17 | Los Alamos National Security, Llc | Efficient convolutional sparse coding |
CN105959015A (zh) * | 2016-04-22 | 2016-09-21 | 西安电子科技大学 | 基于最小多面体模型的ldpc码线性规划译码方法 |
Non-Patent Citations (2)
Title |
---|
JUNLI LIANG ET AL.: "Unimodular Sequence Design Based on Alternating Direction Method of Multipliers", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
王永丽 等: "一种求解单调变分不等式的下降型邻近点交替方向乘子法", 《山东科技大学学报 自然科学版》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107677999A (zh) * | 2017-09-25 | 2018-02-09 | 西北工业大学 | 一种精确控制相关性旁瓣的序列集设计方法 |
CN107728118A (zh) * | 2017-09-25 | 2018-02-23 | 西北工业大学 | 无需拟合协方差矩阵的低旁瓣发射波束图设计方法 |
CN107677999B (zh) * | 2017-09-25 | 2020-09-08 | 西北工业大学 | 一种精确控制相关性旁瓣的序列集设计方法 |
CN107728118B (zh) * | 2017-09-25 | 2020-11-06 | 西北工业大学 | 无需拟合协方差矩阵的低旁瓣发射波束图设计方法 |
CN107703489A (zh) * | 2017-09-26 | 2018-02-16 | 西北工业大学 | 一种mimo雷达恒模波形和接收机的联合设计方法 |
CN107703489B (zh) * | 2017-09-26 | 2021-05-18 | 西北工业大学 | 一种mimo雷达恒模波形和接收机的联合设计方法 |
CN109343006A (zh) * | 2018-10-31 | 2019-02-15 | 中国科学院电子学研究所 | 基于增广拉格朗日遗传算法的nflm信号优化方法及装置 |
CN109446665A (zh) * | 2018-10-31 | 2019-03-08 | 中国科学院电子学研究所 | 一种非线性调频信号优化方法及装置、存储介质 |
CN112630732A (zh) * | 2020-12-14 | 2021-04-09 | 中国电子科技集团公司第五十四研究所 | 一种基于isl约束的抗射频干扰设计方法 |
CN112630732B (zh) * | 2020-12-14 | 2022-07-19 | 中国电子科技集团公司第五十四研究所 | 一种基于isl约束的抗射频干扰设计方法 |
CN113311413A (zh) * | 2021-03-25 | 2021-08-27 | 西北工业大学 | 可控制模糊度函数的声呐波形设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106680776B (zh) | 2019-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106680776A (zh) | 对多普勒信息不敏感的低旁瓣波形设计方法 | |
Cuomo et al. | Ultrawide-band coherent processing | |
Bocquet | Parameter estimation for Pareto and K distributed clutter with noise | |
WO2018045566A1 (zh) | 基于压缩感知的随机脉冲多普勒雷达角度-多普勒成像方法 | |
CN111796272B (zh) | 穿墙雷达人体图像序列的姿态实时识别方法及计算机设备 | |
US10275707B2 (en) | Systems and methods for training multipath filtering systems | |
CN108008389B (zh) | 一种基于gpu的快速频域后向投影三维成像方法 | |
CN110208760B (zh) | 一种基于时域上采样的雷达回波仿真方法 | |
CN104536000A (zh) | 一种实波束扫描雷达角超分辨方法 | |
Pan et al. | Multi-task hidden Markov modeling of spectrogram feature from radar high-resolution range profiles | |
CN108896990B (zh) | 一种利用耦合模式字典学习的建筑物墙体成像方法和装置 | |
CN107843875A (zh) | 基于奇异值分解降噪的贝叶斯压缩感知雷达数据融合方法 | |
CN106802409B (zh) | 基于多cpu处理的外辐射源雷达实时信号处理方法 | |
CN104950297A (zh) | 基于矩阵1范数拟合的阵元误差估计方法 | |
CN104833970B (zh) | 一种适用于存在遮挡效应的旋转目标稀疏多普勒成像方法 | |
Ding et al. | Super‐resolution 3D imaging in MIMO radar using spectrum estimation theory | |
Chang et al. | Applied state estimation and association | |
CN115453528A (zh) | 基于快速sbl算法实现分段观测isar高分辨成像方法及装置 | |
CN111880159A (zh) | 一种基于lstm的雷达序列信号检测方法及系统 | |
Feng et al. | Outlier-robust tri-percentile parameter estimation of compound-Gaussian clutter with lognormal distributed texture | |
CN112904298B (zh) | 一种基于局部网格分裂的网格偏离空时自适应处理方法 | |
Michałowska et al. | Neural operator learning for long-time integration in dynamical systems with recurrent neural networks | |
Liu et al. | Maneuvering extended object tracking based on constrained expectation maximization | |
CN105022025B (zh) | 基于稀疏处理的信号波达方向估计方法 | |
CN105572642B (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 | ||
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: 20190621 Termination date: 20191213 |