CN111859704B - 一种分布式多视角下非刚体目标电磁散射建模方法 - Google Patents
一种分布式多视角下非刚体目标电磁散射建模方法 Download PDFInfo
- Publication number
- CN111859704B CN111859704B CN202010751032.0A CN202010751032A CN111859704B CN 111859704 B CN111859704 B CN 111859704B CN 202010751032 A CN202010751032 A CN 202010751032A CN 111859704 B CN111859704 B CN 111859704B
- Authority
- CN
- China
- Prior art keywords
- electromagnetic
- target
- rigid body
- scattering
- vibration
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000006870 function Effects 0.000 claims abstract description 45
- 238000005070 sampling Methods 0.000 claims abstract description 30
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000005284 excitation Effects 0.000 claims abstract description 10
- 230000005405 multipole Effects 0.000 claims abstract description 10
- 230000001678 irradiating effect Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000006698 induction Effects 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 238000004088 simulation Methods 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 108010015046 cell aggregation factors Proteins 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000005094 computer simulation Methods 0.000 abstract description 17
- 238000004458 analytical method Methods 0.000 abstract description 12
- 238000005516 engineering process Methods 0.000 abstract description 10
- 239000010410 layer Substances 0.000 description 15
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000002592 echocardiography Methods 0.000 description 5
- 239000002184 metal Substances 0.000 description 4
- 229910052751 metal Inorganic materials 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种分布式多视角下非刚体目标电磁散射建模方法,包括以下步骤:步骤1:电磁波照射非刚体目标,基于RWG基函数,获得具有目标振动特性的感应电磁流;步骤2:构建出具有目标振动状态的电磁场描述模型;步骤3:根据电磁场描述模型进行快速多极计算,获得不同采样时刻下的高精度目标电磁散射回波数据;步骤4:根据高精度目标电磁散射回波数据,计算不同激励条件下的具有目标振动状态的电磁场描述模型,从而获得非刚体目标电磁散射特征。此发明解决了常规电磁散射特性仿真建模技术对非刚体目标电磁特性分析逼真度较差的问题,将不同基函数的振动特性映射到电磁场矩量法的空间分组中心,开展迭代计算,提升了仿真建模的适应性和精确度。
Description
技术领域
本发明涉及目标电磁散射特性仿真建模技术领域,具体涉及一种分布式多视角下非刚体目标电磁散射建模方法。
背景技术
大尺寸(飞机)目标在运动过程中固有的振动特性引起各个部件的局部振动,极大的影响目标电磁散射特征。常规电磁散射特性仿真建模技术将目标当做刚体开展电磁求解,目标为精确或全局运动一致,对真实作战环境下的目标电磁特性分析逼真度较差。为了弥补传统电磁散射特性仿真建模的不足,充分发挥电磁散射技术在武器系统设计及评估中的优势,急需深入研究分布式多视角探测下非刚体目标电磁散射特性仿真建模技术。
随着目标仿真建模高精度、高逼真度的要求,常规电磁建模手段逐渐无法满足工程需求。对于大尺寸目标,通过模态分析,可得到机械结构在某一受影响的频率范围内各阶模态的振动特性,以及机械结构在此频段内及在内部或外部各种振源激励下的振动响应结果,用于目标结构设计,然后联合目标电磁散射特性分析即可逼真目标真实作战/飞行中的真实雷达散射回波。因此,急需开展联合目标电磁特性和结构特性的联合仿真建模技术,为分布式下目标电磁特性分析提供有力手段。
2015年周城宏、钱卫平、郭永强等发表在四川兵工学报上的非专利文献“动态目标的电磁散射模型”文章中对动态目标则在准静态原理假设下近似计算,基于Maxwell方程组和入射电磁场边界条件在雷达参考系和目标参考系相互转换规律,在目标参考系中构建了动态电磁散射问题的精确数学形式,非均匀边界条件下的Maxwell方程求解问题,并讨论了准静态近似的理论依据,指出其适用范围局限于非相对论速度的动态目标。该文献为准静态目标电磁特性研究提供了必要的理论依据,但是未考虑目标动态特新引起的电磁特征,对动态问题的研究和论述略显不足。
2009年徐林、袁仕继、李康乐等发表在现代电子技术上的非专利文献“非刚体运动建模与宽带雷达回波仿真”中以空间非刚体运动目标产生的微多普勒特征为空间目标物理特性反演的依据,建立了目标微动数学模型,并基于典型的宽带雷达波形(LFM)推导了微动对回波的调制形式,分析了微多普勒的特点,最后利用仿真实验验证了该模型的有效性。该文献重点考察了非刚体运动目标产生的雷达效应,尚未设计非刚体目标电磁特性。
2013年陈如山、樊振宏等发表的专利文献“薄介质涂敷的金属旋转对称目标电磁散射快速计算方法”,专利公开号CN103177193A,公开了一种薄介质涂敷的金属旋转对称目标电磁散射快速计算方法。针对这种金属介质混合结构,仅需对金属部分建立电场积分方程,介质部分不需要建立方程进行描述。利用目标的旋转对称特性,将3维问题降成2维问题分析,加速求解速度,降低求解内存。该文献针对动态刚体目标为研究对象,为设计非刚体目标的动态过程。
2011年曾涛、胡程等发表的专利文献“一种前向散射雷达地面运动目标信号建模及成像方法,专利公开号为CN102129067A,发明所提出的方法考虑地面反射、大衍射角以及目标运动方向较为随机等实际情况,建立了精确的地面运动目标前向散射信号模型,修正了SISAR成像算法,并对目标像进行了更为精确的估计。修正后的成像算法在大衍射角、目标斜穿基线和多径干涉情况下依然可以获得精确的目标像,这些目标像与目标真实轮廓形状非常接近,可有效地对不同形状的目标进行分类识别,该文献研究对象为刚体目标,为设计非刚体目标引起的目标电磁特性变化。
发明内容
本发明的目的是提供一种分布式多视角下非刚体目标电磁散射建模方法。此方法旨在解决常规电磁散射特性仿真建模技术将目标当做刚体开展电磁求解,对非刚体目标电磁特性分析逼真度较差的问题,将不同基函数的振动特性映射到电磁场矩量法的空间分组中心,开展迭代计算,获取非刚体目标在不同探测角度下的雷达散射回波,提升仿真建模的适应性和精确度,为非刚体目标电磁散射特性分析提供有效技术手段。
为达到上述目的,本发明提供了一种分布式多视角下非刚体目标电磁散射建模方法,包括以下步骤:
步骤1:电磁波照射局部振动的非刚体目标,并基于RWG基函数,获得具有目标振动特性的感应电磁流;
步骤2:根据具有目标振动特性的感应电磁流,构建出具有目标振动状态的电磁场描述模型;
步骤3:基于MLFMA算法,根据具有目标振动状态的电磁场描述模型进行快速多极计算,获得不同采样时刻下的高精度目标电磁散射回波数据;
步骤4:根据不同采样时刻下的高精度目标电磁散射回波数据,计算出不同激励条件下的具有目标振动状态的电磁场描述模型,从而获得非刚体目标局部振动下的电磁散射特征。
最优选的,获得具有目标振动特性的感应电磁流包括以下步骤:
步骤1.1:外来电磁波照射局部振动的非刚体目标,在非刚体目标表面形成感应电磁流J,且感应电磁流J在外界空间产生电磁散射场;
步骤1.2:基于RWG基函数,对感应电磁流J展开,获得具有目标振动特性的感应电磁流。
最优选的,电磁散射场为Es,Hs,且分别满足:
其中,为虚数单位;s为散射场标识;μ为空间背景的磁导率;/>为入射电磁波角频率;τ'为积分标识;k为波数,/>为梯度算子,/>为散度算子;G为自由空间的格林函数,且满足:
G=e-ikR/(4πR)。
最优选的,具有目标振动特性的感应电磁流为J(r),且满足:
其中,r为相邻平面三角形贴片内部的位置矢量;N为基函数的个数;±为RWG基函数对应的两个相邻平面三角形贴片;Ai(t)为相邻平面三角形贴片中第i个顶点的振动幅度信息;s(t)为相邻平面三角形贴片中第i个顶点的相位移动信息;t为时间;Ln为对应基函数的线段长度;an为感应电磁流展开系数;Tn为对应基函数相关的三角形编号;ρi为相邻平面三角形贴片中第i个顶点的位置信息,且满足:
其中,为相邻平面三角形贴片中第i个顶点初始时刻的位置信息;nn为目标的剖分顶点个数。
最优选的,具有目标振动状态的电磁场描述模型为Esct(r(t)),且满足:
其中,r(t)和r'(t)分别对应不同采样时刻下的相邻平面三角形贴片内部的位置矢量。
最优选的,快速多极计算包括以下步骤:
步骤3.1:对电磁场描述模型进行匹配,生成电磁波矩阵方程;
步骤3.2:根据电磁波矩阵方程,计算出不同采样时刻下的目标表面感应电磁流系数;
步骤3.3:基于MLFMA算法,在相邻振动模态采样间隔上,利用上一次仿真得到的目标表面感应电磁流系数作为本次初值开展迭代计算,从而获得不同采样时刻下的高精度目标电磁散射回波数据。
最优选的,电磁波矩阵方程满足:
ZN×N(t0)IN×1(t0)=bN×1(t0)
其中,其中,ZN×N(t0)为t0时刻的阻抗矩阵,IN×1(t0)为t0时刻的感应电磁流系数,bn(t0)为t0时刻的激励矩阵,每个具体元素满足:
I(t)为动态感应系数矩阵,m,n表示ZN×N(t)中对应元素的行号和列号,且满足:
I(t)=[a1(t),a2(t),.....,ap(t)]T
P表示总的采样时刻,满足:
最优选的,迭代计算满足:
其中,Bm表示由m的附近组构成的集合;Gn表示由第n个组中心的基函数构成的集合;为场组内基函数的配置因子,且满足:
为源组内基函数的聚合因子,且满足:
为远组间的耦合矩阵,且满足:
最优选的,目标振动特性包络目标RCS和多普勒回波特征信号。
运用此发明,解决了常规电磁散射特性仿真建模技术将目标作为刚体开展电磁求解,对非刚体目标电磁特性分析逼真度较差的问题,将不同基函数的振动特性映射到电磁场矩量法的空间分组中心,开展迭代计算,获取非刚体目标在不同探测角度下的雷达散射回波,提升了仿真建模的适应性和精确度,为非刚体目标电磁散射特性分析提供了有效技术手段。
相对于现有技术,本发明具有以下有益效果:
本发明提供的分布式多视角下非刚体目标电磁散射建模方法解决了常规电磁散射特性仿真建模技术将目标作为刚体开展电磁求解,对非刚体目标电磁特性分析逼真度较差的问题,将不同基函数的振动特性映射到电磁场矩量法的空间分组中心,开展迭代计算,获取非刚体目标在不同探测角度下的雷达散射回波,提升了仿真建模的适应性和精确度,为非刚体目标电磁散射特性分析提供了有效技术手段。
附图说明
图1为本发明提供的非刚体目标电磁散射建模方法流程图;
图2(a)和图2(b)分别为某机翼模型作为刚体目标和非刚体目标时,在某一振动时刻下的状态示意图;
图3为本发明提供的RWG基函数定义的网格示意图;
图4(a)为本发明提供的MLFMA算法的单层电磁场结构示意图;
图4(b)为本发明提供的MLFMA算法的多层分组结构示意图。
具体实施方式
以下结合附图通过具体实施例对本发明作进一步的描述,这些实施例仅用于说明本发明,并不是对本发明保护范围的限制。
本发明提供了一种分布式多视角下非刚体目标电磁散射建模方法,如图1所示,包括以下步骤:
步骤1:外来电磁波Ei,Hi(i表示激励标识)照射局部振动的非刚体目标,并基于RWG基函数,获得具有目标振动特性的感应电磁流J(r);其中,获得具有目标振动特性的感应电磁流J(r)还包括以下步骤:
形成感应电磁流J和电磁散射场Es,Hs;形成感应电磁流和电磁散射场包括以下步骤:
步骤1.1:外来电磁波Ei,Hi照射局部振动的非刚体目标,在非刚体目标表面形成感应电磁流J;如图2(a)为某机翼模型作为刚体目标,在某一振动时刻下的状态示意图,刚体目标所在平面平行于直角坐标系中x-y的水平面;如图2(b)所示,为某机翼模型作为非刚体目标,在某一振动时刻下的状态示意图,非刚体目标所在平面与直角坐标系中x-y的水平面成一定角度,且其上各点在直角坐标系中z轴上的投影各不相同;且感应电磁流J在外界空间产生电磁散射场Es,Hs,且分别满足:
其中,为虚数单位;s为散射场标识;μ为空间背景的磁导率;/>为入射电磁波角频率;τ'为积分标识;;k为波数,/>为梯度算子,/>为散度算子;G为自由空间的格林函数,且满足:
G=e-ikR/(4πR)。
步骤1.2:采用矩量法,将感应电磁流J添加到矩量法中描述局部感应电流的基函数上,即基于广义的屋脊基函数(RWG基函数),对感应电磁流J展开,获得具有目标振动特性的感应电磁流J(r),且满足:
其中,r为相邻平面三角形贴片内部的位置矢量;N为基函数的个数;±为RWG基函数对应的两个相邻平面三角形贴片;Ai(t)为相邻平面三角形贴片中第i个顶点的振动幅度信息;s(t)为相邻平面三角形贴片中第i个顶点的相位移动信息;t为时间;Ln为对应基函数的线段长度;an为感应电磁流展开系数;Tn为对应基函数相关的三角形编号;;ρi为相邻平面三角形贴片中第i个顶点的位置信息,且满足:
其中,为相邻平面三角形贴片中第i个顶点初始时刻的位置信息;nn为目标的剖分顶点个数。
如图3所示,如果RWG基函数对应的三角形在部件的边界上,则没有对应的三角形对形成互耦,此时RWG仅定义在一个贴片上。在刚体目标电磁仿真建模中,定义基函数的顶点位置ρ不变或全局周期性运动,各个顶点间的相对位移保持不变。然而,在非刚体目标振动中,各个顶点ρ的振动幅度不同,必须添加周期性振动幅度才能完成非刚体目标在振动状态下的电磁建模。
步骤2:根据具有目标振动特性的感应电磁流J(r)和电磁散射场Es,Hs,构建出具有目标振动状态的电磁场描述模型Esct(r(t)),且满足:
其中,r(t)和r'(t)分别对应不同采样时刻下的相邻平面三角形贴片内部的位置矢量。
其中,目标振动状态下的振动特性包络目标RCS、多普勒回波等特征信号。
步骤3:基于多层快速多极(MLFMA)算法,根据具有目标振动状态的电磁场描述模型Esct(r(t))进行快速多极计算,获得不同采样时刻下的高精度目标电磁散射回波数据;其中,快速多极计算包括以下步骤:
步骤3.1:采用伽略金方法对电磁场描述模型Esct(r(t))进行匹配,生成电磁波矩阵方程;电磁波矩阵方程满足:
ZN×N(t0)IN×1(t0)=bN×1(t0)
其中,其中,ZN×N(t0)为t0时刻的阻抗矩阵,IN×1(t0)为t0时刻的感应电磁流系数,bn(t0)为t0时刻的激励矩阵,每个具体元素满足:
I(t)为动态感应系数矩阵,m,n表示ZN×N(t)中对应元素的行号和列号,且满足:
I(t)=[a1(t),a2(t),.....,aP(t)]T
P表示总的采样时刻,满足:
其中,电磁波矩阵方程为时变方程,具有目标振动特性矩量法的阻抗矩阵Zn×n(t)是变化的,则不同目标振动状态下得到的电磁波矩阵方程也不同。
步骤3.2:根据电磁场矩阵方程,计算出不同采样时刻下的目标表面感应电磁流系数ap(t)(p=1,2,......P表示采样时刻)。
步骤3.3:基于多层快速多极(MLFMA)算法,对电磁波矩阵方程进行分组转化,在相邻振动模态时间采样间隔上,利用前一时刻仿真得到的目标表面感应电磁流系数作为本次初值开展迭代计算,从而获得不同采样时刻下的高精度目标电磁散射回波数据;其中图4(a)表示最底层的盒子相互作用,图4(b)表示多层的层级结构,迭代计算满足:
其中,Bm表示由m的附近组构成的集合;Gn表示由第n个组中心的基函数构成的集合;为场组内基函数的配置因子,且满足:
为源组内基函数的聚合因子,且满足:
为远组间的耦合矩阵,且满足:
其中,迭代计算是对目标振动过程进行离散时间采样,不同振动模态下的矩量法矩阵方程的初始解为前一个时刻的解,极大的减小当前采样时刻的迭代计算次数。
非刚体目标振动为不同部件的振动幅度不一致,但振动周期相同,通过对目标进行离散采样(表面或体)后,用各个节点的振动进行描述,并建立与电磁场基函数的关联关系。
不同于刚体目标运动,非刚体目标振动的阻抗矩阵Zn×n(t)在不同采样时刻下是变化的;不同采样时刻下必须重新计算才能确保仿真精度。
在本实施例中,采用MLFMA算法实现非刚体目标的高效快速计算,为保证仿真精度,一个振动周期内选择30个采样时间点;对于不同的采样时刻(除0时刻),每次迭代初值不必从0开始,开始利用前一时刻的采样值。
MLFMA算法是首先将目标区域按照八叉树方式划分为远亲组和近亲组,不同组之间的电磁散射作用耦合到盒子组中心,利用每个盒子中心的电磁波代表整个盒子内电磁波的贡献,不同目标部件上的感应电磁波对其他部件的耦合作用通过相应盒子组中心来转换代替,而不是每个基函数之间的相互作用。MLFMA算法的具体步骤包括:
最细层的多基展开:将最细层盒子组包含的电磁波集中到盒子中心上;
多级聚合:将电磁波源组盒子中心的电磁波集中到父层组盒子的中心位置处,至到第二层结束;
多级转移:从第二层开始,将源组盒子中心的电磁波贡献传递到远亲组盒子的组中心上,完成最高层盒子的耦合;
多级配置:从第二层盒子开始,向子层盒子组中心进行配置,至到最细曾结束,每次配置中既有父层组盒子中心的贡献,也有本层组远亲组盒子的贡献;
部分场展开:将最细层盒子组中心的电磁波贡献传播到盒子内的基函数上,完成电磁波远区耦合;
直接计算的贡献:将上述电磁波贡献与直接计算部分得到场累加完成整个区域的迭代计算。
步骤4:根据不同采样时刻下的高精度目标电磁散射回波数据,计算出不同激励条件下的具有目标振动状态的电磁场描述模型Esct(r(t)),从而获得非刚体目标局部振动下的电磁散射特征。
目标散射回波数据的周期特性与目标自身振动周期是振动幅度等信息是线性映射关系,通过雷达散射回波信号可直接反演非刚体目标振动状态,进而为目标飞行状态(高度、大气密度、速度等)参数确定提供了理论依据。
本发明的工作原理:
外来电磁波照射局部振动的非刚体目标,形成感应电磁流和电磁散射场;基于RWG基函数,对所述感应电磁流展开,获得具有目标振动特性的感应电磁流;根据所述具有目标振动特性的感应电磁流和所述电磁散射场,构建出具有目标振动状态的电磁场描述模型;基于MLFMA算法,根据所述具有目标振动状态的电磁场描述模型进行快速多极计算,获得不同采样时刻下的高精度目标电磁散射回波数据;根据所述不同采样时刻下的高精度目标电磁散射回波数据,计算出不同激励条件下的所述具有目标振动状态的电磁场描述模型,从而获得非刚体目标局部振动下的电磁散射特征。
综上所述,本发明一种分布式多视角下非刚体目标电磁散射建模方法,解决了常规电磁散射特性仿真建模技术将目标当做刚体开展电磁求解,对非刚体目标电磁特性分析逼真度较差的问题,将不同基函数的振动特性映射到电磁场矩量法的空间分组中心,开展迭代计算,获取非刚体目标在不同探测角度下的雷达散射回波,提升了仿真建模的适应性和精确度,为非刚体目标电磁散射特性分析提供了有效技术手段。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (9)
1.一种分布式多视角下非刚体目标电磁散射建模方法,其特征在于,包括以下步骤:
步骤1:电磁波照射局部振动的非刚体目标,并基于RWG基函数,获得具有目标振动特性的感应电磁流;
步骤2:根据所述具有目标振动特性的感应电磁流,构建出具有目标振动状态的电磁场描述模型;
步骤3:基于MLFMA算法,根据所述具有目标振动状态的电磁场描述模型进行快速多极计算,获得不同采样时刻下的高精度目标电磁散射回波数据;
步骤4:根据所述不同采样时刻下的高精度目标电磁散射回波数据,计算出不同激励条件下的所述具有目标振动状态的电磁场描述模型,从而获得非刚体目标局部振动下的电磁散射特征。
2.如权利要求1所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,获得所述具有目标振动特性的感应电磁流包括以下步骤:
步骤1.1:外来电磁波照射局部振动的非刚体目标,在非刚体目标表面形成感应电磁流J,且所述感应电磁流J在外界空间产生电磁散射场;
步骤1.2:基于RWG基函数,对所述感应电磁流J展开,获得具有目标振动特性的感应电磁流。
3.如权利要求2所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述电磁散射场为Es,Hs,且分别满足:
其中,为虚数单位;s为散射场标识;μ为空间背景的磁导率;/>
为入射电磁波角频率;τ'为积分标识,k为波数,为梯度算子,/>为散度算子;G为自由空间的格林函数,且满足:
G=e-ikR/(4πR)。
4.如权利要求3所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述具有目标振动特性的感应电磁流为J(r),且满足:
其中,r为相邻平面三角形贴片内部的位置矢量;N为基函数的个数;±为RWG基函数对应的两个相邻平面三角形贴片;Ai(t)为相邻平面三角形贴片中第i个顶点的振动幅度信息;s(t)为相邻平面三角形贴片中第i个顶点的相位移动信息;t为时间;Ln为对应基函数的线段长度;an为感应电磁流展开系数;Tn为对应基函数相关的三角形编号;ρi为相邻平面三角形贴片中第i个顶点的位置信息,且满足:
其中,为相邻平面三角形贴片中第i个顶点初始时刻的位置信息;nn为目标的剖分顶点个数。
5.如权利要求4所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述具有目标振动状态的电磁场描述模型为Esct(r(t)),且满足:
其中,r(t)和r'(t)分别对应不同采样时刻下的相邻平面三角形贴片内部的位置矢量。
6.如权利要求5所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述快速多极计算包括以下步骤:
步骤3.1:对电磁场描述模型进行匹配,生成电磁波矩阵方程;
步骤3.2:根据电磁波矩阵方程,计算出不同采样时刻下的目标表面感应电磁流系数;
步骤3.3:基于MLFMA算法,在相邻振动模态采样间隔上,利用上一次仿真得到的所述目标表面感应电磁流系数作为本次初值开展迭代计算,从而获得不同采样时刻下的高精度目标电磁散射回波数据。
7.如权利要求6所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述电磁波矩阵方程满足:
ZN×N(t0)IN×1(t0)=bN×1(t0)
其中,ZN×N(t0)为t0时刻的阻抗矩阵,IN×1(t0)为t0时刻的感应电磁流系数,bn(t0)为t0时刻的激励矩阵,每个具体元素满足:
I(t)为动态感应系数矩阵,m,n表示ZN×N(t)中对应元素的行号和列号,且满足:
I(t)=[a1(t),a2(t),.....,ap(t)]T
P表示总的采样时刻,满足:
8.如权利要求6所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述迭代计算满足:
其中,Bm表示由m的附近组构成的集合;Gn表示由第n个组中心的基函数构成的集合;为场组内基函数的配置因子,且满足:
为源组内基函数的聚合因子,且满足:
为远组间的耦合矩阵,且满足:
9.如权利要求1所述的分布式多视角下非刚体目标电磁散射建模方法,其特征在于,所述目标振动特性包络目标RCS和多普勒回波特征信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010751032.0A CN111859704B (zh) | 2020-07-30 | 2020-07-30 | 一种分布式多视角下非刚体目标电磁散射建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010751032.0A CN111859704B (zh) | 2020-07-30 | 2020-07-30 | 一种分布式多视角下非刚体目标电磁散射建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111859704A CN111859704A (zh) | 2020-10-30 |
CN111859704B true CN111859704B (zh) | 2024-06-18 |
Family
ID=72945153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010751032.0A Active CN111859704B (zh) | 2020-07-30 | 2020-07-30 | 一种分布式多视角下非刚体目标电磁散射建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111859704B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112836375B (zh) * | 2021-02-07 | 2023-03-03 | 中国人民解放军空军工程大学 | 一种高效目标电磁散射仿真方法 |
CN114330124A (zh) * | 2021-12-29 | 2022-04-12 | 东南大学 | 一种周期目标电磁散射特性快速仿真方法 |
CN114741880A (zh) * | 2022-04-15 | 2022-07-12 | 北京理工大学 | 一种计算群目标电磁特性的方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105303022A (zh) * | 2014-07-28 | 2016-02-03 | 南京辉宏电子信息有限公司 | 快速获取目标电磁散射特性的高斯波束方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7535228B2 (en) * | 2006-03-21 | 2009-05-19 | Radiation Monitoring Devices, Inc. | Sensor array for nuclear magnetic resonance imaging systems and method |
GB201116518D0 (en) * | 2011-09-23 | 2011-11-09 | Isis Innovation | Investigation of physical properties of an object |
CN106529082A (zh) * | 2016-12-02 | 2017-03-22 | 上海无线电设备研究所 | 一种快速计算电大尺寸目标电磁散射特征的方法 |
CN107271981B (zh) * | 2017-07-31 | 2020-02-21 | 河海大学 | 一种目标电磁散射数据模型的建立方法 |
CN107992676A (zh) * | 2017-11-29 | 2018-05-04 | 上海无线电设备研究所 | 一种运动目标时域散射回波的快速仿真建模方法 |
CN109190149B (zh) * | 2018-07-20 | 2023-04-21 | 北京理工大学 | 一种基于鸟类电磁散射模型提取振翅频率的仿真验证方法 |
CN109472059B (zh) * | 2018-10-17 | 2023-06-06 | 西安电子科技大学 | 基于测量应变的相控阵天线幅相补偿方法 |
CN109376485A (zh) * | 2018-12-03 | 2019-02-22 | 上海无线电设备研究所 | 基于aca-mlfma加速的区域分解非共形网格的快速仿真建模方法 |
CN110502782B (zh) * | 2019-07-04 | 2022-04-26 | 西安电子科技大学 | 基于分区矢量输运的大体量箔条云电磁散射测定方法 |
-
2020
- 2020-07-30 CN CN202010751032.0A patent/CN111859704B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105303022A (zh) * | 2014-07-28 | 2016-02-03 | 南京辉宏电子信息有限公司 | 快速获取目标电磁散射特性的高斯波束方法 |
Non-Patent Citations (1)
Title |
---|
基于有限元方法的粗糙面与目标复合电磁散射快速混合算法研究;何红杰;《中国博士学位论文电子期刊网 基础科学辑》(第2期);A002-134 * |
Also Published As
Publication number | Publication date |
---|---|
CN111859704A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111859704B (zh) | 一种分布式多视角下非刚体目标电磁散射建模方法 | |
CN108318880B (zh) | 一种运动误差参数化的极坐标sar成像方法 | |
CN108872986B (zh) | 一种方位角波数均匀化处理的极坐标sar运动补偿成像方法 | |
CN108594230A (zh) | 一种海船场景的合成孔径雷达图像仿真方法 | |
Peng et al. | Polar format imaging algorithm with wave-front curvature phase error compensation for airborne DLSLA three-dimensional SAR | |
CN108872985A (zh) | 一种近场圆周sar快速三维成像方法 | |
CN109188384B (zh) | 空间目标回波动态观测的电磁仿真方法 | |
CN107271977A (zh) | 基于移动激励源fdtd算法的高精度sar回波仿真方法 | |
CN114740469A (zh) | Isar回波实时精细模拟生成方法、装置和存储介质 | |
CN109991608B (zh) | 一种基于逆合成孔径原理的水下目标声纳成像方法 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN106908789B (zh) | 一种基于空间域支撑区融合的sar成像方法 | |
CN109188436B (zh) | 适用于任意平台轨迹的高效双基sar回波生成方法 | |
CN103630904B (zh) | 机载下视阵列3-d sar成像方法 | |
Zhuravlev et al. | Microwave imaging of concealed objects with linear antenna array and optical tracking of the target for high-performance security screening systems | |
CN104698442A (zh) | 机载下视阵列三维合成孔径雷达分布式三维场景仿真方法 | |
CN109884621B (zh) | 雷达高度表回波相参积累方法 | |
CN108957454B (zh) | 机载干涉合成孔径雷达的干涉相位仿真方法 | |
Elsherbini et al. | Image distortion effects in SAR subsurface imaging and a new iterative approach for refocusing and coregistration | |
JP2003248049A (ja) | レーダ画像模擬装置 | |
Dumont et al. | DIONISOS. An end-to-end SAR simulator | |
Wielgo et al. | Convolutional neural network for 3D ISAR non-cooperative target recognition | |
Zhou et al. | A Novel Method of Three-Dimensional Geometry Reconstruction of Space Targets Based on the ISAR Image Sequence | |
CN111460653A (zh) | 分析多个高速运动旋转对称目标电磁特性的区域分解方法 | |
Feng et al. | Multiview ISAR imaging for complex targets based on improved SBR scattering model |
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 |