CN103425864B - 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 - Google Patents
应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 Download PDFInfo
- Publication number
- CN103425864B CN103425864B CN201310187609.XA CN201310187609A CN103425864B CN 103425864 B CN103425864 B CN 103425864B CN 201310187609 A CN201310187609 A CN 201310187609A CN 103425864 B CN103425864 B CN 103425864B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- group
- 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
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，根据混合结构的散射特性，目标上的总电场等于入射场与散射场之和，均匀平面波通常被用来作为入射电场，散射电场可以用待求的电通密度和感应电流密度来表示；将自由空间中格林函数基于加法定理展开，给出远场部分的聚合因子，转移因子和配置因子的具体表达形式；结合体面积分方程的特性，利用射线追踪的基本原理，合理的选择出转移因子中角谱分量大的部分，将角谱分量小的部分舍弃掉；在更高层采用快速远场近似的方法计算远场作用；矩阵方程求解以及电磁散射参数的计算。本发明对于求解金属复杂非均匀媒质混合目标散射问题需要更少的计算内存以及计算时间。
Description
技术领域
本发明属于目标电磁散射特性的快速计算技术领域，具体涉及一种应用于金属复杂非均匀媒质混合目标的电磁散射分析方法。
背景技术
随着金属复杂非均匀媒质混合目标越来越多的出现在实际应用中，提出一种精确而有效的电磁分析方法就显得极为重要。在计算电磁学领域也相继出现了一些解决此类结构的电磁散射特性问题的数值计算方法。
针对这种金属介质混合结构，金属部分通常被作为理想导电体（PEC）来处理,并且容易被面积分方程方法（SIE）来分析求解，其中RWG基函数（Rao M, Wilton D and GlissonA. Electromagnetic s cattering by surfaces of arbitrary shape. IEEE Transaction on Antennas and Propagation, 1982, 30(3): 409– 418. ）由于其灵活性通常被被用来作为展开未知电流的基函数。介质部分，通常使用体积分方程方法（Schaubert D,Wilton D an d Glisson A. A tetrahedral modeling method for elect romagneticscattering by arbitrarily shaped inhomogeneo us dielectric bodies. IEEETransaction on Antennas an d Propagation, 1984, 32(1): 77–85.）进行分析。
由于复杂非均匀媒质部分采用体剖分造成未知量巨大，即便采用了多层快速多级子（MLFMM） （Tao, S.F. Ding, D.Z.;Chen, R.S. ;Chen, M. Electromagnetic analysisof electrically lar ge and arbitraried shaped FSS using hybrid VSIE comb inedwith parallelized MLFMM. ICMMT 2010, Chengdu, 61 4–616.）加速的矩量法进行求解，依然面临者计算资源消耗高的问题。
发明内容
本发明的目的在于提出本一种多层快速多极子（MLFMM）、快速远场近似（RPFMA）以及射线追踪（FAFFA）相结合的分析金属复杂非均匀媒质混合目标的电磁散射问题的快速分析方法。由于对远场部分采用更高效的计算方式，本发明对于求解金属复杂非均匀媒质混合目标散射问题需要更少的计算内存以及计算时间。
实现本发明目的的技术方案为：
一种应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，具体步骤如下：
第一步，根据混合结构的散射特性，目标上的总电场等于入射场与散射场之和，入射电场为已知激励，均匀平面波通常被用来作为入射电场，散射电场可以用待求的电通密度和感应电流密度来表示，
其中，ZDD代表介质对介质的作用， ZDM表示介质对金属的作用， ZMD都表示介质和金属的相互作用部分，ZMM示金属对金属的作用；
第二步，将自由空间中格林函数基于加法定理展开，结合体面积分方程的表达式，给出远场部分的聚合因子，转移因子和配置因子的具体表达形式；
第三步，结合体面积分方程的特性，利用射线追踪的基本原理，合理的选择出转移因子中角谱分量大的部分，将角谱分量小的部分舍弃掉；
第四步，结合体面积分方程的特性，在更高层采用快速远场近似的方法计算远场作用；
第五步，矩阵方程求解以及电磁散射参数的计算。
优选地,所述体面积分方程中，对转移因子进行加窗处理，采用的窗函数表示如下：
其中，L为球坐标下θ方向L个离散点的个数，l为每个转移因子分量的编号，经过窗函数处理过后的转移因子中较小的各个分量变得更小，设定每层转移因子角谱分量阈值，令θe为设定阈值对应的角谱分量方向与远场作用组中心连线方向的最大夹角，提取并计算角谱分量方向与两作用组中心连线方向夹角不大于θe的角谱分量。
优选地,所述体面积分方程中，把快速远场近似、射线追踪和多层快速多级子方法相结合，基于远场组之间的距离将每一层组与组之间采用快速多极子加速远场作用、射线多极子加速远场作用或者快速远场近似加速远场作用的方法进行处理，从而在整体上实现加速远场计算的目的。
本发明与现有技术相比，其显著优点：本发明是应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，对于求解金属复杂非均匀媒质混合目标散射问题,由于格林函数是自由空间的格林函数，本发明中采用了基于格林函数加法定理展开的快速多级子技术，加速了矩阵求解速度，同时节约了计算内存，并且结合体面积分方程的特性，利用射线追踪的原理，合理的舍去了转移因子较小的分量，对于更高层远场组之间的作用，仅需计算组与组中心连接线方向的转移分量，进一步减少了计算时间和计算内存。
附图说明
图1是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法中的埃瓦耳德球示意图。
图2是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法中多层快速多极子对应的角谱分量示意图。
图3是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法中射线多极子对应的角谱分量示意图。
图4是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法中快速远场近似对应的角谱分量示意图。
图5是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法中双战RCS对比曲线示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明，使本发明的上述及其它目的、特征和优势将更加清晰。
图1是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法中的埃瓦耳德球示意图，图2是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法多层快速多极子对应的角谱分量示意图，图3是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法射线多极子对应的角谱分量示意图，图4是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法快速远场近似对应的角谱分量示意图，如图1至图4所示，本发明提供了一种应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，具体步骤如下：
第一步，根据混合结构的散射特性，目标上的总电场等于入射场与散射场之和，入射电场为已知激励，均匀平面波通常被用来作为入射电场，散射电场可以用待求的电通密度和感应电流密度来表示，
其中，ZDD代表介质对介质的作用， ZDM表示介质对金属的作用， ZMD都表示介质和金属的相互作用部分， ZMM表示金属对金属的作用。表示形式如下，
其中，和分别代表体和面测试基函数， ω为电磁波角频率，和分别为，，和表示待求的电通量密度和金属面电流密度，是自由空间的格林函数。
上式中右边向量是由平面波产生的，可以写成
第二步，对自由空间的格林函数按照加法定理展开，假设在给定的一层中，任一组m中有一个场点ri，任一组n中有一个源点rj，rm、rn分别是场点组和源点组的中心点；两个非空组中场点和源点的空间矢量记为：rij=ri-rj=rim+rmn+rnj，当所分析的场点组和源点组不重合也不相邻时，满足|rim+rnj|＜|rmn|，自由空间的标量格林函数可在角谱空间写成：
整个积分是定义在单位球SE上的，αmn是两个组的转移因子，定义如下：
其中，jl(kd)为球贝塞函数，为第二类球汉克尔函数，为勒让德函数；L是无限级数的截断长度，L=kD+β(kD)1/3，D是分组的尺寸，β为精度参数，β≥2；用（θ，φ）表示单位球SE坐标，积分点数为KL= 2L2，其中在θ方向共采集L个点的一维高斯积分，φ方向共采集2L个点数的梯形法则积分；
给出ZDD、ZMD、ZDM和ZMM聚合因子，转移因子和配置因子的具体表达式分别为：
以上分析可以看出，在体面积分方程中应用快速多极子算法时，所有的单位球上的角谱分量都要用于转移过程，用(θ,φ)来描述单位球( 0≤θ≤π,0≤φ≤2π)，可以在θ方向(0,π)得到L个离散的点，在φ方向(0,2π)得到2L个离散的点，L是与层数以及目标电尺寸有关的参量，且层数与电尺寸越大，L的值越大。因此，当金属介质混合目标的尺寸非常大时，仅ZDD项所需要的转移量就很大，另外还有ZDM、ZMD和 ZMM三项，可知，转移因子所占内存非常大，转移过程是相当耗时的。
第三步，利用射线追踪的基本原理，结合体面积分方程的特性合理的选择出转移因子中角谱分量大的部分，将角谱分量小的部分舍弃掉。转移因子是两个远场作用组中心距离和组的大小的函数，当两个远场作用组中心距离和组的大小确定时，转移因子的特性类似于线性天线阵列的辐射方向图特性，将转移矩阵加上一个窗函数：
其中，L为球坐标下θ方向L个离散点的个数，l为每个转移因子分量的编号。具体的做法是：在转移因子上乘上窗函数的作用，使得两组中心连线方向附近的转移分量的值相对于其它部分转移因子的值有陡峭的变化，在这种情况下，可以省略对转移过程影响较小的角谱分量的作用，从而使得射线多极子获得更高的效率。同时，可以验证这种处理对计算精度影响很小。接着需要设定需要转移的角谱分量的范围。如图1所示，θr为两作用组连线方向与z轴的夹角，φr为两作用组连线方向向量在 XOY面上投影和x轴的夹角，θe为设定的阈值对应的角谱分量方向与两作用组连线方向的夹角。在计算转移因子的过程中，转移因子角谱分量的大小随着与远场作用组中心连线方向的夹角的增大而减小。这样可以对转移因子设定阈值对圆锥区域的大小进行判定。在具体操作中，只需提取并计算角谱分量方向与两作用组中心连线方向夹角不大于θe的角谱分量。在算法验证时可以得到这部分角谱分量可以很好地联系两个远场作用组的作用。
第四步，结合体面积分方程的特性，在更高层采用快速远场近似的方法快速计算远场作用。快速远场近似方法中，当场源之间的距离满足一定的条件时，我们可以只计算场源连线方向上的角谱分量，这样将转移量降到最低。本发明中采用组与组之间采用快速远场近似的判据条件为：
rmn＞3γdlevel （15）
其中dlevel为该层组的电尺寸，γ为结合体面积分方程特点设置的一个经验值，本发明中取为1.5。由此可以将快速远场近似、射线追踪和多层快速多级子方法相结合，基于远场组之间的距离可以将每一层组与组之间的作用采用三种快速方法进行处理：1快速多极子加速远场作用； 2射线多极子加速远场作用；3快速远场近似加速远场作用。从而在整体上实现加速远场计算的目的，本发明中将这种方法称为MLFMA-RPFM A-FAFFA方法，该方法具有进一步简化了转移因子的计算和存储过程的优点，具有重大的工程意义。
第五步，矩阵方程求解以及电磁散射参数的计算。
为了验证方法的效率和精度，下面给出了电大尺寸的金属媒质混合目标的电磁散射的算例，把使用RPFMA+FAFFA+MLFMM方法计算结果与纯M LFMM计算结果相对比，表1中可以看出该方法的高效性。
一个200m*0.2m的金属长条，未知量为10007。用多层快速多极子分成 8层计算，入射平面波频率为300MHZ，垂直入射，TM极化，入射方向沿 Z轴负方向，第四层开始使用RPFMA，其中LC 表示从该层开始使用FA FFA近似处理，图 5是本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法双战RCS对比曲线示意图，给出了其归一化的RCS对比图，并给出了两种计算方法的相关参数对比，如表1所示，
| 方法 | 迭代时间(s) | 远场内存消耗(MB) |
| MLFMM | 44.3 | 1059 |
| 本发明方法（Lc=4） | 3.65 | 70.7 |
| 本发明方法（Lc=5） | 2.03 | 65.8 |
| 本发明方法（Lc=6） | 1.37 | 44.7 |
表1
本算例是在Microsoft Visual Studio 2005编译器release模式调用八个进程下运行的结果。
本发明的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，有益效果如下：实现过程简单，在多层快速多级子的基础上，只需对转移因子做处理即可；所需计算资源少，射线追踪方法中由于只需要转移场源连线方向周围的角谱分量，相较于多层快速多级子中对单位球上的所有角谱分量进行转移，节省了转移因子的计算时间和所需的存储内存。同时当满足快速远场近似的条件下，仅需计算转移因子在场源组中心连线方向上的角谱分量，在满足计算精度的条件下进一步节省了计算资源。
在以上的描述中阐述了很多具体细节以便于充分理解本发明。但是以上描述仅是本发明的较佳实施例而已，本发明能够以很多不同于在此描述的其它方式来实施，因此本发明不受上面公开的具体实施的限制。同时任何熟悉本领域技术人员在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围内。
Claims (4)
1.一种应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，其特征在于，具体步骤如下：
第一步，根据混合结构的散射特性，用待求的电通密度和感应电流密度来表示金属复杂非均匀媒质混合目标上的散射电场，目标上的总电场等于入射场与散射场之和，入射电场为已知激励，均匀平面波被用来作为入射电场，
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>D</mi>
<mi>D</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>M</mi>
<mi>D</mi>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>D</mi>
<mi>M</mi>
</mrow>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>M</mi>
<mi>M</mi>
</mrow>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>v</mi>
<mi>m</mi>
<mi>V</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>v</mi>
<mi>m</mi>
<mi>S</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，ZDD代表介质对介质的作用，ZDM表示介质对金属的作用，ZMD都表示介质和金属的相互作用部分，ZMM表示金属对金属的作用；
第二步，将自由空间中格林函数基于加法定理展开，结合体面积分方程的表达式，给出远场部分的聚合因子、转移因子和配置因子的具体表达形式；
第三步，利用射线追踪的基本原理，选择出转移因子中角谱分量大的部分，舍弃掉角谱分量小的部分；
第四步，结合体面积分方程的特性，在树形结构更高层采用快速远场近似的方法计算远场作用；
第五步，矩阵方程求解以及电磁散射参数的计算。
2.根据权利要求1所述的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法，其特征在于，步骤1中，矩阵方程的具体表达形式如下：
其中，和分别代表体和面测试基函数，ω为电磁波角频率，和分别为和 和表示待求的电通量密度和金属面电流密度，是自由空间的格林函数；
式(1)中右边向量是由平面波产生的，写成
是入射电场。
3.根据权利要求1所述的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法,其特征在于，步骤二中结合体面积分方程的表达式，给出远场部分的聚合因子、转移因子和配置因子的具体表达形式的步骤如下：
步骤2.1，建立金属复杂非均匀媒质混合目标的树形结构；
步骤2.2，设在给定的一层中，任一组m中有一个场点ri，任一组n中有一个源点rj，rm、rn分别是场点组和源点组的中心点；两个非空组中场点和源点的空间矢量记为：rij＝ri-rj＝rim+rmn+rnj，当所分析的场点组和源点组不重合也不相邻时，满足|rim+rnj|<|rmn|，自由空间的标量格林函数可在角谱空间写成：
<mrow>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>k</mi>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</msup>
<mrow>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<munder>
<mo>&Integral;</mo>
<msub>
<mi>S</mi>
<mi>E</mi>
</msub>
</munder>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>k</mi>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>n</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mover>
<mi>k</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
整个积分是定义在单位球SE上的，αmn是两个组的转移因子，表示如下：
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>L</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mi>l</mi>
</msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msubsup>
<mi>h</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mi>r</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
<mrow>
<mo>(</mo>
<mover>
<mi>k</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，jl(kd)为球贝塞函数，为第二类球汉克尔函数，为勒让德函数；L是无限级数的截断长度，L＝kD+β(kD)1/3，D是分组的尺寸，β为精度参数，β≥2；用(θ，φ)表示单位球SE坐标，积分点数为KL＝2L2，其中在θ方向共采集L个点的一维高斯积分，φ方向共采集2L个点数的梯形法则积分，单位球SE上的每个积分点都被称为角谱分量；
步骤2.3，给出ZDD、ZMD、ZDM和ZMM聚合因子，转移因子和配置因子的具体表达式分别为：
步骤2.4，当两个远场作用组中心距离和组的大小确定时，将多个转移因子构成的转移矩阵加上一个窗函数：
<mrow>
<msub>
<mi>w</mi>
<mi>l</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>l</mi>
<mo>&le;</mo>
<mi>L</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>l</mi>
<mi>L</mi>
</mfrac>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>l</mi>
<mo>></mo>
<mi>L</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
在转移因子上乘上窗函数使得两组中心连线方向附近的转移分量的值相对于其它部分转移因子的值有陡峭的变化，舍弃角谱分量小的部分；
步骤2.5，θr为两作用组连线方向与z轴的夹角，φr为两作用组连线方向向量在XOY面上投影和x轴的夹角，θe为设定的阈值对应的角谱分量方向与两作用组连线方向的夹角；转移因子角谱分量的大小随着与远场作用组中心连线方向的夹角的增大而减小，计算角谱分量方向与两作用组中心连线方向夹角不大于θe的角谱分量。
4.根据权利要求1所述的应用于金属复杂非均匀媒质混合目标的电磁散射分析方法,其特征在于，步骤4采用快速远场近似方法计算远场作用中，当场源之间的距离满足式(15)的条件时，仅需计算场源连线方向上的角谱分量；
组与组之间采用快速远场近似的判据条件为：
rmn>3γdlevel (15)
其中dlevel为该层组的电尺寸，γ为设置的一个经验值，本发明中取为1.5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310187609.XA CN103425864B (zh) | 2013-05-20 | 2013-05-20 | 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310187609.XA CN103425864B (zh) | 2013-05-20 | 2013-05-20 | 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103425864A CN103425864A (zh) | 2013-12-04 |
| CN103425864B true CN103425864B (zh) | 2017-11-14 |
Family
ID=49650593
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310187609.XA Active CN103425864B (zh) | 2013-05-20 | 2013-05-20 | 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103425864B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104732050B (zh) * | 2013-12-20 | 2018-04-03 | 南京理工大学 | 雷电脉冲下碳纤维材料飞行目标内的电磁分布预估方法 |
| CN104778151B (zh) * | 2014-01-09 | 2017-11-14 | 南京理工大学 | 基于矩量法和抛物线方程的含腔目标电磁散射分析方法 |
| CN104778293B (zh) * | 2014-01-15 | 2018-06-12 | 南京理工大学 | 非均匀介质目标电磁散射的体积分Nystrom分析方法 |
| CN105095546A (zh) * | 2014-05-16 | 2015-11-25 | 南京理工大学 | 分析多尺度导体目标电磁散射特性的混合阶Nystrom方法 |
| CN105303022B (zh) * | 2014-07-28 | 2018-04-03 | 南京辉宏电子信息有限公司 | 快速获取目标电磁散射特性的高斯波束方法 |
| CN106202594B (zh) * | 2015-05-07 | 2019-12-27 | 南京理工大学 | 分析混合目标瞬态电磁散射特性的时域不连续伽辽金方法 |
| CN106202599B (zh) * | 2015-05-07 | 2019-12-27 | 南京理工大学 | 加速求解混合目标电磁散射特性的多层复点源方法 |
| CN106326659B (zh) * | 2016-08-26 | 2019-04-26 | 上海无线电设备研究所 | 基于射线追踪的金属介质组合目标电磁散射计算方法 |
| CN106597380B (zh) * | 2016-11-28 | 2019-12-06 | 上海无线电设备研究所 | 一种基于非均匀波的角度域电磁散射特性合成方法 |
| CN108920741B (zh) * | 2018-04-28 | 2022-10-21 | 南京航空航天大学 | 目标发生局部变化的电磁散射快速分析方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901965A (zh) * | 2012-09-26 | 2013-01-30 | 同济大学 | 一种海洋杂波环境中电大尺寸目标雷达成像方法 |
-
2013
- 2013-05-20 CN CN201310187609.XA patent/CN103425864B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102901965A (zh) * | 2012-09-26 | 2013-01-30 | 同济大学 | 一种海洋杂波环境中电大尺寸目标雷达成像方法 |
Non-Patent Citations (3)
| Title |
|---|
| Efficient MLFMA,RPFMA and FAFFA algorithms for EM scattering by very large structures;Tie Jun Cui等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20040331;第52卷(第3期);759-769 * |
| 导体介质复合目标的电磁散射特性分析;查丽萍等;《2011年全国微波毫米波会议论文集（下册）》;20110601;1229-1232 * |
| 非均匀复杂结构目标电磁散射理论建模与高效算法研究;何十全;《中国博士学位论文全文数据库信息科技辑》;20111215;I135-3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103425864A (zh) | 2013-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103425864B (zh) | 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 | |
| CN102129523B (zh) | 基于mda和mlssm的分析复杂目标电磁散射的方法 | |
| CN102156764B (zh) | 一种分析天线辐射和电磁散射的多分辨预条件方法 | |
| Eibert et al. | Fast spectral domain algorithm for hybrid finite element/boundary integral modelling of doubly periodic structures | |
| CN106294894A (zh) | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 | |
| CN112329285B (zh) | 泰勒级数展开的混合算法分析多尺度结构瞬态响应方法 | |
| CN104200074A (zh) | 快速获取目标电磁散射特性的多层复波束方法 | |
| You et al. | A new method to analyze broadband antenna-radome interactions in time-domain | |
| CN103425816A (zh) | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 | |
| Wu et al. | Performance enhanced Crank-Nicolson boundary conditions for EM problems | |
| CN103235193B (zh) | 毫米波段内卫星电磁散射特性的数值方法 | |
| CN105303022B (zh) | 快速获取目标电磁散射特性的高斯波束方法 | |
| CN108038313A (zh) | 一种剖分不均匀的目标电磁散射特性的分析方法 | |
| CN104573257A (zh) | 基于面元空间分集的电磁遮挡判断快速算法 | |
| CN112329204A (zh) | 考虑载体平台耦合的重复性结构电磁特性特征模快速分析方法 | |
| CN104915324B (zh) | 腔体含介质目标电磁散射混合分析方法 | |
| Rasool et al. | Radiowave propagation prediction in the presence of multiple knife edges using 3D parabolic equation method | |
| CN112014815A (zh) | 混合型面积分的分析方法、装置、设备及存储介质 | |
| Chen et al. | Multilevel fast dipole method for electromagnetic scattering from perfect electric conducting targets | |
| Chen et al. | Accelerating the multilevel fast multipole method with parallel preconditioner for large-scale scattering problems | |
| Jiang et al. | Preconditioned MDA-SVD-MLFMA for analysis of multi-scale problems | |
| CN111931353A (zh) | 一种应用于仿真fss结构的散射场求解方法 | |
| CN106294898B (zh) | 一种加速分析介质目标电磁散射特性的复点源求解方法 | |
| Gu et al. | Parallel Integral Equation Based Non-overlapping DDM for Fast Solving Electromagnetic Scattering Problems with Changeable Parts | |
| CN106156394A (zh) | 基于显式差分格式的电磁特性提取方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |