加速求解混合目标电磁散射特性的多层复点源方法技术

本发明专利技术公开了一种加速求解混合目标电磁散射特性的多层复点源方法，步骤如下：建立体面积分方程；对介质部分未知量利用SWG基函数离散，金属部分未知量利用RWG基函数离散，并利用伽辽金方法对体面积分方程进行测试，形成待求解的矩阵方程组；引入多层复点源方法，加速矩阵方程组求解；根据求解出的介质体内部和金属表面上的感应电流，确定目标的电磁散射特性。本发明专利技术方法中，互为近场的组对之间相互作用利用传统矩量法求解，互为远场的组对之间相互作用可以引入多层复点源方法加速求解，基于多层复点源加速的体面积分方程可以更加高效的求解混合目标的电磁散射问题，降低计算复杂度。

【技术实现步骤摘要】
一
本专利技术属于电磁仿真
，特别是一种加速求解混合目标电磁散射特性的多层复点源方法。二
技术介绍
目标电磁散射特性的获取与分析是电磁问题中的一个非常重要研究领域，目标的电磁散射波是雷达探测、遥感观测以及地质勘测邓众多应用的信息来源，散射特性的定量分析是这些应用系统在设计和工作时的主要依据。金属介质混合目标在现实中有较广泛的应用领域，例如，金属目标涂覆隐身材料，FSS等等。因此，对其电磁散射特性的分析具有特别重要的现实意义。体面积分方程(VSIE)(Lu CC,Chew WC,“A coupled surface-volume integral equation approach for the calculation of electromagnetic scattering from composite metallic and material targets,”IEEE Transactions on Antennas Propagation,vol.48,pp.1866–1868,2000)作为一种灵活有效的方法，可以用来求解金属介质混合目标的电磁特性，特别是对于分析含有不均匀介质的混合目标。但是，随着求解目标的未知量越来越大，传统的体面积分方程将无法高效地进行求解。此时，需要利用某些高效的快速方法来打破瓶颈。复点源方法(CSB)(K.Tap,P.H.Pathak,and R.J.Burkholder,“Complex source beam-moment method procedure for accelerating numerical integral equation solutions of radiation and scattering problems,”IEEE Transactions on Antennas Propagation,vol.62,pp.2052–2062,2014)是近年来兴起的一种加速求解电磁散射问题的快速方法。利用八叉树分层分组思想，将待求目标上基函数之间的作用分为近场组对之间的作用和远场组对之间的作用，近场组对之间仍采用矩量法进行求解，远场组对之间利用复点源方法加速求解。将基函数等效到等效球面的复点源上，利用复点源之间的互作用代替远场组对之间基函数的相互作用。但是目前的复点源方法只应用于分析金属目标的电磁散射特性，对于介质目标或混合目标尚未研究。三
技术实现思路
本专利技术的目的在于提供一种能够用于介质目标或混合目标的加速求解混合目标电
磁散射特性的多层复点源方法。实现本专利技术目的的技术解决方案为：一种加速求解混合目标电磁散射特性的多层复点源方法，步骤如下：步骤1，建立体面积分方程；步骤2，对介质部分未知量利用SWG基函数离散，金属部分未知量利用RWG基函数离散，并利用伽辽金方法对体面积分方程进行测试，形成待求解的矩阵方程组；步骤3，引入多层复点源方法，加速矩阵方程组求解；步骤4，根据求解出的介质体内部和金属表面上的感应电流，确定目标的电磁散射特性。本专利技术与现有技术相比，其显著优点是：(1)将复点源思想推广到体面积分方程，将复点源方法与体面积分方程相结合，加速了体面积分方程的求解；(2)拓展了复点源思想的应用范围，将复点源思想推广到体积分方程，分析含介质的目标，将介质体上的SWG基函数用复点源进行展开，可以更加高效地求解含介质材料的混合目标的电磁散射特性；(3)利用多层的结构以及格林函数平移不变性的思想进一步扩大了复点源方法的效率。四附图说明图1为本专利技术中等效面上复点源分布的示意图。图2为本专利技术中互为远场组的复点源截取角度示意图。图3为本专利技术中介质涂覆金属球目标的示意图。图4为本专利技术中双站RCS与解析方法结果对比示意图。五具体实施方式下面结合附图对本专利技术作进一步详细描述。本专利技术加速求解混合目标电磁散射特性的多层复点源方法，步骤如下：步骤1，建立体面积分方程，具体如下：令电磁波照射到金属介质混合目标上，在介质体内产生感应体电流金属表面产生感应面电流根据电场边界条件，得到体面积分方程VSIE的表达式，如下： E → i = D → ( r → ) ϵ ( r → ) + jω μ 0 ∫ G ‾ ( r → , r → ′ ) J → v ( r → ′ ) d r ′ → + jω μ 0 ∫ G ‾ ( r → , r → ′ ) J → s ( r → ′ ) d r → ′ 本文档来自技高网...

【技术保护点】
一种加速求解混合目标电磁散射特性的多层复点源方法，其特征在于，步骤如下：步骤1，建立体面积分方程；步骤2，对介质部分未知量利用SWG基函数离散，金属部分未知量利用RWG基函数离散，并利用伽辽金方法对体面积分方程进行测试，形成待求解的矩阵方程组；步骤3，引入多层复点源方法，加速矩阵方程组求解；步骤4，根据求解出的介质体内部和金属表面上的感应电流，确定目标的电磁散射特性。

【技术特征摘要】
1.一种加速求解混合目标电磁散射特性的多层复点源方法，其特征在于，步骤如下：步骤1，建立体面积分方程；步骤2，对介质部分未知量利用SWG基函数离散，金属部分未知量利用RWG基函数离散，并利用伽辽金方法对体面积分方程进行测试，形成待求解的矩阵方程组；步骤3，引入多层复点源方法，加速矩阵方程组求解；步骤4，根据求解出的介质体内部和金属表面上的感应电流，确定目标的电磁散射特性。2.根据权利要求1所述的加速求解混合目标电磁散射特性的多层复点源方法，其特征在于，步骤1中所述建立体面积分方程，具体如下：令电磁波照射到金属介质混合目标上，在介质体内产生感应体电流金属表面产生感应面电流根据电场边界条件，得到体面积分方程VSIE的表达式，如下： E → i = D → ( r → ) ϵ ( r → ) + jω μ 0 ∫ G ‾ ( r → , r → ′ ) J → v ( r → ′ ) d r → ′ + jω μ 0 ∫ G ‾ ( r → , r → ′ ) J → s ( r → ′ ) d r → ′ ( r → ∈ V ) - - - ( 1 ) ]]> E → tan i = [ jω μ 0 ∫ G ‾ ( r → , r → ′ ) J → v ( r → ′ ) d r → ′ + jω μ 0 ∫ G ‾ ( r → , r → ′ ) J → s ( r → ′ ) d r → ′ ] tan ( r → ∈ S ) - - - ( 2 ) ]]>其中，表示照射在混合目标上的入射电场，表示并矢格林函数，表示电通量密度，μ0表示自由空间的磁导率，ε0表示自由空间的介电常数，ε表示介质的介电常数，ω表示角频率，为场的位置坐标，为源的位置坐标，c表示真空中的光速，tan表示金属表面的切向分量，V表示介质区域，S表示金属区域。3.根据权利要求1所述的加速求解混合目标电磁散射特性的多层复点源方法，其特征在于，步骤2中所述对介质部分未知量利用SWG基函数离散，金属部分未知量利用RWG基函数离散，并利用伽辽金方法对体面积分方程进行测试，形成待求解的矩阵方程组，具体如下：介质目标的体电流离散表示如下： J → v ( r → ) = jω Σ n = 1 N D K ( r → ) f → n V ( r → ) D n - - - ( 3 ) ]]>金属目标的面电流离散表示如下： J → s ( r → ) = Σ n = 1 N S f → n s ( r → ) I n - - - ( 4 ) ]]>其中，称为对比度，为SWG基函数，为RWG基函数,ND为介质部分未知量的个数，NS为金属部分未知量的个数，Dn为第n个介质未知量的待求系数，In为第n个金属未知量的待求系数；将式(3)、(4)带入式(1)、(2)，并利用伽辽金方法对对体面积分方程进行测试，形成了如下式所示的矩阵方程形式： Z DD Z MD Z DM Z MM D I = v V ...

【专利技术属性】
技术研发人员：陈如山，陈培林，丁大志，樊振宏，呼延龙，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏;32