关键词:
电磁逆散射问题
变分玻恩迭代方法
分层范围限制策略
高斯混合模型
自适应阻尼安德森加速方法
摘要:
电磁逆散射问题通过接收器接收的电磁场信号来反演未知目标区域内的散射体位置、形状以及相对介电常数、电导率等信息,也称为电磁成像问题,电磁逆问题,电磁反演问题等。由于其无损、轻量化等特性,该问题为电磁学领域近三四十年来的一个热门研究课题。电磁逆散射问题具有广泛的应用前景的同时,由于其非线性与病态性,求解的准确性与稳定性也面临着极大的挑战。本文聚焦电磁逆散射问题中的介质反演问题,主要在场型反演方法玻恩迭代方法及变分玻恩迭代方法的框架下,为提高其反演的准确性与稳定性,针对几个应用场景,提出了几项改进算法。改进算法表现在对反演算法的目标函数解空间的降维、压缩、以及对迭代解的修正上,使得改进后迭代反演过程更加稳定快速。本文的主要创新内容包括:
1.基于分层范围限制策略的反演方法
待反演目标通常不占据初始划分的整个反演区域。因此,本文提出了分层范围限制策略(Hierarchical Domain Constraint,HDC),该策略将整个反演过程分为若干层,每层包含若干次迭代。HDC在层首迭代反演后通过蒙特卡洛方法及点估计,将整个反演区域的网格分类为散射体网格与背景网格,每层层中迭代反演限制在散射体网格区域内。此举减少了反演中未知数的数目,加快了计算速度,同时因为减少了解的维度,降低了逆散射问题的病态性,从而提高反演精度。HDC的分层策略使得散射体区域的判断是可逆的,避免了因误判将真实待反演散射体所在的网格排除在散射体区域之外而导致的迭代发散问题。同时,HDC引入虚拟背景,在高对比度目标与低对比度目标共存时,增强了对低对比度目标所在网格的识别能力,减少了传统范围限制方法在不同场景下对经验依赖的分类阈值选择。
2.基于高斯混合模型与假设检验的上下限先验融合方法
当已知待反演目标的电磁参数上下限先验知识时,为将此先验知识稳定地融入反演过程中以提升反演精度,本文提出了二次规划-变分玻恩迭代方法-高斯混合模型-分层范围限制(Quadratic Programming-Variational Born Iterative Method-Gaussian Mixture Model-Hierarchical Domain C onstraint,QP-VBIM-GMM-HDC)。首先,QP-VBIM通过二次规划的不等式约束,在反演过程中将待反演参数限制在上下限内,以提升反演结果的精度。然而,直接强制约束待反演参数在其真实上下限内,具有潜在的迭代反演过程发散风险。因此,本文通过基于HDC的GMM-HDC分层自适应地缩放上限解决这一稳定性问题。GMM-HDC同样将反演过程分为多层,在每层反演时,通过GMM聚类及假设检验方法自适应确定反演参数上限及反演范围,Q P-VBIM则在给定的动态参数上限及反演范围内进行反演。由于GMM-HDC对反演参数上限的动态调整,有效规避了因固定上限引发的难以察觉的迭代发散风险。此外,基于GMM的聚类方法能够更精准地识别未知待反演目标所在网格,且无需人为设定阈值。借助假设检验方法,即便在缺乏待反演目标电磁参数上下限先验知识的情况下,GMM-HDC所给出的上限依然能够依概率收敛至真实上限。
3.基于自适应阻尼安德森加速方法的反演
为在无先验知识及外源数据的情况下提升逆散射方法的收敛速度及稳定性,本文基于计算数学中加速不动点迭代的安德森加速方法,提出了应用于逆散射问题的自适应阻尼安德森加速(Adaptive Damping Anderson Acceleration,ADAA)方法,并将其与场型迭代方法变分玻恩迭代方法(Born Iterative Method,BIM)及玻恩迭代方法(Born Iterative Method,BIM)结合验证,形成 VBIM-ADAA 与 BIM-ADAA方法。由于逆散射问题的病态性,数据方程的求解需在正则化条件下进行,这导致所求迭代解无法精确满足数据方程。因此,ADAA在每次迭代后通过计算前几次迭代历史解及当前解在当前数据方程下的最小残差线性组合系数,并以此系数组合历史解以获取更精确满足当前数据方程的解。从而使得VBIM-ADAA及BIM-ADAA实现了更高的迭代稳定性与更快的收敛速度,且无需外源信息及先验知识。同时,本文借助QR分解方法进一步说明了 ADAA在不同正则化参数下性能变化的原因。ADAA可在当前迭代发散时回溯至前几次迭代的解以保证迭代稳定性,这使得ADAA允许VBIM及BIM在更大范围内选择正则化参数。理论分析及数值实验表明,ADAA方法具有极小的计算负担。并且作为一种通用的数学方法,ADAA可与其他逆散射方法结合使用。
本文针对电磁逆散射问题中介质反演问题的不同场景,分析了现有方法的局限性,并提出了几种提升反演精度与稳定性的改进策略。本文将这些策略嵌入到传统场型