CN103777186A - 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 - Google Patents
基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 Download PDFInfo
- Publication number
- CN103777186A CN103777186A CN201310544994.9A CN201310544994A CN103777186A CN 103777186 A CN103777186 A CN 103777186A CN 201310544994 A CN201310544994 A CN 201310544994A CN 103777186 A CN103777186 A CN 103777186A
- Authority
- CN
- China
- Prior art keywords
- radar
- target
- field
- near field
- bin
- 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/41—Details 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
-
- 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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/883—Radar 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)中波束分解公式如下:
其中,为具有较大波束宽度的宽波束在空间上的场值,pB(α′)为宽波束在谱域上的角分布,pb(α′)为具有较小波束宽度的窄波束在谱域上的角分布,u(X)窄波束的加权幅值。
优选的,步骤4)中窄波束照射面元满足远场条件公式为:
其中,D表示窄波束照射面元的最大线尺寸,Rn表示被照射面元与第n个子波束间的距离,λ为入射波的波长。
优选的,步骤5)计算近场后向回波功率公式为:
其中,Pr为接收机接收到的回波功率,Pt为发射机发射电磁波的功率,Gt为雷达发射天线的增益,Gr为雷达接受天线的增益,若采用收发共用天线则Gt与Gr相同,r1和r2分别是发射机到目标和目标到接收机的距离,λ为入射波波长。
采用了上述技术方案,本发明的有益效果为:本发明采用波束分解的方法将宽波束分解为众多窄波束并利用物理光学局部近似的特点将远场计算方法推广到近场计算,将远场雷达截面带入雷达方程获得雷达回波功率时间序列。本发明所述的方法既有清晰的物理意义又适合工程测量的物理量,并结合物理光学方法解决了计算电大尺寸运动目标的近场雷达回波重大技术难题,对运动目标近场电磁回波仿真模拟问题的研究在导弹引信设计方面具有重要的意义。
图1是本发明实施例的运动目标近场回波计算流程图;
图2是本发明实施例脱靶量和脱靶方位示意图;
其中:●表示探测器;○表示目标;ρ:脱靶量;α:脱靶方位
具体实施方式
下面结合实施例对本发明进一步说明。
参照图1,本发明的具体实现步骤如下:
步骤1,确定雷达(探测器)坐标系、运动目标坐标系及两者间的变换关系
运动目标和雷达(探测器)往往都定义在以各自为中心的局部坐标系下,需要依据不同的需求对其进行变换:如当确定脱靶量和脱靶角时需将探测器变换到以目标为中心的坐标系中,即探测器相对目标运动;当要确定波束照射目标时目标个部件的照射情况时,需将目标变换到以探测器为中心的坐标系中,目标相对探测器运动。
步骤2,探测器和目标的相对速度坐标系,由此获得脱靶量
探测器和目标都以一定的速度进行运动,为了方便计算脱靶量确定一个以目标为中心,探测器相对目标运动的相对速度坐标系。通过目标的几何中心并与目标中心轴垂直的平面称为脱靶平面。在脱靶平面内定义脱靶量和脱靶方位。探测器相对运动轨迹与脱靶平面的交点P称为脱靶点,脱靶点P与目标中心O的连线OP称为脱靶量,用ρ表示。OP连线在脱靶平面上的方位角称为脱靶方位,用α表示。脱靶量和脱靶方位见图2.
步骤3,给出雷达波束描述对入射波束进行分解,对目标进行面元剖分,目标局部照射区域确定和面元遮挡判断
波束分解公式如下:
其中,为具有较大波束宽度的宽波束在空间上的场值,PB(α′)为宽波束在谱域上的角分布,pb(α′)为具有较小波束宽度的窄波束在谱域上的角分 布,u(X)窄波束的加权幅值。
确定每一个窄波束在目标上的局部照射区域,并对面元进行遮挡判断,完成散射场计算前的预处理过程。
步骤4,子波束照射面元满足远场条件,应用物理光学法叠加获得后向雷达截面
窄波束照射面元满足远场条件公式为:
其中,D表示窄波束照射面元的最大线尺寸,Rn表示被照射面元与第n个子波束间的距离,λ为入射波的波长。
将物理光学计算的每个区域的散射场叠加得到总散射场。
步骤5,将后向雷达截面代入雷达方程计算出近场回波功率
计算近场雷达回波功率的公式如下
其中,Pr为接收机接收到的回波功率,Pt为发射机发射电磁波的功率,Gt为雷达发射天线的增益,Gr为雷达接受天线的增益,若采用收发共用天线则Gt与Gr相同,r1和r2分别是发射机到目标和目标到接收机的距离,λ为入射波波长。
步骤6,重复步骤4)和步骤5)计算出不同时刻下由于目标与探测器相对位置的变化而得到的不同回波功率,即回波功率时间序列。对回波功率时间序列进行傅里叶变换得到相应的频谱。
本发明不局限于上述具体的实施方式,本领域的普通技术人员从上述构思出发,不经过创造性的劳动,所作出的种种变换,均落在本发明的保护范围之内。
Claims (4)
1.一种基于波束分解和局部照射下运动目标近场雷达回波特性计算方法,其特征在于,包括以下步骤:
1)确定雷达(探测器)坐标系、运动目标坐标系及两者间的变换关系;
2)探测器和目标的相对速度坐标系,由此获得脱靶量;
3)给出雷达波束描述对入射波束进行分解,对目标进行面元剖分,目标局部照射区域确定和面元遮挡判断;
4)子波束照射面元满足远场条件,应用物理光学法叠加获得后向雷达截面;
5)将后向雷达截面代入雷达方程计算出近场回波功率;
6)重复步骤4)和步骤5)计算出不同时刻下由于目标与探测器相对位置的变化而得到的不同回波功率,即回波功率时间序列以及对应的频谱特性。
3.如权利要求1所述的基于波束分解和局部照射下运动目标近场雷达回波特性计算方法,其特征在于,步骤4)窄波束照射面元满足远场条件公式如下:
其中,D表示窄波束照射面元的最大线尺寸,Rn表示被照射面元与第n个子波束间的距离,λ为入射波的波长。
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 true CN103777186A (zh) | 2014-05-07 |
CN103777186B 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) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104076353A (zh) * | 2014-06-30 | 2014-10-01 | 西安空间无线电技术研究所 | 一种面目标回波波束中心速度测量方法 |
CN104992064A (zh) * | 2015-07-09 | 2015-10-21 | 西安电子科技大学 | 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法 |
CN105486180A (zh) * | 2015-12-15 | 2016-04-13 | 西安电子科技大学 | 基于波束分解和局部照射的激光引信近场回波功率计算 |
CN108061883A (zh) * | 2017-11-28 | 2018-05-22 | 上海无线电设备研究所 | 局部散射源反演的近场散射函数转换弹目交会回波的方法 |
CN109520383A (zh) * | 2017-09-20 | 2019-03-26 | 南京理工大学 | 基于matlab的体目标回波模拟方法 |
CN111896816A (zh) * | 2020-08-10 | 2020-11-06 | 西安电子科技大学 | 基于远场平面扫描的毫米波辐射器近场辐射剂量测量方法 |
CN112069713A (zh) * | 2020-09-10 | 2020-12-11 | 上海无线电设备研究所 | 一种近场散射特性建模方法、电子设备及存储介质 |
CN112946610A (zh) * | 2021-01-26 | 2021-06-11 | 中国科学院半导体研究所 | 多参数模拟的激光目标模拟器 |
CN113687446A (zh) * | 2021-07-19 | 2021-11-23 | 西安空间无线电技术研究所 | 一种基于星载综合孔径微波辐射计的近场外定标方法 |
Citations (1)
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) |
-
2013
- 2013-11-02 CN CN201310544994.9A patent/CN103777186B/zh not_active Expired - Fee Related
Patent Citations (1)
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)
Title |
---|
伍振兴 等: "使用波束分解算法分析海洋表面电磁散射", 《电波科学学报》 * |
吴振森 等: "激光引信目标近场散射特性研究", 《激光与红外》 * |
李铁 等: "无线电引信弹目交会信号的数字仿真", 《兵工学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104076353A (zh) * | 2014-06-30 | 2014-10-01 | 西安空间无线电技术研究所 | 一种面目标回波波束中心速度测量方法 |
CN104076353B (zh) * | 2014-06-30 | 2016-08-17 | 西安空间无线电技术研究所 | 一种面目标回波波束中心速度测量方法 |
CN104992064A (zh) * | 2015-07-09 | 2015-10-21 | 西安电子科技大学 | 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法 |
CN105486180A (zh) * | 2015-12-15 | 2016-04-13 | 西安电子科技大学 | 基于波束分解和局部照射的激光引信近场回波功率计算 |
CN105486180B (zh) * | 2015-12-15 | 2017-03-22 | 西安电子科技大学 | 基于波束分解和局部照射的激光引信近场回波功率计算方法 |
CN109520383A (zh) * | 2017-09-20 | 2019-03-26 | 南京理工大学 | 基于matlab的体目标回波模拟方法 |
CN108061883A (zh) * | 2017-11-28 | 2018-05-22 | 上海无线电设备研究所 | 局部散射源反演的近场散射函数转换弹目交会回波的方法 |
CN108061883B (zh) * | 2017-11-28 | 2021-06-04 | 上海无线电设备研究所 | 局部散射源反演的近场散射函数转换弹目交会回波的方法 |
CN111896816A (zh) * | 2020-08-10 | 2020-11-06 | 西安电子科技大学 | 基于远场平面扫描的毫米波辐射器近场辐射剂量测量方法 |
CN112069713A (zh) * | 2020-09-10 | 2020-12-11 | 上海无线电设备研究所 | 一种近场散射特性建模方法、电子设备及存储介质 |
CN112069713B (zh) * | 2020-09-10 | 2022-12-23 | 上海无线电设备研究所 | 一种近场散射特性建模方法、电子设备及存储介质 |
CN112946610A (zh) * | 2021-01-26 | 2021-06-11 | 中国科学院半导体研究所 | 多参数模拟的激光目标模拟器 |
CN112946610B (zh) * | 2021-01-26 | 2023-10-03 | 中国科学院半导体研究所 | 多参数模拟的激光目标模拟器 |
CN113687446A (zh) * | 2021-07-19 | 2021-11-23 | 西安空间无线电技术研究所 | 一种基于星载综合孔径微波辐射计的近场外定标方法 |
CN113687446B (zh) * | 2021-07-19 | 2024-02-09 | 西安空间无线电技术研究所 | 一种基于星载综合孔径微波辐射计的近场外定标方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103777186B (zh) | 2016-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103777186A (zh) | 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 | |
CN103176170B (zh) | 一种基于gpu并行计算的sar回波模拟方法 | |
CN103530469B (zh) | 一种角反射器雷达截面积的计算方法 | |
CN104865562A (zh) | 基于混合模型的雷达非合作目标的识别方法 | |
CN104007440A (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN106886021A (zh) | 高分辨率星载sar成像质量提升方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN107783092A (zh) | 基于链条关系式的近场后向rcs测量***及方法 | |
CN104992064B (zh) | 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法 | |
Ding et al. | The Qujing incoherent scatter radar: system description and preliminary measurements | |
CN106842162A (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 | |
CN106772284A (zh) | 考虑近海面大气环境影响的编队内雷达远场预测方法 | |
CN104101869B (zh) | 一种极坐标下的地波雷达运动目标仿真建模方法 | |
Li et al. | Cascaded method for ionospheric decontamination and sea clutter suppression for high‐frequency hybrid sky‐surface wave radar | |
CN103729485B (zh) | 一种基于dem数据的宽带雷达相干杂波仿真方法 | |
CN103699810A (zh) | 一种粗糙面微波段双向反射分布函数的建模方法 | |
CN103954969B (zh) | 基于激光一维距离像反演动态锥形目标姿态角和几何尺寸 | |
CN101846739A (zh) | Sar扩展场景原始数据混合域仿真方法 | |
CN103675777B (zh) | 基于拟合法的机载雷达杂波模拟方法及装置 | |
Knepp | Multiple phase screen calculation of two‐way spherical wave propagation in the ionosphere | |
Borzov et al. | Mathematical modeling and simulation of the input signals of short-range radar systems | |
Weinmann | Accurate prediction of EM scattering by wind turbines | |
CN104914425A (zh) | 一种超电大尺寸的强电磁脉冲环境时频空多维分析模型 | |
Wu et al. | A semi-empirical model of sea clutter based on zero memory nonlinearity |
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 |
Addressee: Patent of Xi'an University of Electronic Science and technology The person in charge Document name: payment instructions |
|
DD01 | Delivery of document by public notice | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161005 Termination date: 20211102 |
|
CF01 | Termination of patent right due to non-payment of annual fee |