预条件技术和迭代算法在高效实现MLFMA中的对比研讨.pdfVIP

2005’全国微波毫米波会议论文集 预条件技术和迭代算法在高效实现MLFMA中 的对比研究 董健，柴舜连，毛钧杰 (国防科学技术大学电子科学与工程学院，湖南长沙，410073) 摘 要：本文使用了不同的迭代算法和预条件技术结合MLFMA来加速计算电磁场积 空问迭代算法。使用的预条件算法包括主对角预条件、块对角预条件、ILuT预条件和基于 MLFMA近场矩阵的完全LU分解预条件。文中通过不同的算例比较了不同迭代算法和预条 件技术与MLFMA算法结合后对整个计算过程的收敛速度和内存需求的影响。 关键词：积分方程；矩量法；多层快速多极算法；迭代算法；预条件方法 1引言 近年来，快速多极方法(FMM)被广泛应用于电磁场积分方程迭代求解过程中[1]。作为 20世纪十大算法发明之一，其多层实现算法(MLFMA)可以把迭代法中的矩阵向量乘积运算 N1，存储量需求亦降为 降为O(NLogN)，进而整个求解过程的计算量降为O(N“”N log 本身的一些参数选择影响整个计算过程的精度、效率、复杂度和存储需求外，还有一些重要 的因素会对整个计算过程产生很大的影响，包括： (1)选择哪种类型的积分方程来分析电磁问题； (2)使用何种迭代算法来求解方程； (3)迭代过程中是否使用预条件技术，使用何种预条件算法； (4)迭代时初值的选择； 在笔者实现MLFMA算法的过程中，为验证不同迭代算法和预处理技术对整个计算过程的 主对角预条件、块对角预条件、ILUT预条件和对MLFMA近场矩阵作完全Lu分解做预条件。 我们使用了大量的计算实例对这些迭代算法和预条件技术做了验证和比较，这些实例包括标 准的验证模型(球、NASA杏仁体等)、复杂的验证模型(Slicy模型)、真实目标模型(飞机、 导弹模型)等。 2电磁场积分方程与多层快速多极算法 对于理想导体目标的散射问题，可以采用求解电磁场积分方程的方法来解决。使用矩量 法，可把积分方程离散化，得到一复系数稠密线性方程组。可以使用迭代法求解这一方程组， 存储量和计算量需求为研酽1。快速多极算法是一种用以快速计算多个点源的场的算法[2]。 其可被用来加速用迭代法求解积分方程时的矩阵向量乘积的运算。矩量法离散积分方程得到 的系数矩阵与待求电流向量的乘积可认为是计算由N个源点产生的场在N个场点的值。直接 计算N个场点的值需要的计算量为研Ⅳz1。通过把场点和源点进行空间分组，把场源点之间 s1。把空间多层分组，采用层迭算 的作用转化为组与组之间的作用，可把计算量降为O(N1 法，可把计算量进一步降为研Nlog(N))。 4Krylov子空间迭代算法和预条件技术 在求解大规模线性方程组的迭代算法中，Krylov子空问迭代算法是目前最为有效的一 类算法[3]。已经有很多Krylov子空间迭代算法被用于求解大型线性方程组，如CG、BiCG、 1017 2005’全国微波毫米波会议论文集 收敛性能要好一些，但其代价为存储量需求较大，其参数m正是为控制存储量而设置的。 使用迭代法求解线性方程组，当其系数矩阵是病态矩阵或近似病态矩阵，即条件数比较 大时，求解过程中迭代算法的收敛速度很慢，甚至不能收敛到准确的结果。可通过预条件方 法来解决这一问题。预条件矩阵一般由原线性方程组的系数矩阵产生。在用MLFMA来加速分 析电磁问题时，系数矩阵z被分为近场矩阵刀一和远场矩阵z加两部分，只有近场矩阵被显 式的计算和存储，远场矩阵与向量的乘积运算借助于MLFMA来完成。因而只能由稀松矩阵 z～来构造预条件矩阵M。我们实现了几种预条件算法，包括： (1)主对角预条件：M的主对角元素与z””相同，其余元素全为0。。 (2)块对角预条件：M为分块对角矩阵，各予块内部与Z”“对应元素相同，其余为O。其 存储量需求为lxm2，l为分块数目，m为子块维数，且有关系m=N／l。 (3)ILUT预条件：M通过对近场矩阵z””做ILU分解来得到，在分解过程中，通过两