CN109738883B - 栅瓣抑制的宽带多阶频率步进线性调频信号设计方法 - Google Patents
栅瓣抑制的宽带多阶频率步进线性调频信号设计方法 Download PDFInfo
- Publication number
- CN109738883B CN109738883B CN201811535679.9A CN201811535679A CN109738883B CN 109738883 B CN109738883 B CN 109738883B CN 201811535679 A CN201811535679 A CN 201811535679A CN 109738883 B CN109738883 B CN 109738883B
- Authority
- CN
- China
- Prior art keywords
- frequency
- order
- stepping
- modulation signal
- waveform
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种栅瓣抑制的宽带多阶频率步进线性调频信号设计方法,对频率步进线性调频信号波形进行周期延拓,同时整体引入一个新的步进频率,得到多阶频率步进调频信号,然后推导得出其距离自相关函数,利用零极点对消的方法设计波形参数,实现对其距离自相关函数栅瓣的抵消,从而消除模糊峰。本发明不仅能够有效地抑制雷达波形的距离栅瓣,而且可以极大地提高信号的整体带宽,满足雷达高距离分辨力的要求。
Description
技术领域
本发明属于雷达波形技术领域,特别是一种栅瓣抑制的宽带多阶频率步进线性调频(Multi-Stage Stepped-Frequency LFM,MS-SFLFM)信号设计方法。
背景技术
宽带线性调频(Linear Frequency Modulated,LFM)信号是高距离分辨力成像雷达常用的发射波形,由于该波形通过在一个脉冲持续时间内频率的连续线性变化来获取较大的信号带宽,对于后端的模数转换芯片(Analog Digital Converter,ADC)以及信号处理速率压力较大。为了降低对ADC芯片和信号处理的要求,一般使用频率步进信号合成宽带信号,频率步进信号是一串脉间载频线性跳变的脉冲串,能够在不增加信号瞬时带宽的情况下,通过相参合成以获得高的距离分辨力,降低了对采样与处理的硬件要求。
频率步进信号主要分为频率步进点频信号和频率步进线性调频(Stepped-frequency LFM,SFLFM)信号。频率步进点频信号是在点频脉冲串信号上添加脉间步进频率形成的,能够获得较大的整体带宽,但是其距离自相关函数会引入严重的距离栅瓣。SFLFM信号则是在前者基础上用LFM脉冲替换固定频率脉冲形成的。当脉间步进频率不小于脉冲内调频带宽时,SFLFM信号就是雷达系统中常用的离散频率SFLFM信号(DiscreteFrequency SFLFM,DF-SFLFM),这种信号由于子脉冲频谱的不重叠特性能够降低距离自相关函数的旁瓣,也能获得很大的整体带宽,但同样地也会引入较严重的距离栅瓣;当脉间步进频率小于脉冲内调频带宽时,通过合理设计步进频率,就能利用距离自相关函数中sinc函数的性质去消除栅瓣,但这种方法会使信号的整体带宽降低。
发明内容
本发明的目的在于提供一种栅瓣抑制的宽带多阶频率步进线性调频信号设计的方法。
实现本发明的技术解决方案具体为:一种栅瓣抑制的宽带多阶频率步进线性调频信号设计方法,包括以下步骤:
步骤1、建立宽带多阶频率步进线性调频信号波形数学模型|uMS_SF(t)|;
步骤2、根据宽带多阶频率步进线性调频信号波形数学模型得到宽带多阶频率步进线性调频信号波形中心带条模糊函数表达式用ξ=0的切割面对宽带多阶频率步进线性调频信号波形模糊图进行切割,得到宽带多阶频率步进线性调频信号的距离自相关函数|RMS_SF(τ)|;
步骤3、根据宽带多阶频率步进线性调频信号的距离自相关函数特征,利用零极点对消的方法,使一阶步进频率引入的梳齿包络和二阶步进频率引入的梳齿包络互相用零点消除对方的不重合的极大值点,用线性调频子脉冲引入的近似sinc包络的零点消除两个梳齿包络重合的极大值点,由此确定抑制栅瓣的宽带多阶频率步进线性调频信号波形参数的约束条件;
步骤4、根据宽带多阶频率步进线性调频信号波形数学模型建立整体带宽表达式,选择一组满足栅瓣抑制的波形约束条件的参数,得到宽带多阶频率步进线性调频的整体带宽,完成宽带多阶频率步进线性调频信号设计。
优选地,步骤1中建立的宽带多阶频率步进线性调频信号波形数学模型为:
其中:
式中,u1(t)为LFM子脉冲复包络,tr为脉冲重复周期,tp为脉冲宽度,k为脉宽内的调频斜率,若单个子脉冲带宽为B,则k=B/tp,N2表示频率步进线性调频脉冲组个数,Δf2表示频率步进线性调频脉冲组之间的步进频率,N1表示频率步进线性调频脉冲组内子脉冲间的个数,Δf1表示子脉冲间步进频率,宽带多阶频率步进线性调频波形中共包含N1·N2个子脉冲。
优选地,步骤2中得到的MS-SFLFM信号的距离自相关函数具体为:
|RMS_SF(τ)|=|RMS_SF1(τ)||RMS_SF2(τ)||RMS_SF3(τ)|
其中:
优选地,步骤3中利用零极点对消的方法确定抑制栅瓣的宽带多阶频率步进线性调频信号波形参数的约束条件的具体步骤为:
步骤3-1、确定消除一阶步进频率引入的梳齿包络|RMS_SF2(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的不重合极大值点的条件,具体为:
|RMS_SF2(τ)|的零点位置为极值点位置为|RMS_SF3(τ)|的零点位置为极值点位置为其中k1,k2,k3,k4均为正整数,使|RMS_SF2(τ)|和|RMS_SF3(τ)|的零点位置分别与|RMS_SF3(τ)|和|RMS_SF2(τ)|的极值点位置重合,则有:
令
其中M1,M2表示互质且大于1的正整数,得到:
为了使得任意一个极值点都存在零点相抵消,则N1,N2满足:
其中k5,k6均为正整数;
步骤3-2、确定消除一阶步进频率引入的梳齿包络|RMS_SF(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的的重合极大值点的条件,具体为:
|RMS_SF(τ)|和|RMS_SF3(τ)|两者重合的极值点位置为:
其中,k7为正整数,第q和第r重合极值点位置分别对应|RMS_SF(τ)|的第m个和第n个零点(q>r,m>n),则有:
令q=2,r=1,则有:
步骤3-3、根据步骤3-1和步骤3-2得到宽带多阶频率步进线性调频信号为消除其距离自相关函数的栅瓣,参数设计需要满足的总约束条件为:
优选地,步骤4中满足Δf1=B条件下,宽带多阶频率步进线性调频的整体带宽为BMS_SF1具体为:
本发明与现有技术相比,其显著优点为:(1)本发明通过多阶频率步进信号合成宽带信号,有效降低雷达系统对后级ADC芯片和信号处理的要求;(2)本发明能有效抑制雷达波形距离自相关函数的栅瓣,从而消除模糊峰;(3)本发明能有效增大雷达波形的整体带宽,满足雷达高距离分辨力的要求。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明的多阶频率步进线性调频信号设计的流程图。
图2是MS-SFLFM信号波形的频率-时间关系图。
图3是MS-SFLFM信号距离自相关函数各部分包络波形图。
图4是MS-SFLFM信号距离自相关函数波形图。
具体实施方式
一种栅瓣抑制的宽带多阶频率步进线性调频信号设计方法,包括以下步骤:
步骤1、建立MS-SFLFM信号波形数学模型,如图2所示,MS-SFLFM信号由N2个SFLFM脉冲组组成,SFLFM脉冲组之间的步进频率为Δf2,单个SFLFM脉冲组则由N1个脉间步进频率为Δf1的LFM脉冲组成,MS-SFLFM波形中共由N1·N2个脉冲组成。因此MS-SFLFM信号波形数学模型即复包络归一化表达式|uMS_SF(t)|为:
其中:
u1(t)为LFM子脉冲复包络,tr为脉冲重复周期,tp为脉冲宽度,k为脉宽内的调频斜率,若单个子脉冲带宽为B,则k=B/tp。
令
|RMS_SF(τ)|=|RMS_SF(τ)||RMS_SF(τ)||RMS_SF(τ)|
其中:
步骤3、确定抑制栅瓣的MS-SFLFM信号波形参数。根据MS-SFLFM信号的距离自相关函数特征,利用零极点对消的方法,使一阶步进频率引入的梳齿包络|RMS_SF2(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的零点分别对其极大值点相消,而对于两者极大值重合的点,则利用近似sinc包络的|RMS_SF1(τ)|的零点对其进行消除,sinc包络是信号方向基础的波形,具体消除方法为:
步骤3-1、确定消除一阶步进频率引入的梳齿包络|RMS_SF2(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的不重合极大值点的条件,具体为:
|RMS_SF2(τ)|的零点位置为极值点位置为同样地,|RMS_SF3(τ)|的零点位置为极值点位置为其中k1,k2,k3,k4均为正整数。使|RMS_SF2(τ)|和|RMS_SF3(τ)|的零点位置分别与|RMS_SF3(τ)|和|RMS_SF2(τ)|的极值点位置重合,则有
为了更好地描述N1,Δf1,N2,Δf2之间的关系,令
其中M1,M2表示互质且大于1的正整数,可以得到
为了使得任意一个极值点都存在零点相抵消,则N1,N2需满足
其中k5,k6均为正整数。
步骤3-2、确定消除一阶步进频率引入的梳齿包络|RMS_SF2(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的的重合极大值点的条件,具体为:
|RMS_SF2(τ)|和|RMS_SF3(τ)|两者重合的极值点位置为
其中k7为正整数。假设第q和第r重合极值点位置分别对应|RMS_SF1(τ)|的第m个和第n个零点(q>r,m>n),则有
一般要求消除前两个重合的栅瓣,而这种情况下往往能消除所有重合的栅瓣,即令q=2,r=1,则有
步骤3-3、综合步骤3-1和步骤3-2得到MS-SFLFM信号为消除其距离自相关函数的栅瓣,参数设计需要满足的总约束条件为:
步骤4、上述设计的MS-SFLFM信号,设其整体带宽为BMS_SF,则
BMS_SF=(N2-1)Δf2+(N1-1)Δf1+B
根据步骤3得到的波形约束条件,选取一组参数m,n,M1,M2,使得Δf1=B,得到MS-SFLFM的整体带宽为:
实施例
如图1所示,本发明抑制栅瓣的宽带多阶频率步进线性调频信号设计的方法,步骤如下:
步骤1、为MS-SFLFM信号波形建立数学模型,MS-SFLFM信号由N2个SFLFM脉冲组组成,SFLFM脉冲组之间的步进频率为Δf2,单个SFLFM脉冲组则由N1个脉间步进频率为Δf1的LFM脉冲组成,MS-SFLFM波形中共由N1·N2个脉冲组成。由此得到MS-SFLFM信号的复包络归一化表达式|uMS_SF(t)|。
步骤3、确定抑制栅瓣的MS-SFLFM信号波形参数。根据MS-SFLFM信号的距离自相关函数特征,利用零极点对消的方法,使梳齿包络|RMS_SF2(τ)|和|RMS_SF3(τ)|的零点分别对其极大值点相消,而对于两者极大值重合的点,则利用近似sinc包络的|RMS_SF1(τ)|的零点对其进行消除。得到MS-SFLFM信号为消除其距离自相关函数的栅瓣,参数设计需要满足的总约束条件为
根据约束条件,选择一组系数n=3,m=4,M1=4,M2=3,得到满足约束条件的波形参数N1=3,N2=4,tpΔf1=tpB=16,tpΔf2=12。如图3是MS-SFLFM信号距离自相关函数三部分包络在一个tp内的波形图。
步骤4、上述设计的MS-SFLFM信号,设其整体带宽为BMS_SF,选取信号的脉宽tp=5*10-6s,则得到MS-SFLFM的整体带宽为BMS_SF=16.8MHz。
通过本发明设计波形参数N1=3,N2=4,tpΔf1=tpB=16,tpΔf2=12,其距离自相关函数如图4所示,栅瓣得到了很好的抑制,其整体带宽为16.8MHz。
由上可知,本发明设计的MS-SFLFM在满足波形约束条件下,能有效地抑制距离栅瓣,同时具有较大的整体带宽,从而提高了雷达的距离分辨力。
Claims (5)
1.一种栅瓣抑制的宽带多阶频率步进线性调频信号设计方法,其特征在于,包括以下步骤:
步骤1、建立宽带多阶频率步进线性调频信号波形数学模型|uMS_SF(t)|;
步骤2、根据宽带多阶频率步进线性调频信号波形数学模型得到宽带多阶频率步进线性调频信号波形中心带条模糊函数表达式用ξ=0的切割面对宽带多阶频率步进线性调频信号波形模糊图进行切割,得到宽带多阶频率步进线性调频信号的距离自相关函数|RMS_SF(τ)|;
步骤3、根据宽带多阶频率步进线性调频信号的距离自相关函数特征,利用零极点对消的方法,使一阶步进频率引入的梳齿包络和二阶步进频率引入的梳齿包络互相用零点消除对方的不重合的极大值点,用线性调频子脉冲引入的近似sinc包络的零点消除两个梳齿包络重合的极大值点,由此确定抑制栅瓣的宽带多阶频率步进线性调频信号波形参数的约束条件;
步骤4、根据宽带多阶频率步进线性调频信号波形数学模型建立整体带宽表达式,选择一组满足栅瓣抑制的波形约束条件的参数,得到宽带多阶频率步进线性调频的整体带宽,完成宽带多阶频率步进线性调频信号设计。
4.根据权利要求1所述的栅瓣抑制的宽带多阶频率步进线性调频信号设计方法,其特征在于,步骤3中利用零极点对消的方法确定抑制栅瓣的宽带多阶频率步进线性调频信号波形参数的约束条件的具体步骤为:
步骤3-1、确定消除一阶步进频率引入的梳齿包络|RMS_SF2(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的不重合极大值点的条件,具体为:
|RMS_SF2(τ)|的零点位置为极值点位置为|RMS_SF3(τ)|的零点位置为极值点位置为其中k1,k2,k3,k4均为正整数,使|RMS_SF2(τ)|和|RMS_SF3(τ)|的零点位置分别与|RMS_SF3(τ)|和|RMS_SF2(τ)|的极值点位置重合,则有:
令
其中M1,M2表示互质且大于1的正整数,得到:
为了使得任意一个极值点都存在零点相抵消,则N1,N2满足:
其中k5,k6均为正整数;
步骤3-2、确定消除一阶步进频率引入的梳齿包络|RMS_SF2(τ)|和二阶步进频率引入的梳齿包络|RMS_SF3(τ)|的重合极大值点的条件,具体为:
|RMS_SF2(τ)|和|RMS_SF3(τ)|两者重合的极值点位置为:
其中,k7为正整数,第q和第r重合极值点位置分别对应|RMS_SF1(τ)|的第m个和第n个零点, q>r,m>n,则有:
令q=2,r=1,则有:
步骤3-3、根据步骤3-1和步骤3-2得到宽带多阶频率步进线性调频信号为消除其距离自相关函数的栅瓣,参数设计需要满足的总约束条件为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811535679.9A CN109738883B (zh) | 2018-12-14 | 2018-12-14 | 栅瓣抑制的宽带多阶频率步进线性调频信号设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811535679.9A CN109738883B (zh) | 2018-12-14 | 2018-12-14 | 栅瓣抑制的宽带多阶频率步进线性调频信号设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109738883A CN109738883A (zh) | 2019-05-10 |
CN109738883B true CN109738883B (zh) | 2022-11-01 |
Family
ID=66359449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811535679.9A Active CN109738883B (zh) | 2018-12-14 | 2018-12-14 | 栅瓣抑制的宽带多阶频率步进线性调频信号设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109738883B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5044358B2 (ja) * | 2007-10-22 | 2012-10-10 | 株式会社東芝 | レーダ装置 |
US7843382B2 (en) * | 2008-12-15 | 2010-11-30 | Adly T. Fam | Mismatched filter |
CN102608594B (zh) * | 2012-03-20 | 2013-07-31 | 北京理工大学 | 一种调频步进信号合成高分辨像栅瓣抑制方法 |
CN103138799B (zh) * | 2013-03-21 | 2014-12-17 | 哈尔滨工业大学 | 一种低旁瓣随机跳频脉冲信号的调制方法 |
CN103323820B (zh) * | 2013-06-19 | 2015-07-15 | 北京理工大学 | 一种抑制频率步进sar一维距离图像栅瓣的方法 |
CN104765033B (zh) * | 2015-04-28 | 2017-03-08 | 西北工业大学 | 利用互相关函数抑制步进频率成像中距离旁瓣的方法 |
EP3109939B1 (en) * | 2015-06-26 | 2024-01-03 | HENSOLDT Sensors GmbH | Dual-band phased array antenna with built-in grating lobe mitigation |
-
2018
- 2018-12-14 CN CN201811535679.9A patent/CN109738883B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109738883A (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105259534A (zh) | 一种非线性调频信号的生成方法和装置 | |
US9645228B1 (en) | Shaping the spectrum of random-phase radar waveforms | |
US4379295A (en) | Low sidelobe pulse compressor | |
CN112748403B (zh) | 一种基于频谱调制捷变的非线性调频脉冲串波形设计方法 | |
DE102008033988A1 (de) | Rampenlinearisierung bei FMCW-Radar mit digitaler Abwärtswandlung eines abgetasteten VCO-Signals | |
CN109613507B (zh) | 一种针对高阶机动目标雷达回波的检测方法 | |
CN110888131B (zh) | 一种基于频率步进sar雷达的多通道相位补偿方法 | |
CN109738883B (zh) | 栅瓣抑制的宽带多阶频率步进线性调频信号设计方法 | |
WO2024179616A1 (zh) | 一种高自由度参数化调频编码波形的设计方法 | |
US8159280B2 (en) | Noise generator | |
US9983295B2 (en) | Direct sampling of received signals in radar | |
US5757848A (en) | Method and apparatus for a decimating digital PN correlator | |
EP0782259A2 (en) | Digital center line filter | |
ITRM20080511A1 (it) | Procedimento di elaborazione di un segnale di eco radar, prodotto da un bersaglio, per compensare effetti di degradazione di detto segnale dovuti al moto di detto bersaglio. | |
RU2439609C2 (ru) | Многоканальное корреляционно-фильтровое приемное устройство с селекцией движущихся целей | |
CN107769815A (zh) | 线性调频近程探测系统噪声调幅干扰抑制方法 | |
CA2152854C (en) | Ramp weighted correction with oversampling | |
RU2319170C1 (ru) | Цифровое многоканальное корреляционно-фильтровое приемное устройство с селекцией движущихся целей | |
CN108169725A (zh) | 一种用于距离副瓣抑制的自适应恒虚警方法 | |
RU2685972C1 (ru) | Способ и устройство фильтрации частотно-модулированных сигналов | |
Doerry et al. | Random-phase radar waveforms with shaped spectrum | |
RU2503971C1 (ru) | Способ подавления боковых лепестков автокорреляционной функции широкополосного сигнала | |
Liao et al. | A low side-lobe waveform design method for random hopping frequency signal | |
Galati et al. | Optimal chirp waveforms for long/medium range surveillance radar | |
RU2628475C1 (ru) | Устройство подавления боковых лепестков при импульсном сжатии симметрично усеченных многофазных кодов (варианты) |
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 |