CN106597389A - 步进频信号固定主瓣宽度的低旁瓣加权系数设计方法 - Google Patents

步进频信号固定主瓣宽度的低旁瓣加权系数设计方法 Download PDF

Info

Publication number
CN106597389A
CN106597389A CN201611014141.4A CN201611014141A CN106597389A CN 106597389 A CN106597389 A CN 106597389A CN 201611014141 A CN201611014141 A CN 201611014141A CN 106597389 A CN106597389 A CN 106597389A
Authority
CN
China
Prior art keywords
frequency signal
stepped frequency
main lobe
sidelobe
weight coefficient
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
CN201611014141.4A
Other languages
English (en)
Other versions
CN106597389B (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.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research Institute
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 Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN201611014141.4A priority Critical patent/CN106597389B/zh
Publication of CN106597389A publication Critical patent/CN106597389A/zh
Application granted granted Critical
Publication of CN106597389B publication Critical patent/CN106597389B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

Abstract

本发明涉及一种步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,包含:S1、设定步进频信号的约束的主瓣宽度;S2、利用距离自相关函数,将步进频信号的最大旁瓣最小化约束条件转化为凸优化问题;S3、利用凸优化进行加权系数的求解;S4、对步进频信号的回波基带数据进行频域加窗,并进行逆傅里叶变换,求得固定主瓣宽度后的低旁瓣一维高分辨距离像。本发明改善了步进频信号的旁瓣抑制问题,同时可适应于固定分辨力情况下的步进频信号的旁瓣抑制问题;并且可改善强杂波干扰环境中降低杂波旁瓣对目标信号的影响。

Description

步进频信号固定主瓣宽度的低旁瓣加权系数设计方法
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种步进频信号固定主瓣宽度的低旁瓣加权系数的设计方法。
背景技术
现有技术中,对于高分辨率波形信号,如线性调频信号、多相码信号、步进频信号等,其固有缺陷都是通过匹配滤波器之后的输出有较高的距离旁瓣。在多目标环境中,强目标回波的输出旁瓣很可能混淆或者掩盖弱目标回波的主峰,从而影响对弱目标的探测和分辨。为了提高分辨目标的能力,因此必须采用旁瓣抑制或加权技术。
在对信号进行匹配滤波或者步进频信号合成高分辨率距离像的过程中,如果旁瓣过高会影响对目标的检测。传统的旁瓣抑制方法通常采用时域或频域加窗方法,如Hamming窗函数或Kaiser窗函数等。但是这些加窗方法侧重于得到高的主旁瓣比,加窗后的旁瓣得到抑制,但输出主瓣有不同程度的展宽,并且对主瓣展宽的程度不可控,从而导致分辨力降低,且降低量未知。在目标的要害部位定位、目标精确成像、多目标分辨等应用场合中,对距离分辨力要求苛刻。
因此,为了兼顾旁瓣抑制和高距离分辨力需求,本发明提出一种步进频信号固定主瓣宽度的低旁瓣加权系数设计方法。
发明内容
本发明的目的在于提供一种步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,改善了步进频信号的旁瓣抑制问题,同时可适应于固定分辨力情况下的步进频信号的旁瓣抑制问题;并且可改善强杂波干扰环境中降低杂波旁瓣对目标信号的影响。
为了达到上述目的,本发明提供一种步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,包含以下步骤:
S1、设定步进频信号的约束的主瓣宽度;
S2、利用距离自相关函数,将步进频信号的最大旁瓣最小化约束条件转化为凸优化问题;
S3、利用凸优化进行加权系数的求解;
S4、对步进频信号的回波基带数据进行频域加窗,并进行逆傅里叶变换,求得固定主瓣宽度后的低旁瓣一维高分辨距离像。
所述的S1中,步进频信号的一个波形周期由N个脉冲组成,每个脉冲宽度为Tp,并进行N点频域采样,则步进频信号对应的时间分辨率为Δt=Tp/N,主瓣宽度为TD=hΔt=hTp/N,其中,h=2,3,…。
所述的S2中,具体包含以下步骤:
步进频信号对点目标的回波延时为τ,步进量为Δf,第n个脉冲步进量为fn=(n-1)Δf,n=0,1,2,…,N-1,且频点对应的加权系数为ω,则加权的距离自相关函数为:
加权系数ω的向量表示为:
ω=[ω0 ω1 … ωN-1]T
步进频信号一个波形周期的基带信号表示为:
a(τ)=[exp(j2πf0τ),exp(j2πf1τ),…,exp(j2πfN-1τ)]T
根据加权的距离自相关函数给出的关系R(τ),可得步进频信号的主瓣以外的旁瓣部分为|τ|>TD的部分,将其最大值最小化的过程转化为凸优化问题,表示为:
其中,||ω||1为1范数,即所有系数绝对值求和。
所述的S3中,求解得到的加权系数可按一固定倍数进行缩放。
所述的S4中,具体包含以下步骤:
步进频信号的回波基带数据为:
s=[s(0),s(1),…,s(N-1)]T
则步进频信号在固定主瓣宽度后的低旁瓣一维距离像为:
y=IFFT{ωTs}。
综上所述,本发明提供的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,改善了步进频信号的旁瓣抑制问题,同时可适应于固定分辨力情况下的步进频信号的旁瓣抑制问题;并且可改善强杂波干扰环境中降低杂波旁瓣对目标信号的影响。本发明适用于一般步进频信号的旁瓣抑制、固定分辨力情况下的步进频信号的加窗处理、以及强杂波干扰环境中降低杂波旁瓣对目标信号的影响等。
附图说明
图1为本发明中步进频信号固定主瓣宽度的低旁瓣加权系数设计方法的流程图;
图2为本发明中步进频信号未加窗的距离自相关输出示意图;
图3为本发明中步进频信号未加窗的距离像示意图;
图4为本发明中加权系数求解结果示意图;
图5为本发明中步进频信号进行系数加权后的距离自相关输出示意图;
图6为本发明中步进频信号进行系数加权后的距离像示意图。
具体实施方式
以下结合图1~图6,详细说明本发明的优选实施例。
如图1所示,为本发明所述的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,包含以下步骤:
S1、设定步进频信号的约束的主瓣宽度;
S2、利用距离自相关函数,将步进频信号的最大旁瓣最小化约束条件转化为凸优化问题;
S3、利用凸优化进行加权系数的求解;
S4、对步进频信号的回波基带数据进行频域加窗,并进行逆傅里叶变换(IFFT),求得固定主瓣宽度后的低旁瓣一维高分辨距离像。
所述的S1中,步进频信号的一个波形周期由N个脉冲组成,每个脉冲宽度为Tp,并进行N点频域采样,则每个频点表示原始距离分辨率的1/N,即步进频信号对应的时间分辨率为Δt=Tp/N;因此主瓣宽度TD=hΔt=hTp/N,其中,h=2,3,…。
所述的S2中,具体包含以下步骤:
步进频信号对点目标的回波延时为τ,步进量为Δf,第n个脉冲步进量为fn=(n-1)Δf,n=0,1,2,…,N-1,且频点对应的加权系数为ω,则加权的距离自相关函数为:
加权系数ω的向量表示为:
ω=[ω0 ω1 … ωN-1]T
步进频信号一个波形周期的基带信号表示为:
a(τ)=[exp(j2πf0τ),exp(j2πf1τ),…,exp(j2πfN-1τ)]T
根据加权的距离自相关函数给出的关系R(τ),可得步进频信号的主瓣以外的旁瓣部分为|τ|>TD的部分,将其最大值最小化的过程转化为凸优化问题,表示为:
其中,||ω||1为1范数,即所有系数绝对值求和。
所述的S3中,实际使用中,加权系数的求解结果ω可按某一固定倍数进行缩放,以适应数字信号处理的要求。本实施例中,可以借助MATLAB软件中的凸优化(CVX)工具箱进行加权系数的求解。
所述的S4中,具体包含以下步骤:
步进频信号的回波基带数据为:
s=[s(0),s(1),…,s(N-1)]T
则步进频信号在固定主瓣宽度后的低旁瓣一维距离像为:
y=IFFT{ωTs}。
以下通过一个具体实施例,详细说明本发明。
本实施实例中,设定模拟雷达回波信号的参数为:载频为35GHz,步进频信号一个波形周期内的脉冲数为256,步进量为5MHz,脉冲宽度为0.1μs;设置3个点目标距离分别为10m、15m和20m,且目标散射截面积都为1;并且确定步进频信号的主瓣宽度为TD=3Tp/N。
图2和图3分别给出了步进频信号未加窗进行旁瓣抑制的距离自相关函数及3个点目标的距离像结果。从图2中可以看出步进频信号在未加窗优化时,主旁瓣比为13dB。
采用本发明方法的结果如图4~图6所示。其中,图4给出了固定主瓣宽度TD=3Tp/N后,并利用凸优化所求得的每个频点对应的加权系数结果;图5给出了利用本发明方法进行系数加权后的步进频信号的距离自相关函数,可以看出固定主瓣宽度TD=3Tp/N后,得到的主旁瓣比为52dB;从图6的距离像结果可以看出,利用本发明方法后,主瓣宽度固定,即分辨力不变,同时旁瓣能够得到明显抑制。
综上所述,本发明所提供的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,通过固定主瓣宽度,并利用距离自相关函数对最大旁瓣最小化准则进行优化,并将该优化过程转化为凸优化过程,利用凸优化快速求解加权系数,从而使得步进频信号在固定主瓣宽度的情况下,能最大限度的降低旁瓣。
本发明所提供的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,能够解决步进频信号在旁瓣抑制过程中,加窗导致的主瓣宽度不能定量控制的问题。本发明改善了步进频信号的旁瓣抑制问题,同时可适应于固定分辨力情况下的步进频信号的旁瓣抑制问题;并且可改善强杂波干扰环境中降低杂波旁瓣对目标信号的影响。本发明适用于一般步进频信号的旁瓣抑制、固定分辨力情况下的步进频信号的加窗处理、以及强杂波干扰环境中降低杂波旁瓣对目标信号的影响等。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。

Claims (5)

1.一种步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,其特征在于,包含以下步骤:
S1、设定步进频信号的约束的主瓣宽度;
S2、利用距离自相关函数,将步进频信号的最大旁瓣最小化约束条件转化为凸优化问题;
S3、利用凸优化进行加权系数的求解;
S4、对步进频信号的回波基带数据进行频域加窗,并进行逆傅里叶变换,求得固定主瓣宽度后的低旁瓣一维高分辨距离像。
2.如权利要求1所述的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,其特征在于,所述的S1中,步进频信号的一个波形周期由N个脉冲组成,每个脉冲宽度为Tp,并进行N点频域采样,则步进频信号对应的时间分辨率为Δt=Tp/N,主瓣宽度为TD=hΔt=hTp/N,其中,h=2,3,…。
3.如权利要求2所述的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,其特征在于,所述的S2中,具体包含以下步骤:
步进频信号对点目标的回波延时为τ,步进量为Δf,第n个脉冲步进量为fn=(n-1)Δf,n=0,1,2,…,N-1,且频点对应的加权系数为ω,则加权的距离自相关函数为:
加权系数ω的向量表示为:
ω=[ω0 ω1 … ωN-1]T
步进频信号一个波形周期的基带信号表示为:
a(τ)=[exp(j2πf0τ),exp(j2πf1τ),…,exp(j2πfN-1τ)]T
根据加权的距离自相关函数给出的关系R(τ),可得步进频信号的主瓣以外的旁瓣部分为|τ|>TD的部分,将其最大值最小化的过程转化为凸优化问题,表示为:
其中,||ω||1为1范数,即所有系数绝对值求和。
4.如权利要求3所述的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,其特征在于,所述的S3中,求解得到的加权系数可按一固定倍数进行缩放。
5.如权利要求4所述的步进频信号固定主瓣宽度的低旁瓣加权系数设计方法,其特征在于,所述的S4中,具体包含以下步骤:
步进频信号的回波基带数据为:
s=[s(0),s(1),…,s(N-1)]T
则步进频信号在固定主瓣宽度后的低旁瓣一维距离像为:
y=IFFT{ωTs}。
CN201611014141.4A 2016-11-17 2016-11-17 步进频信号固定主瓣宽度的低旁瓣加权系数设计方法 Active CN106597389B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611014141.4A CN106597389B (zh) 2016-11-17 2016-11-17 步进频信号固定主瓣宽度的低旁瓣加权系数设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611014141.4A CN106597389B (zh) 2016-11-17 2016-11-17 步进频信号固定主瓣宽度的低旁瓣加权系数设计方法

Publications (2)

Publication Number Publication Date
CN106597389A true CN106597389A (zh) 2017-04-26
CN106597389B CN106597389B (zh) 2019-07-19

Family

ID=58592118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611014141.4A Active CN106597389B (zh) 2016-11-17 2016-11-17 步进频信号固定主瓣宽度的低旁瓣加权系数设计方法

Country Status (1)

Country Link
CN (1) CN106597389B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108874743A (zh) * 2018-04-16 2018-11-23 北京理工大学 一种控制频谱主副瓣特性的加窗方法
CN110618411A (zh) * 2019-10-23 2019-12-27 电子科技大学 一种机载雷达杂波实时信号产生方法
CN111628954A (zh) * 2020-05-22 2020-09-04 锐迪科微电子(上海)有限公司 基于sc-fdma的数据处理方法及装置、存储介质、终端
CN113589250A (zh) * 2021-07-29 2021-11-02 上海无线电设备研究所 一种稀疏学习的相参捷变频雷达距离高分辨处理方法
CN113691479A (zh) * 2021-07-15 2021-11-23 中国电子科技集团公司第二十九研究所 一种信号消除旁瓣影响的方法、装置及可读存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104198996A (zh) * 2014-08-30 2014-12-10 西安电子科技大学 高码率低距离分辨率的低旁瓣相位编码信号设计方法
CN104199052A (zh) * 2014-09-22 2014-12-10 哈尔滨工程大学 一种基于范数约束的波束旁瓣抑制方法
CN104635230A (zh) * 2015-03-05 2015-05-20 北京航空航天大学 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法
CN105044685A (zh) * 2015-06-19 2015-11-11 南京理工大学 线性调频脉冲信号的超低旁瓣脉压器设计方法
CN105842702A (zh) * 2016-05-13 2016-08-10 南京航空航天大学 多波束成像声呐的旁瓣抑制方法、阵列稀疏方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104198996A (zh) * 2014-08-30 2014-12-10 西安电子科技大学 高码率低距离分辨率的低旁瓣相位编码信号设计方法
CN104199052A (zh) * 2014-09-22 2014-12-10 哈尔滨工程大学 一种基于范数约束的波束旁瓣抑制方法
CN104635230A (zh) * 2015-03-05 2015-05-20 北京航空航天大学 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法
CN105044685A (zh) * 2015-06-19 2015-11-11 南京理工大学 线性调频脉冲信号的超低旁瓣脉压器设计方法
CN105842702A (zh) * 2016-05-13 2016-08-10 南京航空航天大学 多波束成像声呐的旁瓣抑制方法、阵列稀疏方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵德华: "调频脉冲信号旁瓣抑制技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108874743A (zh) * 2018-04-16 2018-11-23 北京理工大学 一种控制频谱主副瓣特性的加窗方法
CN110618411A (zh) * 2019-10-23 2019-12-27 电子科技大学 一种机载雷达杂波实时信号产生方法
CN110618411B (zh) * 2019-10-23 2022-11-04 电子科技大学 一种机载雷达杂波实时信号产生方法
CN111628954A (zh) * 2020-05-22 2020-09-04 锐迪科微电子(上海)有限公司 基于sc-fdma的数据处理方法及装置、存储介质、终端
CN113691479A (zh) * 2021-07-15 2021-11-23 中国电子科技集团公司第二十九研究所 一种信号消除旁瓣影响的方法、装置及可读存储介质
CN113589250A (zh) * 2021-07-29 2021-11-02 上海无线电设备研究所 一种稀疏学习的相参捷变频雷达距离高分辨处理方法
CN113589250B (zh) * 2021-07-29 2024-02-23 上海无线电设备研究所 一种稀疏学习的相参捷变频雷达距离高分辨处理方法

Also Published As

Publication number Publication date
CN106597389B (zh) 2019-07-19

Similar Documents

Publication Publication Date Title
CN106597389A (zh) 步进频信号固定主瓣宽度的低旁瓣加权系数设计方法
CN103616679B (zh) 基于差波束调制和波形分析的pd雷达测距测角方法
CN103364765B (zh) 基于模拟解调/调制的数字射频存储器
CN114594428B (zh) 基于脉间频率编码脉内线性调频的间歇采样干扰抑制方法
CN109257051A (zh) 一种数字信道化射频存储系统
CN104991239B (zh) 一种基于脉冲压缩雷达的距离旁瓣抑制方法
JP2014132250A (ja) レーダ装置
CN111610503B (zh) 基于改进的lvd的线性调频信号参数估计方法
CN107346017A (zh) 基于频域匹配滤波的脉冲压缩雷达密集目标模拟方法
CN104076369B (zh) 基于自适应门限判决的频域抗干扰方法及装置
CN110109089B (zh) 一种线性调频连续波探测系统距离测量精度的改进方法
JP3740434B2 (ja) パルス受信分析装置及びパルス受信分析方法
CN107167777B (zh) 锯齿波线性调频信号参数提取方法
CN108020834A (zh) 基于改进edpca的运动目标检测方法、装置及电子设备
CN114545342A (zh) 利用多通道侦察接收机的雷达脉冲信号参数测量方法
CN116113843A (zh) 调频连续波雷达系统中的多普勒处理
CN107255809B (zh) 一种基于宽带聚焦矩阵的阻塞阵波束形成方法
CN112630740B (zh) 密集多目标的分辨方法及装置
CN107656268A (zh) 一种基于漏泄同轴电缆传感器的多目标定位系统
CN107769815A (zh) 线性调频近程探测系统噪声调幅干扰抑制方法
CN111427035A (zh) 一种基于子带微多普勒差异的雷达实测数据集扩充方法
CN113325382A (zh) 基于全局最优化相位调制的ddma mimo雷达速度解模糊方法
Lei et al. A TD-CF preprocessing method of FMCW radar for Dynamic Hand Gesture Recognition
CN111999725B (zh) 窄带信号引导下的宽带多项式相位信号去斜方法及装置
WO2021210480A1 (ja) 通信装置、通信方法

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