内容简介
《电磁分析中的预条件方法》主要介绍了预条件方法的基本理论及其在电磁分析中的应用，包括计算电磁学中的主要数值方法、Krylov子空间迭代方法、预条件技术、迭代算法的自适应加速技术、预条件技术的优化措施、基于物理模型的预条件技术、基于特征谱信息的快速迭代算法及预条件技术、高阶有限元及多重网格迭代法、高阶矩量法及多重网格方法、块迭代算法、并行预条件技术等，重点介绍了多种预条件技术在矩量法和有限元方法中的应用。
目录
目录
前言
第1章 绪论 1
1.1 计算电磁学发展现状 2
1.2 迭代解法和预条件技术 5
1.3 内容安排 9
参考文献 10
第2章 计算电磁学中的主要数值方法 26
2.1 有限元法 26
2.1.1 电磁场边值问题 27
2.1.2 伽辽金加权余量法与里茨变分法 27
2.1.3 有限元法的步骤 28
2.1.4 数值结果 31
2.2 矩量法 35
2.2.1 矩量法的离散化过程 35
2.2.2 积分方程的选取 36
2.2.3 散射场的计算 39
2.2.4 多层快速多极子方法 40
2.2.5 并行多层快速多极子方法 44
2.2.6 数值结果 47
参考文献 49
第3章 Krylov子空间迭代方法 53
3.1 直接解法和迭代解法简介 53
3.2 迭代方法的分类 54
3.3 共轭梯度类迭代方法 55
3.4 广义*小余量迭代算法 56
3.5 常用Krylov子空间迭代算法的比较 58
3.6 常用迭代算法在体积分方程中的应用 58
3.7 常用迭代算法在表面积分方程中的应用 64
参考文献 68
第4章 预条件技术 70
4.1 预条件技术概述 70
4.2 稠密矩阵的稀疏化 71
4.3 预条件广义*小余量迭代算法 72
4.4 对角预条件技术 73
4.5 对称超松弛预条件技术 74
4.6 不完全LU分解预条件技术 74
4.7 稀疏近似逆预条件技术 75
4.8 几种常用预条件技术性能的比较 77
参考文献 82
第5章 迭代算法的自适应加速技术 84
5.1 GMRES迭代算法收敛性分析 84
5.2 基于GMRES迭代算法的自适应加速技术概述 86
5.3 Krylov子空间扩大技术 87
5.3.1 扩大子空间的广义*小余量迭代算法 87
5.3.2 松散的广义*小余量迭代算法 90
5.4 特征谱重复循环技术 92
5.4.1 隐式循环的广义*小余量迭代算法 92
5.4.2 显式循环的广义*小余量迭代算法 95
5.5 特征谱预条件的广义*小余量迭代算法 97
5.6 内外迭代技术 99
5.6.1 灵活的广义*小余量迭代算法 99
5.6.2 嵌套的广义*小余量迭代算法 102
5.7 几种加速技术性能的比较 103
5.8 其他迭代加速技术 107
参考文献 111
第6章 预条件技术的优化措施 114
6.1 对称超松弛预条件技术的有效实现 114
6.2 不完全LU分解预条件技术中的扰动技术 117
6.2.1 对角线扰动技术 117
6.2.2 MFIE主值项扰动技术 120
6.3 多层快速多极子方法中一种有效的稀疏近似逆预条件技术 125
6.4 混合预条件技术 129
6.4.1 双步混合预条件技术 129
6.4.2 SSOR预条件技术与GMRESR及FGMRES结合算法 134
6.5 多重预条件技术 138
6.5.1 多重预条件共轭梯度算法 138
6.5.2 多重预条件广义*小余量算法 139
6.6 预条件矩阵插值 142
6.6.1 基于有理函数模型的阻抗矩阵插值技术 142
6.6.2 基于有理函数模型的稀疏近似逆预条件矩阵插值技术 145
参考文献 149
第7章 基于物理模型的预条件技术 151
7.1 电场矢量有限元方程的病态特性 151
7.2 基于A-V场的预条件技术 153
7.2.1 A-V场有限元公式 153
7.2.2 数值结果与分析 155
7.3 基于转移Laplace算子的预条件技术 158
7.3.1 转移Laplace算子的预条件 158
7.3.2 数值结果与分析 160
7.4 基于吸收边界条件的预条件技术 167
7.4.1 快速多极子结合有限元方法理论及公式 167
7.4.2 利用吸收边界条件构造预条件矩阵 170
7.4.3 数值结果与分析 173
参考文献 177
第8章 基于特征谱信息的快速迭代算法及预条件技术 180
8.1 改进的扩大子空间广义*小余量迭代算法 180
8.1.1 GMRESE迭代算法基本原理 180
8.1.2 GMRESE迭代算法的收敛性能 182
8.1.3 GMRESE迭代算法的性能随参数变化情况 186
8.1.4 GMRESE迭代算法在单站RCS计算中的应用 189
8.2 基于特征谱信息的代数多重网格迭代算法 192
8.2.1 基于特征谱信息的代数多重网格迭代算法基本原理 192
8.2.2 SMG迭代算法的收敛性能 195
8.2.3 SMG迭代算法的性能随参数变化情况 197
8.2.4 SMG迭代算法在单站RCS计算中的应用 199
8.2.5 SMG性能随未知量变化情况 201
8.3 基于特征谱信息的多步混合预条件技术 203
8.3.1 基于特征谱信息的双步混合预条件技术的基本思想 203
8.3.2 基于特征谱信息的双步混合预条件技术的性能 205
8.3.3 基于特征谱信息的多步混合预条件 208
8.3.4 多步混合预条件技术在单站RCS计算中的应用 211
8.3.5 基于等级基函数的双步谱预条件技术 216
参考文献 222
第9章 高阶有限元及多重
试读
第1章 绪论
1864年，英国物理学家麦克斯韦(Maxwell)在“A dynamical theory of the electromagnetic field”一文中提出了位移电流的概念［1］，认为不仅变化的磁场能产生电场，而且变化的电场也能产生磁场，电磁场由此交替地向外传播，同时提出了麦克斯韦方程组，并首次预言了电磁波的存在。1887年，德国物理学家赫兹用火花放电的实验证实了麦克斯韦的预言，得出了电磁能量可以越过空间传播的结论。赫兹的发现轰动了全世界科学界，开启了人们研究电磁场理论及其应用的大门。随着研究的不断深入，电磁场的理论已经在航空航天、无线电通信、雷达、半导体芯片及封装、遥测遥感、生物医疗以及光电子等行业得到了广泛的应用。
目前，利用电磁波工作的通信、雷达设备已经在国防军事和民用领域发挥着重要的作用，与人们的日常生活息息相关。随着社会的发展，人们希望这些电磁设备的功能越来越齐全、结构越来越小型化、电磁波频谱的利用率越来越高。为了降低电磁设备的设计成本，提高设备性能，需要分析电磁设备的电磁特性，因此各种电磁设备的抗电磁干扰性及彼此间的电磁耦合研究日益引起人们的重视。理论上，以上需求都可归结为复杂的电磁分析问题，即研究电磁波和具有复杂材料构成及几何结构的物质之间的相互作用。实际设备模型的复杂性，使经典电磁学解析求解时存在较大的局限性，甚至根本无法实现。另一个有效的途径是，通过现场实验来研究实际的电磁问题，但在实际环境中开展实验的人力、物力成本巨大，研究周期长，有时因为实验条件的缺失而无法进行研究。
伴随着计算机技术水平的不断提高，利用计算机进行仿真的理论手段越来越受到重视。相应地，在电磁学领域内出现了多种分析电磁问题的数值算法，从而诞生了一门解决复杂电磁场理论与微波工程问题的新兴学科——计算电磁学［2，4］。在电磁学领域的科研工作者的不懈努力下，计算电磁学已经取得了丰硕的成果。随着多种新型算法的提出和改进，计算电磁学能够解决的问题也越来越复杂。计算电磁学可以看成数学方法、电磁场理论和计算机技术结合的产物。使用计算电磁方法，可避免高昂的现场测试费用，同时其数值解对电磁理论的发展和现场实验也有一定的指导意义。随着计算电磁学的发展，利用电磁仿真进行研究已经成为设计电路、天线、电磁兼容和电磁散射等实际工程应用领域内不可或缺的、甚至主要的研究手段。
1.1 计算电磁学发展现状
电磁场数值分析是根据麦克斯韦方程，利用适当的边界条件确定所关心区域内的电磁场或电流分布，进而求出所需要的物理参量。回顾计算电磁学发展的历史，早在1864年，麦克斯韦已用偏微分方程的形式给出了电磁波现象中电场和磁场的统一表达式，他的研究成果被誉为19世纪*显著的科学成就之一。而求解电磁场问题的方法归纳起来可分为三大类，第一类是解析法，第二类是数值法，第三类是半解析数值法，其中每一类又包含若干种方法。
经典的数学分析方法是近百年来电磁学学科发展中一个极为重要的手段。解析法包括建立和求解偏微分方程或积分方程。严格求解偏微分方程的经典方法是分离变量法；严格求解积分方程的方法主要是变换数学法。解析法的优点是：
(1)可将解答表示为已知函数的显式，从而可计算出精确的结果；
(2)可以作为近似解和数值解的检验标准；
(3)在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对结果所起的作用。
用解析法求解电磁场的边值问题可以得到精确的数值结果，并能根据参量的变化推断出解的变化趋势。但是这种方法所能解决的问题不多，满足不了不断增加的工程方面的需要。于是，人们开始致力于研究求解复杂边值问题的近似方法和数值方法。早在1897年，麦克斯韦就尝试用积分方程的数值解来计算矩形金属板间的电容量。在计算机出现以前，应用数值法求解更复杂的边值问题并非易事，但随着高速度大存储量计算机的发展和数值方法应用的日益广泛，过去看来难解的问题现在已能比较容易解决，并能得到足够精确的数值解。
根据问题求解域的不同，数值计算方法可分为时域和频域两大类。时域方法通常直接离散时域麦克斯韦方程，模拟电磁波的传播，随时间迭代计算，适合于宽频带特征瞬态电磁场的数值仿真，如典型的时域有限差分(finite difference time domain， FDTD)法［5，10］、时域有限元(finite element time domain， FETD)法［11，13］、时域谱元法［14，23］和时域积分方法［24，31］。时域方法的一个突出优点是可以给出关于问题空间丰富的时域瞬态信息，能更直观地反映问题的物理现象。频域方法以不包含时间变量或时间变量导数的频域麦克斯韦方程为出发点，求解系统的相位和频率响应。
频域方法可分为高频近似方法和低频数值方法。在高频方法分析时，电磁波和目标之间的相互作用是一种局部现象，只与作用点附近的几何结构、材质参数和入射波性质有关
