CN113408122B - 一种多自由度调频信号的生成方法、装置及存储介质 - Google Patents
一种多自由度调频信号的生成方法、装置及存储介质 Download PDFInfo
- Publication number
- CN113408122B CN113408122B CN202110662568.XA CN202110662568A CN113408122B CN 113408122 B CN113408122 B CN 113408122B CN 202110662568 A CN202110662568 A CN 202110662568A CN 113408122 B CN113408122 B CN 113408122B
- Authority
- CN
- China
- Prior art keywords
- frequency modulation
- degree
- function
- freedom
- modulation signal
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/086—Learning methods using evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Computational Linguistics (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请公开了一种多自由度调频信号的生成方法,所述方法包括:基于分段非线性函数,构建多自由度调频信号的时域函数;根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的冲激响应宽度IRW和峰值旁瓣比PSLR,基于所述IRW和所述PSLR建立信号优化模型;基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型;根据选取的窗函数和驻定相位原理POSP,确定多自由度调频信号的初始化调频率控制点;利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
Description
技术领域
本申请涉及雷达发射信号领域,涉及但不限于一种多自由度调频信号生成方法、装置及存储介质。
背景技术
随着高分辨率雷达卫星的发展,大量的星载合成孔径雷达(Synthetic ApertureRadar,SAR)被应用于遥感领域,如TerraSAR-X、Sentinel-1和COSMO-SkyMed等。传统SAR系统使用线性调频(Linear Frequency Modulation,LFM)信号作为发射波形时会产生高的旁瓣电平,而高的旁瓣电平会干扰甚至覆盖目标附近的散射体。为了抑制目标产生的高旁瓣,通常对LFM信号进行加窗处理。然而加窗的操作会使滤波过程不再与信号精确匹配,从而降低雷达系统信噪比(Signal-Noise Ratio,SNR)。为了克服旁瓣抑制造成的SNR损失的难题,多自由度调频信号因其能同时保持雷达的信噪比又能达到通过加窗处理抑制旁瓣的效果,得到了广大雷达发射信号处理工作者的青睐。
然而现有技术中,在使用多自由度调频信号的过程中,会在降低旁瓣水平的同时不可避免地造成了主瓣的展宽,主瓣的展宽会导致多自由度调频信号的分辨率下降,进而降低SAR图像的成像质量。
多自由度调频信号因其非线性的时频关系使其具有更大的设计自由度,因此,如何设计和优化多自由度调频信号,使其能够最大限度降低旁瓣的同时保持较窄的主瓣是当前亟需解决的问题。
发明内容
本申请实施例提供一种多自由度调频信号的生成方法、装置及存储介质,生成的多自由度调频信号能够在最大限度降低旁瓣的同时最低限度展宽主瓣。
本申请的技术方案是这样实现的:
第一方面,本申请实施例提供一种多自由度调频信号的生成方法,所述方法包括:
基于分段非线性函数,构建多自由度调频信号的时域函数;
根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的冲激响应宽度(Impulse Response Width,IRW)和峰值旁瓣比(Peak Side LobeRatio,PSLR),基于所述IRW和所述PSLR建立信号优化模型;
基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型;
根据选取的窗函数和驻定相位原理(Principle of Stationary Phase,POSP),确定多自由度调频信号的初始化调频率控制点;
利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
第二方面,本申请实施例提供一种多自由度调频信号的生成装置,所述装置包括:函数建立单元、确定单元、模型建立单元、初始化单元和迭代单元;
所述函数建立单元,用于基于分段非线性函数,构建多自由度调频信号的时域函数;
所述确定单元,用于根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的IRW和PSLR,基于所述IRW和所述PSLR建立信号优化模型;
所述模型建立单元,用于基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型;
所述初始化单元,用于根据选取的窗函数和POSP,确定多自由度调频信号的初始化调频率控制点;
所述迭代单元,用于利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
第三方面,本申请实施例提供一种存储介质,所述存储介质上存储有计算机可执行指令,所述计算机可执行指令被处理器执行时,实现上述多自由度调频信号的生成方法。
本申请实施例提供了一种基于分段线性函数,建立分段非线性函数的多自由度调频信号的时域函数;根据所述分段非线性函数的多自由度调频信号的时域函数,确定所述分段非线性函数的多自由度调频信号的自相关函数性能中的IRW和PSLR,基于所述IRW和所述PSLR建立信号优化模型;基于所述信号优化模型的约束条件和目标函数,确定增广拉格朗日遗传算法的算法模型;根据选取的窗函数和POSP,确定多自由度调频信号的初始化调频率控制点;利用所述分段非线性函数的多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行优化,直到所述增广拉格朗日遗传算法收敛;如此,得到的多自由度调频信号能够在最大限度降低旁瓣的同时,最低限度展宽主瓣。
附图说明
图1为本申请实施例提供的一种多自由度调频信号的生成方法的流程图;
图2为本申请实施例提供的一种多自由度调频信号的调频率函数;
图3为本申请实施例提供的一种分段非线性的多自由度调频信号的时频函数;
图4为本申请实施例提供的一种多自由度调频信号的时频函数;
图5为本申请实施例提供的一种多自由度调频信号的时域波形;
图6为本申请实施例提供的一种多自由度调频信号的自相关输出波形;
图7为本申请实施例提供的一种多自由度调频信号的生成装置的可选地结构示意图;
图8为本申请实施例提供的一种多自由度调频信号的生成装置的可选地结构示意图。
具体实施方式
在本申请实施例中,基于分段非线性函数,构建多自由度调频信号的时域函数;根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的IRW和PSLR,基于所述IRW和所述PSLR建立信号优化模型;基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型;根据选取的窗函数和POSP,确定多自由度调频信号的初始化调频率控制点;利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行优化,直到所述增广拉格朗日遗传算法收敛。
下面结合附图对技术方案的实施作进一步的详细描述。
本申请实施例提供一种多自由度调频信号生成方法,如图1所示,所述方法包括S101至S105:
S101、基于分段非线性函数,构建多自由度调频信号的时域函数。
在一些实施例中,在上述S101之前,所述方法还包括:
基于分段线性函数的多自由度调频信号的瞬时调频率函数,确定分段非线性函数的多自由度调频信号的时频函数。
这里,时频函数为瞬时频率函数,自变量为t,变量为f。
在一些实施例中,上述S101包括:根据所述分段非线性函数的多自由度调频信号的时频函数,确定所述分段非线性函数的多自由度调频信号的相位函数;根据所述相位函数和所述多自由度调频信号对应的幅度,确定所述多自由度调频信号的时域函数。
在实际应用中,通常将所述多自由度调频信号对应的幅度A进行归一化处理,即将A取为1。
在一些实施例中,上述基于分段线性函数的多自由度调频信号的瞬时调频率函数,确定所述分段非线性函数的多自由度调频信号的时频函数之前,所述方法还包括:在调频率关系平面内确定时间分段点t和调频率控制点r;基于所述时间分段点t和所述调频率控制点r,确定所述分段线性函数的多自由度调频信号的瞬时调频率函数。
这里,在笛卡尔坐标系中定义调频率关系(t,r)平面,其中,t为横坐标,用于表示分段线性函数的多自由度调频信号的时间,r为纵坐标,用于表示分段线性函数的多自由度调频信号的瞬时调频率。
在调频率关系平面内,预先设置多自由度调频信号的脉冲宽度为T,若将脉冲宽度为T的信号均匀分成2n+2段,则横坐标上共有2n+3个均匀分布的时间分段点各时间分段点在横坐标上是均匀分布的,其中,T1i为在调频率关系平面内构成的一段线性函数的横坐标,T2i为在在调频率关系平面内构成的另一段线性函数的横坐标;纵坐标对应横坐标共有4n+4个调频率控制点r=[r2n2,r2n1,...,r202,r201,r101,r102,...,r1n1,r1n2],其中,n为正整数。
在确定了时间分段点t和调频率控制点r后,如图2所示,每两个调频率控制点r为一组,构成2n+2段不连续的分段线性函数,其中所述分段线性函数的多自由度调频信号对应的瞬时调频率函数r(t)为公式(1):
在一些实施例中,所述基于所述分段线性函数的多自由度调频信号的瞬时调频率函数,确定所述分段非线性函数的多自由度调频信号的时频函数,包括:
对所述分段线性函数的多自由度调频信号的瞬时调频率函数进行积分,确定所述分段非线性函数的多自由度调频信号的时频函数。
这里,瞬时调频率函数用于表示某一时刻的瞬时频率的变化率。时频函数是瞬时频率函数,自变量为t,变量为f。
其中,如图3所示,为确定的所述分段非线性函数的多自由度调频信号的时频函数。
所述分段非线性函数的多自由度调频信号的时频函数f(t)为公式(2):
其中,上述系数为公式(3)至公式(5):
其中,在公式(3)中,r102、r101、r1i2、r1i1、r1n2、r1n1、r202、r201、r2i2、r2i1、r2n2和r2n1分别为在调频率关系平面内构成的线性函数的纵坐标,T11、T1(i+1)、T1i、T1n、T21、T2(i+1)、T2i、和T2n分别为在调频率关系平面内构成的线性函数的横坐标。
在一些实施例中,上述根据所述分段非线性函数的多自由度调频信号的时频函数,确定所述分段非线性函数的多自由度调频信号的相位函数,包括:
对所述分段非线性的多自由度调频信号的时频函数进行积分,确定所述分段非线性的多自由度调频信号的相位函数θ(t),其中,θ(t)可通过公式(6)表示:
θ(t)=2π·∫f(t)dt 公式(6)。
S102、根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的冲激响应宽度IRW和PSLR,基于所述IRW和所述PSLR建立信号优化模型。
这里,冲激响应宽度IRW为冲激响应的3dB主瓣宽度。
多自由度调频信号的时域函数为分段非线性函数的多自由度调频信号的时域函数。
根据多自由度调频信号的时域函数,确定多自由度调频信号的自相关函数为:将所述多自由度调频信号的幅度转变为dB形式,计算多自由度调频信号的IRW和PSLR。
其中,时域函数s(t)为公式(7):
s(t)=Aexp{-jθ(t)} 公式(7);
自相关函数R(t)为公式(8):
其中,在公式(8)中,τ为参数;
IRW为主瓣3dB处对应的宽度,PSLR是除主瓣之外的旁瓣最大值,IRW和PSLR的定义如下:
IRW:3dB主瓣宽度的大小,一般归一化为采样点;
PSLR:最高旁瓣与主瓣峰值高度的比值,单位为dB;
在一些实施例中,基于所述IRW和PSLR建立信号优化模型,包括:
所述信号优化模型为:
min PSLR(r),such that IRW(r)-β≤0,
其中,min PSLR(r),such that IRW(r)-β≤0用于表征在满足不展宽主瓣的情况下降低旁瓣,PSLR(r)为根据调频控制点求得的多自由度调频信号的PSLR;IRW(r)为根据调频控制点求得的多自由度调频信号的IRW;β为主瓣宽度控制因子,该控制因子用于确保主瓣宽度在优化过程中不超过该值;r1i1和r1i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;r2i1和r2i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;为脉冲宽度T平均分成2n+2段的时间间隔;B为预设的带宽;若和则可控制多自由度调频信号的频率在的范围内。
本申请实施例中,基于IRW和PSLR建立的信号优化模型,可以在不展宽主瓣的情况下尽可能地降低旁瓣。
S103、基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型。
这里,所述基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型包括:
基于所述信号优化模型,利用增广拉格朗日算法确定目标函数Ψ(r,λ,s),其中,目标函数Ψ(r,λ,s)为公式(9):
Ψ(r,λ,s)=PSLR(r)-λlog(s-(IRW(r)-β)): 公式(9);
其中,r为调频率控制点,λ为拉格朗日算子,s为偏移量,λ和s均为非负数,约束条件为IRW(r)-β,β为常数。
利用增广拉格朗日算法确定参数函数,其中,参数函数λk+1为公式(10):
其中,μ为保证s-(IRW(r)-β)大于0的数值;k为非负数。
所述目标函数和所述参数函数组成所述增广拉格朗日遗传算法的算法模型。
根据所述目标函数计算当前多自由度调频信号集合中各个多自由度调频信号的适应度;根据各个多自由度调频信号的适应度数值的大小,从所述多自由度调频信号集合中选择父类多自由度调频信号;对所述父类多自由度调频信号进行交叉处理、变异处理,得到下一次迭代的父类多自由度调频信号。
根据所述参数函数计算当前拉格朗日算子和偏移量对应的下一次迭代的拉格朗日算子和偏移量。
增广拉格朗日算法把具体的求解问题分成两部分:一部分为传统的遗传算法,一部分为增广拉格朗日算法。在传统的遗传算法中,适应度函数为PSLR(r),约束条件为IRW(r)-β。本申请实施例中,由于了引入拉格朗日算子λ,因此,可以将约束条件,即IRW(r)-β转换为惩罚项添加到目标函数中,即PSLR(r)-λlog(s-(IRW(r)-β)),从而可以将约束问题转换为一个无约束问题。
增广拉格朗日算法用于解决约束问题,可以根据上述参数函数,不断更新λ和s。
遗传算法是模拟“适者生存”的自然选择的动态优化过程。首先,基于增广拉格朗日遗传算法模型,计算初始多自由度调频信号集合中的适应度,根据各个多自由度调频信号集合中的适应度从中选择父类,对所述父类进行交叉处理和变异处理,从而可以得到下一次迭代的多自由度调频信号;再根据所述参数函数计算当前拉格朗日算子和偏移量对应的下一次迭代的拉格朗日算子和偏移量,将适应度好的,即较优的多自由度调频信号选择出来作为父类,经过交叉处理和变异处理提高子代的多样性,并再次评估适应度的优劣。在整个迭代更新优化过程中,适应度高的个体会被保留下来,适应度差的个体则将被淘汰,这样,经过多次的循环可以得到最优的多自由度调频信号。
S104、根据选取的窗函数和POSP,确定多自由度调频信号的初始化调频率控制点。
在一些实施例中,上述S104包括:根据所述窗函数和驻定相位原理POSP,计算多自由度调频信号的群时延函数;根据所述群时延函数,确定所述多自由度调频信号对应的初始时频函数;根据所述多自由度调频信号对应的初始时频函数,确定所述多自由度调频信号的初始调频率函数;根据所述多自由度调频信号的初始调频率函数,确定所述多自由度调频信号的初始调频率控制点。
这里,根据选取的窗函数和POSP,计算多自由度调频信号对应的群时延函数,其中,群时延函数T(f)为公式(11):
本申请实施例中,驻定相位原理POSP指出,在预设一个窗函数,例如,泰勒窗后,设计一个多自由度调频信号的功率谱,使该功率谱接近预设的窗函数,从而可以生成低旁瓣的多自由度调频信号。这样,可以通过窗函数的性能特性来很好地预测生成的多自由度调频信号的性能特性,从而可以根据不同的需求选择相应的窗函数。其中,选取窗函数为sll=-35的泰勒窗P(f)作为期望的功率谱,指的是希望设计的多自由度调频信号的功率谱接近泰勒窗。
群时延函数与初始时频函数是一对反函数,在确定群时延函数后,可以根据f(t)=T-1(f),确定初始时频函数。
在一示例中,可以根据f(t)=T-1(f),确定一曲线形状近似为“S”型的初始时频函数。
本申请实施例中,可以对所述多自由度调频信号对应的初始时频函数进行求导,从而确定所述多自由度调频信号的初始调频率函数。
在一示例中,可以通过对曲线形状近似为“S”型的初始时频函数进行求导,以得到曲线形状近似为“U”型的初始调频率函数。
在构建的调频率关系平面内,由于横坐标上有2n+3个时间分段点t,纵坐标上有4n+4个调频率控制点r,因此,在得到多自由度调频信号的初始调频率函数时,在点处对应r2n2,在t=T2n处对应r2n1和r2(n-1)2,…,在t=0处对应r201和r101,t=T11处对应r102和r111,…,在t=T1n处对应r1(n-1)2和r1n1,在处对应r1n2。在后续优化过程中,在t=T2n处对应的r2n1和r2(n-1)2会逐渐分开,…,在t=T1n处对应r1(n-1)2和r1n1也会逐渐分开,其中,逐渐分开代表着调频率控制点在优化的过程中性能会变好。根据所述多自由度调频信号的初始调频率函数,获得所述多自由度调频信号的初始调频率控制点。
S105、利用所述多自由度调频信号时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
这里,当在S104中确定多自由度调频信号的初始化调频率控制点后,通过S101确定的分段非线性函数的多自由度调频信号的时域函数,以及通过S103确定的增广拉格朗日遗传算法的算法模型对所述多自由度调频信号的初始化调频率控制点进行迭代,直到所述增广拉格朗日遗传算法收敛。
在一些实施例中,上述S105包括:利用所述分段非线性函数的多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,确定第一多自由度调频信号的适应度和所述第一多自由度调频信号的迭代参数,基于所述第一多自由度调频信号的适应度和所述第一多自由度调频信号的迭代参数,对各所述第一多自由度调频信号进行交叉处理和变异处理,得到第二多自由度调频信号;其中,在第一次迭代中,所述第一多自由度调频信号为初始多自由度调频信号;根据所述增广拉格朗日遗传算法,判断所述第二多自由度调频信号的适应度是否满足预设值,若不满足所述预设值,则继续对所述第二多自由度调频信号进行所述交叉处理和所述变异处理,直到进行所述交叉处理和所述变异处理的多自由度调频信号的适应度满足所述预设值。
需要说明的是,本申请实施例不对预设值进行限定,预设值可以为1*10-6,也可以为其他数值。
这里,在对多自由度调频信号进行迭代之前,可先根据参数函数,确定迭代参数,该迭代参数包括拉格朗日算子λ和偏移量s。初始化的迭代参数可称为初始迭代参数,用于迭代过程中的第一次迭代。
在第一次迭代时,利用分段非线性函数的多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,计算各多自由度调频信号的适应度,根据各多自由度调频信号的适应度从初始多自由度调频信号中选择多个父类多自由度调频信号,对选取的多个父类多自由度调频信号进行交叉处理和变异处理,得到第二次迭代的多自由度调频信号,并通过参数函数计算初始迭代参数对应的第二次迭代参数。
在第二次迭代时,利用目标函数和第二次迭代参数计算各第二次迭代的多自由度调频信号的适应度,根据各第二次迭代的多自由度调频信号的适应度从第二次迭代的多自由度调频信号中选取多个父类多自由度调频信号,对选取的多个父类多自由度调频信号进行交叉处理和变异处理,得到第三次迭代的多自由度调频信号,并通过参数函数计算第二次迭代参数对应的第三次迭代参数。依次类推,直到进行交叉处理和变异处理的多自由度调频信号的适应度满足预设值。
本申请实施例中,根据所述增广拉格朗日遗传算法的算法模型,对第二次迭代的多自由度调频信号和所述第二次迭代参数进行迭代,得到第三次迭代的多自由度调频信号和第三次迭代参数;基于所述增广拉格朗日遗传算法的算法模型,持续对迭代的多自由度调频信号和迭代参数进行迭代,直到所述增广拉格朗日遗传算法的算法模型收敛为止,其中,增广拉格朗日遗传算法收敛表明已完成了对多自由度调频信号的设计与优化。
下面结合一具体示例,例如,合成孔径雷达SAR系统常用的大时宽信号设计参数,对多自由度调频信号的生成方法进行具体说明。
大时宽信号设计参数包括:
脉宽为10us,带宽为100MHz,采样频率为130MHz,时频信号被分成130段(即n=64),主瓣宽度控制因子β=1.18,以及初始的多自由度调频信号是基于驻定相位原理以sll=-35的泰勒窗函数为功率谱生成的。
基于上述设计参数,对于多自由度调频信号进行优化前和优化后的结果如表1所示:
表1
信号 | PSLR(dB) | IRW |
初始信号 | -34.9 | 1.19 |
优化后的信号 | -40.6 | 1.17 |
根据表1可知,通过本申请实施例提供的多自由度调频信号的生成方法,可以在不展宽主瓣的条件下,使旁瓣电平降低了5.7dB。此外,由于泰勒窗是一种性能较优的窗函数,所以此时初始多自由度调频信号的PSLR和IRW指标已经很好,通过本申请实施例提供的多自由度调频信号的生成能在不展宽主瓣的情况使旁瓣电平继续降低5.7dB,充分证实了本申请实施例提供的多自由度调频信号的生成方法的有效性。
如图4所示,为通过本申请实施例提供的多自由度调频信号的生成方法所生成的时频函数图,其中,横坐标为时间:脉宽为10us,纵坐标为频率:带宽为100MHz,在图中全部归一化。从图4中可以看出,优化后的信号也大致具有“S”型的瞬时频率函数,但是其存在微小的抖动,然而,正是这种抖动带来了信号性能的提升。
如图5所示,为优化后的多自由度调频信号的时域波形图,其中,横坐标为时间,纵坐标为幅度,幅度为1,脉宽为10us;实线表征该信号时域波形中的实部信息,虚线表征该时域波形中的虚部信息。
如图6所示,为优化后的多自由度调频信号自相关输出波形图,从该图可以看出,通过本申请实施例得到的多自由度调频信号与加窗产生等效旁瓣水平的相同时宽LFM信号相比,具有1.2dB的系统信噪比保持能力,具有很大的实际工程意义。
如图7所示,本申请实施例提供一种多自由度调频信号的生成装置700,所述装置包括:函数建立单元701、确定单元702、模型建立单元703、初始化单元704和迭代单元705;
所述函数建立单元701,用于基于分段非线性函数,构建多自由度调频信号的时域函数;
所述确定单元702,用于根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的冲激响应宽度IRW和峰值旁瓣比PSLR,基于所述IRW和所述PSLR建立信号优化模型;
所述模型建立单元703,用于基于所述信号优化模型,确定增广拉格朗日遗传算法的算法模型;
所述初始化单元704,用于根据选取的窗函数和驻定相位原理POSP,确定多自由度调频信号的初始化调频率控制点;
所述迭代单元705,用于利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
在一些实施例中,所述函数建立单元701,还用于基于分段线性函数的多自由度调频信号的瞬时调频率函数,确定多自由度调频信号的时频函数;
所述函数建立单元701,还用于根据所述分段非线性函数的多自由度调频信号的时频函数,确定所述分段非线性函数的多自由度调频信号的相位函数;
所述函数建立单元701,还用于根据所述相位函数和所述多自由度调频信号对应的幅度,确定所述多自由度调频信号的时域函数。
在一些实施例中,所述函数建立单元701,还用于在调频率关系平面内确定时间分段点t和调频率控制点r;
所述函数建立单元701,还用于基于所述时间分段点和所述调频率控制点,确定所述分段线性函数的多自由度调频信号的瞬时调频率函数;
所述函数建立单元701,还用于对所述分段线性函数的多自由度调频信号的瞬时调频率函数进行积分,确定分段非线性函数的多自由度调频信号的时频函数。
在一些实施例中,所述函数建立单元701,还用于对所述分段非线性函数的多自由度调频信号的时频函数进行积分,确定所述分段非线性函数的多自由度调频信号的相位函数。
在一些实施例中,确定单元702,还用于确定信号优化模型,所述信号优化模型为:
min PSLR(r),such that IRW(r)-β≤0,
其中,min PSLR(r),such that IRW(r)-β≤0用于表征在满足不展宽主瓣的情况下降低旁瓣,PSLR(r)为根据调频控制点r向量求得的多自由度调频信号的PSLR;IRW(r)为根据调频控制点r向量求得的多自由度调频信号的IRW;β为主瓣宽度控制因子;r1i1和r1i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;r2i1和r2i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;ΔT为脉冲宽度T平均分成2n+2段的时间间隔;B为预设的带宽。
在一些实施例中,所述模型建立单元703,还用于基于所述信号优化模型,利用增广拉格朗日算法确定目标函数;
Ψ(r,λ,s)=PSLR(r)-λlog(s-(IRW(r)-β));
其中,λ为拉格朗日算子,s为偏移量,λ和s均为非负数,β为参数;
所述模型建立单元703,还用于利用增广拉格朗日算法确定参数函数;
其中,μ为保证s-(IRW(r)-β)大于0的数值,k为非负数;
所述目标函数和所述参数函数组成所述增广拉格朗日遗传算法的算法模型。
所述初始单元704,还用于根据所述窗函数和驻定相位原理POSP,计算多自由度调频信号的群时延函数;根据所述群时延函数,确定所述多自由度调频信号对应的初始时频函数;根据所述多自由度调频信号对应的初始时频函数,确定所述多自由度调频信号的初始调频率函数;根据所述多自由度调频信号的初始调频率函数,确定所述多自由度调频信号的初始调频率控制点。
在一些实施例中,所述迭代单元705,还用于利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,确定第一多自由度调频信号的适应度和所述第一多自由度调频信号的迭代参数,基于所述第一多自由度调频信号的适应度和所述第一多自由度调频信号的迭代参数,对各所述第一多自由度调频信号进行交叉处理和变异处理,得到第二多自由度调频信号;其中,在第一次迭代中,所述第一多自由度调频信号为初始多自由度调频信号;
所述迭代单元705,还用于根据所述增广拉格朗日遗传算法,判断所述第二多自由度调频信号的适应度是否满足预设值,若不满足所述预设值,则继续对所述第二多自由度调频信号进行所述交叉处理和所述变异处理,直到进行所述交叉处理和所述变异处理的多自由度调频信号的适应度满足所述预设值。
需要说明的是,以上装置实施例的描述,与上述方法实施例的描述是类似的,具有同方法实施例相似的有益效果。对于本申请装置实施例中未披露的技术细节,请参照本申请方法实施例的描述而理解。
需要说明的是,本申请实施例中,如果以软件功能模块的形式实现上述的即时通讯方法,并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台即时通讯设备(可以是终端、服务器等)执行本申请各个实施例所述方法的全部或部分。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read Only Memory,ROM)、磁碟或者光盘等各种可以存储程序代码的介质。这样,本申请实施例不限制于任何特定的硬件和软件结合。
本申请实施例提供一种多自由度调频信号的生成装置,图8为本申请实施例多自由度调频信号的生成装置的组成结构示意图,如图8所示,所述装置800包括:一个处理器801、至少一个通信总线802、用户接口803、至少一个外部通信接口804和存储器805。其中,通信总线802配置为实现这些组件之间的连接通信。其中,用户接口803可以包括显示屏,外部通信接口804可以包括标准的有线接口和无线接口。其中所述处理器801,配置为执行存储器中存储的多自由度调频信号的生成程序,以实现上述实施例提供的多自由度调频信号的生成方法的步骤。
相应地,本申请实施例再提供一种存储介质(即计算机存储介质),所述计算机存储介质上存储有计算机可执行指令,所述计算机可执行指令被处理器执行时实现上述实施例提供的多自由度调频信号的生成方法的步骤。
以上多自由度调频信号的生成装置和计算机存储介质实施例的描述,与上述方法实施例的描述是类似的,具有同方法实施例相似的有益效果。对于本申请多自由度调频信号的生成装置和计算机存储介质实施例中未披露的技术细节,请参照本申请方法实施例的描述而理解。
应理解,说明书通篇中提到的“一个实施例”或“一实施例”意味着与实施例有关的特定特征、结构或特性包括在本申请的至少一个实施例中。因此,在整个说明书各处出现的“在一个实施例中”或“在一实施例中”未必一定指相同的实施例。此外,这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。应理解,在本申请的各种实施例中,上述各过程的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本申请实施例的实施过程构成任何限定。上述本申请实施例序号仅仅为了描述,不代表实施例的优劣。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
在本申请所提供的几个实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的。另外,所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口,设备或单元的间接耦合或通信连接,可以是电性的、机械的或其它形式的。
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的,作为单元显示的部件可以是、或也可以不是物理单元;既可以位于一个地方,也可以分布到多个网络单元上;可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
另外,在本申请各实施例中的各功能单元可以全部集成在一个处理单元中,也可以是各单元分别单独作为一个单元,也可以两个或两个以上单元集成在一个单元中;上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(Read Only Memory,ROM)、磁碟或者光盘等各种可以存储程序代码的介质。
或者,本申请上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本申请各个实施例所述方法的全部或部分。而前述的存储介质包括:移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种多自由度调频信号的生成方法,其特征在于,所述方法包括:
基于分段非线性函数,构建多自由度调频信号的时域函数;
根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的冲激响应宽度IRW和峰值旁瓣比PSLR,基于所述IRW和所述PSLR建立信号优化模型,所述信号优化模型为:
min PSLR(r),such that IRW(r)-β≤0,
其中,min PSLR(r),such that IRW(r)-β≤0用于表征在满足不展宽主瓣的情况下降低旁瓣,PSLR(r)为根据调频控制点求得的多自由度调频信号的PSLR;IRW(r)为根据调频控制点求得的多自由度调频信号的IRW;β为主瓣宽度控制因子;r1i1和r1i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;r2i1和r2i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;ΔT为脉冲宽度T平均分成2n+2段的时间间隔;B为预设的带宽;
基于所述信号优化模型,利用增广拉格朗日算法确定目标函数Ψ(r,λ,s):Ψ(r,λ,s)=PSLR(r)-λlog(s-(IRW(r)-β));
其中,r为调频率控制点,λ为拉格朗日算子,s为偏移量,λ和s均为非负数,β为常数;
利用增广拉格朗日算法确定参数函数λk+1:
其中,μ为保证s-(IRW(r)-β)大于0的数值,k为非负数;
所述目标函数和所述参数函数组成所述增广拉格朗日遗传算法的算法模型;
根据选取的窗函数和驻定相位原理POSP,确定多自由度调频信号的初始化调频率控制点;
利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
2.根据权利要求1所述的方法,其特征在于,所述基于分段非线性函数,构建多自由度调频信号的时域函数之前,所述方法还包括:
基于分段线性函数的多自由度调频信号的瞬时调频率函数,确定分段非线性函数的多自由度调频信号的时频函数;
所述基于分段非线性函数,构建多自由度调频信号的时域函数,包括:
根据所述分段非线性函数的多自由度调频信号的时频函数,确定所述分段非线性函数的多自由度调频信号的相位函数;
根据所述相位函数和所述多自由度调频信号对应的幅度,确定所述多自由度调频信号的时域函数。
3.根据权利要求2所述的方法,其特征在于,所述基于分段线性函数的多自由度调频信号的瞬时调频率函数,确定分段非线性函数的多自由度调频信号的时频函数之前,所述方法还包括:
在调频率关系平面内确定时间分段点和调频率控制点;
基于所述时间分段点和所述调频率控制点,确定所述分段线性函数的多自由度调频信号的瞬时调频率函数;
所述基于分段线性函数的多自由度调频信号的瞬时调频率函数,确定分段非线性函数的多自由度调频信号的时频函数,包括:
对所述分段线性函数的多自由度调频信号的瞬时调频率函数进行积分,确定所述分段非线性函数的多自由度调频信号的时频函数。
4.根据权利要求2所述的方法,其特征在于,所述根据所述分段非线性函数的多自由度调频信号的时频函数,确定所述分段非线性函数的多自由度调频信号的相位函数,包括:
对所述分段非线性函数的多自由度调频信号的时频函数进行积分,确定所述分段非线性函数的多自由度调频信号的相位函数。
5.根据权利要求1所述的方法,其特征在于,所述根据选取的窗函数和驻定相位原理POSP,确定多自由度调频信号的初始化调频率控制点,包括:
根据所述窗函数和驻定相位原理POSP,计算多自由度调频信号的群时延函数;
根据所述群时延函数,确定所述多自由度调频信号对应的初始时频函数;
根据所述多自由度调频信号对应的初始时频函数,确定所述多自由度调频信号的初始调频率函数;
根据所述多自由度调频信号的初始调频率函数,确定所述多自由度调频信号的初始调频率控制点。
6.根据权利要求1所述的方法,其特征在于,所述利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛,包括:
利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,确定第一多自由度调频信号的适应度和所述第一多自由度调频信号的迭代参数,基于所述第一多自由度调频信号的适应度和所述第一多自由度调频信号的迭代参数,对各所述第一多自由度调频信号进行交叉处理和变异处理,得到第二多自由度调频信号;其中,在第一次迭代中,所述第一多自由度调频信号为初始多自由度调频信号;
根据所述增广拉格朗日遗传算法,判断所述第二多自由度调频信号的适应度是否满足预设值,若不满足所述预设值,则继续对所述第二多自由度调频信号进行所述交叉处理和所述变异处理,直到进行所述交叉处理和所述变异处理的多自由度调频信号的适应度满足所述预设值。
7.一种多自由度调频信号的生成装置,其特征在于,所述装置包括:函数建立单元、确定单元、模型建立单元、初始单元和迭代单元;
所述函数建立单元,用于基于分段非线性函数,构建多自由度调频信号的时域函数;
所述确定单元,用于根据所述多自由度调频信号的时域函数,确定所述多自由度调频信号的自相关函数性能中的冲激响应宽度IRW和峰值旁瓣比PSLR,基于所述IRW和所述PSLR建立信号优化模型,所述信号优化模型为:
min PSLR(r),such that IRW(r)-β≤0,
其中,min PSLR(r),such that IRW(r)-β≤0用于表征在满足不展宽主瓣的情况下降低旁瓣,PSLR(r)为根据调频控制点求得的多自由度调频信号的PSLR;IRW(r)为根据调频控制点求得的多自由度调频信号的IRW;β为主瓣宽度控制因子;r1i1和r1i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;r2i1和r2i2为在所述调频率关系平面内构成的一段线性函数的纵坐标;ΔT为脉冲宽度T平均分成2n+2段的时间间隔;B为预设的带宽;
所述模型建立单元,用于基于所述信号优化模型,利用增广拉格朗日算法确定目标函数Ψ(r,λ,s):Ψ(r,λ,s)=PSLR(r)-λlog(s-(IRW(r)-β));
其中,r为调频率控制点,λ为拉格朗日算子,s为偏移量,λ和s均为非负数,β为常数;
利用增广拉格朗日算法确定参数函数λk+1:
其中,μ为保证s-(IRW(r)-β)大于0的数值,k为非负数;
所述目标函数和所述参数函数组成所述增广拉格朗日遗传算法的算法模型;
所述初始单元,用于根据选取的窗函数和驻定相位原理POSP,确定多自由度调频信号的初始化调频率控制点;
所述迭代单元,用于利用所述多自由度调频信号的时域函数和所述增广拉格朗日遗传算法的算法模型,对所述多自由度调频信号的初始化调频率控制点进行持续迭代,直到所述增广拉格朗日遗传算法收敛。
8.一种存储介质,其特征在于,所述存储介质上存储有计算机可执行指令,所述计算机可执行指令被处理器执行时,实现权利要求1至6任一项所述的多自由度调频信号的生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110662568.XA CN113408122B (zh) | 2021-06-15 | 2021-06-15 | 一种多自由度调频信号的生成方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110662568.XA CN113408122B (zh) | 2021-06-15 | 2021-06-15 | 一种多自由度调频信号的生成方法、装置及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113408122A CN113408122A (zh) | 2021-09-17 |
CN113408122B true CN113408122B (zh) | 2022-06-24 |
Family
ID=77683949
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110662568.XA Active CN113408122B (zh) | 2021-06-15 | 2021-06-15 | 一种多自由度调频信号的生成方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113408122B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108804736A (zh) * | 2018-03-20 | 2018-11-13 | 中国科学院电子学研究所 | 一种多自由度调频信号设计及优化的方法和装置 |
CN109343006A (zh) * | 2018-10-31 | 2019-02-15 | 中国科学院电子学研究所 | 基于增广拉格朗日遗传算法的nflm信号优化方法及装置 |
-
2021
- 2021-06-15 CN CN202110662568.XA patent/CN113408122B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108804736A (zh) * | 2018-03-20 | 2018-11-13 | 中国科学院电子学研究所 | 一种多自由度调频信号设计及优化的方法和装置 |
CN109343006A (zh) * | 2018-10-31 | 2019-02-15 | 中国科学院电子学研究所 | 基于增广拉格朗日遗传算法的nflm信号优化方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113408122A (zh) | 2021-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108804736B (zh) | 一种多自由度调频信号设计及优化的方法和装置 | |
CN107561484B (zh) | 基于内插互质阵列协方差矩阵重建的波达方向估计方法 | |
CN109765528B (zh) | 一种非线性调频信号的实时生成方法及装置 | |
CN108872929A (zh) | 基于内插虚拟阵列协方差矩阵子空间旋转不变性的互质阵列波达方向估计方法 | |
CN103984676A (zh) | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 | |
JP2017166918A (ja) | 信号処理装置、レーダ装置、信号処理方法およびプログラム | |
CN104079520B (zh) | 一种ofdm系统的脉冲干扰抑制方法 | |
US11041946B2 (en) | Model-based protection algorithms | |
CN109343059B (zh) | 一种正交非线性调频信号生成方法及装置 | |
CN113408122B (zh) | 一种多自由度调频信号的生成方法、装置及存储介质 | |
CN109358327B (zh) | 一种非线性调频信号的迭代生成方法、终端和存储介质 | |
CN109343006B (zh) | 基于增广拉格朗日遗传算法的nflm信号优化方法及装置 | |
CN107193002A (zh) | 一种可抑制宽带相位噪声的一维距离像高分辨成像方法 | |
Rosenberg et al. | Continuous sea clutter models for the mean backscatter and K-distribution shape | |
Wang et al. | Comb waveform optimisation with low peak‐to‐average power ratio via alternating projection | |
CN109471073B (zh) | 基于增广拉格朗日粒子群算法的nlfm信号生成方法及装置 | |
CN111025253A (zh) | 用于雷达脉冲压缩处理的窗函数波形联合设计方法 | |
CN110967686B (zh) | 一种相位编码信号生成方法及系统 | |
CN109492291B (zh) | 基于增广拉格朗日粒子群算法的nlfm信号优化方法及装置 | |
Zhao et al. | Over-sampled polyphase code design for physical implementation with spectral and correlation consideration | |
Kodituwakku et al. | Strong clutter suppression in non‐uniform PRF radar: techniques based on interpolation and adaptive processing | |
Xie et al. | Low sidelobe quasi-orthogonal NLFM waveforms with reciprocating frequency modulation | |
CN109446665B (zh) | 一种非线性调频信号优化方法及装置、存储介质 | |
CN109655793B (zh) | 高频雷达从瞬态干扰中提取回波的方法 | |
CN112162240A (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 |