CN109520383A - 基于matlab的体目标回波模拟方法 - Google Patents

基于matlab的体目标回波模拟方法 Download PDF

Info

Publication number
CN109520383A
CN109520383A CN201710853067.3A CN201710853067A CN109520383A CN 109520383 A CN109520383 A CN 109520383A CN 201710853067 A CN201710853067 A CN 201710853067A CN 109520383 A CN109520383 A CN 109520383A
Authority
CN
China
Prior art keywords
face element
echo
body target
signal
moment
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
Application number
CN201710853067.3A
Other languages
English (en)
Other versions
CN109520383B (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201710853067.3A priority Critical patent/CN109520383B/zh
Publication of CN109520383A publication Critical patent/CN109520383A/zh
Application granted granted Critical
Publication of CN109520383B publication Critical patent/CN109520383B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42—AMMUNITION; BLASTING
    • F42C—AMMUNITION FUZES; ARMING OR SAFETY MEANS THEREFOR
    • F42C21/00—Checking fuzes; Testing fuzes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种基于matlab的体目标回波模拟方法,其特征在于，包括如下步骤：(10)体目标3D建模：在matlab中将体目标划分成多个面元，并对每一个面元在追踪过程中是否有效进行判定，得到有效面元；(20)面元RCS计算：将每一有效面元等效为金属平板，通过物理光学法得到每一个有效面元的RCS；(30)体目标回波生成：在各个仿真时刻将所有有效面元产生的回波信号进行叠加后得到该仿真时刻的回波信号，再将所有仿真时刻的信号时域叠加，得到所述体目标回波模拟信号。本发明的基于matlab的体目标回波模拟方法,模拟信号更加接近真实回波信号。

Description

基于matlab的体目标回波模拟方法
技术领域
本发明属于引信目标回波模拟技术领域，特别是一种模拟信号更加接近真实回波信号的基于matlab的体目标回波模拟方法。
背景技术
随着无线电引信的不断进步与发展，尤其是在其研制调试阶段，对引信的全面测试就显得非常重要。而采用回波模拟器进行测试可以大大减少对靶场实验的依赖性，极大的节省了人力物力和财力。
但以往在模拟体目标回波时，通常是将体目标近似成N个点，每个点看成是点目标并产生回波信号，最后进行矢量叠加。在系统仿真学报的引信目标回波模拟器的建模与实现的论文中，详细的讲述了该方法的具体过程。
但该方法有以下问题：在弹目交汇过程中，自始至终都通过体目标近似的N个点形成回波，并不准确。因为在不同的姿态角之下，目标部分部位会被遮挡，并且目标离导弹距离很近的情况下，引信很有可能只是照射到目标的某一部分。因此该方法产生的回波模拟器可能会造成信号功率过大，影响到最后的起爆判断，不能形成最大的杀伤力。
总之，现有技术存在的问题是：对体目标的回波模拟误差大、不真实。
发明内容
本发明的目的在于提供一种基于matlab的体目标回波模拟方法,模拟信号更加接近真实回波信号。
实现本发明目的的技术解决方案为：
一种基于matlab的体目标回波模拟方法，包括如下步骤：
(10)体目标3D建模：在matlab中将体目标划分成多个面元，并对每一个面元在追踪过程中是否有效进行判定，得到有效面元；
(20)面元RCS计算：将每一有效面元等效为金属平板，通过物理光学法得到每一个有效面元的RCS；
(30)体目标回波生成：在各个仿真时刻将所有有效面元产生的回波信号进行矢量叠加后得到该仿真时刻的回波信号，再将所有仿真时刻的信号时域叠加，得到所述体目标回波模拟信号。
本发明与现有技术相比，其显著优点为：
在本发明中，加入了面元的判定，可以有效的选取出在交汇过程中真正被引信照射到的部分。并且，在对每一个面元的RCS计算方法中，将每一个面元近似为一个金属平板，通过物理光学法计算得出的RCS，不仅不需考虑引信的极化方式，并且不需要考虑因面元法而产生的虚构劈的处理，计算更为简单，但对RCS计算的准确度稍微不足。
本发明的回波模拟方法，可较准确的模拟出，在弹目逐渐接近的过程中，照射区域由大到小而产生的回波，更为真实准确。
下面结合附图和具体实施方式对本发明作进一步的详细描述。
附图说明
图1是本发明基于matlab的体目标回波模拟方法的主流程图。
图2是图1中体目标3D建模步骤的流程图。
图3是图1中体目标回波生成步骤的流程图。
图4是基于面元划分原则构建的PAC-3导弹模型。
图5是基于有效面元判定的仿真实例。
图6是基于比例导引法仿真出的追踪过程图示。
图7是本发明在追踪过程中的回波模型。
具体实施方式
如图1所示，本发明基于matlab的体目标回波模拟方法,在matlab上对体目标进行建模，求得第k个仿真时刻每一有效面元的RCS，之后根据延时信号，求得该仿真时刻的总的回波信号，最后将这些仿真时刻的回波信号时域叠加得到最终的回波信号。全过程通过matlab实现。包括如下步骤：
(10)体目标3D建模：在matlab中将体目标划分成多个面元，并对每一个面元在追踪过程中是否有效进行判定，得到有效面元；
如图2所示，所述(10)体目标3D建模步骤包括：
(11)面元划分：在满足远场条件的前提下，将属于电大尺寸复杂目标的体目标剖分成若干电小尺寸的面元；
根据面元划分原则，在微波段的雷达散射截面理论定义式中，目标必须满足远场条件，即：
“式(1)”中R为弹目距离，D为目标的最大线性尺寸，λ为引信的工作波长，2D2/λ为引信与目标之间的最短距离Rmtn的值。根据“式(1)”计算出满足远场条件时目标面元的最大线性尺寸D的值。再根据“式(2)”可以得到目标表面面元划分数目N的大致范围。
参考图4，PAC-3导弹全长5.2m，弹翼展0.5m，尾翼展0.6m。根据面元划分准则，将其分成1438个点。
(12)面元有效判定：对面元划分得到的每一个面元，首先去除后向面面元，之后去除波束照射外的面元，最后一步去除被遮挡面元，得到有效面元。
在每一个仿真时刻，都要对目标进行面元判断，得到第k个仿真时刻能被引信照射到的Mk个面元。判断方法如下：
后向面判别：
在每一个仿真时刻，目标与雷达引信间的角度关系是确定的，因此能够根据目标表面上每一个面元的外法向量与每一个面元相对雷达引信的角度进行判断。若面元的外法向量与入射波方向矢量夹角大于90°时，该面元能被入射光束照射到；当面元的外法线单位向量与入射波方向矢量夹角小于或等于90°时，说明该面元面向相对引信的反方向，该面元则不能被入射光束照射到。
波束内判别：
在近距离追踪过程中，目标不一定整个都在雷达引信的波束照射下，因此，根据引信的照射中心方向矢量与每一个面元相对引信的角度进行判断。若面元相对引信的角度与引信的照射中心方向矢量夹角小于3dB波束宽度时，该面元能被入射光束照射到；当面元相对引信的角度与引信的照射中心方向矢量夹角大于3dB波束宽度时，说明该面元没有被引信所探测到。
遮挡判断：
在对目标进行上述处理后，需要对处理过后剩下的n个面元进行互遮挡判断。首先将这n个面元依次计算到引信中心的距离值并按距离值从小到大进行排序，即第n个面元的距离值最大。判断第n个面元中心点到引信的连线与其他n-1个面元是否相交。若存在交点则第n个面元判定被遮挡，若不存在交点则判定第n个面元无遮挡，n＝n-1，跳转上一步。当n＝0时遮挡判断完成。
根据上述方法，得到有效面元示意图，如图5所示。由图可知，当弹目之间距离较近时，引信真正能照射到的部位并不多。如果按照原有方法，则得到的回波模拟信号功率将比真实信号要大得多。因此，该方法能更加准确的模拟出回波信号的真实功率。
(20)面元RCS计算：将每一有效面元等效为金属平板，通过物理光学法得到每一个有效面元的RCS；
得到第k个仿真时刻的有效面元Mk后，将每一个面元模拟成一个金属平板，通过物理光学法计算，则可以得到每一个面元的RCS.
计算每一面元的RCS时将面元等效为金属平板，所述(20)面元RCS计算，按下式计算各有效面元的RCS：
式中，σk(m)为第k个仿真时刻第m个面元的RCS，a为面元边长，λ为波长，θ为引信波束中心线与面元法向量之间的夹角。
(30)体目标回波生成：在各仿真时刻将所有有效面元产生的回波信号进行叠加后得到该仿真时刻的回波信号，再将所有仿真时刻的信号时域叠加，得到所述体目标回波模拟信号。
如图3所示，所述(30)体目标回波生成步骤包括：
(31)面元回波信号产生：按下式生成第k个仿真时刻第m个面元的回波信号：
式中，ULMk(m)为第k个仿真时刻第m个面元的回波幅度,dt为仿真步长，τk(m为第k个仿真时刻第m个面元的回波延时，ω0为载波角速度，θ0为初始相位，为矩形脉冲信号，NT为脉冲数，T为脉冲重复周期。
其中，
式中，Rk(m)为第m个面元第k个仿真时刻的弹目距离,c为光速，Pt为发射功率，Ls为系统损耗，gvigvr为天线的发射增益和接受增益，λ为波长，σk(m)为第m个面元第k个仿真时刻的RCS，k1为信号在空气中传播衰减系数。
(32)第k个仿真时刻的回波信号生成：按下式为第k个仿真时刻所有有效面元产生的信号叠加后的回波信号，
式中Mk为第k个仿真时刻有效面元数。
(33)体目标回波模拟信号生成：按下式得到体目标回波模拟信号，
U＝Uk|k＝1：N，
式中，N为仿真时刻总数。
实施例
追踪参数：
目标为PAC-3导弹，全长5.2m，弹翼展0.5m，尾翼展0.6m。
追踪雷达运动速度设为600，初始位置为(-100，0，-100)，目标PAC-3导弹运动速度设为420，初始位置为(0，0，0)，波长1cm，脉冲重复周期为0.01s，脉冲宽度为0.003s。图6所示为比例导引法追踪过程的示意图。图7为在追踪过程中的回波模拟信号。
从图6、图7可以看出，在图6所示比例导引法追踪的过程中，弹目逐渐交汇，但图7所示回波信号幅度与距离的三次方并非严格按照反比例函数的关系，在大概0.3s前后信号幅度有凹陷，在最后接近的过程中，信号幅度的增大也趋于缓和。这说明，在追踪过程中，前后时刻照射部位不同，有可能会导致RCS较大的变化，最终造成回波信号功率的变化。而以往的回波模拟方法并未考虑这一点。因此，本发明所述的体目标回波模拟方法较以往方法更加真实准确。

Claims (4)

1.一种基于matlab的体目标回波模拟方法,其特征在于，包括如下步骤：
(10)体目标3D建模：在matlab中将体目标划分成多个面元，并对每一个面元在追踪过程中是否有效进行判定，得到有效面元；
(20)面元RCS计算：将每一有效面元等效为金属平板，通过物理光学法得到每一个有效面元的RCS；
(30)体目标回波生成：在各个仿真时刻将所有有效面元产生的回波信号进行叠加后得到该仿真时刻的回波信号，再将所有仿真时刻的信号时域叠加，得到所述体目标回波模拟信号。
2.根据权利要求1所述的体目标回波模拟方法,其特征在于在matlab上对体目标进行建模和判定，所述(10)体目标3D建模步骤包括：
(11)面元划分：在满足远场条件的前提下，将属于电大尺寸复杂目标的体目标剖分成若干电小尺寸的面元；
(12)面元有效判定：对面元划分得到的每一个面元，首先去除后向面面元，之后去除波束照射外的面元，最后一步去除被遮挡面元，得到有效面元。
3.根据权利要求1所述的体目标回波模拟方法,其特征在于，所述(20)面元RCS计算，按下式计算各有效面元的RCS：
式中，σk(m)为第k个仿真时刻第m个面元的RCS，a为面元边长，λ为波长，θ为引信波束中心线与面元法向量之间的夹角。
4.根据权利要求3所述的体目标回波模拟方法,所述(30)体目标回波生成步骤包括：
(31)面元回波信号产生：按下式生成第k个仿真时刻第m个面元的回波信号：
式中，ULMk(m)为第k个仿真时刻第m个面元的回波幅度,dt为仿真步长，τk(m)为第k个仿真时刻第m个面元的回波延时，ω0为载波角速度，θ0为初始相位，为矩形脉冲信号，NT为脉冲数，T为脉冲重复周期；
其中，
式中，Rk(m)为第m个面元第k个仿真时刻的弹目距离,c为光速，Pt为发射功率，Ls为系统损耗，gvigvr为天线的发射增益和接受增益，λ为波长，σk(m)为第m个面元第k个仿真时刻的RCS，k1为信号在空气中传播衰减系数；
(32)第k个仿真时刻的回波信号生成：按下式为第k个仿真时刻所有有效面元产生的信号叠加后的回波信号，
式中Mk为第k个仿真时刻有效面元数；
(33)体目标回波模拟信号生成：按下式得到体目标回波模拟信号，
U＝Uk|k＝1：N，
式中，N为仿真时刻总数。
CN201710853067.3A 2017-09-20 2017-09-20 基于matlab的体目标回波模拟方法 Active CN109520383B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710853067.3A CN109520383B (zh) 2017-09-20 2017-09-20 基于matlab的体目标回波模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710853067.3A CN109520383B (zh) 2017-09-20 2017-09-20 基于matlab的体目标回波模拟方法

Publications (2)

Publication Number Publication Date
CN109520383A true CN109520383A (zh) 2019-03-26
CN109520383B CN109520383B (zh) 2020-07-31

Family

ID=65767808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710853067.3A Active CN109520383B (zh) 2017-09-20 2017-09-20 基于matlab的体目标回波模拟方法

Country Status (1)

Country Link
CN (1) CN109520383B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110095035A (zh) * 2019-04-22 2019-08-06 南京理工大学 导弹引信探测仿真方法
CN110991030A (zh) * 2019-11-28 2020-04-10 北京电子工程总体研究所 一种通用复杂引战系统仿真的方法及系统
CN110991031A (zh) * 2019-11-28 2020-04-10 北京电子工程总体研究所 一种多支路引信目标探测建模的方法
CN113030900A (zh) * 2021-03-26 2021-06-25 中国人民解放军国防科技大学 基于面元分布的动态匹配反射系数缩比测量方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879768A (zh) * 2012-09-14 2013-01-16 北京航空航天大学 一种基于稳态rcs的星载sar高保真回波仿真方法
CN103176179A (zh) * 2013-02-22 2013-06-26 清华大学 Sar雷达回波实时模拟方法
CN103513232A (zh) * 2012-06-26 2014-01-15 北京华清瑞达科技有限公司 一种引信体目标模拟方法及系统
CN103777186A (zh) * 2013-11-02 2014-05-07 西安电子科技大学 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法
CN105486180A (zh) * 2015-12-15 2016-04-13 西安电子科技大学 基于波束分解和局部照射的激光引信近场回波功率计算

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513232A (zh) * 2012-06-26 2014-01-15 北京华清瑞达科技有限公司 一种引信体目标模拟方法及系统
CN102879768A (zh) * 2012-09-14 2013-01-16 北京航空航天大学 一种基于稳态rcs的星载sar高保真回波仿真方法
CN103176179A (zh) * 2013-02-22 2013-06-26 清华大学 Sar雷达回波实时模拟方法
CN103777186A (zh) * 2013-11-02 2014-05-07 西安电子科技大学 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法
CN105486180A (zh) * 2015-12-15 2016-04-13 西安电子科技大学 基于波束分解和局部照射的激光引信近场回波功率计算

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110095035A (zh) * 2019-04-22 2019-08-06 南京理工大学 导弹引信探测仿真方法
CN110095035B (zh) * 2019-04-22 2021-06-29 南京理工大学 导弹引信探测仿真方法
CN110991030A (zh) * 2019-11-28 2020-04-10 北京电子工程总体研究所 一种通用复杂引战系统仿真的方法及系统
CN110991031A (zh) * 2019-11-28 2020-04-10 北京电子工程总体研究所 一种多支路引信目标探测建模的方法
CN110991030B (zh) * 2019-11-28 2023-10-13 北京电子工程总体研究所 一种通用复杂引战系统仿真的方法及系统
CN110991031B (zh) * 2019-11-28 2023-11-14 北京电子工程总体研究所 一种多支路引信目标探测建模的方法
CN113030900A (zh) * 2021-03-26 2021-06-25 中国人民解放军国防科技大学 基于面元分布的动态匹配反射系数缩比测量方法和装置

Also Published As

Publication number Publication date
CN109520383B (zh) 2020-07-31

Similar Documents

Publication Publication Date Title
CN109520383A (zh) 基于matlab的体目标回波模拟方法
CN105093215B (zh) 基于多普勒信息的雷达对低空慢速小目标的跟踪方法
CN107765226A (zh) 一种sar卫星雷达回波模拟方法、系统和介质
CN104597430B (zh) 机载气象雷达低空风切变三维回波仿真方法
US11899130B2 (en) Method and device for determining a radar cross section, method for training an interaction model, and radar target emulator and test facility
CN103941233B (zh) 一种多平台主被动传感器协同跟踪的雷达间歇交替辐射控制方法
CN104614713A (zh) 一种适合于艇载雷达系统的雷达回波信号模拟器
CN105467369B (zh) 一种目标回波仿真方法和装置
CN104076355A (zh) 基于动态规划的强杂波环境中弱小目标检测前跟踪方法
CN104076338B (zh) 基于数字高程和数字地表覆盖的机载雷达杂波仿真方法
CN107607943A (zh) 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法
CN103093057A (zh) 一种船舶导航雷达信号仿真方法
CN103675769A (zh) 基于分布式接收机的斜视sar欺骗干扰方法
CN105260524B (zh) 一种扫描状态下船舶导航雷达二维回波序列像仿真方法
CN108459308A (zh) 一种基于时变rcs数据的回波模拟方法及装置
CN107255805A (zh) 基于加权最小二乘的雷达目标rcs的预测方法
CN105738887A (zh) 基于多普勒通道划分的机载雷达杂波功率谱的优化方法
CN109444837A (zh) 一种雷达回波仿真方法和装置
CN105372653B (zh) 面向岸基空管雷达系统中一种高效转弯机动目标跟踪方法
CN103954949B (zh) 一种隐身滑跃式机动目标的多模型椭圆Hough变换积累检测方法
CN106556833A (zh) 基于时域弹跳射线法快速近场计算的isar成像仿真方法
CN112098999A (zh) 一种高动态雷达导引头掠海目标电磁信号建模方法
CN107356906A (zh) 一种利用布儒斯特效应的超低空目标镜像干扰抑制方法
CN104101869B (zh) 一种极坐标下的地波雷达运动目标仿真建模方法
CN111965619A (zh) 一种基于射线追踪的复杂目标时域散射信号仿真方法

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