圆柱周期特性的介质目标的频域电磁散射特性分析方法
【技术领域】
[0001] 本发明属于目标电磁散射特性数值计算技术，特别是一种分析圆柱周期结构介质 目标的频域电磁散射特性的方法。
【背景技术】
[0002] 在现代科学工程应用中，由于各种电子设备的工作频率的不断升高，很多物体的 电尺寸（物理尺寸与电磁波波长的比值）也越来越大，从而导致了采用数值计算方法处理 问题时产生的未知量也越来越大，消耗的计算机内存也成倍增加，造成计算资源的极大消 耗。通常处理的许多电大尺寸复杂目标的电磁散射求解问题中，有许多研究目标在结构上 具有某些特性，如旋转对称体、圆柱周期结构（即离散的旋转对称体）等，这种特性可以在 数值计算的算法中充分利用，可以解决未知量较大的问题。传统的分析介质目标电磁散射 特性方法有体积分方法、面面积分方法。它们都没有利用这种圆柱周期结构特性，都需要对 整个物体进行建模剖分。也就意味着需要更多的计算资源。
【发明内容】
[0003] 本发明的目的在于提供一种分析圆柱周期结构介质目标的频域电磁散射特性的 数值方法。
[0004] 实现本发明目的的技术解决方案为：一种圆柱周期特性的介质目标的频域电磁散 射特性分析方法，步骤如下：
[0005] 第一步，介质体目标表面PMCHff积分方程的建立。利用等效原理，将介质问题分成 对介质体外部场等效和对介质体内部场等效两个独立的问题。通过分界面上的等效电流和 磁流可得这些等效电流和磁流在内域或外域中所产生的散射场。根据等效原理，由于仅在 外域存在着入射场（相当于外加激励源），而内域中并不存在，因此在分界面S上强迫使用 切向连续性条件，则可以得到介质体目标表面PMCHff积分方程。
[0006] 第二步，根据圆柱周期结构目标的旋转对称性，将目标分成若干个相同的区域，其 他区域点的位置就可以由第一个区域的点的位置旋转得到。所以只需要对第一个区域进 行建模就可以了。
[0007] 第三步，阻抗矩阵计算，对积分方程进行伽辽金测试，得到矩阵方程。由于旋转对 称特性，只需要填充第1个区域的自作用以及与其它区域的互作用矩阵元素即可。为了加 快矩阵的迭代求解，可以充分利用阻抗矩阵的特性，因此利用快速傅里叶变换技术来来加 速矩阵迭代求解时矢量乘操作。
[0008] 第四步，矩阵方程求解以及电磁散射参数的计算。
[0009] 本发明与现有技术相比，其显著优点：（1)建模简单。只需要对周期结构中的一个 完整周期区域进行建模即可，避免了建立繁复的模型。（2)减少未知量。因为只对周期结构 中的一个完整周期区域进行建模，所以未知量大大减少。（3)节省内存。矩阵元素填充时， 由于旋转对称特性，只需要填充第1个区域的自作用以及与其它区域的互作用矩阵元素即 可，因而节省了内存。
【附图说明】
[0010] 图1是本发明圆柱周期结构的几何示意图。
[0011] 图2是本发明介质环双站RCS曲线示意图。
【具体实施方式】
[0012] 下面结合附图对本发明作进一步详细描述。
[0013] 结合图1，本发明圆柱周期特性的介质目标的频域电磁散射特性分析方法，步骤如 下：
[0014] 第一步，根据等效原理，介质目标PMCHff方程的建立。
[0026] 第二步，旋转对称特性。
[0027] 圆柱周期结构如图1所示，整个结构是由Ns个小区域组成的，每个区域的旋转角 度为灼=21/假设第1个周期区域为VO部分，整个结构的各个区域可以由VO区域绕z 轴旋转一定的角度 < 得到。
[0028] 在建模时只需要对VO区域建模剖分即可，其它部分可以有几何关系推导出来，这 样就节省了大量未知量和内存消耗。VO区域上点的向量为r'，那么第η个区域上对应点 的向量为：
[0053] 循环矩阵C与向量x的矩矢乘运算可以写为循环卷积形式，从而可以进一步利用 快速傅里叶变换（FFT)来加速运算，表示如下：
[0054] Cx = c*x = FFT 1 [FFT (c) XFFT (x) ] (25)其中 c = [c0cv.. cn J，x = [X0X1... Xn J，星号*表示循环卷积运算。
[0055] 相应的阻抗矩阵调整为：
[0056]
[0062] 循环矩阵矢量乘：
[0063] {V} = [Z] {1} (29)
[0064] 可以转化为循环卷积运算：
[0065] GO)
[0066] 其中表TK循环卷积算符，其它变量与原变量的关系如下：
[0067] T' = T'1 I ^ n ^ Ns (31)
[0068] 由式（28)可以看出循环矩阵[Ζ]只需要存储第一列即可，内存需求大大减少了， 而且这里的循环卷积可以采用快速傅里叶变换来加速求解。
[0069] \V\(32)
[0070] 式中FFT表示的是快速傅里叶变换处理，FFT 1表示的是快速傅里叶逆变换处理。
[0071] 第四步，求解矩阵方程，得到电流系数，再根据互易定理由电流系数计算电磁散射 参量。
[0072]为了验证本发明的正确性与有效性，下面分析了圆柱周期结构介质目标电磁散射 的典型算例。算例在主频2. 83GHz、内存3. 37GB的个人计算机上实现。
[0073] 目标为外半径0. 8m，内半径0. 6m，厚度0. 2m的介质圆环状结构，介电常数为4,入 射波采用均匀平面波，入射频率为150MHz，θ Θ极化，入射角度为Θ i = 〇°，Cj5i = 〇°，散 射角为〇°彡9s<180°，Φ3 = 0°。采用平面三角形剖分，剖分尺度0. 1波长。采用离 散旋转对称体电磁散射的频域矩量法建模：介质环等分成64份相同区域，每个区域旋转角 度为5. 625度时，DBOR建模只需要对一个区域进行建模并处理好周期性边界处的三角形即 可。DBOR建模得到未知量为48。图2给出了双站RCS值，可以看出本发明的结果与FEKO 吻合的很好。
【主权项】
1. 一种圆柱周期特性的介质目标的频域电磁散射特性分析方法，其特征在于步骤如 下： 第一步，根据等效原理，得到介质体目标表面PMCHW积分方程为：￡中e和 d分别表示介质体目标外部与内部，和h分别为真空中的磁导率以及介电参数，1^和ed分别为介质的磁导率以及介电参数，r和r'分别为场和源的位置坐标； 1表示e或者d，A: (r)和Fi(r)分别为矢量磁位与矢量电位函数，?彳分别表 示标量电位与标量磁位函数，如下式：第二步，设圆柱周期结构整个结构由队个子区域组成，每个子区域的旋转角度为 % =ZT/iV，设第1个子区域为V0部分，整个结构的其它子区域由V0区域绕z轴旋转角度 乾得到，每个区域被离散为Q个面片，则表面电磁流表示如下：第三步，阻抗矩阵计算；将式（3) (4)代入式（1) (2)，通过伽辽金测试， 得到矩阵方程：由i的数学特性得到：(?) 把式（7)代入式（6)得到：对于阻抗矩阵Z'ZEM，Z'Z?，分别填充第一行的子矩阵中的元素； 将式（8)的阻抗矩阵调整为：式（9)中的每个NSXNS的[Z]pq都是一个循环矩阵，即第i行的第1至第Ns-1个元素 右移一列，就得到了第i+1行的第2至第Ns个元素，第i+1行的第一个元素等于第i行的 第Ns个元素；在存储时只需要每个小矩阵[Z]pq的第一行的全部元素的值即可； 第四步，求解矩阵方程，得到电磁流系数，再根据互易定理由电流系数计算电磁散射参 量。2.根据权利要求1所述的分析圆柱周期结构的介质目标的频域电磁散射特性分析方 法，其特征在于：所述步骤2和3中的旋转对称特性如下： 阻抗矩阵ZE\ZEM，Z'Z?都具有块循环特性，阻抗计算时只需要计算填充第一行的子矩 阵元素，进而利用快速傅里叶变换技术来加速矩阵迭代求解时矢量乘操作。
【专利摘要】本发明公开了一种圆柱周期特性的介质目标频域电磁散射特性分析方法。采用面面积积分方程（PMCHW）计算圆柱周期结构介质目标的电磁散射特性时，需要对整个的周期结构目标进行电磁建模，尤其是对有较多周期的复杂结构，建模方式比较繁琐。针对具有圆柱周期特性的目标，采用离散旋转对称法（DBOR）计算，只需要对周期结构中的一个完整周期区域进行建模即可，避免了建立繁复的模型，降低剖分未知量，进而大大减少了矩阵填充时间并节省内存。
【IPC分类】G06F17/50
【公开号】CN105095541
【申请号】CN201410201521
【发明人】陈如山, 樊振宏, 丁大志, 程光尚
【申请人】南京理工大学
【公开日】2015年11月25日
【申请日】2014年5月13日