CN113742937A - 一种特征结构散射特性计算方法 - Google Patents

一种特征结构散射特性计算方法 Download PDF

Info

Publication number
CN113742937A
CN113742937A CN202111075284.7A CN202111075284A CN113742937A CN 113742937 A CN113742937 A CN 113742937A CN 202111075284 A CN202111075284 A CN 202111075284A CN 113742937 A CN113742937 A CN 113742937A
Authority
CN
China
Prior art keywords
reflection coefficient
electric field
feature
reflection
calculating
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.)
Pending
Application number
CN202111075284.7A
Other languages
English (en)
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.)
Shanghai Radio Equipment Research Institute
Original Assignee
Shanghai Radio Equipment Research Institute
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 Shanghai Radio Equipment Research Institute filed Critical Shanghai Radio Equipment Research Institute
Priority to CN202111075284.7A priority Critical patent/CN113742937A/zh
Publication of CN113742937A publication Critical patent/CN113742937A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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

Abstract

一种特征结构散射特性计算方法，包含：步骤S1、建立特征结构的有限大阵列，根据计算规模选择低频算法；步骤S2、获取与入射方向关于特征结构法线方向对称的反射方向上距离特征结构表面固定倍数波长的位置上的反射电场；步骤S3、根据反射系数的定义，结合步骤S2中获取的反射电场计算得出特征结构的近似反射系数；步骤S4、利用步骤S3中得到的近似反射系数，采用基于特征反射系数的弹跳射线法仿真计算特征结构的散射特性。本发明可以高效计算并精确仿真难以求解的电大尺寸特征结构单元的反射系数问题。

Description

一种特征结构散射特性计算方法
技术领域
本发明涉及雷达目标特性仿真技术，尤其涉及一种特征结构散射特性计算方法。
背景技术
采用直接建模方式获取特征结构散射特性在时间复杂度和空间复杂度上均难以满足现实需求。传统解决此类问题的方法是针对周期结构的单元采用Floquet周期边界求解获取特征结构的反射特性，然后借助传统高频方法开展求解。然而，对于特征结构单元为电大尺寸的情形，上述两类方法都无法有效开展计算。
发明内容
本发明的目的在于提供一种特征结构散射特性计算方法，可以高效计算并精确仿真难以求解的电大尺寸特征结构单元的反射系数问题。
为了达到上述目的，本发明提供一种特征结构散射特性计算方法，包含以下步骤：
步骤S1、建立特征结构的有限大阵列，根据计算规模选择低频算法；
步骤S2、获取与入射方向关于特征结构法线方向对称的反射方向上距离特征结构表面固定倍数波长的位置上的反射电场；
步骤S3、根据反射系数的定义，结合步骤S2中获取的反射电场计算得出特征结构的近似反射系数；
步骤S4、利用步骤S3中得到的近似反射系数，采用基于特征反射系数的弹跳射线法仿真计算特征结构的散射特性。
所述的有限大阵列是指周期数为10个×10个或以上。
所述低频算法采用多层快速多极子MLFMM。
当有限大阵列的周期数大于等于10个且小于20个时，选择距离特征结构表面5-10倍波长位置上的反射电场；当有限大阵列的周期数为20个×20个或以上时，选择距离特征结构表面10-15倍波长位置上的反射电场。
所述步骤S3中，反射系数Γ定义表达为：
Figure BDA0003262052820000021
式中，Er为反射电场标量值，Ei为入射电场标量值，反射系数相位可以用反射电场相位确定，Γ为频率、极化和入射方向的函数，在计算时，形成扫频、扫角、极化中的一种或几种的反射系数列表。
所述步骤S4中，反射场强表达为：
Figure BDA0003262052820000022
式中，
Figure BDA0003262052820000023
分别为入射电场的垂直和水平分量，
Figure BDA0003262052820000024
分别为入射平面的两个基，Γ⊥、Γ||分别为垂直反射系数和水平反射系数，由步骤S3中所提供，入射电场和反射电场同时设为垂直或水平。
本发明借助传统积分方程方法和反射系数定义，提出了一种求解电大尺寸、有耗介质特征结构单元的反射系数的高效方法，进而可以实现特征结构散射特性的快速仿真。
附图说明
图1是本发明提供的一种特征结构散射特性计算方法的流程图。
图2是场强位置示意图。
具体实施方式
以下根据图1～图2，具体说明本发明的较佳实施例。
如图1所示，本发明提供一种特征结构散射特性计算方法，包含以下步骤：
步骤S1、建立特征结构的有限大阵列(10个周期×10个周期或以上)，根据计算规模选择合适的低频算法；
对于大规模计算问题，低频算法一般选择多层快速多极子(MLFMM)，该MLFMM求解电大尺寸的时间复杂度和空间复杂度均为O(NlogN)，效率极高；
步骤S2、如图2所示，获取与入射方向关于特征结构法线方向对称的反射方向上距离特征结构表面5-10倍波长的位置上的反射电场；
在获取入射方向关于特征结构法线方向对称的反射方向上的场强时，选择计算所有特征结构单元往往因为边界效应会引入误差，一般选择不包含边界的特征结构单元部分的散射贡献，进而更加接近理论真值，在距离选择上，一般选择5-10倍波长，这样可以减少特征结构单元的绕射贡献。当有限大阵列周期数为20个×20个或以上时，距离一般选择10-15倍波长。
步骤S3、根据反射系数的定义，结合步骤S2中获取的反射电场计算得出特征结构的近似反射系数；
由经典电磁场理论，反射系数Γ定义表达为：
Figure BDA0003262052820000031
式中，Er为反射电场标量值，Ei为入射电场标量值，反射系数相位可以用反射电场相位确定，Γ为频率、极化和入射方向的函数，在计算时，形成扫频、扫角、极化(视问题需要计算一种或组合形式)反射系数列表，以便调用；
步骤S4、采用基于特征反射系数的弹跳射线法仿真计算特征结构的散射特性；
弹跳射线法在场强追踪时，反射场强表达为：
Figure BDA0003262052820000032
式中，
Figure BDA0003262052820000033
分别为入射电场的垂直和水平分量，
Figure BDA0003262052820000034
分别为入射平面的两个基，Γ⊥、Γ||分别为垂直反射系数和水平反射系数，由步骤S3中所提供，入射电场和反射电场同时设为垂直或水平。
特征结构反射系数的计算一般需要计算无限大周期，这样会导致求解效率的大幅降低。本发明选择有限大阵列作为求解对象，计算场强分布，结合反射系数的定义，获取反射系数。
本发明具有以下优点：
1、采用有限大周期阵列作为计算对象，结合反射系数的定义，可以高效计算精确仿真难以求解的电大尺寸特征结构单元反射系数问题。
2、基于反射系数文件替代特征结构本体，可以极为方便地嵌入到传统弹跳射线法，实现场强的射线追踪。
3、当特征结构透波时，对照步骤S3，结合透射系数的定义，只需形成透射系数文件，即可方便拓展求解透射问题。
本发明特别适用于特征结构单元为电大尺寸、有耗介质的情形，这种情形采用Floquet周期边界只是理论上可以求解，在内存、时间上无法满足工程需求。
需要说明的是，在本发明的实施例中，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述实施例，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。

Claims (6)

1.一种特征结构散射特性计算方法，其特征在于，包含以下步骤：
步骤S1、建立特征结构的有限大阵列，根据计算规模选择低频算法；
步骤S2、获取与入射方向关于特征结构法线方向对称的反射方向上距离特征结构表面固定倍数波长的位置上的反射电场；
步骤S3、根据反射系数的定义，结合步骤S2中获取的反射电场计算得出特征结构的近似反射系数；
步骤S4、利用步骤S3中得到的近似反射系数，采用基于特征反射系数的弹跳射线法仿真计算特征结构的散射特性。
2.如权利要求1所述的特征结构散射特性计算方法，其特征在于，所述的有限大阵列是指周期数为10个×10个或以上。
3.如权利要求1所述的特征结构散射特性计算方法，其特征在于，所述低频算法采用多层快速多极子MLFMM。
4.如权利要求2所述的特征结构散射特性计算方法，其特征在于，当有限大阵列的周期数大于等于10个且小于20个时，选择距离特征结构表面5-10倍波长位置上的反射电场；当有限大阵列的周期数为20个×20个或以上时，选择距离特征结构表面10-15倍波长位置上的反射电场。
5.如权利要求4所述的特征结构散射特性计算方法，其特征在于，所述步骤S3中，反射系数Γ定义表达为：
Figure FDA0003262052810000011
式中，Er为反射电场标量值，Ei为入射电场标量值，反射系数相位可以用反射电场相位确定，Γ为频率、极化和入射方向的函数，在计算时，形成扫频、扫角、极化中的一种或几种的反射系数列表。
6.如权利要求5所述的特征结构散射特性计算方法，其特征在于，所述步骤S4中，反射场强表达为：
Figure FDA0003262052810000012
式中，
Figure FDA0003262052810000015
分别为入射电场的垂直和水平分量，
Figure FDA0003262052810000014
分别为入射平面的两个基，Γ⊥、Γ||分别为垂直反射系数和水平反射系数，由步骤S3中所提供，入射电场和反射电场同时设为垂直或水平。
CN202111075284.7A 2021-09-14 2021-09-14 一种特征结构散射特性计算方法 Pending CN113742937A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111075284.7A CN113742937A (zh) 2021-09-14 2021-09-14 一种特征结构散射特性计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111075284.7A CN113742937A (zh) 2021-09-14 2021-09-14 一种特征结构散射特性计算方法

Publications (1)

Publication Number Publication Date
CN113742937A true CN113742937A (zh) 2021-12-03

Family

ID=78738842

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111075284.7A Pending CN113742937A (zh) 2021-09-14 2021-09-14 一种特征结构散射特性计算方法

Country Status (1)

Country Link
CN (1) CN113742937A (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003078326A (ja) * 2001-09-04 2003-03-14 Katsuo Tanaka 電磁波透過構造体の評価方法
US20110137625A1 (en) * 2009-12-08 2011-06-09 Asml Netherlands B.V. Methods and Apparatus for Determining Electromagnetic Scattering Properties and Structural Parameters of Periodic Structures
CN106291496A (zh) * 2016-07-26 2017-01-04 上海无线电设备研究所 基于四路径模型的运动目标与地面复合散射建模方法
CN106772301A (zh) * 2016-12-02 2017-05-31 上海无线电设备研究所 一种基于介质射线追踪的多层非平行界面介质电磁散射仿真算法
CN108763153A (zh) * 2018-04-11 2018-11-06 杭州电子科技大学 一种地下各向同性介质球电磁散射的计算方法
CN109581340A (zh) * 2018-10-11 2019-04-05 上海无线电设备研究所 一种基于时域弹跳射线法的等离子体电磁散射建模方法
US20200025911A1 (en) * 2018-05-24 2020-01-23 New York University System, method and computer-accessible medium for real time imaging using a portable device

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003078326A (ja) * 2001-09-04 2003-03-14 Katsuo Tanaka 電磁波透過構造体の評価方法
US20110137625A1 (en) * 2009-12-08 2011-06-09 Asml Netherlands B.V. Methods and Apparatus for Determining Electromagnetic Scattering Properties and Structural Parameters of Periodic Structures
CN106291496A (zh) * 2016-07-26 2017-01-04 上海无线电设备研究所 基于四路径模型的运动目标与地面复合散射建模方法
CN106772301A (zh) * 2016-12-02 2017-05-31 上海无线电设备研究所 一种基于介质射线追踪的多层非平行界面介质电磁散射仿真算法
CN108763153A (zh) * 2018-04-11 2018-11-06 杭州电子科技大学 一种地下各向同性介质球电磁散射的计算方法
US20200025911A1 (en) * 2018-05-24 2020-01-23 New York University System, method and computer-accessible medium for real time imaging using a portable device
CN109581340A (zh) * 2018-10-11 2019-04-05 上海无线电设备研究所 一种基于时域弹跳射线法的等离子体电磁散射建模方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
XIAOBING WANG等: "Efficient GPU implementation of SBR for fast computation of composite scattering from electrically large target over a randomly rough surface", 《IEEE》, pages 1666 - 1667 *
丁建军;陈磊;刘志伟;陈如山;: "基于时域弹跳射线法分析电大尺寸目标的散射", 系统工程与电子技术, no. 09, pages 60 - 63 *
刘健等: "飞行器太赫兹散射特性仿真及抑制方法研究", 《制导与引信》, vol. 42, no. 1, pages 42 - 48 *
方志成;: "多层介质板散射特性的分析", 武汉交通职业学院学报, no. 04, pages 89 - 91 *
高火涛等: "部件分解法/高频混合法在复杂目标散射近场预估中的应用", 《微波学报》, no. 2, pages 54 - 59 *

Similar Documents

Publication Publication Date Title
US9735899B2 (en) Device and method for calibrating antenna array systems
US9335359B2 (en) Far electromagnetic field estimation method and apparatus, and near electromagnetic field measurement apparatus
US20110193566A1 (en) Multichannel absorberless near field measurement system
US10928432B2 (en) Reverberation chamber
CN108562897A (zh) 一种mimo穿墙雷达的结构稀疏成像方法和装置
US10394204B1 (en) Methods and systems for synthesis of a waveguide array antenna
WO2009046516A1 (en) Multichannel absorberless near field measurement system
CN107703478B (zh) 基于互相关矩阵的扩展孔径二维doa估计方法
CN102708229A (zh) 复杂分层媒质结构的矩阵分解结合新奇异值分解方法
Zhang et al. The low-frequency sound power measuring technique for an underwater source in a non-anechoic tank
US10725146B2 (en) Wideband radio-frequency antenna
CN113742937A (zh) 一种特征结构散射特性计算方法
CN109376422A (zh) 一种均匀圆阵优化设计评估方法和装置
CN103760416B (zh) 一种定向耦合器测量高功率微波相位的误差分析方法
Park et al. An improved technique for the evaluation of transverse slot discontinuities in rectangular waveguide
Chen et al. Pyramid-shaped grid for elastic wave propagation
CN115169174A (zh) 一种同构蜂群目标电磁散射特性快速计算方法
TWI783713B (zh) 偶極矩模型建立方法與系統
CN102749538A (zh) 大功率发射源主波束照射下抛物面天线耦合电平数值仿真的校模方法
CN115542315B (zh) 基于adi-fdtd的探地雷达正演模拟方法及系统
CN108268747B (zh) 基于力密度法的低栅瓣索网天线机电综合设计方法
CN117332533A (zh) 一种计算导体超宽带散射的网格稳定、高阶稳定的阶数步进时域积分方程方法
CN115081204A (zh) 超电大平台超宽谱电磁脉冲时域瞬态场预测方法
Bekasiewicz et al. Low-cost EM-simulation-driven multi-objective optimization of antennas
Okada et al. RCS Measurement and PO Simulation of a Scale Model Rocket

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