CN103777186B - 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 - Google Patents

基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 Download PDF

Info

Publication number
CN103777186B
CN103777186B CN201310544994.9A CN201310544994A CN103777186B CN 103777186 B CN103777186 B CN 103777186B CN 201310544994 A CN201310544994 A CN 201310544994A CN 103777186 B CN103777186 B CN 103777186B
Authority
CN
China
Prior art keywords
radar
target
prime
field
alpha
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.)
Expired - Fee Related
Application number
CN201310544994.9A
Other languages
English (en)
Other versions
CN103777186A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201310544994.9A priority Critical patent/CN103777186B/zh
Publication of CN103777186A publication Critical patent/CN103777186A/zh
Application granted granted Critical
Publication of CN103777186B publication Critical patent/CN103777186B/zh
Expired - Fee Related 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/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/883Radar or analogous systems specially adapted for specific applications for missile homing, autodirectors

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于波束分解和局部照射下运动目标近场雷达回波特性的计算方法,包括以下步骤:确定雷达探测器坐标系、运动目标坐标系及两者间的变换关系;雷达探测器和目标的相对速度坐标系,由此获得脱靶量;给出雷达波束描述对入射波束进行分解,对目标进行面元剖分,目标局部照射区域确定和面元遮挡判断;子波束照射面元满足远场条件,应用物理光学法叠加获得后向雷达截面;利用雷达方程获得目标近场雷达回波功率时间序列,以及对应的频谱特性。本发明利用波束分解和物理光学局部近似的特点将远场计算方法推广到近场计算,由雷达方程计算出近场回波功率的时间序列和相应的频谱。

Description

基于波束分解和局部照射下运动目标近场雷达回波特性计算方法
技术领域
本发明涉及一种雷达回波计算方法,采用基于波束分解和局部照射下运动目标近场雷达回波特性计算方法,波束分解能有效减少宽波束散射计算对内存的需求。局部照射使得目标在局部区域上满足远场条件,可用远场物理光学法计算近场回波,降低了计算复杂度。
背景技术
一直以来大量的文献资料都关注于雷达目标回波的远场特性,然而雷达目标的近场特性也有着广泛的用途。例如,在暗室测量时,由于目标与测量雷达间的距离受暗室尺寸的限制,目标于雷达间的距离有时并不能严格满足远场条件;具有近炸引信的导弹在对飞机、导弹等飞行目标进行攻击时,引信与目标间的相互作用也属于近场作用;另外,许多电磁兼容问题也都属于近场电磁问题的范畴。
电大(几何尺寸远大于雷达波长)尺寸运动目标在波束局部照射下的回波计算也是一个有待解决的问题,虽然解析方法和数值方法能得到较精确的结果,但是它们受限于正交坐标系和计算机内存的限制,只能计算规则目标或者电小尺寸复杂目标的问题。物理光学方法是一种近似方法能较快的计算具有复杂形状的目标在高频区的散射问题,但常用的物理光学法一般是基于远场条件的,近场问题往往会变得更为复杂。宽波束入射会使得计算所需的内存剧增。
发明内容
本发明所要解决的技术问题是:提供一种基于波束分解和局部照射下运动目标近场雷达回波特性计算方法,采用波束分解的方法将宽波束分解为众多窄波束,利用物理光学局部近似的特点将远场计算方法推广到近场计算,结合雷达方程获得雷达回波的时间序列。
为解决上述技术问题,本发明的技术方案是:一种基于波束分解和局部照射下运动目标近场雷达回波特性计算方法,包括以下步骤:
1)确定雷达探测器坐标系、运动目标坐标系及两者间的变换关系;
2)雷达探测器和目标的相对速度坐标系,由此获得脱靶量;
3)给出雷达波束描述对入射波束进行分解,对目标进行面元剖分,目标局部照射区域确定和面元遮挡判断;
4)子波束照射面元满足远场条件,应用物理光学法叠加获得后向雷达截面;
5)将后向雷达截面代入雷达方程计算出近场回波功率;
6)重复步骤4)和步骤5)计算出不同时刻下由于目标与雷达探测器相对位置的变化而得到的不同回波功率,即回波功率时间序列以及对应的频谱特性;
优选的,步骤3)中波束分解公式如下:
E B i ( x , y ) = ∫ - k k p B ( α ′ ) e i ( α ′ x - β ′ y ) dα ′ = ∫ - ∞ + ∞ u ( X ) [ ∫ - k k p b ( α ′ ) e i ( α ′ ( x - X ) - β ′ y ) dα ′ ] d X
其中,为具有较大波束宽度的宽波束在空间上的场值,pB(α′)为宽波束在谱域上的角分布,pb(α′)为具有较小波束宽度的窄波束在谱域上的角分布,u(X)窄波束的加权幅值。
优选的,步骤4)中窄波束照射面元满足远场条件公式为:
R n ≥ 2 D 2 λ
其中,D表示窄波束照射面元的最大线尺寸,Rn表示被照射面元与第n个子波束间的距离,λ为入射波的波长。
优选的,步骤5)计算近场后向回波功率公式为:
P r = G t P t 4 πr 1 2 σ 4 πr 2 2 G r λ 2 4 π
其中,Pr为接收机接收到的回波功率,Pt为发射机发射电磁波的功率,Gt为雷达发射天线的增益,Gr为雷达接受天线的增益,若采用收发共用天线则Gt与Gr相同,r1和r2分别是发射机到目标和目标到接收机的距离,λ为入射波波长。
采用了上述技术方案,本发明的有益效果为:本发明采用波束分解的方法将宽波束分解为众多窄波束并利用物理光学局部近似的特点将远场计算方法推广到近场计算,将远场雷达截面带入雷达方程获得雷达回波功率时间序列。本发明所述的方法既有清晰的物理意义又适合工程测量的物理量,并结合物理光学方法解决了计算电大尺寸运动目标的近场雷达回波重大技术难题,对运动目标近场电磁回波仿真模拟问题的研究在导弹引信设计方面具有重要的意义。
具体实施方式
下面结合实施例对本发明进一步说明。
参照图1,本发明的具体实现步骤如下:
步骤1,确定雷达探测器坐标系、运动目标坐标系及两者间的变换关系
运动目标和雷达探测器往往都定义在以各自为中心的局部坐标系下,需要依据不同的需求对其进行变换:如当确定脱靶量和脱靶角时需将雷达探测器变换到以目标为中心的坐标系中,即雷达探测器相对目标运动;当要确定波束照射目标时目标个部件的照射情况时,需将目标变换到以雷达探测器为中心的坐标系中,目标相对雷达探测器运动。
步骤2,雷达探测器和目标的相对速度坐标系,由此获得脱靶量
雷达探测器和目标都以一定的速度进行运动,为了方便计算脱靶量确定一个以目标为中心,雷达探测器相对目标运动的相对速度坐标系。通过目标的几何中心并与目标中心轴垂直的平面称为脱靶平面。在脱靶平面内定义脱靶量和脱靶方位。雷达探测器相对运动轨迹与脱靶平面的交点P称为脱靶点,脱靶点P与目标中心O的连线OP称为脱靶量,用ρ表示。OP连线在脱靶平面上的方位角称为脱靶方位,用α表示。脱靶量和脱靶方位见图2.
步骤3,给出雷达波束描述对入射波束进行分解,对目标进行面元剖分,目标局部照射区域确定和面元遮挡判断
波束分解公式如下:
E B i ( x , y ) = ∫ - k k p B ( α ′ ) e i ( α ′ x - β ′ y ) dα ′ = ∫ - ∞ + ∞ u ( X ) [ ∫ - k k p b ( α ′ ) e i ( α ′ ( x - X ) - β ′ y ) dα ′ ] d X
其中,为具有较大波束宽度的宽波束在空间上的场值,pB(α′)为宽波束在谱域上的角分布,pb(α′)为具有较小波束宽度的窄波束在谱域上的角分布,u(X)窄波束的加权幅值。
确定每一个窄波束在目标上的局部照射区域,并对面元进行遮挡判断,完成散射场计算前的预处理过程。
步骤4,子波束照射面元满足远场条件,应用物理光学法叠加获得后向雷达截面
窄波束照射面元满足远场条件公式为:
R n ≥ 2 D 2 λ
其中,D表示窄波束照射面元的最大线尺寸,Rn表示被照射面元与第n个子波束间的距离,λ为入射波的波长。
将物理光学计算的每个区域的散射场叠加得到总散射场。
步骤5,将后向雷达截面代入雷达方程计算出近场回波功率
计算近场雷达回波功率的公式如下
P r = G t P t 4 πr 1 2 σ 4 πr 2 2 G r λ 2 4 π
其中,Pr为接收机接收到的回波功率,Pt为发射机发射电磁波的功率,Gt为雷达发射天线的增益,Gr为雷达接受天线的增益,若采用收发共用天线则Gt与Gr相同,r1和r2分别是发射机到目标和目标到接收机的距离,λ为入射波波长。
步骤6,重复步骤4)和步骤5)计算出不同时刻下由于目标与雷达探测器相对位置的变化而得到的不同回波功率,即回波功率时间序列。对回波功率时间序列进行傅里叶变换得到相应的频谱。
本发明不局限于上述具体的实施方式,本领域的普通技术人员从上述构思出发,不经过创造性的劳动,所作出的种种变换,均落在本发明的保护范围之内。

Claims (1)

1.一种基于波束分解和局部照射下运动目标近场雷达回波特性计算方法,其特征在于,包括以下步骤:
1)确定雷达探测器坐标系、运动目标坐标系及两者间的变换关系;
2)雷达探测器和目标的相对速度坐标系,由此获得脱靶量;
3)给出雷达波束描述对入射波束进行分解,对目标进行面元剖分,目标局部照射区域确定和面元遮挡判断;
4)子波束照射面元满足远场条件,应用物理光学法叠加获得后向雷达截面;
5)将后向雷达截面代入雷达方程计算出近场回波功率;
6)重复步骤4)和步骤5)计算出不同时刻下由于目标与雷达探测器相对位置的变化而得到的不同回波功率,即回波功率时间序列以及对应的频谱特性;
步骤3)中波束分解公式如下:
E B i ( x , y ) = ∫ - k k p B ( α ′ ) e i ( α ′ x - β ′ y ) dα ′ = ∫ - ∞ + ∞ u ( X ) [ ∫ - k k p b ( α ′ ) e i ( α ′ ( x - X ) - β ′ y ) dα ′ ] d X \ * M E R G E F O R M A T - - - ( 1 )
其中,为具有较大波束宽度的宽波束在空间上的场值,pB(α′)为宽波束在谱域上的角分布,pb(α′)为具有较小波束宽度的窄波束在谱域上的角分布,u(X)窄波束的加权幅值;
步骤4)窄波束照射面元满足远场条件公式如下:
R n ≥ 2 D 2 λ \ * M E R G E F O R M A T - - - ( 2 )
其中,D表示窄波束照射面元的最大线尺寸,Rn表示被照射面元与第n个子波束间的距离,λ为入射波的波长;
将步骤4),步骤5)和步骤6)的计算运动目标近场后向回波功率的时间序列,依据下式计算近场后向回波功率:
P r = G t P t 4 πr 1 2 σ 4 πr 2 2 G r λ 2 4 π \ * M E R G E F O R M A T - - - ( 3 )
其中,Pr为接收机接收到的回波功率,Pt为发射机发射电磁波的功率,Gt为雷达发射天线的增益,Gr为雷达接受天线的增益,若采用收发共用天线则Gt与Gr相同,r1和r2分别是发射机到目标和目标到接收机的距离,λ为入射波波长。
CN201310544994.9A 2013-11-02 2013-11-02 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 Expired - Fee Related CN103777186B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310544994.9A CN103777186B (zh) 2013-11-02 2013-11-02 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310544994.9A CN103777186B (zh) 2013-11-02 2013-11-02 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法

Publications (2)

Publication Number Publication Date
CN103777186A CN103777186A (zh) 2014-05-07
CN103777186B true CN103777186B (zh) 2016-10-05

Family

ID=50569652

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310544994.9A Expired - Fee Related CN103777186B (zh) 2013-11-02 2013-11-02 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法

Country Status (1)

Country Link
CN (1) CN103777186B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104076353B (zh) * 2014-06-30 2016-08-17 西安空间无线电技术研究所 一种面目标回波波束中心速度测量方法
CN104992064B (zh) * 2015-07-09 2017-11-10 西安电子科技大学 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法
CN105486180B (zh) * 2015-12-15 2017-03-22 西安电子科技大学 基于波束分解和局部照射的激光引信近场回波功率计算方法
CN109520383B (zh) * 2017-09-20 2020-07-31 南京理工大学 基于matlab的体目标回波模拟方法
CN108061883B (zh) * 2017-11-28 2021-06-04 上海无线电设备研究所 局部散射源反演的近场散射函数转换弹目交会回波的方法
CN111896816A (zh) * 2020-08-10 2020-11-06 西安电子科技大学 基于远场平面扫描的毫米波辐射器近场辐射剂量测量方法
CN112069713B (zh) * 2020-09-10 2022-12-23 上海无线电设备研究所 一种近场散射特性建模方法、电子设备及存储介质
CN112946610B (zh) * 2021-01-26 2023-10-03 中国科学院半导体研究所 多参数模拟的激光目标模拟器
CN113687446B (zh) * 2021-07-19 2024-02-09 西安空间无线电技术研究所 一种基于星载综合孔径微波辐射计的近场外定标方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8325417B1 (en) * 2008-09-23 2012-12-04 The United States Of America As Represented By The Secretary Of The Navy Apparatus and method for producing an annular far-field patterned beam(s)

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8325417B1 (en) * 2008-09-23 2012-12-04 The United States Of America As Represented By The Secretary Of The Navy Apparatus and method for producing an annular far-field patterned beam(s)

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
使用波束分解算法分析海洋表面电磁散射;伍振兴 等;《电波科学学报》;20050630;第20卷(第3期);第325-329页 *
无线电引信弹目交会信号的数字仿真;李铁 等;《兵工学报》;20050930;第26卷(第5期);第634-637页 *
激光引信目标近场散射特性研究;吴振森 等;《激光与红外》;20011031;第31卷(第5期);第294-297页 *

Also Published As

Publication number Publication date
CN103777186A (zh) 2014-05-07

Similar Documents

Publication Publication Date Title
CN103777186B (zh) 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法
CN103235292B (zh) 平面相控阵调零保形校正的全维和差测角方法
CN102253127A (zh) 用于变压器局部放电定位的l形超声阵列传感器及其方法
CN104391183A (zh) 近场测量天线远区场特性快速计算方法
CN107976660A (zh) 弹载多通道雷达超低空目标分析与多径回波建模方法
CN104992064B (zh) 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法
CN106872981A (zh) 雨量雷达的降水强中心跟踪与预报方法
CN106886021A (zh) 高分辨率星载sar成像质量提升方法
CN109932719A (zh) 基于sar成像的rcs高精度测量方法
CN107783092A (zh) 基于链条关系式的近场后向rcs测量系统及方法
CN104502911A (zh) 一种穿墙成像雷达的墙壁参数估计方法
Gill et al. On the development of a second-order bistatic radar cross section of the ocean surface: A high-frequency result for a finite scattering patch
CN113376597A (zh) 基于数字高程地图和gpu的复杂地形电磁散射快速仿真方法
CN110244273A (zh) 一种基于均匀分布式阵列的目标角度估计方法
Ya‐jun et al. Analysis and simulation for broadening first‐order sea clutter spectrum in high frequency hybrid sky‐surface wave propagation mode
CN104101869B (zh) 一种极坐标下的地波雷达运动目标仿真建模方法
CN107329117A (zh) 一种基于改进omp的双基机载雷达自适应杂波谱补偿方法
CN103729485B (zh) 一种基于dem数据的宽带雷达相干杂波仿真方法
CN103954969B (zh) 基于激光一维距离像反演动态锥形目标姿态角和几何尺寸
CN106646418A (zh) 一种基于改进omp的机载雷达杂波空时谱快速估计方法
Knepp Multiple phase screen calculation of two‐way spherical wave propagation in the ionosphere
CN103675777B (zh) 基于拟合法的机载雷达杂波模拟方法及装置
Li et al. Clutter suppression approach for end-fire array airborne radar based on adaptive segmentation
CN112816794A (zh) 天线波束照射下近场动态海面电磁波强度分布可视化方法
Weinmann Accurate prediction of EM scattering by wind turbines

Legal Events

Date Code Title Description
C06 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
DD01 Delivery of document by public notice
DD01 Delivery of document by public notice

Addressee: Patent of Xi'an University of Electronic Science and technology The person in charge

Document name: payment instructions

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: 20161005

Termination date: 20211102