CN113030900B - 基于面元分布的动态匹配反射系数缩比测量方法和装置 - Google Patents

基于面元分布的动态匹配反射系数缩比测量方法和装置 Download PDF

Info

Publication number
CN113030900B
CN113030900B CN202110325787.9A CN202110325787A CN113030900B CN 113030900 B CN113030900 B CN 113030900B CN 202110325787 A CN202110325787 A CN 202110325787A CN 113030900 B CN113030900 B CN 113030900B
Authority
CN
China
Prior art keywords
model
coating
rcs
reflection coefficient
scaling
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
Application number
CN202110325787.9A
Other languages
English (en)
Other versions
CN113030900A (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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202110325787.9A priority Critical patent/CN113030900B/zh
Publication of CN113030900A publication Critical patent/CN113030900A/zh
Application granted granted Critical
Publication of CN113030900B publication Critical patent/CN113030900B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S7/411Identification of targets based on measurements of radar reflectivity
    • G01S7/412Identification of targets based on measurements of radar reflectivity based on a comparison between measured values and known or stored values

Landscapes

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

Abstract

本申请涉及一种基于面元分布的动态匹配反射系数缩比测量方法、装置、计算机设备和存储介质。所述方法包括:通过获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值,根据金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及涂覆缩比模型的涂覆材料的第一反射系数信息,得到金属缩比模型在入射波角度下的RCS面元分布信息;获取目标原型涂覆材料的第二反射系数信息,根据第二反射系数信息、面元分布信息和预设的几何缩比因子,得到目标原型的RCS值。本发明有利于提高反演精度,扩大了缩比测量的适用范围,为太赫兹频段缩比测量的材料色散问题提供了有效的解决方法。

Description

基于面元分布的动态匹配反射系数缩比测量方法和装置
技术领域
本申请涉及雷达技术领域,特别是涉及一种面元分布的动态匹配反射系数缩比测量方法、装置、计算机设备和存储介质。
背景技术
近年来,随着对超大尺寸目标缩比测量需求的增加,对缩比测量的频段和缩比反演方法提出了更高的要求。为了实现大缩比因子的缩比测量,缩比频率上升到太赫兹频段,材料的色散特性不容忽略,经典电磁相似率的条件无法满足,如何将太赫兹频段的缩比测量结果反演得到目标原型的RCS成为亟待解决的问题。
针对缩比测量中的材料色散问题,现有的解决方法主要包括:放宽相似法则法、待定物理缩比因子法、寻找替代材料法。现有方法中,寻找替代材料法涉及到替代材料的选取或者合成,实现起来周期长、成功率不高,比较困难;放宽相似法则法的有效性局限于高频、远场和后向散射,且用于复杂目标时误差明显增大;待定物理缩比因子法用于反射系数变化趋势相似的两种材料效果较好,当两种材料的反射系数变化趋势不同时,该方法出现较大误差。放宽相似法则法和待定物理缩比因子法的缺点在于没有充分考虑涂覆材料在诸多影响因素下出现的复杂变化趋势,导致缩比反演方法失效。现有方法忽略了材料反射系数的诸多影响因素(如厚度、损耗正切、角度等)对缩比关系的影响,对反射系数变化复杂的材料涂覆目标以及结构复杂的目标缩比反演精度低。
发明内容
基于此,有必要针对上述技术问题,提供一种能够提高由缩比模型RCS测量值反演出目标原型RCS值的反演精度的基于面元分布的动态匹配反射系数缩比测量方法、装置、计算机设备和存储介质。
一种基于面元分布的动态匹配反射系数缩比测量方法,所述方法包括:
获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同;
根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;
根据所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述第一反射系数信息,得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息;
获取目标原型涂覆材料的第二反射系数信息,根据所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,得到所述目标原型的RCS值。
在其中一个实施例中,还包括:获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型的个数为面元的类别数减1,每个所述涂覆缩比模型的涂覆材料不同;所述面元根据法向方向分类,法向方向相同的面元为一类。
在其中一个实施例中,还包括:根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;所述第一反射系数信息随着雷达入射波与所述涂覆缩比模型的不同面元形成的夹角大小变化而变化。
在其中一个实施例中,还包括:获取所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述第一反射系数信息;
将所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值和所述第一反射系数信息代入以下公式:
Figure BDA0002994611560000031
其中,
Figure BDA0002994611560000032
表示第1…m-1个涂覆缩比模型的RCS测量值,
Figure BDA0002994611560000033
表示雷达入射波的入射角;
Figure BDA0002994611560000034
表示所述金属缩比模型RCS测量值;m表示面元的类别数;θj,j=1…m表示入射波与第j类面元的法向方向形成的夹角,Rkj),k=1…m-1,j=1…m表示入射波与第k个涂覆缩比模型第j类面元的法向方向形成的夹角为θj时的反射系数;
Figure BDA0002994611560000035
表示面元法向为
Figure BDA0002994611560000036
的所有面元复RCS的叠加;
求解得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息
Figure BDA0002994611560000037
在其中一个实施例中,还包括:获取目标原型涂覆材料的第二反射系数信息;
将所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,代入以下公式:
Figure BDA0002994611560000038
Figure BDA0002994611560000039
其中,
Figure BDA00029946115600000310
表示所述目标原型的RCS值,Rprotj),j=1…m表示入射波与目标原型第j类面元的法向方向形成的夹角为θj时的反射系数,s为预设的几何缩比因子;
根据公式计算结果得到所述目标原型的RCS值。
在其中一个实施例中,还包括:当所述涂覆缩比模型的面元分为两类时,获取所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述涂覆缩比模型的涂覆材料的第一反射系数信息;
将所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值和所述第一反射系数信息代入以下公式:
Figure BDA0002994611560000041
求解得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息
Figure BDA0002994611560000042
在其中一个实施例中,还包括:获取目标原型涂覆材料的第二反射系数信息;
将所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,代入以下公式:
Figure BDA0002994611560000043
Figure BDA0002994611560000044
根据公式计算结果得到所述目标原型的RCS值。
一种基于面元分布的动态匹配反射系数缩比测量装置,所述装置包括:
缩比模型RCS测量值获取模块,用于获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同;
第一反射系数信息获取模块,用于根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;
RCS面元分布信息获取模块,用于根据所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述第一反射系数信息,得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息;
目标原型的RCS值获取模块,用于获取目标原型涂覆材料的第二反射系数信息,根据所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,得到所述目标原型的RCS值。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同;
根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;
根据所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述第一反射系数信息,得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息;
获取目标原型涂覆材料的第二反射系数信息,根据所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,得到所述目标原型的RCS值。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同;
根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;
根据所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述第一反射系数信息,得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息;
获取目标原型涂覆材料的第二反射系数信息,根据所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,得到所述目标原型的RCS值。
上述基于面元分布的动态匹配反射系数缩比测量方法、装置、计算机设备和存储介质,通过获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值,根据金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及涂覆缩比模型的涂覆材料的第一反射系数信息,得到金属缩比模型在入射波角度下的RCS面元分布信息;获取目标原型涂覆材料的第二反射系数信息,根据第二反射系数信息、面元分布信息和预设的几何缩比因子,得到目标原型的RCS值。本发明基于入射波与面元法向的夹角动态匹配反射系数,实现不同角度下反射系数贡献的权重分配,再根据缩比模型与原型的面元分布的相似性,对反射系数进行替换,最终实现缩比模型到目标原型的反演。对于反射系数变化趋势不相同的两种材料,有利于提高反演精度,扩大了缩比测量的适用范围,为太赫兹频段缩比测量的材料色散问题提供了有效的解决方法,同时不增加***和算法的复杂度,稳定性高、普适性好,具有较好的实用性。
附图说明
图1为一个实施例中基于面元分布的动态匹配反射系数缩比测量方法的应用场景图;
图2为一个实施例中发明的原理图;
图3为一个实施例中基于面元分布的动态匹配反射系数缩比测量方法;
图4为一个实施例中SLICY模型俯视图;
图5为一个实施例中SLICY模型观测角度为φ=270°,θ=45°时示意图;
图6为一个实施例中缩比模型A和原型涂敷材料反射系数示意图;
图7为一个实施例中由缩比模型A反演得到的原型的RCS结果显示图;
图8为一个实施例中缩比模型B和原型涂敷材料反射系数示意图;
图9为一个实施例中由缩比模型B反演得到的原型的RCS结果显示图;
图10为一个实施例中基于面元分布的动态匹配反射系数缩比测量装置的结构框图;
图11为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
本申请提供的基于面元分布的动态匹配反射系数缩比测量方法,可以应用于以下应用环境中。通过终端执行一种基于面元分布的动态匹配反射系数缩比测量方法,通过获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值,根据金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及涂覆缩比模型的涂覆材料的第一反射系数信息,得到金属缩比模型在入射波角度下的RCS面元分布信息;获取目标原型涂覆材料的第二反射系数信息,根据第二反射系数信息、面元分布信息和预设的几何缩比因子,得到目标原型的RCS值。其中,终端可以但不限于是各种个人计算机、笔记本电脑、智能手机、平板电脑和便携式可穿戴设备。
在一个实施例中,如图1所示,提供了一种基于面元分布的动态匹配反射系数缩比测量方法,包括以下步骤:
步骤102,获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值。
本发明主要解决雷达散射截面(radar-cross section,RCS)缩比测量中,当材料参数在测量频段发生色散(或者缩比模型材料参数与原型材料参数不一致)时,如何由缩比模型的RCS测量值反演出目标原型RCS。
缩比测量是获取目标RCS的有效方法之一。当对某个全尺寸目标的RCS测量无法实现时,可以对全尺寸目标的缩比模型进行测量,进而根据缩比关系反演出目标原型的RCS。缩比测量的前提条件是构建与原型测量***相似的模型测量***,两个相似的电磁***中对应的物理量需要满足一定的关系,电磁相似率就是对这一系列关系的总结。当材料参数不随频率发生变化时,电磁相似率的条件得以简化,仅需保证模型电磁***的几何缩比因子与波长缩比因子相同即可实现缩比测量。随着待测目标原型尺寸的增大(如超大尺寸的航母等),缩比因子变得更大,缩比测量的频率将上升到太赫兹频段。太赫兹(Terahertz,THz)波通常指频率在0.1THz到10THz之间的电磁波,在该频段材料的色散特性逐步显现,同一种材料在微波频段与太赫兹频段呈现出不同的材料参数(如介电常数),电磁相似率的条件无法满足。对于上升到太赫兹频段的超大尺寸目标的缩比测量,为了使太赫兹频段的缩比测量结果能够高精度反演出目标原型的RCS,需要考虑材料参数发生变化时对RCS缩比关系的影响,并进一步提出缩比模型和原型材料参数不相同时的RCS缩比反演方法。
图2所示为本发明的原理图(这里的目标模型以圆柱为例),
Figure BDA0002994611560000082
表示不同面元的法向,θ1、θ2分别表示入射角度为
Figure BDA0002994611560000081
时,入射波与不同指向的面元的夹角。不同入射角度对应的θ1、θ2不同,相应的反射系数取值也不同。每一涂覆面元的复RCS可以视为涂覆材料的反射系数对相应的金属面元复RCS的调制结果,因此,通过涂覆缩比模型和金属缩比模型的RCS测量值以及材料反射系数,可以获取不同角度下的面元分布。
涂覆缩比模型是对金属缩比模型的表面进行了材料涂覆,涂覆缩比模型可为一个或多个,每个涂覆缩比模型的涂覆材料不同。
步骤104,根据雷达入射波与涂覆缩比模型面元法向形成的夹角获取涂覆材料的第一反射系数信息。
不同入射角度对应的θ1、θ2不同,相应的反射系数取值也不同。
步骤106,根据金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及第一反射系数信息,得到金属缩比模型在入射波角度下的RCS面元分布信息。
在相同的观测角度下,形状相似的缩比模型和原型不同法向面元的分布比例相同,这一相似性是反演涂覆目标原型的重要理论依据。因此,得到金属缩比模型在入射波角度下的RCS面元分布信息,可用于涂覆目标原型的反演。
步骤108,获取目标原型涂覆材料的第二反射系数信息,根据第二反射系数信息、面元分布信息和预设的几何缩比因子,得到目标原型的RCS值。
利用金属缩比模型在不同角度下的面元分布、动态匹配原型的材料反射系数从而获得中间模型的RCS值,该中间模型为假想模型,它与目标原型满足电磁相似率的缩比条件,通过该中间模型以及电磁相似率的RCS缩比关系即可反演出目标原型的RCS。
上述基于面元分布的动态匹配反射系数缩比测量方法中,通过获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值,根据金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及涂覆缩比模型的涂覆材料的第一反射系数信息,得到金属缩比模型在入射波角度下的RCS面元分布信息;获取目标原型涂覆材料的第二反射系数信息,根据第二反射系数信息、面元分布信息和预设的几何缩比因子,得到目标原型的RCS值。本发明基于入射波与面元法向的夹角动态匹配反射系数,实现不同角度下反射系数贡献的权重分配,再根据缩比模型与原型的面元分布的相似性,对反射系数进行替换,最终实现缩比模型到目标原型的反演。对于反射系数变化趋势不相同的两种材料,有利于提高反演精度,扩大了缩比测量的适用范围,为太赫兹频段缩比测量的材料色散问题提供了有效的解决方法,同时不增加***和算法的复杂度,稳定性高、普适性好,具有较好的实用性。
在其中一个实施例中,还包括:获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;涂覆缩比模型是对金属缩比模型的表面进行了材料涂覆,涂覆缩比模型的个数为面元的类别数减1,每个涂覆缩比模型的涂覆材料不同;面元根据法向方向分类,法向方向相同的面元为一类。
多面元的情况下,若面元类别数为m,则需要m-1个缩比模型的RCS值。
在其中一个实施例中,还包括:根据雷达入射波与涂覆缩比模型面元法向形成的夹角获取涂覆材料的第一反射系数信息;第一反射系数信息随着雷达入射波与涂覆缩比模型的不同面元形成的夹角大小变化而变化。
在其中一个实施例中,还包括:获取金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及第一反射系数信息;将金属缩比模型RCS测量值、涂覆缩比模型RCS测量值和第一反射系数信息代入以下公式:
Figure BDA0002994611560000101
其中,
Figure BDA0002994611560000102
表示第1…m-1个涂覆缩比模型的RCS测量值,
Figure BDA0002994611560000103
表示雷达入射波的入射角;
Figure BDA0002994611560000104
表示金属缩比模型RCS测量值;m表示面元的类别数;θj,j=1…m表示入射波与第j类面元的法向方向形成的夹角,
Rkj),k=1…m-1,j=1…m表示入射波与第k个涂覆缩比模型第j类面元的法向方向形成的夹角为θj时的反射系数;
Figure BDA0002994611560000105
表示面元法向为
Figure BDA0002994611560000106
的所有面元复RCS的叠加;
求解得到金属缩比模型在入射波角度下的RCS面元分布信息
Figure BDA0002994611560000107
在其中一个实施例中,还包括:获取目标原型涂覆材料的第二反射系数信息;将第二反射系数信息、面元分布信息和预设的几何缩比因子,代入以下公式:
Figure BDA0002994611560000108
Figure BDA0002994611560000109
其中,
Figure BDA00029946115600001010
表示目标原型的RCS值,Rprotj),j=1…m表示入射波与目标原型第j类面元的法向方向形成的夹角为θj时的反射系数,s为预设的几何缩比因子;
根据公式计算结果得到目标原型的RCS值。
在其中一个实施例中,还包括:当涂覆缩比模型的面元分为两类时,获取金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及涂覆缩比模型的涂覆材料的第一反射系数信息;
将金属缩比模型RCS测量值、涂覆缩比模型RCS测量值和第一反射系数信息代入以下公式:
Figure BDA0002994611560000111
求解得到金属缩比模型在入射波角度下的RCS面元分布信息
Figure BDA0002994611560000112
在其中一个实施例中,还包括:获取目标原型涂覆材料的第二反射系数信息;将第二反射系数信息、面元分布信息和预设的几何缩比因子,代入以下公式:
Figure BDA0002994611560000113
Figure BDA0002994611560000114
根据公式计算结果得到目标原型的RCS值。
应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个具体实施例中,如图3所示,提供一种基于面元分布的动态匹配反射系数缩比测量方法,模型为两面元分布模型,处理过程可分为三部分:第一步主要进行RCS测量,对涂覆缩比模型和金属缩比模型分别进行测量,获取相应模型的复RCS值;第二步结合模型材料的反射系数动态获取不同观测角度下法向面元的分布;第三步引入原型材料反射系数以及几何缩比因子反演得到涂覆原型的RCS。
在另一个具体实施例中,仿真实验设计了涂覆材料参数不相同、材料反射系数变化趋势相同的缩比模型A和目标原型。
仿真实验选取的目标原型为SLICY模型,模型俯视图如图4所示,观测角度为φ=270°,θ=45°时的模型如图5所示。缩比反演的观测角度为φ=270°,θ=0°-180°。缩比模型A和目标原型的参数详见表1。
表1缩比模型A和目标原型的仿真参数
Figure BDA0002994611560000121
缩比模型A和原型涂敷材料反射系数如图6所示。从图中可以看出,两条曲线变化趋势相同,对应于材料参数不同、但反射系数变化趋势相同的情况。采用动态匹配反射系数法由缩比模型A反演得到的原型的RCS结果如图7所示,反演结果与原型的仿真计算结果一致,证明了本发明方法的有效性。
在另一个具体实施例中,仿真实验设计了涂覆材料参数不相同、材料反射系数变化趋势不相同的缩比模型B和目标原型。
仿真模型仍为SLICY模型,缩比模型B与原型的涂覆材料反射系数变化趋势不相同,缩比模型B和目标原型的参数详见表2。
表2缩比模型B和目标原型的仿真参数
Figure BDA0002994611560000122
缩比模型B和原型涂敷材料反射系数如图8所示。从图中可以看出,两种材料的反射系数变化趋势明显不同。由此可以看出反射系数在诸多影响因素的作用下具有复杂的变化趋势,这种复杂的变化趋势也是影响缩比反演的关键因素。由缩比模型B反演得到的原型结果如图9所示,图中给出了采用所提方法和待定物理缩比因子法反演的原型RCS结果,通过与仿真计算的原型RCS进行比对,证明了本发明方法的有效性。
此外,通过与待定物理缩比因子法的反演结果进行比对,可以证明基于面元分布的动态匹配反射系数方法的缩比反演精度比待定物理缩比因子法有明显提升。
在一个实施例中,如图10所示,提供了一种基于面元分布的动态匹配反射系数缩比测量装置,包括:缩比模型RCS测量值获取模块1002、第一反射系数信息获取模块1004、RCS面元分布信息获取模块1006和目标原型的RCS值获取模块1008,其中:
缩比模型RCS测量值获取模块1002,用于获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;涂覆缩比模型是对金属缩比模型的表面进行了材料涂覆,可为一个或多个,每个涂覆缩比模型的涂覆材料不同;
第一反射系数信息获取模块1004,用于根据雷达入射波与涂覆缩比模型面元法向形成的夹角获取涂覆材料的第一反射系数信息;
RCS面元分布信息获取模块1006,用于根据金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及第一反射系数信息,得到金属缩比模型在入射波角度下的RCS面元分布信息;
目标原型的RCS值获取模块1008,用于获取目标原型涂覆材料的第二反射系数信息,根据第二反射系数信息、面元分布信息和预设的几何缩比因子,得到目标原型的RCS值。
RCS面元分布信息获取模块1006还用于获取金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及第一反射系数信息;
将金属缩比模型RCS测量值、涂覆缩比模型RCS测量值和第一反射系数信息代入以下公式:
Figure BDA0002994611560000141
其中,
Figure BDA0002994611560000142
表示第1…m-1个涂覆缩比模型的RCS测量值,
Figure BDA0002994611560000143
表示雷达入射波的入射角;
Figure BDA0002994611560000144
表示金属缩比模型RCS测量值;m表示面元的类别数;θj,j=1…m表示入射波与第j类面元的法向方向形成的夹角,Rkj),k=1…m-1,j=1…m表示入射波与第k个涂覆缩比模型第j类面元的法向方向形成的夹角为θj时的反射系数;
Figure BDA0002994611560000145
表示面元法向为
Figure BDA0002994611560000146
的所有面元复RCS的叠加;
求解得到金属缩比模型在入射波角度下的RCS面元分布信息
Figure BDA0002994611560000147
目标原型的RCS值获取模块1008还用于获取目标原型涂覆材料的第二反射系数信息;
将第二反射系数信息、面元分布信息和预设的几何缩比因子,代入以下公式:
Figure BDA0002994611560000148
Figure BDA0002994611560000149
其中,
Figure BDA00029946115600001410
表示目标原型的RCS值,Rprotj),j=1…m表示入射波与目标原型第j类面元的法向方向形成的夹角为θj时的反射系数,s为预设的几何缩比因子;
根据公式计算结果得到目标原型的RCS值。
RCS面元分布信息获取模块1006还用于当涂覆缩比模型的面元分为两类时,获取金属缩比模型RCS测量值、涂覆缩比模型RCS测量值,以及涂覆缩比模型的涂覆材料的第一反射系数信息;
将金属缩比模型RCS测量值、涂覆缩比模型RCS测量值和第一反射系数信息代入以下公式:
Figure BDA0002994611560000151
求解得到金属缩比模型在入射波角度下的RCS面元分布信息
Figure BDA0002994611560000152
目标原型的RCS值获取模块1008还用于获取目标原型涂覆材料的第二反射系数信息;
将第二反射系数信息、面元分布信息和预设的几何缩比因子,代入以下公式:
Figure BDA0002994611560000153
Figure BDA0002994611560000154
根据公式计算结果得到目标原型的RCS值。
关于基于面元分布的动态匹配反射系数缩比测量装置的具体限定可以参见上文中对于基于面元分布的动态匹配反射系数缩比测量方法的限定,在此不再赘述。上述基于面元分布的动态匹配反射系数缩比测量装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是终端,其内部结构图可以如图11所示。该计算机设备包括通过***总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***和计算机程序。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于面元分布的动态匹配反射系数缩比测量方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图11中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,该存储器存储有计算机程序,该处理器执行计算机程序时实现上述方法实施例中的步骤。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现上述方法实施例中的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (8)

1.一种基于面元分布的动态匹配反射系数缩比测量方法,其特征在于,所述方法包括:
获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同;
根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;
将所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值和所述第一反射系数信息代入以下公式:
Figure FDA0003654869830000011
其中,
Figure FDA0003654869830000012
表示第1…m-1个涂覆缩比模型的RCS测量值,
Figure FDA0003654869830000013
表示雷达入射波的入射角;
Figure FDA0003654869830000014
表示所述金属缩比模型RCS测量值;m表示面元的类别数;θj,j=1…m表示入射波与第j类面元的法向方向形成的夹角,Rkj),k=1…m-1,j=1…m表示入射波与第k个涂覆缩比模型第j类面元的法向方向形成的夹角为θj时的反射系数;
Figure FDA0003654869830000015
表示面元法向为
Figure FDA0003654869830000016
的所有面元复RCS的叠加;
求解得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息
Figure FDA0003654869830000017
获取目标原型涂覆材料的第二反射系数信息,将所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,代入以下公式:
Figure FDA0003654869830000021
Figure FDA0003654869830000022
其中,
Figure FDA0003654869830000023
表示所述目标原型的RCS值,Rprotj),j=1…m表示入射波与目标原型第j类面元的法向方向形成的夹角为θj时的反射系数,s为预设的几何缩比因子,根据公式计算结果得到所述目标原型的RCS值。
2.根据权利要求1所述的方法,其特征在于,所述获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同,包括:
获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,所述涂覆缩比模型的个数为面元的类别数减1,每个所述涂覆缩比模型的涂覆材料不同;所述面元根据法向方向分类,法向方向相同的面元为一类。
3.根据权利要求2所述的方法,其特征在于,根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息,包括:
根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;所述第一反射系数信息随着雷达入射波与所述涂覆缩比模型的不同面元形成的夹角大小变化而变化。
4.根据权利要求3所述的方法,其特征在于,当所述涂覆缩比模型的面元分为两类时,根据所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述第一反射系数信息,得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息的步骤包括:
获取所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值,以及所述涂覆缩比模型的涂覆材料的第一反射系数信息;
将所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值和所述第一反射系数信息代入以下公式:
Figure FDA0003654869830000031
求解得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息
Figure FDA0003654869830000032
5.根据权利要求4所述的方法,其特征在于,当所述涂覆缩比模型的面元分为两类时,根据所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,得到所述目标原型的RCS值的步骤包括:
获取目标原型涂覆材料的第二反射系数信息;
将所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,代入以下公式:
Figure FDA0003654869830000033
Figure FDA0003654869830000034
根据公式计算结果得到所述目标原型的RCS值。
6.一种基于面元分布的动态匹配反射系数缩比测量装置,其特征在于,所述装置包括:
缩比模型RCS测量值获取模块,用于获取金属缩比模型的金属缩比模型RCS测量值和涂覆缩比模型的涂覆缩比模型RCS测量值;所述涂覆缩比模型是对所述金属缩比模型的表面进行了材料涂覆,可为一个或多个,每个所述涂覆缩比模型的涂覆材料不同;
第一反射系数信息获取模块,用于根据雷达入射波与所述涂覆缩比模型面元法向形成的夹角获取所述涂覆材料的第一反射系数信息;
RCS面元分布信息获取模块,用于将所述金属缩比模型RCS测量值、所述涂覆缩比模型RCS测量值和所述第一反射系数信息代入以下公式:
Figure FDA0003654869830000041
其中,
Figure FDA0003654869830000042
表示第1…m-1个涂覆缩比模型的RCS测量值,
Figure FDA0003654869830000043
表示雷达入射波的入射角;
Figure FDA0003654869830000044
表示所述金属缩比模型RCS测量值;m表示面元的类别数;θj,j=1…m表示入射波与第j类面元的法向方向形成的夹角,Rkj),k=1…m-1,j=1…m表示入射波与第k个涂覆缩比模型第j类面元的法向方向形成的夹角为θj时的反射系数;
Figure FDA0003654869830000045
表示面元法向为
Figure FDA0003654869830000046
的所有面元复RCS的叠加;
求解得到所述金属缩比模型在所述入射波角度下的RCS面元分布信息
Figure FDA0003654869830000047
目标原型的RCS值获取模块,用于获取目标原型涂覆材料的第二反射系数信息,将所述第二反射系数信息、所述面元分布信息和预设的几何缩比因子,代入以下公式:
Figure FDA0003654869830000051
Figure FDA0003654869830000052
其中,
Figure FDA0003654869830000053
表示所述目标原型的RCS值,Rprotj),j=1…m表示入射波与目标原型第j类面元的法向方向形成的夹角为θj时的反射系数,s为预设的几何缩比因子,根据公式计算结果得到所述目标原型的RCS值。
7.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法的步骤。
CN202110325787.9A 2021-03-26 2021-03-26 基于面元分布的动态匹配反射系数缩比测量方法和装置 Active CN113030900B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110325787.9A CN113030900B (zh) 2021-03-26 2021-03-26 基于面元分布的动态匹配反射系数缩比测量方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110325787.9A CN113030900B (zh) 2021-03-26 2021-03-26 基于面元分布的动态匹配反射系数缩比测量方法和装置

Publications (2)

Publication Number Publication Date
CN113030900A CN113030900A (zh) 2021-06-25
CN113030900B true CN113030900B (zh) 2022-08-09

Family

ID=76474367

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110325787.9A Active CN113030900B (zh) 2021-03-26 2021-03-26 基于面元分布的动态匹配反射系数缩比测量方法和装置

Country Status (1)

Country Link
CN (1) CN113030900B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115356703B (zh) * 2022-10-17 2023-02-03 中国人民解放军国防科技大学 基于面元分布的粗糙目标rcs缩比测量方法及装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1983608B1 (en) * 2007-04-20 2013-02-27 Saab AB Airborne vehicle integrated antenna
JP5058838B2 (ja) * 2008-02-01 2012-10-24 キヤノン株式会社 情報処理装置および方法
CN104573368B (zh) * 2015-01-13 2018-02-23 北京航空航天大学 基于面元投影的三角形截面射线管电磁射线追踪方法
CN105823756B (zh) * 2016-03-24 2018-11-13 西安电子科技大学 一种金属太赫兹至远红外复折射率的联合反演方法
CN109520383B (zh) * 2017-09-20 2020-07-31 南京理工大学 基于matlab的体目标回波模拟方法
CN109614637B (zh) * 2018-10-29 2023-05-23 上海无线电设备研究所 一种非金属结构体电磁缩比设计方法
CN109614652B (zh) * 2018-11-13 2023-05-23 上海无线电设备研究所 一种含抑制缝隙散射涂层的缩比目标构造方法

Also Published As

Publication number Publication date
CN113030900A (zh) 2021-06-25

Similar Documents

Publication Publication Date Title
Strobl et al. Approximate kernel-based conditional independence tests for fast non-parametric causal discovery
Chávez-Hurtado et al. Polynomial-based surrogate modeling of RF and microwave circuits in frequency domain exploiting the multinomial theorem
CN113030900B (zh) 基于面元分布的动态匹配反射系数缩比测量方法和装置
Jacobs et al. Computationally efficient multi-fidelity Bayesian support vector regression modeling of planar antenna input characteristics
Pérez et al. Analysis of antennas on board arbitrary structures modeled by NURBS surfaces
CN113821878A (zh) 一种改善高超声速气动热热流分布异常的计算方法和装置
Koziel et al. Design-oriented modeling of antenna structures by means of two-level kriging with explicit dimensionality reduction
Gu et al. A trimmed moving total least-squares method for curve and surface fitting
Hafver et al. Environmental contours as Voronoi cells
Koziel et al. Tolerance-aware optimization of microwave circuits by means of principal directions and domain-restricted metamodels
Chen et al. Are you all normal? It depends!
CN115356703B (zh) 基于面元分布的粗糙目标rcs缩比测量方法及装置
Koziel et al. Uniform sampling in constrained domains for low-cost surrogate modeling of antenna input characteristics
Liu et al. Analysis of electromagnetic scattering with higher-order moment method and NURBS model
Pietrenko-Dabrowska et al. Cost-efficient surrogate modeling of high-frequency structures using nested kriging with automated adjustment of model domain lateral dimensions
CN111948463A (zh) 一种植入式天线工作现场性能测试方法及***
CN116381467A (zh) 双端口近场探头非对称补偿方法、装置、设备和存储介质
Kaz’Min et al. Evaluating the accuracy of reconstruction of the electrical and geometric parameters of multilayer dielectric coatings by a multifrequency radio-wave method for slow surface electromagnetic waves
CN112634425B (zh) 模型的渲染方法、装置、存储介质及计算机设备
CN114942233A (zh) 近红外光谱特征波长选择方法、装置、设备和存储介质
CN111044983B (zh) 双波段气象雷达外部标定验证方法、装置和计算机设备
Hurisse On the use of Glimm‐like schemes for transport equations on multidimensional domain
Roohani Ghehsareh et al. Numerical investigation of electromagnetic scattering problems based on the compactly supported radial basis functions
Jarem Rigorous coupled wave analysis of bipolar cylindrical systems: Scattering from inhomogeneous dielectric material, eccentric, composite circular cylinders
Sari et al. Valuation of the American put option as a free boundary problem through a high-order difference scheme

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zeng Yang

Inventor after: Pang Shuang

Inventor after: Yang Qi

Inventor after: Deng Bin

Inventor after: Wang Hongqiang

Inventor before: Zeng Yang

Inventor before: Feng Shuang

Inventor before: Yang Qi

Inventor before: Deng Bin

Inventor before: Wang Hongqiang