CN103630882A - 一种利用多斜率调制实现分布式雷达波形的方法 - Google Patents
一种利用多斜率调制实现分布式雷达波形的方法 Download PDFInfo
- Publication number
- CN103630882A CN103630882A CN201310597526.8A CN201310597526A CN103630882A CN 103630882 A CN103630882 A CN 103630882A CN 201310597526 A CN201310597526 A CN 201310597526A CN 103630882 A CN103630882 A CN 103630882A
- Authority
- CN
- China
- Prior art keywords
- radar
- chirp rate
- centerdot
- waveform
- value
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/282—Transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种利用多斜率调制实现分布式雷达波形的方法,该方法克服了现有的分布式雷达发射的正交波形多普勒敏感性高、占用频率资源大的缺点。与现有技术中的分布式雷达波形在信号时宽、带宽相同情况下,与相位编码相比,本发明设计的波形具有更低的多普勒敏感性;并且与OFDM-LFM波形相比,本发明设计的波形可使分布式雷达系统中各发射波形共享同一带宽,实现对信号频率资源的充分利用。
Description
技术领域
本发明属于分布式雷达信号处理技术领域,具体涉及一种利用多斜率调制实现分布式雷达波形的方法的设计。
背景技术
分布式雷达系统由多个分布在四周的雷达组成,雷达间协同工作,并从多个不同的角度观测目标,将所得到的目标回波数据进行融合处理,从而提高整个系统的目标检测性能。由于雷达间距离较近,雷达间信号会相互干扰,多个雷达必须发射互不相干的正交波形才能使各雷达接收的目标回波信号相互独立,从而实现各发射信号的区分与联合信号处理。
现有用于分布式雷达的正交波形主要有:正交相位编码波形,OFDM-LFM波形等。虽然这两种波形在分布式雷达领域均得到了应用,但它们还是暴露出一些缺点。由于正交相位编码波形对多普勒频率较为敏感,当分布式雷达探测移动目标或雷达本身在运动时,目标回波信号的相关峰会出现较大程度的衰减,从而增大了目标的检测难度。而采用OFDM-LFM波形需要对每部雷达分配独立的带宽资源,当分布式雷达系统中的雷达数量较大时,将会占据巨大的频率资源,造成频率资源的浪费。
发明内容
本发明的目的是为了解决现有技术中用于分布式雷达的正交波形多普勒敏感性高、占用频率资源大的缺点,而提出一种利用多斜率调制实现分布式雷达波形的方法。
本发明的技术方案是:一种利用多斜率调制实现分布式雷达波形的方法,具体包括:
步骤1:根据分布式雷达系统对目标的分辨率ρr的要求,确定信号的总带宽B,其中C为光速;
步骤2:根据系统对目标的速度变化范围[vmin,vmax]的探测要求,得到多普勒频移f的上限值fmax,其中λ为雷达的载波波长;
步骤5:由步骤3和步骤4得到的T2·Δk的取值上限和下限,计算调频斜率差Δk的取值范围,其表达式为
步骤7:取公式得到的最小整数作为正交波形所需要的调频斜率总数N,将N代入步骤6的各表达式可得到信号时宽T、子信号时宽Tb、平均调频斜率k和调频斜率差Δk的取值,同时确定N个调频斜率的取值为:n=1,2,…,N,则第m个雷达的调频斜率组合可表示为{km,km+1,…,kM+K-1,k0,…,km-1};
步骤8:将步骤7得到的参数N、Tb和调频斜率组合代入式(1),即得到分布式雷达各雷达Sm(t)的发射波形;
其中kmn(n=1,2,…,N)为第m个雷达的N个斜率,km0=0,rect(·)为矩形窗函数,j为虚数单位,t为时间变量。
本发明的有益效果:本发明的一种利用多斜率调制实现分布式雷达波形的方法与现有技术中的分布式雷达波形在信号时宽、带宽相同情况下,与相位编码相比,本发明设计的波形具有更低的多普勒敏感性;并且与OFDM-LFM波形相比,本发明设计的波形可使分布式雷达系统中各发射波形共享同一带宽,实现对信号频率资源的充分利用。
附图说明
图1是本发明实施例一种利用多斜率调制实现分布式雷达波形的方法的流程框图。
具体实施方式
下面结合附图和具体的实施例对本发明作进一步的阐述。
如图1所示为是本发明实施例一种利用多斜率调制实现分布式雷达波形的方法的流程框图,具体包括:
步骤2:根据系统对目标的速度变化范围[vmin,vmax]的探测要求,可得到多普勒频移f的上限值fmax,其中λ为雷达的载波波长;
步骤4:根据系统指标对不同波形的互相值的要求|C(τ)|≤α,根据不等式计算出T2·Δk的取值下限,其中,C(τ)表示任意一部雷达的发射波形与其它雷达的发射波形的互相关值,M为雷达数量,α为系统对不同波形互相关值得要求,根据指标确定,在同一个系统,一般为一常数。
步骤5:由步骤3和步骤4得到调频斜率差Δk的取值范围表达式为
步骤6:考虑到小斜率信号对自相关主瓣展宽的影响,Δk的取值应尽可能小,因此我们通常取步骤5中Δk取值范围的最小值作为Δk的值,即同时我们将最小调频斜率与平均调频斜率之差控制在k/2以内,则调频斜率差需满足其中N为调频斜率总数,k为平均调频斜率k=B/T,结合步骤5,可知求解该关系式可确定信号时宽T的取值范围,进而得到各子信号的时宽
步骤7:取公式得最小整数作为正交波形所需要的调频斜率总数N,将N代入步骤6的各表达式可得到信号时宽T、子信号时宽Tb、平均调频斜率k和调频斜率差Δk的取值,同时我们可以确定N个调频斜率的取值为:n=1,2,…,N,则第m个雷达的调频斜率组合可表示为{km,km+1,…,kM+K-1,k0,…,km-1};
通过上述设计,能得到多斜率调制法设计的分布式雷达波形的各个参数,上述多数参数都是以不等式的形式出现,但在给出的范围内都能满足系统的要求,因此在实际波形设计中,各参数可在分布式雷达系统允许的范围内进行取值。
步骤8:将步骤7得到的参数N、Tb和调频斜率组合代入式(1),便可得到分布式雷达各雷达Sm(t)的发射波形。
其中kmn(n=1,2,…,N)为第m个雷达的N个斜率,km0=0,rect(·)为矩形窗函数,j为虚数单位,t为时间变量。
为了本领域技术人员能够理解并且实施本发明技术方案,将以雷达数M为3的分布式雷达系统为例,对本发明进行详细说明,设系统分辨率ρr为1.5m,雷达的载波波长λ为3cm,目标相对于雷达的最大速度vmax为150m/s,不同波形的互相关值|C(τ)|≤0.1,则分布式雷达的波形设计方法如下:
步骤1:根据分布式雷达系统对目标的分辨率ρr=1.5m的要求,可确定信号的总带宽
步骤2:根据系统对目标的速度变化范围[0,150]的探测要求,得到多普勒频移f的上限值fmax=10kHz;
步骤4:根据系统指标对不同波形的互相值的要求|C(τ)|≤0.1,根据不等式 计算出T2·Δk≥1.6×103;
步骤5:由步骤3和步骤4得到调频斜率差Δk的取值范围表达式为
步骤6:取同时我们将最小调频斜率与平均调频斜率之差控制在k/2以内,则调频斜率差需满足其中N为调频斜率总数,k为平均调频斜率k=108/T,结合步骤5,可知求解该关系式可确定信号时宽T的取值范围,T≥1.6×10-5(N-1),进而得到各子信号的时宽
步骤7:计算取大于4.18的最小整数5作为正交波形所需要的调频斜率总数,即取N=5,将N=5代入步骤6的各表达式可得到信号时宽T≥6.4×10-5、子信号时宽Tb≥1.28×10-5、平均调频斜率k≤1.56×1012和调频斜率差Δk≤3.9×1011,这里,我们取T≥64μs、Tb=12.8μs、k≤1.56×1012Hz/s和Δk=3.9×1011Hz/s作为波形参数的取值,同时我们可以确定N个调频斜率的取值为:k1=2.34×1012Hz/s,k2=1.95×1012Hz/s,k3=1.56×1012Hz/s,k4=1.17×1012Hz/s,k5=0.78×1012Hz/s,则3部雷达的调频斜率组合可表示为:
第一部雷达:{2.34×1012,1.95×1012,1.56×1012,1.17×1012,0.78×1012}
第二部雷达:{1.95×1012,1.56×1012,1.17×1012,0.78×1012,2.34×1012}
第三部雷达:{1.56×1012,1.17×1012,0.78×1012,2.34×1012,1.95×1012}
步骤8:将步骤7求出的参数值代入波形表达式(1)中,即可得到分布式雷达的最终波形。
其中kmn(n=1,2,…,N)为步骤7中第m(m=1,2,3)部雷达的N个斜率,km0=0。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为发明的保护范围并不局限于这样的特别陈述和实施例。凡是根据上述描述做出各种可能的等同替换或改变,均被认为属于本发明的权利要求的保护范围。
Claims (1)
1.一种利用多斜率调制实现分布式雷达波形的方法,其特征在于,具体包括:
步骤5:由步骤3和步骤4得到的T2·Δk的取值上限和下限,计算调频斜率差Δk的取值范围,其表达式为
步骤7:取公式得到的最小整数作为正交波形所需要的调频斜率总数N,将N代入步骤6的各表达式可得到信号时宽T、子信号时宽Tb、平均调频斜率k和调频斜率差Δk的取值,同时确定N个调频斜率的取值为:n=1,2,...,N,则第m个雷达的调频斜率组合可表示为{km,km+1,…,kM+K-1,k0,…,km-1};
步骤8:将步骤7得到的参数N、Tb和调频斜率组合代入式(1),即得到分布式雷达各雷达Sm(t)的发射波形;
其中kmn(n=1,2,…,N)为第m个雷达的N个斜率,km0=0,rect(·)为矩形窗函数,j为虚数单位,t为时间变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310597526.8A CN103630882B (zh) | 2013-11-22 | 2013-11-22 | 一种利用多斜率调制实现分布式雷达波形的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310597526.8A CN103630882B (zh) | 2013-11-22 | 2013-11-22 | 一种利用多斜率调制实现分布式雷达波形的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630882A true CN103630882A (zh) | 2014-03-12 |
CN103630882B CN103630882B (zh) | 2015-05-27 |
Family
ID=50212133
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310597526.8A Expired - Fee Related CN103630882B (zh) | 2013-11-22 | 2013-11-22 | 一种利用多斜率调制实现分布式雷达波形的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630882B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158736A (zh) * | 2015-06-03 | 2015-12-16 | 中国电子科技集团公司第二十八研究所 | 一种mimo雷达发射方向图与波形设计方法 |
CN109239675A (zh) * | 2018-10-23 | 2019-01-18 | 西安空间无线电技术研究所 | 一种基于同带宽不同调频斜率的分布式星载雷达lfm波形确定方法 |
CN113567935A (zh) * | 2021-07-20 | 2021-10-29 | 电子科技大学 | 一种基于多子脉冲处理的高速目标探测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6492938B1 (en) * | 2002-02-11 | 2002-12-10 | Delphi Technologies, Inc. | Method of associating target data in a multi-slope FMCW radar system |
CN1610840A (zh) * | 2002-03-26 | 2005-04-27 | 雷斯昂公司 | 高效宽带波形生成和信号处理设计 |
CN1645163A (zh) * | 2005-01-28 | 2005-07-27 | 武汉大学 | 线性数字调频信号的产生方法 |
CN2840070Y (zh) * | 2004-12-10 | 2006-11-22 | 武汉大学 | 微波雷达系统线性调频信号发生器 |
US7619556B1 (en) * | 2008-02-28 | 2009-11-17 | Rockwell Collins, Inc. | System and method for synthesizing localizer and glide slope deviations from weather radar |
US20100045513A1 (en) * | 2008-08-22 | 2010-02-25 | Microsoft Corporation | Stability monitoring using synthetic aperture radar |
CN102073032A (zh) * | 2010-11-02 | 2011-05-25 | 中国兵器工业第二○六研究所 | 多波形雷达信号的模块化产生方法 |
CN102540187A (zh) * | 2010-12-13 | 2012-07-04 | 电子科技大学 | 编队卫星sar的一种正交波形设计方法 |
CN102955155A (zh) * | 2011-08-26 | 2013-03-06 | 中国科学院空间科学与应用研究中心 | 一种分布式有源相控阵雷达及其波束形成方法 |
CN102998656A (zh) * | 2012-10-11 | 2013-03-27 | 北京理工大学 | 一种基于频率步进的宽带分布式雷达时间同步方法 |
-
2013
- 2013-11-22 CN CN201310597526.8A patent/CN103630882B/zh not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6492938B1 (en) * | 2002-02-11 | 2002-12-10 | Delphi Technologies, Inc. | Method of associating target data in a multi-slope FMCW radar system |
CN1610840A (zh) * | 2002-03-26 | 2005-04-27 | 雷斯昂公司 | 高效宽带波形生成和信号处理设计 |
CN2840070Y (zh) * | 2004-12-10 | 2006-11-22 | 武汉大学 | 微波雷达系统线性调频信号发生器 |
CN1645163A (zh) * | 2005-01-28 | 2005-07-27 | 武汉大学 | 线性数字调频信号的产生方法 |
US7619556B1 (en) * | 2008-02-28 | 2009-11-17 | Rockwell Collins, Inc. | System and method for synthesizing localizer and glide slope deviations from weather radar |
US20100045513A1 (en) * | 2008-08-22 | 2010-02-25 | Microsoft Corporation | Stability monitoring using synthetic aperture radar |
CN102073032A (zh) * | 2010-11-02 | 2011-05-25 | 中国兵器工业第二○六研究所 | 多波形雷达信号的模块化产生方法 |
CN102540187A (zh) * | 2010-12-13 | 2012-07-04 | 电子科技大学 | 编队卫星sar的一种正交波形设计方法 |
CN102955155A (zh) * | 2011-08-26 | 2013-03-06 | 中国科学院空间科学与应用研究中心 | 一种分布式有源相控阵雷达及其波束形成方法 |
CN102998656A (zh) * | 2012-10-11 | 2013-03-27 | 北京理工大学 | 一种基于频率步进的宽带分布式雷达时间同步方法 |
Non-Patent Citations (2)
Title |
---|
宗竹林等: "编队卫星雷达正交码分LFM信号设计及性能分析", 《电子科技大学学报》 * |
鲍坤超等: "多发射体制下小卫星分布式雷达系统的波形设计", 《电子与信息学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158736A (zh) * | 2015-06-03 | 2015-12-16 | 中国电子科技集团公司第二十八研究所 | 一种mimo雷达发射方向图与波形设计方法 |
CN105158736B (zh) * | 2015-06-03 | 2017-06-06 | 中国电子科技集团公司第二十八研究所 | 一种mimo雷达发射方向图与波形设计方法 |
CN109239675A (zh) * | 2018-10-23 | 2019-01-18 | 西安空间无线电技术研究所 | 一种基于同带宽不同调频斜率的分布式星载雷达lfm波形确定方法 |
CN113567935A (zh) * | 2021-07-20 | 2021-10-29 | 电子科技大学 | 一种基于多子脉冲处理的高速目标探测方法 |
CN113567935B (zh) * | 2021-07-20 | 2023-09-19 | 电子科技大学 | 一种基于多子脉冲处理的高速目标探测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103630882B (zh) | 2015-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2628566C1 (ru) | Способ работы радиолокационной станции с повышенными допплеровскими характеристиками | |
US10048353B2 (en) | Interference cancellation in an FMCW radar | |
US11582305B2 (en) | Vehicle radar system with a shared radar and communication system | |
EP3436845B1 (en) | Direct detection lidar system and method with frequency modulation (fm) transmitter and quadrature receiver | |
CN102798866B (zh) | 激光雷达系统及采用该系统的正弦调幅-脉冲相位编码调制的复合测距和测速方法 | |
CN108415010B (zh) | 一种基于梯形lfmcw调制的雷达多目标检测方法 | |
Shah et al. | Step-frequency radar with compressive sampling (SFR-CS) | |
US20160131742A1 (en) | Angle-resolving fmcw radar sensor | |
CN105103003A (zh) | 快速线性调频fmcw雷达 | |
CN105738889B (zh) | 一种调频连续波测速测距方法 | |
AU2013227470A1 (en) | Determining position of underwater node | |
CN103605112B (zh) | 多发多收干涉合成孔径雷达时频二维信号波形设计方法 | |
CN104931973A (zh) | 用于激光雷达系统的非对称相位编码测距方法 | |
CN103605113B (zh) | 多发多收干涉合成孔径雷达空时二维信号波形设计方法 | |
CN104931975A (zh) | 基于微波光子信号处理的相位编码激光成像雷达 | |
Capria et al. | DVB-T passive radar for vehicles detection in urban environment | |
CN103630882B (zh) | 一种利用多斜率调制实现分布式雷达波形的方法 | |
CN105022037A (zh) | 一种基于超混沌编码的汽车雷达交叉干扰抑制方法 | |
CN109239675A (zh) | 一种基于同带宽不同调频斜率的分布式星载雷达lfm波形确定方法 | |
KR20200085903A (ko) | 상호 지연된 직교 코드에 의해 개선된 레인지 해상도를 갖는 레이더 시스템의 작동 방법 및 레이더 시스템 | |
CN102998659A (zh) | 基于脉间调制的多普勒频谱赋形方法及系统 | |
US20220236403A1 (en) | Ofdm radar sensor system having an actively retransmitting repeater | |
RU2699240C1 (ru) | Способ определения координат цели в рлс с непрерывным излучением | |
RU2580507C2 (ru) | Способ радиолокации и радиолокатор с доплеровским передатчиком для его реализации | |
US3242490A (en) | Apparatus for and a method of determining the position of electromagnetic noise source |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150527 Termination date: 20161122 |