CN103675781B - 一种准确的近场回波获取方法 - Google Patents
一种准确的近场回波获取方法 Download PDFInfo
- Publication number
- CN103675781B CN103675781B CN201210362109.0A CN201210362109A CN103675781B CN 103675781 B CN103675781 B CN 103675781B CN 201210362109 A CN201210362109 A CN 201210362109A CN 103675781 B CN103675781 B CN 103675781B
- Authority
- CN
- China
- Prior art keywords
- field
- rho
- scattering
- scattering unit
- antenna
- 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
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/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
-
- 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/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
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
本发明属于信号特征控制技术领域,具体涉及一种准确的近场回波获取方法。该方法包括如下步骤:步骤一:建立目标坐标系O-XtYtZt,天线坐标系O-XaYaZa,散射单元坐标系O-XmYmZm;步骤二:在散射体表面,获得入射到第m个散射单元上的场;步骤三:获得散射单元的散射场;步骤四:获得接收天线接收到的第m个散射单元的散射场;步骤五:将接收天线收到的每个散射单元的散射场进行叠加,得到总的散射场;步骤六:计算出总散射场,进而得到近场回波强度和功率特性。该方法利用PO+PTD+GO法建立了复杂目标近场散射模型,并引入近场极化散射矩阵的概念,将不同机理的散射特性用统一表达式表示,简化了计算过程。
Description
技术领域
本发明属于信号特征控制技术领域,具体涉及一种准确的近场回波获取方法。
背景技术
防空导弹武器系统定型及作战经常受到无线电引信系统不启动或误启动的困扰,引信要确定最佳启动点和启动区,需要研究目标近场回波特性。而近场回波特性的获取需要充分考虑近场散射的特点,如天线方向图、球面波、距离、局部照射、极化等影响,而其理论模型的准确性必须经过复杂目标近场测试数据的校验才能得到确认。近年来一些单位相继开展了复杂目标近场回波建模方法的研究,并取得了一定的成果,但是未能够全面考虑这些方面的影响,也未能给出令人信服的理论模型校验。尤其是针对天线方向图和局部照射的影响,国内外相关近场特性相关文献很少给出定量的计算描述。
发明内容
本发明的目的在于提供一种准确的近场回波获取方法,以克服现有技术存在的上述不足。
为达到上述目的,本发明所采取的技术方案为:
一种准确的近场回波获取方法,该方法包括如下步骤:
步骤一:建立目标坐标系O-XtYtZt,天线坐标系O-XaYaZa,散射单元坐标系O-XmYmZm;
步骤二:在散射体表面,获得入射到第m个散射单元上的场;
步骤三:获得散射单元的散射场;
步骤四:获得接收天线接收到的第m个散射单元的散射场;
步骤五:将接收天线收到的每个散射单元的散射场进行叠加,得到总的散射场;
步骤六:计算出总散射场,进而得到近场回波强度和功率特性。
所述步骤二的具体实现方法如下:在散射体表面,将近场入射波表示成球面波形式,入射到第m个散射单元的场强为:
式中:i表示发射;Pi为发射的信号功率;Di为发射天线的增益;Z0为空气的波阻抗,且Z0=120πΩ;Rim是发射天线相位中心与散射单元中心之间的距离;ω=2πf为电磁振荡的角频率;fim为第m个散射单元处的天线方向性函数;φ为场的相位特性;E0m为第m个散射单元处电场的极化特性;k表示电磁波在自由空间内的波数;t表示时间参数;f为载频频率;
入射到第m个散射单元的场强写成以下矩阵形式:
式中表示发射天线的极化分量。
所述步骤三中获得散射单元的散射场,表示为:
其中:
式中:Esm为散射电场强度,Esm1、Esm2为Esm的两个分量;Rim为散射单元中心到发射天线相位中心的距离;Rsm为散射单元中心到接收天线相位中心的距离;为散射单元的极化散射矩阵,其四个分量是四种反射系数,其中对应于入射和反射电场的极化都在上进行时的反射系数;对应于入射电场的极化在上进行,而反射电场的极化在上进行时的反射系数;应于入射电场的极化在上进行,而反射电场的极化在上进行时的反射系数;对应于入射和反射电场的极化都在上进行时的反射系数;
反射系数的计算公式为:
式中:Eim,Esm为当发射天线和接收天线均在线性极化上进行时,散射单元的散射场,为了表示方便,用表示对应的
所述步骤四:获得接收天线接收到的第m个散射单元的散射场,表示为:
式中表示接收天线的极化分量;
令 则:
所述步骤五:将接收天线收到的每个散射单元的散射场进行叠加,得到总的散射场:
所述步骤六:根据公式(11)计算出总散射场,进而得到近场回波强度和功率特性;由于划分的散射单元的尺寸足够小,引信处于面元的散射远场区,则散射单元的回波强度可以用远场RCS表式为:
引信天线接收的回波功率为:
式中Pi为发射天线的辐射功率,Pr为接收天线处的接收功率,R为天线相位中心到场点的距离;λ为雷达发射信号频率;Di(Ds)为发射(接收)天线的增益;fi(fs)为发射(接收)天线的方向性函数;σ为公式(12)计算得到的近场回波强度。
本发明所取得的有益效果为:
本发明所述一种准确的近场回波获取方法,充分考虑到近场的特点,如天线方向图、球面波、距离、极化等影响,利用PO+PTD+GO法建立了复杂目标近场散射模型,并引入近场极化散射矩阵的概念,将不同机理的散射特性用统一表达式表示,简化了计算过程;其中将入射波表达式表示成球面波形式,解决了球面波照射和距离等问题;考虑天线方向图的影响,通过建立天线坐标系和天线局部坐标系描述天线方向图影响下的近场电磁散射,解决了天线方向图、局部照射和极化等问题,并给出复杂目标近场电磁散射理论模型与实测数据的校验,通过仿真能够得到较为准确的近场回波数据,为引信半实物仿真和数字仿真提供可靠的数据输入。
附图说明
图1为本发明所述一种准确的近场回波获取方法流程图;
图2为目标坐标系、天线坐标系及散射单元坐标系示意图;
图3为某目标的近场回波强度结果对比图;
图4为采用不同方向图时平板的散射强度随距离变化曲线图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
如图1所示,本发明所述一种准确的近场回波获取方法包括如下步骤:
步骤一:建立目标坐标系O-XtYtZt,天线坐标系O-XaYaZa,散射单元坐标系O-XmYmZm,如图2所示;
步骤二:在散射体表面,获得入射到第m个散射单元上的场;
考虑到近场的特点,在散射体表面,将近场入射波表示成球面波形式,图2中入射到第m个散射单元的场强为:
式中:i表示发射;Pi为发射的信号功率;Di为发射天线的增益;Z0为空气的波阻抗,且Z0=120πΩ;Rim是发射天线相位中心与散射单元中心之间的距离;ω=2πf为电磁振荡的角频率;fim为第m个散射单元处的天线方向性函数;φ为场的相位特性;E0m为第m个散射单元处电场的极化特性;k表示电磁波在自由空间内的波数;t表示时间参数;f为载频频率。
入射到第m个散射单元的场强可以写成以下矩阵形式:
式中表示发射天线的极化分量。
步骤三:获得散射单元的散射场,表示为:
其中:
式中:Esm为散射电场强度,Esm1、Esm2为Esm的两个分量;Rim为散射单元中心到发射天线相位中心的距离;Rsm为散射单元中心到接收天线相位中心的距离;为散射单元的极化散射矩阵,其四个分量是四种反射系数,其中对应于入射和反射电场的极化都在上进行时的反射系数;对应于入射电场的极化在上进行,而反射电场的极化在上进行时的反射系数;应于入射电场的极化在上进行,而反射电场的极化在上进行时的反射系数;对应于入射和反射电场的极化都在上进行时的反射系数;而公式(6)所定义的反射系数可以将物理光学反射系数、边缘绕射场和多次反射场的反射系数统一表示,这样极大简化了电磁散射的计算过程,并能从理论上将不同机理的散射特性用一种统一的观点去描述,这是本发明的一个特点和贡献。
反射系数的计算公式为:
式中:Eim,Esm为当发射天线和接收天线均在线性极化上进行时,散射单元的散射场,为了表示方便,本文都用表示对应的
步骤四:获得接收天线接收到的第m个散射单元的散射场,表示为:
式中表示接收天线的极化分量;
令 则:
步骤五:将接收天线收到的每个散射单元的散射场进行叠加,便得到总的散射场:
步骤六:根据公式(11)计算出总散射场,进而得到近场回波强度和功率特性。
由于划分的散射单元的尺寸足够小,引信处于面元的散射远场区,则散射单元的回波强度可以用远场RCS表式为:
引信天线接收的回波功率为:
式中Pi为发射天线的辐射功率,Pr为接收天线处的接收功率,R为天线相位中心到场点的距离;λ为雷达发射信号频率;Di(Ds)为发射(接收)天线的增益;fi(fs)为发射(接收)天线的方向性函数;σ为公式(12)计算得到的近场回波强度。
如图3所示,为根据本发明计算出的某目标的近场回波测试结果仿真曲线,图中实线为实际测量结果,虚线为采用本发明计算的结果,从图中可以看出两种结果符合较好。
如图4所示,为采用不同方向图时平板的近场回波强度随距离变化曲线,图中虚线为使用刀状波束天线方向图的结果,实线为不考虑天线方向图的结果,从图中可以看出使用刀状波束天线方向图的结果明显小于不考虑天线方向图,即平面波完全照射的结果,这是因为在近场条件下考虑方向图时目标处于局部照射,不考虑方向图时目标处于完全照射。从计算结果上看,天线方向图对近场散射的影响是不容忽视的。
Claims (1)
1.一种准确的近场回波获取方法,其特征在于:该方法包括如下步骤:
步骤一:建立目标坐标系O-XtYtZt,天线坐标系O-XaYaZa,散射单元坐标系O-XmYmZm;
步骤二:在散射体表面,获得入射到第m个散射单元上的场;
步骤三:获得散射单元的散射场;
步骤四:获得接收天线接收到的第m个散射单元的散射场;
步骤五:将接收天线收到的每个散射单元的散射场进行叠加,得到总的散射场;
步骤六:计算出总散射场,进而得到近场回波强度和功率特性;
所述步骤二的具体实现方法如下:在散射体表面,将近场入射波表示成球面波形式,入射到第m个散射单元的场强为:
式中:i表示发射;Pi为发射的信号功率;Di为发射天线的增益;Z0为空气的波阻抗,且Z0=120πΩ;Rim是发射天线相位中心与散射单元中心之间的距离;ω=2πf为电磁振荡的角频率;fim为第m个散射单元处的天线方向性函数;φ为场的相位特性;E0m为第m个散射单元处电场的极化特性;k表示电磁波在自由空间内的波数;t表示时间参数;f为载频频率;
入射到第m个散射单元的场强写成以下矩阵形式:
式中表示发射天线的极化分量;
所述步骤三中获得散射单元的散射场,表示为:
其中:
式中:Esm为散射电场强度,Esm1、Esm2为Esm的两个分量;Rim为散射单元中心到发射天线相位中心的距离;Rsm为散射单元中心到接收天线相位中心的距离;为散射单元的极化散射矩阵,其四个分量是四种反射系数,其中对应于入射和反射电场的极化都在上进行时的反射系数;对应于入射电场的极化在上进行,而反射电场的极化在上进行时的反射系数;应于入射电场的极化在上进行,而反射电场的极化在上进行时的反射系数;对应于入射和反射电场的极化都在上进行时的反射系数;
反射系数的计算公式为:
式中:Eim,Esm为当发射天线和接收天线均在线性极化上进行时,散射单元的散射场,为了表示方便,用表示对应的
所述步骤四:获得接收天线接收到的第m个散射单元的散射场,表示为:
式中表示接收天线的极化分量;
令 则:
所述步骤五:将接收天线收到的每个散射单元的散射场进行叠加,得到总的散射场:
所述步骤六:根据公式(11)计算出总散射场,进而得到近场回波强度和功率特性;由于划分的散射单元的尺寸足够小,引信处于面元的散射远场区,则散射单元的回波强度可以用远场RCS表式为:
引信天线接收的回波功率为:
式中Pi为发射天线的辐射功率,Pr为接收天线处的接收功率,R为天线相位中心到场点的距离;λ为雷达发射信号频率;Di(Ds)为发射(接收)天线的增益;fi(fs)为发射(接收)天线的方向性函数;σ为公式(12)计算得到的近场回波强度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210362109.0A CN103675781B (zh) | 2012-09-25 | 2012-09-25 | 一种准确的近场回波获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210362109.0A CN103675781B (zh) | 2012-09-25 | 2012-09-25 | 一种准确的近场回波获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675781A CN103675781A (zh) | 2014-03-26 |
CN103675781B true CN103675781B (zh) | 2016-05-18 |
Family
ID=50313922
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210362109.0A Active CN103675781B (zh) | 2012-09-25 | 2012-09-25 | 一种准确的近场回波获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675781B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106526547A (zh) * | 2016-12-09 | 2017-03-22 | 北京环境特性研究所 | 基于InSAR技术的直线扫描近场RCS测试杂波抑制方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166126B (zh) * | 2014-07-21 | 2016-08-24 | 西安空间无线电技术研究所 | 一种用于连续波雷达的回波信号模拟方法 |
CN105403881A (zh) * | 2015-12-14 | 2016-03-16 | 北京环境特性研究所 | 基于大面元模型的引信多普勒频率获取方法 |
CN106093872B (zh) * | 2016-06-02 | 2018-09-18 | 中国舰船研究设计中心 | 一种考虑近海面面环境影响的编队内雷达远场预测方法 |
CN106772290B (zh) * | 2016-12-09 | 2019-01-15 | 中国人民解放军国防科学技术大学 | 一种全极化双基地雷达目标动态回波模拟方法 |
CN108061883B (zh) * | 2017-11-28 | 2021-06-04 | 上海无线电设备研究所 | 局部散射源反演的近场散射函数转换弹目交会回波的方法 |
CN108445303B (zh) * | 2018-03-08 | 2020-06-26 | 北京环境特性研究所 | 一种近场电磁散射特性模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4172255A (en) * | 1977-08-08 | 1979-10-23 | The United States Of America As Represented By The Secretary Of The Interior | HF coastal current mapping radar system |
CN102401893B (zh) * | 2010-09-10 | 2015-03-04 | 上海无线电设备研究所 | 一种近场目标对称极化rcs测试方法 |
-
2012
- 2012-09-25 CN CN201210362109.0A patent/CN103675781B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106526547A (zh) * | 2016-12-09 | 2017-03-22 | 北京环境特性研究所 | 基于InSAR技术的直线扫描近场RCS测试杂波抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103675781A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675781B (zh) | 一种准确的近场回波获取方法 | |
CN108445303B (zh) | 一种近场电磁散射特性模拟方法 | |
CN104063544B (zh) | 近场效应误差分析方法 | |
CN104733850B (zh) | 基于全息调制的人工张量阻抗表面天线及其实现方法 | |
CN106886020A (zh) | 一种单天线接收条件下的电磁涡旋成像方法 | |
CN102401893A (zh) | 一种近场目标对称极化rcs测试方法 | |
CN109031226B (zh) | 基于机电耦合的变形有源相控阵雷达探测性能快速评估方法 | |
CN103064072B (zh) | 雷达散射截面测量中的背景信号提取方法 | |
CN105738891B (zh) | 一种机载数字阵列雷达对弱机动目标角度进行跟踪的方法 | |
CN105372633B (zh) | 一种相控阵雷达降维四通道抗主副瓣干扰的方法 | |
CN102608599B (zh) | 一种超宽带多发多收阵列配置方法 | |
CN104391183A (zh) | 近场测量天线远区场特性快速计算方法 | |
CN103033812B (zh) | 机载双基地sar的波束追赶空间同步方法 | |
CN106707255A (zh) | 一种相控阵雷达仿真系统及方法 | |
CN104535973A (zh) | 一种机载预警雷达目标检测方法 | |
CN107333290A (zh) | 一种基站天线电性能参数监测方法 | |
CN105044688B (zh) | 基于迭代子空间跟踪算法的雷达稳健空时自适应处理方法 | |
CN110376560A (zh) | 一种基于单距离门的机载双基地mimo雷达幅相误差校正方法 | |
CN109520383B (zh) | 基于matlab的体目标回波模拟方法 | |
Pola et al. | Model comparison of bistatic radar clutter | |
CN104992035A (zh) | 一种太赫兹频段表面粗糙目标电磁散射的快速计算方法 | |
CN104237867A (zh) | 引信接收机的有效多普勒带宽的确定方法和装置 | |
Borzov et al. | Mathematical modeling and simulation of the input signals of short-range radar systems | |
CN103729485B (zh) | 一种基于dem数据的宽带雷达相干杂波仿真方法 | |
Bae et al. | Automatic target recognition with unknown orientation and adaptive waveforms |
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 |