抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法技术

本发明专利技术公开了一种抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法，包括：利用雷达发射端和接收端自由度，建立关于非周期发射波形集和非匹配滤波器组的最小化距离旁瓣电平的目标函数；结合发射波形的恒模特性、非匹配滤波器序列的恒能量特性和信噪比损失情况，建立无需模板的分式优化模型；基于ADMM算法框架实现变量解耦，对分式优化模型进行迭代求解，获取最优的波形集和非匹配滤波器组序列。本发明专利技术可以有效避免微弱目标淹没、虚警率升高等不良情况，可以满足更严格的距离旁瓣要求，提高雷达探测和参数估计性能。提高雷达探测和参数估计性能。提高雷达探测和参数估计性能。

【技术实现步骤摘要】
抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法


[0001]本专利技术属于雷达信号
，涉及一种抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法。

技术介绍

[0002]在多输入多输出雷达信号处理过程中，发射波形和接收滤波器的相关性对目标探测效果有着决定性作用。然而传统的一些方法仅开展发射端波形设计，并在接收端应用匹配滤波器作为接收滤波器，难以满足更加严格的距离旁瓣要求。现阶段用于MIMO雷达多波形设计的WeCAN算法以最小化积分旁瓣电平为优化目标，在多个距离区间上设置加权因子，获得局部低旁瓣效果；李风等在WeCAN算法的基础上，借助0/1加权因子构造ISL准则下的优化目标函数，并基于空间信赖域算法开展波形优化；针对固定的发射波形，Sun等人设计出对应的接收滤波器，在特定区间上实现距离旁瓣的有效抑制；Zhou等人假设当前电磁环境中的干扰情况已知并构造出间歇转发式干扰信号，以最大化新干噪比为准则，开展发射波形和非匹配滤波器的联合优化。
[0003]现有方法中为获得更好的距离旁瓣抑制效果，常设置接收端非匹配滤波器序列的长度为信号长度的2
‑
3倍，无疑会增加优化计算复杂度，消耗计算资源，不利于工程实现。

技术实现思路

[0004]本专利技术的目的在于解决现有技术中为获得更好的距离旁瓣抑制效果，常将接收端非匹配滤波器序列的长度设置为信号长度的2
‑
3倍，增加优化计算复杂度，消耗计算资源，不利于工程实现的问题，提供一种抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法。
[0005]为达到上述目的，本专利技术采用以下技术方案予以实现：
[0006]抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法，包括：
[0007]利用雷达发射端和接收端自由度，建立关于非周期发射波形集和非匹配滤波器组的最小化距离旁瓣电平的目标函数；
[0008]结合发射波形的恒模特性、非匹配滤波器序列的恒能量特性和信噪比损失情况，建立无需模板的分式优化模型；
[0009]基于ADMM算法框架实现变量解耦，对分式优化模型进行迭代求解，获取最优的波形集和非匹配滤波器组序列。
[0010]本专利技术的进一步改进在于：
[0011]进一步的，建立关于非周期发射波形集和非匹配滤波器组的最小化距离旁瓣电平的目标函数，具体为：
[0012]令和分别表示待设计的M个非周期发射波形和非匹配滤波器组成的序列集合，集合中每个序列的长度为N，即x
m
＝[x
m
(1),x
m
(2),...,x
m
(N)]T
,h
m
＝[h
m
(1),h
m
(2),...,h
m
(N)]T
；定义维度为N
×
N的选择矩阵S
d
中的第g行h列的元素如下：
[0013]S
d
(g,h)＝δ(g
‑
h
‑
d+N),d∈Ω (1)
[0014]其中，δ(
·
)表示狄拉克函数；当函数输入变量值为0时，函数值为1；当函数输入变量值不为0时，函数值为0；Ω表示距离位移d的索引，Ω＝[1,2,...,2N
‑
1]；
[0015]根据前述变量设置情况，建立目标函数模型：
[0016][0017]其中，Ω
m,p
表示根据波形序列x
m
和非匹配滤波器序列h
p
的组合情况的低旁瓣区域；α
m,p,d
表示在旁瓣水平的加权因子；目标函数中分子项表示波形x
p
经过滤波器h
m
滤波后，在距离位移k处的输出情况，分母项表示第m种发射波形和第m种非匹配滤波器输出的主瓣峰值电平；x
m,n
表示第m个波形中的第n个元素，其需要满足恒模约束，而第m个滤波器序列需要满足恒能量约束；引入发射波形和对应非匹配滤波器序列相似性约束用相似性因子β调控接收端非匹配滤波器输出的信噪比损失。
[0018]进一步的，结合发射波形的恒模特性、非匹配滤波器序列的恒能量特性和信噪比损失情况，建立无需模板的分式优化模型，具体为：
[0019]将非周期发射波形集和非匹配滤波器组多个短序列变量整合成一个长序列变量w：
[0020][0021]定义选择矩阵U
m
、V
m
以及辅助矩阵Ξ：
[0022][0023]借助选择矩阵将x
m
和h
m
重新表示为：
[0024][0025]利用公式(4)和公式(5)，将公式(2)中的优化模型转化为：
[0026][0027]引入辅助变量：
[0028][0029]将公式(6)转化为：
[0030][0031]进一步的，基于ADMM算法框架实现变量解耦，对分式优化模型进行迭代求解具体为：
[0032]构建拉格朗日增广函数：
[0033][0034]其中，ρ表示步长，{λ
m,p,d
,κ
m,p,d
,ξ}为拉格朗日乘子，表示取括号内部元素的实部；
[0035]迭代跳出条件设置为：
[0036][0037]其中，ε＞0。
[0038]进一步的，获取最优的波形集和非匹配滤波器组序列，包括：
[0039]对分式优化模型进行拆分，获取待求解的子问题；
[0040]对待求解的子问题进行处理，获取最优化的波形集序和非匹配滤波器组。
[0041]进一步的，对分式优化模型进行拆分，获取待求解的子问题，具体为：
[0042]在第(t+1)次迭代过程中，{η(t+1),y
m,p,d
(t+1),z
m,p,d
(t+1),k(t+1)}的求解步骤为：
[0043]对公式(9)进行拆分得到两个待优化的子问题：
[0044][0045][0046]其中：
[0047][0048]w(t),λ
m,p,d
(t),κ
m,p,d
(t),ξ(t)表示变量w,λ
m,p,d
,κ
m,p,d
,ξ在第t次优化得到的结果。
[0049]进一步的，对待求解的子问题进行处理，获取最优化的波形集序和非匹配滤波器组，具体为：
[0050]当η确定时，得到y
m,p,d
(t+1)和z
m,p,d
(t+1)：
[0051][0052][0053]定义阶跃函数
[0054][0055]将公式(14)至公式(16)带入到公式(11)得到
[0056][0057]公式(17)转化为关于η的单变量问题，针对此单变量问题求出最优值η(t+本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法，其特征在于，包括：利用雷达发射端和接收端自由度，建立关于非周期发射波形集和非匹配滤波器组的最小化距离旁瓣电平的目标函数；结合发射波形的恒模特性、非匹配滤波器序列的恒能量特性和信噪比损失情况，建立无需模板的分式优化模型；基于ADMM算法框架实现变量解耦，对分式优化模型进行迭代求解，获取最优的波形集和非匹配滤波器组序列。2.根据权利要求1所述的抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法，其特征在于，所述建立关于非周期发射波形集和非匹配滤波器组的最小化距离旁瓣电平的目标函数，具体为：令和分别表示待设计的M个非周期发射波形和非匹配滤波器组成的序列集合，集合中每个序列的长度为N，即x
m
＝[x
m
(1),x
m
(2),...,x
m
(N)]
T
,h
m
＝[h
m
(1),h
m
(2),...,h
m
(N)]
T
；定义维度为N
×
N的选择矩阵S
d
中的第g行h列的元素如下：S
d
(g,h)＝δ(g
‑
h
‑
d+N),d∈Ω (1)其中，δ(
·
)表示狄拉克函数；当函数输入变量值为0时，函数值为1；当函数输入变量值不为0时，函数值为0；Ω表示距离位移d的索引，Ω＝[1,2,...,2N
‑
1]；根据前述变量设置情况，建立目标函数模型：其中，Ω
m,p
表示根据波形序列x
m
和非匹配滤波器序列h
p
的组合情况的低旁瓣区域；α
m,p,d
表示在旁瓣水平的加权因子；目标函数中分子项表示波形x
p
经过滤波器h
m
滤波后，在距离位移k处的输出情况，分母项表示第m种发射波形和第m种非匹配滤波器输出的主瓣峰值电平；x
m,n
表示第m个波形中的第n个元素，其需要满足恒模约束，而第m个滤波器序列需要满足恒能量约束；引入发射波形和对应非匹配滤波器序列相似性约束用相似性因子β调控接收端非匹配滤波器输出的信噪比损失。3.根据权利要求2所述的抑制距离旁瓣的发射波形集和接收滤波器组联合设计方法，其特征在于，所述结合发射波形的恒模特性、非匹配滤波器序列的恒能量特性和信噪比损失情况，建立无需模板的分式优化模型，具体为：将非周期发射波形集和非匹配滤波器组多个短序列变量整合成一个长序列变量w：
定义选择矩阵U
m
、V
m
以及辅助矩阵Ξ：U
m
＝[0
N
×
(m
‑
1)N
,I
N
×
N
,0
N
×
(2M
‑
m)N
],m＝1,2,...,MV
m
＝[0
N
×
(m
‑
1+M)N
,I
N
×

【专利技术属性】
技术研发人员：梁军利，汪涛，王冠阳，宋珂满，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：