CN109343006B - 基于增广拉格朗日遗传算法的nflm信号优化方法及装置 - Google Patents

Abstract

本发明实施例公开了一种基于增广拉格朗日遗传算法的NLFM信号优化方法，所述方法包括：基于贝塞尔曲线建立NLFM信号的时域函数；根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW和所述PSLR建立优化数学模型；基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；基于初始化函数初始化NLFM信号集合中各NLFM信号，得到相应的初始NLFM信号；利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛；本发明实施例还公开了一种基于增广拉格朗日遗传算法的NLFM信号优化装置。

Description

基于增广拉格朗日遗传算法的NFLM信号优化方法及装置

技术领域

本发明实施例涉及雷达发射信号领域，涉及但不限于一种基于增广拉格朗日遗传算法的非线性调频(Non-linear frequency modulation，NLFM)信号优化方法及装置。

背景技术

线性调频(Linear Frequency Modulation，LFM)信号因其对目标具有高分辨率成像能力，被常用于现代先进的雷达中，比如合成孔径雷达(Synthetic Aperture Radar,SAR)。雷达使用LFM信号产生的旁瓣水平高。为了抑制目标产生的旁瓣虚假信息并保持相位信息，通常对LFM信号进行加窗处理。然而，加窗抑制旁瓣会导致雷达系统信噪比(SNR)下降的问题，比如，现有德国宇航局的TerraSAR-X及意大利的Cosmo-SkyMed等SAR雷达信号加窗处理后，系统SNR下降1～2dB左右，这相当于发射机的能量下降了20％～37％。为了克服旁瓣抑制导致的信噪比丢失的难题，NLFM信号被雷达信号设计工作者所青睐。NLFM信号能保持雷达系统的信噪比的同时，又能达到加窗抑制旁瓣的效果。

此外，LFM加窗处理抑制旁瓣不可避免地会展宽主瓣，损失分辨率，其分辨率损失地大小依赖于窗函数的性能。而NLFM信号具有更多的自由度，可以通过优化的方法，找到最优的NLFM信号。相比于LFM信号加窗，NLFM信号具有更低的旁瓣以及更小的主瓣。

相关技术中，NLFM信号的研究主要集中在信号的设计及其应用。在脉冲多普勒雷达体制上进行NLFM信号的设计，主要包括三种设计方法：

1)基于驻定相位原理，通过设计特定的窗函数，进而求得完整的信号；

2)基于一些最优化方法如最小二乘逼近特定窗函数的功率谱；

3)为克服NLFM对多普勒频域敏感的问题，通过幅度加窗结合第1)种方法设计NLFM信号。

然而，上述方法在降低旁瓣的同时不可避免地展宽主瓣，并没有对主瓣宽度进行制约，从而导致不可避免的降低分辨率。此外一些现有的一些优化方案，如将时频空间分割成多个网格控制点，进行遍历搜索的方法，一次优化往往需要几天或者几个月，浪费大量的计算资源和时间，所以提升优化算法速度也是迫在眉睫的。贝塞尔函数，可以通过少量的控制点往往不超过10个，即可描述时频关系，从而大大降低了搜索时间。所以本发明专利结合贝塞尔函数和增广拉格朗日遗传算法优化设计NLFM信号，将大大减少优化时间，且能得到较好的优化结果。

发明内容

有鉴于此，本发明实施例希望提供一种基于增广拉格朗日遗传算法的 NLFM信号优化方法及装置，生成的NLFM信号能够在最大限度降低旁瓣的同时最低限度展宽主瓣。

本发明实施例的技术方案是这样实现的：

本发明实施例提供一种基于增广拉格朗日遗传算法的NLFM信号优化方法，所述方法包括：

基于贝塞尔曲线建立NLFM信号的时域函数；

根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW和所述PSLR建立优化数学模型；

基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；

基于初始化函数初始化NLFM信号集合中各NLFM信号，得到相应的初始 NLFM信号；

利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛。

本发明实施例还提供一种基于增广拉格朗日遗传算法的NLFM信号优化装置，所述装置包括：函数建立单元、条件确定单元、模型建立单元、初始单元和迭代单元；其中，

所述函数建立单元，用于基于贝塞尔曲线建立NLFM信号的时域函数；

所述条件确定单元，用于根据所述NLFM信号的时域函数确定所述NLFM 信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW 和所述PSLR建立优化数学模型；

所述模型建立单元，用于基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；

所述初始单元，用于基于初始化函数初始化NLFM信号集合中各NLFM信号的频率控制点，得到相应的初始NLFM信号；

所述迭代单元，用于利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛。

本发明实施例提供一种基于增广拉格朗日遗传算法的NLFM信号优化方法及装置，其中，所述方法包括：基于贝塞尔曲线建立NLFM信号的时域函数；根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW和所述PSLR建立优化数学模型；基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；基于初始化函数初始化NLFM信号集合中各NLFM信号的频率控制点，得到相应的初始NLFM信号；利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛；如此，得到的NLFM 信号在最大限度降低旁瓣的同时，最低限度展宽主瓣，适用于高分辨率合成孔径雷达高灵敏度低旁瓣成像的发射大时宽信号。

附图说明

图1为本发明实施例一提供的基于增广拉格朗日遗传算法的NLFM信号优化方法的流程示意图；

图2为本发明实施例二提供的基于增广拉格朗日遗传算法的NLFM信号优化方法的流程示意图；

图3为贝塞尔曲线表示的NLFM信号的时频关系示意图；

图4为本发明实施例三提供的NLFM信号时频关系曲线图；

图5为本发明实施例三提供的NLFM信号时域波形图；

图6为本发明实施例三提供的NLFM信号的自相关函数图；

图7为本发明实施例四提供的基于增广拉格朗日遗传算法的NLFM信号优化装置的结构示意图；

图8为本发明实施例五提供的基于增广拉格朗日遗传算法的NLFM信号优化装置的结构示意图。

具体实施方式

在本发明实施例中，基于贝塞尔曲线建立NLFM信号的时域函数；根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的MW和 PSLR，基于所述MW和所述PSLR建立优化数学模型；基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；基于初始化函数初始化NLFM信号集合中各NLFM信号的频率控制点，得到相应的初始NLFM信号；利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛。

下面结合附图对技术方案的实施作进一步的详细描述。

实施例一

本发明实施例一提供一种基于增广拉格朗日遗传算法的NLFM信号优化方法，如图1所示，所述方法包括：

S101、基于贝塞尔曲线建立NLFM信号的时域函数；

在一实施例中，所述基于贝塞尔曲线建立NLFM信号的时域函数包括：基于贝塞尔曲线建立NLFM信号的时频关系，根据所述NLFM信号的时频关系确定NLFM信号的相位函数，根据所述NLFM信号的相位函数确定所述NLFM 信号的时域函数。

在一实施例中，所述基于贝塞尔曲线建立NLFM信号的时频关系，根据所述NLFM信号的时频关系确定NLFM信号的相位函数，根据所述NLFM信号的相位函数确定所述NLFM信号的时域函数，包括：

在时频坐标关系平面内给定n+1个控制点，则NLFM信号为包括n+1个控制点的第一向量P，P＝{(x_i,y_i)|i＝0,1...,n}，x_i和y_i分别是第i个控制点P_i的横坐标和纵坐标，各控制点的横坐标x_i在时间轴上是均匀分布的，

为已知量；

采用贝塞尔曲线建立所述NLFM信号的时频关系为：

其中，

是一个二次项的系数，T_r为所述NLFM信号的脉冲宽度，x(t) 为控制点在时间轴上的向量，

f(t)为控制点在频率轴的向量，

根据所述时频关系确定所述NLFM信号的相位函数θ(t)为：

θ(t)＝2π·∫f(t)·x'(t)dt；

根据所述相位函数确定幅度为A的NLFM信号的时域函数为：

s(t)＝Aexp{-jθ(t)}。

在笛卡儿坐标系中，定义NLFM信号的时频关系坐标为(t,f)，NLFM 信号的脉冲宽度为T_r，信号带宽为B_r，采样频率为f_s，则NLFM信号时间间隔向量为

各控制点的横坐标x_i在时间轴上是均匀分布的，

为已知量，为了生成对称的波形，可由含

个频率控制点的第二向量p来定义NLFM信号，第二向量p为第一向量P的n+1个控制点中的前

个控制点在频率轴上的向量，p＝{y₁,y₂,...,y_(n/2-1)},

则y_(n-i)＝-y_i,

这里，将第二向量p称为频率控制点向量。由此，可以求得NLFM信号的调频率函数为：

这里，调频率函数表征NLFM信号的调频斜率。

其中，

在实际应用中，第二向量p也可为第一向量P的n+1个控制点在频率轴上的向量，此时，p＝{y₀,y₁,...,y_n}。

需要说明的是，P是在时间/频率平面上的控制点向量，控制点向量是由平面上的多个控制点构成的向量，一个控制点包含横坐标和纵坐标，横坐标在这里表示时间，是均匀分布的；纵坐标表示频率，一般不是均匀分布的，纵坐标的位置构成一个第二向量p，将第二向量p中各纵坐标称为频率控制点。

S102、根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度(Main lobe Width，MW)和峰值旁瓣比(Peak Side Lobe Ratio，PSLR)，基于所述MW和所述PSLR建立优化数学模型；

NLFM信号的时频信号为：s(t)＝Aexp{-jθ(t)}时，则NLFM信号的自相关函数为：

将NLFM的信号幅度转变为 dB形式，计算NLFM信号的PSLR和MW。PSLR是除了主瓣之外的旁瓣最大值，MW可为主瓣-3dB处对应宽度，PSLR和MW的定义如下：

PSLR：最高旁瓣与主瓣峰值高度的比值，单位为dB；

3dB MW：3dB主瓣宽度的大小，一般归一化为采样点。

基于计算的MW和PSLR建立优化数学模型，其中，建立的优化数学模型包括两种优化数学模型：

模型一、所述基于所述MW和所述PSLR建立优化数学模型包括：满足不展宽主瓣的情况下尽可能地降低旁瓣的条件，所述优化数学模型为：

模型二、所述基于所述MW和所述PSLR建立优化数学模型包括：满足不升高旁瓣的情况下可能地降低主瓣的条件，所述优化数学模型为：

其中，PSLR(p)为根据第二向量p求得的NLFM信号的所述PSLR，MW(p) 为根据第二向量p求得的NLFM信号的MW；第二向量p为第一向量P的n+1 个控制点中的前

个控制点在频率轴上的向量， p＝{y₁,y₂,...,y_(n/2-1)},

y_i表示第i个控制点的纵坐标， i＝0,1,...,n，B_r为NLFM信号的带宽。

其中，在模型一中，PSLR(p)为优化数学模型的目标函数，MW(p)代表非线性约束；在模型二中，MW(p)为优化数学模型的目标函数，PSLR(p)代表非线性约束。

需要说明的是，MW(p)也可表示为MW(p)-a，a≥0，a表征NLFM信号的主瓣的放宽程度。

NLFM信号的自相关函数的理想性能为：尽量窄的主瓣，尽可能低的峰值旁瓣比与快速下降的旁瓣波动包络。然而这三个理想的性能是不能同时满足的。通常意义下，自相关函数的主瓣宽度定义为辐射强度低于峰值的3 dB的两点间的夹角间的波瓣宽度。这里，控制点决定了NLFM信号的具体形式，一旦控制点的个数被确定后，控制点的横坐标x＝{x₀,x₁,...,x_n}在时间轴上均匀分布，为已知量，因此，NLFM信号能由含

个频率控制点的第二向量p来定义。

S103、基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；

这里，所述基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型包括：

基于所述优化数学模型，利用增广拉格朗日算法确定目标函数：

Θ(p,λ,s)＝g(p)-λslog(s-c(p))；

其中，λ为拉格朗日算子，s为偏移量，λ、s均为非负数；g(p)为第二向量p对应的PSLR或MW；c(p)为第二向量对应的MW或PSLR；

利用增广拉格朗日算法确定参数函数：

μ_k+1＝μ_k，s_k+1＝μ_k+1λ_k+1；

其中，μ为保证s-c(p)大于0的数值；

根据所述目标函数计算当前NLFM信号集合中各NLFM信号的适应度；根据各NLFM信号的适应度从NLFM信号集合中选择父体NLFM信号；对所述父体NLFM信号进行交叉处理、变异处理，得到下一次迭代的NLFM信号；

根据所述参数函数计算当前拉格朗日算子和偏移量对应的下一次迭代的拉格朗日算子和偏移量。

这里，当优化数学模型为模型一时，g(p)为第二向量p对应的PSLR；c(p) 为第二向量对应的MW；当优化数学模型为模型二时，g(p)为第二向量p对应的MW；c(p)为第二向量对应的PSLR。

其中，增广拉格朗日遗传算法是遗传算法的推广形式，是遗传算法和增广拉格朗日算法结合的解决复杂约束优化的算法。

增广拉格朗日算法的数学描述为：

其中，λ_i为拉格朗日乘子，为一个非负的数，s_i为一个非负的数，表征偏移量，来保证对数的真数非0，ρ为惩戒因子，ceq_i(x)和c_i(x)分别代表等式约束和非线性不等式约束，g(x)为适应度函数，m代表非线性约束的个数，mt代表总约束个数。

当引入优化数学模型时，增广拉格朗日算法即目标函数可以描述成：

Θ(p,λ,s)＝g(p)-λslog(s-c(p))；

其中，优化数学模型为模型一时，g(p)＝PSLR(p)为根据第二向量p求得的信号的PSLR，c(p)＝MW(p)为根据第二向量p求得的信号的MW。优化数学模型为模型二时，g(p)＝MW(p)为根据第二向量p求得的信号的,MW， c(p)＝PSLR(p)为根据第二向量p求得的信号的PSLR。

增广拉格朗日遗传算法把具体的求解问题分成两部分：一部分为传统的遗传算法，一部分为增广拉格朗日算法。

增广拉格朗日算法用于解决约束问题，具体地，它根据以下公式不断的更新λ和s：

μ_k+1＝μ_k，s_k+1＝μ_k+1λ_k+1。

遗传算法将优化问题模拟成“适者生存”的自然选择的动态优化过程。首先，依据目标函数值进行染色体适度按一定的规则进行选择染色体。其次，被选中染色体对依据交配概率r_c进行交叉产生后代；最后，依据一定的变异概率 r_m对染色体的基因进行变异操作，在搜索变量空间中产生新的个体。在整个迭代更新优化过程中，适应度高的染色体个体被选中用于产生后代的概率大，适应度差的个体被更优的后代取代。

在实际应用中，根据各NLFM信号的适应度从NLFM信号集合中选择父体 NLFM信号时，各个NLFM信号被选中的概率与其适应度大小成正比，这里，可采用轮盘选择法(Roulette Wheel Selection，RWS)或其他的算法从NLFM信号集合中选择父体NLFM信号。

S104、基于初始化函数初始化NLFM信号集合中各NLFM信号，得到相应的初始NLFM信号；

针对NLFM信号，把每个含

个频率控制点的第二向量p看成一个染色体,把第二向量p的每个频率控制点看成基因，对群体中各染色体进行初始化，即为对NLFM信号集合中各NLFM信号进行初始化。

在一实施例中，所述基于初始化函数初始化NLFM信号集合中各NLFM信号，包括：

基于初始化函数对各NLFM信号的频率控制点进行初始化；其中，所述 NLFM信号由包括

个频率控制点的第二向量p获得；

对第k个第二向量p的第j个频率控制点进行初始化的所述初始化函数为：

其中，B_r为NLFM信号的带宽，rand为随机数。

S105、利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛。

当在S104中对NLFM信号集合中各NLFM信号初始化得到初始NLFM 信号后，通过S103确定的算法模型对初始NLFM信号进行迭代，直至增广拉格朗日遗传算法收敛，以对初始NLFM信号进行优化。

在一实施例中，所述利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛包括：

利用所述算法模型，根据所述初始迭代参数计算所述初始NLFM信号的适应度和第二次迭代参数，以及基于各初始NLFM信号的适应度，对各初始NLFM 信号进行选择、交叉处理和变异处理得到第二次迭代的NLFM信号；

利用所述算法模型，对所述第二次迭代的NLFM信号和所述第二次迭代参数进行迭代，得到第三次迭代的NLFM信号和第三次迭代参数；并利用所述算法模型，持续对迭代的NLFM信号和迭代参数进行迭代，直到所述增广拉格朗日遗传算法收敛。

在对初始NLFM信号进行迭代之前，可先初始化算法参数，算法参数包括迭代参数和遗传参数，其中，迭代参数包括拉格朗日算子λ和偏移量s，遗传参数包括交叉概率r_c和变异概率r_m，遗传参数还可包括染色体个数popsize，染色体个数popsize表征NLFM信号集合中NLFM信号的个数。初始化的迭代参数可称为初始迭代参数，以用于迭代过程中的第一次迭代。

在第一次迭代时，利用所述目标函数和初始化迭代参数计算各初始NLFM 信号的适应度，根据各初始NLFM信号的适应度从初始NLFM信号中选取多个父体NLFM信号，对选取的多个父体NLFM信号进行交叉处理和变异处理，得到第二次迭代的NLFM信号，并通过参数函数计算初始化迭代参数对应的第二次迭代参数。

在第二次迭代时，利用所述目标函数和第二次迭代参数计算各第二次迭代的NLFM信号的适应度，根据各第二次迭代的NLFM信号的适应度从第二次迭代的NLFM信号中选取多个父体NLFM信号，对选取的多个父体NLFM信号进行交叉处理和变异处理，得到第三次迭代的NLFM信号，并通过参数函数计算第二次迭代参数对应的第三次迭代参数。

其中，第一次迭代、第二次迭代都是基于算法模型对当前的NLFM信号和迭代参数进行迭代更新，在第二次迭代完成后，基于算法模型持续进行迭代，直到增广拉格朗日遗传算法收敛。

这里，可当两次迭代中的最优染色体对应的适应度不再改变时，确定增广拉格朗日遗传算法收敛；也可以设定迭代阈值，当迭代次数达到迭代阈值时，确定增广拉格朗日遗传算法收敛。

这里，RWS的一种算法为：当从popsize个NLFM信号中选择父体NLFM 信号时，个体p的适应度为Θ(p,λ,s)，即第j个个体的适应度为Θ(p_j) j＝(1,2,...,popsize)，p_j的选择概率为

RWS的过程为：

S1、在[0,1]内产生一个均匀分布的随机数r；

S2、若r≤q₁，该NLFM信号p₁被选中；

S3、若q_m-1＜r≤q_m(2≤m≤popsize)，则NLFM信号p_m被选中。

其中，q_c被称为NLFM信号p_c(1≤c≤popsize)的积累概率，积累概率的公式为：

在本发明实施例中，对RWS的具体的算法不进行任何限制。

在迭代过程中，定义选中的父体NLFM信号的总个数为K，在交叉过程中，先为这K个父体NLFM信号产生K个随机数，如果父体NLFM信号对应的随机数低于交叉概率r_c，则表明这些父体NLFM信号被选中用于交叉操作。这里，采取1个交叉点交叉操作，交叉点的位置由随机产生。父体 NLFM信号在交叉点互换频率控制点产生第二次迭代的NLFM信号。变异是指染色体中的基因发生变化的操作。变异的频率控制点被随机选中。因 NLFM信号优化，每个父体NLFM信号的频率控制点个数为

则频率控制点的总个数L为

变异的个数M由变异概率决定，具体为 M＝r_mL。在L个频率控制点中随机选中M个频率控制点进行变异操作，当第 k个NLFM信号的第i个频率控制点变异时，变异运算为：

p_k(i)＝p_k(i)*(1+rand)。

在本发明实施例中，基于贝塞尔曲线建立NLFM信号的时域函数；根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW和所述PSLR建立优化数学模型；基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；基于初始化函数初始化NLFM信号集合中各NLFM信号，得到相应的初始 NLFM信号；利用所述算法模型，对设定的初始迭代参数和所述初始NLFM 信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛。如此，能够在不展宽主瓣的情况下尽可能地降低旁瓣，或者不升高旁瓣的情况下尽可能缩小主瓣，得到的NLFM信号在最大限度降低旁瓣的同时，最低限度展宽主瓣，适用于高分辨率合成孔径雷达高灵敏度低旁瓣成像的发射大时宽信号。

实施例二

在本发明实施例中，对本发明实施例一提供的基于增广拉格朗日遗传算法的NLFM信号优化方法进一步进行描述，如图2所示，所述方法包括：

S201、利用贝塞尔曲线描述NLFM信号的时频关系，求解NLFM信号的相位函数，进而定义信号时域函数；

在笛卡儿坐标系中，把NLFM信号的时频关系坐标定义为(t,f)。假设信号的脉冲宽度为T_r，信号带宽为B_r，采样频率为f_s，则信号时间间隔向量为

在时频坐标关系平面内给定n+1个控制点，定义为 P＝{(x_i,y_i)|i＝0,1...,n}，采用贝塞尔曲线描述NLFM信号的时频关系，可以描述为：

其中，

是一个二次项的系数，控制点在时间轴上的向量x(t)为：

控制点在频率轴的向量f(t)为

其中，贝塞尔曲线表示的NLFM信号的时频关系如图3所示。

式(2)和(3)中的x_i和y_i分别是控制点P_i的横坐标和纵坐标。控制点的横坐标在时间轴上是均匀分布的，

为了生成对称的波形，用向量p表示前

个控制点在频率轴上的分量， p＝{y₁,y₂,...,y_(n/2-1)},

则y_(n-i)＝-y_i,

由此，可以求得NLFM信号的调频率函数为：

式(4)中的

和

的表达式分别为：

设NLFM信号的相位为θ(t)，对相位求导可得瞬时频率：

由此，可以得到NLFM信号的相位的表达式：

θ(t)＝2π·∫f(t)·x'(t)dt 式(8)

则幅度为A的NLFM信号的时域函数可以表示为：

s(t)＝Aexp{-jθ(t)} 式(9)

S202、根据信号的自相关函数的性能，定义了信号优化数学模型；

NLFM信号的自相关函数的理想性能为：尽量窄的主瓣，尽可能低的峰 PSLR与快速下降的旁瓣波动包络。然而这三个理想的性能是不能同时满足的。通常意义下，自相关函数的主瓣宽度定义为低于峰值的3dB的两点间的夹角间的波瓣宽度。由式(1)至式(9)可知，控制点决定了NLFM信号具体形式，NLFM信号的控制点的个数被确定后，各控制点的横坐标 x＝{x₀,x₁,...,x_n}在时间轴上均匀分布，且为已知量。NLFM信号即能由含

个频率控制点的向量p来定义。这里，选取SAR图像点目标冲击响应函数的两个关键性能PSLR和MW来提出NLFM信号优化数学模型，即优化数学模型，其为一个非线性约束优化问题。

当要求在不展宽主瓣的情况下尽可能地降低旁瓣，该问题可以描述成：

当要求在不升高旁瓣的情况下可能地降低主瓣，该问题可以描述成：

其中，式(10)中，PSLR(p)为优化数学模型的目标函数，MW(p)代表非线性约束；式(11)中，MW(p)为优化数学模型的目标函数，PSLR(p)代表非线性约束。

上述两个问题本质上为一个问题，在之后讨论中，以式(10)为例。

S203、根据优化数学模型，利用增广拉格朗日遗传算法进行优化，得到最优NLFM信号。

增广拉格朗日遗传算法是遗传算法的推广形式，是遗传算法和广义拉格朗日算法结合的解决复杂约束优化的先进算法。

广义拉格朗日算法的数学描述为：

其中，λ_i为拉格朗日乘子，为一个非负的数；s_i为一个非负的数，代表整体的偏移量来保证对数的真数非0；ρ为惩戒因子；ceq_i(x)和c_i(x)分别代表等式约束和非线性不等式约束；g(x)为适应度函数，m代表非线性约束的个数，mt 代表总约束个数。

要求在不展宽主瓣的情况下尽可能地降低旁瓣时，该问题即目标函数可以描述成：

Θ(p,λ,s)＝g(p)-λslog(s-c(p)) 式(13)

其中，g(p)＝PSLR(p)为根据向量p求得的NLFM信号的PSLR， c(p)＝MW(p)为根据向量p求得的NLFM信号的MW。

其中，在式(13)中，在不加约束的时候，第二向量p的目标函数即PSLR (p)，加了约束，则目标函数应为PSLR(p)加上约束部分，目标函数值即为适应度。

增广拉格朗日遗传算法把具体的求解问题分成两部分：一部分为传统的遗传算法，一部分为增广拉格朗日算法。增广拉格朗日算法用于解决约束问题，具体地，它根据式(14)不断的更新λ和s，

遗传算法将优化问题模拟成“适者生存”的自然选择的动态优化过程。在搜索空间中，染色体代表具体求解问题的所要确定的变量，遗传算法通常包括染色体的选择、交叉和变异操作。首先，根据求解问题对变量进行编码，依据目标函数值进行染色体适度计算，并按一定的规则进行选择染色体；其次，被选中染色体对依据交配概率r_c进行交叉产生后代；最后，依据一定的变异概率r_m对染色体的基因进行变异操作，在搜索变量空间中产生新的个体。在整个迭代更新优化过程中，适应度高的染色体个体被选中用于产生后代的概率大，适应度差的个体被更优的后代取代。

针对NLFM信号优化，在问题编码过程中，把每个含

个频率控制点的向量p看成一个染色体,把向量p的每个频率控制点看成基因，按照一定的规则对p进行初始化，对第k个染色体的第j个基因按式(15)进行初始化；

该发明采用RWS进行选择父体进行交叉过程。假设选中的染色体的总个数为。在交叉过程中，先为这K个染色体产生K个随机数。如果染色体对应的随机数低于交叉概率r_c，则表明这些染色体被选中用于交叉操作。这里，采取1 个交叉点交叉操作，交叉点的位置由随机产生。父体在交叉点互换基因产生新的染色体。变异是指染色体中的基因发生变化的操作。变异的基因被随机选中。因NLFM信号优化，每个染色体的基因个数为

则基因的总个数L为

变异的个数M由变异概率决定，具体为M＝r_mL。在L随机选中M个进行变异操作，变异运算为：

p_k(i)＝p_k(i)*(1+rand) 式(16)

综上，基于遗传算法的NLFM信号优化过程可总结成算法如下所示：

步骤1、初始化算法参数：设定初始化算法参数：染色体个数popsize，交叉概率r_c，变异概率r_m，拉格朗日算子λ，偏移量s。

其中，染色体个数popsize为正数，popsize的值越大，优化结果越好，收敛速度越慢；0<r_c<1，r_c的值越大，全局优化结果越好，收敛速度越慢；0<r_m<1， r_m的值越大，全局优化结果越好，收敛速度越慢。

步骤2、根据式(15)初始化向量p，即初始化染色体。

步骤3、对每个染色体所表示的频率控制点，依据(1)至(9)式确定对应的NLFM信号，进行自相关运算,然后进行计算输出信号的主瓣宽度MW及峰值旁瓣比PSLR，最后依据式(13)计算染色体的适应度。并依据式(14) 利用增广拉格朗日算法计算下一次优化的λ和偏移量s。

步骤4、根据轮盘选择法选择染色体。

步骤5、根据交叉概率r_c进行一个交叉点交叉操作。

步骤6、根据变异概率r_m进行式(16)变异操作。

步骤7、循环步骤2至步骤6，直至算法收敛。

算法收敛表明已完成了NLFM信号设计和优化。

本发明实施例提供的基于增广拉格朗日遗传算法的NLFM信号优化方法，基于这样一个事实：利用贝塞尔曲线描述NLFM信号的时频关系，求解NLFM 信号的相位函数，进而定义时域函数；根据信号的自相关函数的性能，定义了信号优化数学模型；根据优化数学模型对信号进行初始化，利用增广拉格朗日遗传算法进行快速优化，即可得到最优信号。采用本发明的方法，能设计得到一种适用于高分辨率合成孔径雷达高灵敏度低旁瓣成像的发射大时宽信号。

实施例三

在本发明实施例中，结合具体的合成孔径雷达系统常用的大时宽信号设计参数对本发明实施例提供基于增广拉格朗日遗传算法的NLFM信号优化方法得到的NLFM信号进行说明。

大时宽信号设计参数包括：

脉宽为10us，带宽为100MHz，采样频率为130MHz，其中，初始化NLFM 信号是根据式(1)-(9)和(15)生成的。

基于上述设计参数分别设计包括LFM波形、泰勒窗设计波形、优化波形一和优化波形二，其中，优化波形一和优化波形二都为根据本发明实施例提供的基于增广拉格朗日遗传算法的NLFM信号优化方法生成的NFLM信号，优化波形一和优化波形二的约束MW的放宽程度不一样。各波形的性能比较如表1所示。

表1、不同波形信号的性能比较

从表1可以看出，使用本发明实施例提供的基于增广拉格朗日遗传算法的NLFM信号优化方法生成的NLFM信号相比于LFM信号其在旁瓣轻微展宽的情况下，旁瓣能够降低10dB以上。此外，值得一提的是泰勒窗是一种性能优秀的窗函数，相比于泰勒窗这一性能较好的窗，优化的结果也有轻微的提高。此外本次优化时间仅仅为一个小时，优化时间大大缩减，提升了效率。

图4至图6分别为优化波形一对应的NLFM信号的时频关系曲线、时域波形及自相关函数输出波形。图4为优化信号的时频关系曲线图，与加窗产生等效旁瓣水平的相同时宽带宽的LFM信号相比，具有0.5～1.2dB的系统信噪比保持能力，具有很大的实际工程意义。

实施例四

本发明实施例提供一种基于增广拉格朗日遗传算法的NLFM信号优化装置，所述装置包括：函数建立单元701、条件确定单元702、模型建立单元703、初始单元704和迭代单元705；其中，

函数建立单元701，用于基于贝塞尔曲线建立NLFM信号的时域函数；

条件确定单元702，用于根据所述NLFM信号的时域函数确定所述NLFM 信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW 和所述PSLR建立优化数学模型；

模型建立单元703，用于基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；

初始单元704，用于基于初始化函数初始化NLFM信号集合中各NLFM信号，得到相应的初始NLFM信号；

迭代单元705，用于利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛。

在一实施例中，函数建立单元701，用于：

基于贝塞尔曲线建立NLFM信号的时频关系，根据所述NLFM信号的时频关系确定NLFM信号的相位函数，根据所述NLFM信号的相位函数确定所述 NLFM信号的时域函数。

在一实施例中，函数建立单元701，用于：

在时频坐标关系平面内给定n+1个控制点，由n+1个控制点构成所述NLFM 信号的第一向量P，其中，P＝{(x_i,y_i)|i＝0,1...,n}，x_i和y_i分别是第i个控制点P_i的横坐标和纵坐标，各控制点的横坐标x_i在时间轴上是均匀分布的，

为已知量；

采用贝塞尔曲线建立所述NLFM信号的时频关系为：

其中，

是一个二次项的系数，T_r为所述NLFM信号的脉冲宽度，x(t) 为控制点在时间轴上的向量，

f(t)为控制点在频率轴的向量，

根据所述时频关系确定所述NLFM信号的相位函数θ(t)为：

θ(t)＝2π·∫f(t)·x'(t)dt；

根据所述相位函数确定幅度为A的NLFM信号的时域函数为：

s(t)＝Aexp{-jθ(t)}。

在一实施例中，条件确定单元702，用于：

满足不展宽主瓣的情况下尽可能地降低旁瓣的条件，所述优化数学模型为：

其中，PSLR(p)为根据第二向量p求得的NLFM信号的所述PSLR，MW(p) 为根据第二向量p求得的NLFM信号的MW；第二向量p为第一向量P的n+1 个控制点中的前

个控制点在频率轴上的向量，p＝{y₁,y₂,...,y_(n/2-1)},

y_i表示第i个控制点的纵坐标， i＝0,1,...,n，B_r为NLFM信号的带宽。

在一实施例中，条件确定单元702，用于：

满足不升高旁瓣的情况下可能地降低主瓣的条件，所述优化数学模型为：

其中，PSLR(p)为根据第二向量p求得的NLFM信号的所述PSLR，MW(p) 为根据第二向量p求得的NLFM信号的MW；第二向量p为第一向量P的n+1 个控制点中的前

个控制点在频率轴上的向量， p＝{y₁,y₂,...,y_(n/2-1)},

y_i表示第i个控制点的纵坐标， i＝0,1,...,n，B_r为NLFM信号的带宽。

在一实施例中，模型建立单元703，用于：

基于所述优化数学模型，利用增广拉格朗日算法确定目标函数：

Θ(p,λ,s)＝g(p)-λslog(s-c(p))；

其中，λ为拉格朗日算子，s为偏移量，λ、s均为非负数；g(p)为第二向量p对应的PSLR或MW；c(p)为第二向量对应的MW或PSLR；

利用增广拉格朗日算法确定参数函数：

μ_k+1＝μ_k，s_k+1＝μ_k+1λ_k+1；

其中，μ为保证s-c(p)大于0的数值；

根据所述目标函数计算当前NLFM信号集合中各NLFM信号的适应度；根据各NLFM信号的适应度从NLFM信号集合中选择父体NLFM信号；对所述父体NLFM信号进行交叉处理、变异处理，得到下一次迭代的NLFM信号；

根据所述参数函数计算当前拉格朗日算子和偏移量对应的下一次迭代的拉格朗日算子和偏移量。

在一实施例中，初始单元704，用于：

基于初始化函数对各NLFM信号的频率控制点进行初始化；其中，所述 NLFM信号由包括

个频率控制点的第二向量p获得；

对第k个第二向量p的第j个频率控制点进行初始化的所述初始化函数为：

其中，B_r为NLFM信号的带宽，rand为随机数。

在一实施例中，迭代单元705，用于：

所述利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛包括：

利用所述算法模型，根据所述初始迭代参数计算所述初始NLFM信号的适应度和第二次迭代参数，以及基于各初始NLFM信号的适应度，对各初始NLFM 信号进行选择、交叉处理和变异处理得到第二次迭代的NLFM信号；

利用所述算法模型，对所述第二次迭代的NLFM信号和所述第二次迭代参数进行迭代，得到第三次迭代的NLFM信号和第三次迭代参数；并利用所述算法模型，持续对迭代的NLFM信号和迭代参数进行迭代，直到所述增广拉格朗日遗传算法收敛。

需要说明的是，以上装置实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本发明装置实施例中未披露的技术细节，请参照本发明方法实施例的描述而理解。

需要说明的是，本发明实施例中，如果以软件功能模块的形式实现上述的即时通讯方法，并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台即时通讯设备(可以是终端、服务器等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read Only Memory，ROM)、磁碟或者光盘等各种可以存储程序代码的介质。这样，本发明实施例不限制于任何特定的硬件和软件结合。

实施例五

本发明实施例提供一种基于增广拉格朗日遗传算法的NLFM信号优化装置，图8为本发明实施例基于增广拉格朗日遗传算法的NLFM信号优化装置的组成结构示意图，如图8所示，所述装置80包括：一个处理器801、至少一个通信总线802、用户接口803、至少一个外部通信接口804和存储器805。其中，通信总线802配置为实现这些组件之间的连接通信。其中，用户接口803可以包括显示屏，外部通信接口804可以包括标准的有线接口和无线接口。其中所述处理器801，配置为执行存储器中存储的计算机程序，以实现上述实施例提供的基于增广拉格朗日遗传算法的NLFM信号优化方法的步骤。

相应地，本发明实施例再提供一种存储介质(即计算机存储介质)，所述计算机存储介质上存储有计算机可执行指令，所述该计算机可执行指令被处理器执行时实现上述实施例提供的基于增广拉格朗日遗传算法的NLFM信号优化方法的步骤。

以上基于增广拉格朗日遗传算法的NLFM信号优化装置和计算机存储介质实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本发明基于增广拉格朗日遗传算法的NLFM信号优化装置和计算机存储介质实施例中未披露的技术细节，请参照本发明方法实施例的描述而理解。

应理解，说明书通篇中提到的“一个实施例”或“一实施例”意味着与实施例有关的特定特征、结构或特性包括在本发明的至少一个实施例中。因此，在整个说明书各处出现的“在一个实施例中”或“在一实施例中”未必一定指相同的实施例。此外，这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。应理解，在本发明的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。

需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。

在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其它形式的。

上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元；既可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。

另外，在本发明各实施例中的各功能单元可以全部集成在一个处理单元中，也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一个单元中；上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。

本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器(Read Only Memory，ROM)、磁碟或者光盘等各种可以存储程序代码的介质。

或者，本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括：移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (8)

1.一种基于增广拉格朗日遗传算法的非线性调频NLFM信号优化方法，其特征在于，所述方法包括：

基于贝塞尔曲线建立NLFM信号的时域函数；

根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW和所述PSLR建立优化数学模型；其中，根据不等式约束的非线性规划模型，在满足不展宽主瓣的情况下降低旁瓣的条件时，将PLSR(p)作为优化数学模型的目标函数，将MW(p)作为约束函数，建立所述优化数学模型，或，在满足不升高旁瓣的情况下降低主瓣的条件时，将MW(p)作为优化数学模型的目标函数，将PLSR(p)作为约束函数，建立所述优化数学模型；其中，在时频坐标关系平面内给定n+1个控制点，由n+1个控制点构成所述NLFM信号的第一向量P，第二向量p为第一向量P的n+1个控制点中的前

个控制点在频率轴上的向量，所述PSLR(p)为根据第二向量p求得的NLFM信号的所述PSLR，MW(p)为根据第二向量p求得的NLFM信号的MW；

基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；

基于初始化函数初始化NLFM信号集合中各NLFM信号，得到相应的初始NLFM信号；

利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛；

所述基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型包括：

基于所述优化数学模型，利用增广拉格朗日算法确定目标函数：

Θ(p,λ,s)＝g(p)-λslog(s-c(p))；

其中，λ为拉格朗日算子，s为偏移量，λ、s均为非负数；g(p)为第二向量p对应的PSLR或MW，c(p)为第二向量对应的MW或PSLR；

利用增广拉格朗日算法确定参数函数：

μ_k+1＝μ_k，s_k+1＝μ_k+1λ_k+1；

其中，μ为保证s-c(p)大于0的数值；

根据所述目标函数计算当前NLFM信号集合中各NLFM信号的适应度；根据各NLFM信号的适应度从NLFM信号集合中选择父体NLFM信号；对所述父体NLFM信号进行交叉处理、变异处理，得到下一次迭代的NLFM信号；

根据所述参数函数计算当前拉格朗日算子和偏移量对应的下一次迭代的拉格朗日算子和偏移量。

2.根据权利要求1所述的方法，其特征在于，所述基于贝塞尔曲线建立NLFM信号的时域函数包括：

基于贝塞尔曲线建立NLFM信号的时频关系，根据所述NLFM信号的时频关系确定NLFM信号的相位函数，根据所述NLFM信号的相位函数确定所述NLFM信号的时域函数。

3.根据权利要求2所述的方法，其特征在于，所述基于贝塞尔曲线建立NLFM信号的时频关系，根据所述NLFM信号的时频关系确定NLFM信号的相位函数，根据所述NLFM信号的相位函数确定所述NLFM信号的时域函数，包括：

所述第一向量P＝{(x_i,y_i)|i＝0,1...,n}，x_i和y_i分别是第i个控制点P_i的横坐标和纵坐标，各控制点的横坐标x_i在时间轴上是均匀分布的，

为已知量；

采用贝塞尔曲线建立所述NLFM信号的时频关系为：

其中，

是一个二次项的系数，T_r为所述NLFM信号的脉冲宽度，x(t)为控制点在时间轴上的向量，

f(t)为控制点在频率轴的向量，

根据所述时频关系确定所述NLFM信号的相位函数θ(t)为：

θ(t)＝2π·∫f(t)·x'(t)dt；

根据所述相位函数确定幅度为A的NLFM信号的时域函数为：

s(t)＝Aexp{-jθ(t)}。

4.根据权利要求3所述的方法，其特征在于，所述基于所述MW和所述PSLR建立优化数学模型包括：

满足不展宽主瓣的情况下降低旁瓣的条件，所述优化数学模型为：

其中，所述第二向量p＝{y₁,y₂,...,y_(n/2-1)},

y_i表示第i个控制点的纵坐标，i＝0,1,...,n，B_r为NLFM信号的带宽。

5.根据权利要求3所述的方法，其特征在于，所述基于所述MW和所述PSLR建立优化数学模型包括：

满足不升高旁瓣的情况下降低主瓣的条件，所述优化数学模型为：

其中，所述第二向量p＝{y₁,y₂,...,y_(n/2-1)},

y_i表示第i个控制点的纵坐标，i＝0,1,...,n，B_r为NLFM信号的带宽。

6.根据权利要求3所述的方法，其特征在于，所述基于初始化函数初始化NLFM信号集合中各NLFM信号，包括：

基于初始化函数对各NLFM信号的频率控制点进行初始化；其中，所述NLFM信号由包括

个频率控制点的第二向量p获得；

对第k个第二向量p的第j个频率控制点进行初始化的所述初始化函数为：

其中，B_r为NLFM信号的带宽，rand为随机数。

7.根据权利要求1所述的方法，其特征在于，所述利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛包括：

利用所述算法模型，根据所述初始迭代参数计算所述初始NLFM信号的适应度和第二次迭代参数，以及基于各初始NLFM信号的适应度，对各初始NLFM信号进行选择、交叉处理和变异处理得到第二次迭代的NLFM信号；

利用所述算法模型，对所述第二次迭代的NLFM信号和所述第二次迭代参数进行迭代，得到第三次迭代的NLFM信号和第三次迭代参数；并利用所述算法模型，持续对迭代的NLFM信号和迭代参数进行迭代，直到所述增广拉格朗日遗传算法收敛。

8.一种基于增广拉格朗日遗传算法的非线性调频NLFM信号优化装置，其特征在于，所述装置包括：函数建立单元、条件确定单元、模型建立单元、初始单元和迭代单元；其中，

所述函数建立单元，用于基于贝塞尔曲线建立NLFM信号的时域函数；

所述条件确定单元，用于根据所述NLFM信号的时域函数确定所述NLFM信号的自相关函数性能中的主瓣宽度MW和峰值旁瓣比PSLR，基于所述MW和所述PSLR建立优化数学模型；其中，根据不等式约束的非线性规划模型，在满足不展宽主瓣的情况下降低旁瓣的条件时，将PLSR(p)作为优化数学模型的目标函数，将MW(p)作为约束函数，建立所述优化数学模型，或，在满足不升高旁瓣的情况下降低主瓣的条件时，将MW(p)作为优化数学模型的目标函数，将PLSR(p)作为约束函数，建立所述优化数学模型；其中，在时频坐标关系平面内给定n+1个控制点，由n+1个控制点构成所述NLFM信号的第一向量P，第二向量p为第一向量P的n+1个控制点中的前

个控制点在频率轴上的向量，所述PSLR(p)为根据第二向量p求得的NLFM信号的所述PSLR，MW(p)为根据第二向量p求得的NLFM信号的MW；

所述模型建立单元，用于基于所述优化数学模型确定增广拉格朗日遗传算法的算法模型；

所述初始单元，用于基于初始化函数初始化NLFM信号集合中各NLFM信号的频率控制点，得到相应的初始NLFM信号；

所述迭代单元，用于利用所述算法模型，对设定的初始迭代参数和所述初始NLFM信号进行持续迭代，直到所述增广拉格朗日遗传算法收敛；

所述模型建立单元，还用于基于所述优化数学模型，利用增广拉格朗日算法确定目标函数：

Θ(p,λ,s)＝g(p)-λslog(s-c(p))；

其中，λ为拉格朗日算子，s为偏移量，λ、s均为非负数；g(p)为第二向量p对应的PSLR或MW，c(p)为第二向量对应的MW或PSLR；

利用增广拉格朗日算法确定参数函数：

μ_k+1＝μ_k，s_k+1＝μ_k+1λ_k+1；

其中，μ为保证s-c(p)大于0的数值；

根据所述目标函数计算当前NLFM信号集合中各NLFM信号的适应度；根据各NLFM信号的适应度从NLFM信号集合中选择父体NLFM信号；对所述父体NLFM信号进行交叉处理、变异处理，得到下一次迭代的NLFM信号；

根据所述参数函数计算当前拉格朗日算子和偏移量对应的下一次迭代的拉格朗日算子和偏移量。

Images