基于统计计算方法的反散射问题数值解研究

来源 :电子科技大学 | 被引量 : 0次 | 上传用户:hunterxjtu
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着计算机技术的革新和工程应用的实际需要,人们对于反散射问题的研究不再仅仅限于理论上的分析,更多的时候人们希望能够有效的数值模拟反散射问题的解。然而,与正散射问题相比,反散射问题往往是不适定的,即问题解的存在性、唯一性和稳定性总是被破坏。另外,在数值求解算法上,反散射问题与正散射问题也有很大不同。本文以统计计算方法作为基本反演手段,对声波和弹性波中几类反散射问题的求解进行了研究。研究的主要内容和成果归纳如下:1.针对有限孔径的声波反散射问题,提出了一种将扩展采样法(Extended Sampling Method,ESM)与Bayes方法相结合的新技术来重建声软障碍物的位置和形状。通过Bayes公式将反问题重新表述为一个统计模型进行求解,并证明了转化后模型后验概率测度的适定性。在使用Markov chain Monte Carlo(MCMC)方法对后验分布进行采样的过程中,事先知道目标位置对于加速MCMC方法的收敛速度具有重要的作用。为此,我们事先采用了ESM对目标的位置进行了反演。ESM-Bayes方法分为两个步骤:i)ESM获取目标位置,ii)Bayes方法重构边界。这两个步骤都基于相同的物理模型,并使用相同的测量数据。2.提出了一种结合ESM与集合Kalman滤波(En KF)的两步方法来求解弹性波障碍反散射问题。作为算法的第一步,我们首先使用ESM来反演障碍物的近似位置。为了能使算法适用于有限孔径的数据,我们在已有文献的基础上对ESM进行了改进。第二步,利用En KF对ESM获得的位置进行优化,并利用相同的数据集重建未知障碍物的形状。ESM得到的近似位置可用于En KF初始粒子集的构建。数值实验表明,该方法继承了两种方法各自的优点,且能够利用有限孔径数据对目标障碍物进行有效的恢复。3.鉴于前面两章的研究思路,我们将直接采样法(Direct Sampling Method,DSM)与Bayes方法相结合,提出了一种求解声波反源问题的高质量反演方法。第一步,采用DSM来反演未知源的数目和位置。第二步,对于源的细节信息,采用Bayes方法进行反演。在Bayes统计中,模型中的参数被视为随机变量。第一步中获得的源的数目和位置被编码在Bayes方法的先验分布中。然后,利用Bayes公式,得到了未知参数的后验分布,证明了该方法的适定性,并采用preconditioned Crank Nicolson Metropolis Hastings(p CN-MH)算法对后验密度函数进行求解。数值算例表明该方法具有很强的稳健性,且能够在部分观测数据下对未知源的数目、位置及强度进行反演。4.考虑了一个具有重要应用背景的介质反散射问题,即从均匀背景介质中识别可穿透层状散射体界面的几何形状。我们将集合Kalman方法作为反演求解器。为了能够灵活的处理散射体界面边界的拓扑变化,我们在En KF算法中引入了水平集技术。在整个反演过程中,不需要计算正演算子及其伴随算子的Fréchet导数。这与传统的梯度型算法相比明显是一个优势。通过数值实验说明了该算法的有效性和灵活性。该算法不仅可以处理一个散射体,而且还可以对多个散射体及环状散射体进行恢复。5.考虑了一个利用散射场数据来恢复背景弹性介质中腔体形状的反散射问题。由于与该反问题相对应的正散射问题计算成本较高,若采用MCMC算法进行求解,势必会导致Bayes公式中似然势函数的计算量较大,从而使整个反演过程变的十分耗时。为了提高反演效率,我们采用Kriging代理模型对Bayes公式中的似然势函数进行替代。数值实验表明,Kriging代理模型不仅可以提供合理的计算精度,而且能够极大的缩短MCMC算法的反演时间。
其他文献
相依序列的收敛性质是近代概率极限理论的研究热点之一,它在概率统计、金融与保险、可靠性理论、复杂系统以及计量经济学等领域有着十分广泛的应用。本文致力于研究包括NSD序列、AANA序列、ANA序列以及END序列在内的4类相依序列的极限收敛性质,利用相依序列的矩不等式和一些概率不等式,进一步研究了相依序列的完全收敛性、完全矩收敛性、强大数律以及含随机系数线性过程的强收敛定理,并获得了一些新的结果。本论文
由于能够很好地解释物理、化学、生物等领域的某些重要现象和规律,偏微分方程的理论与应用已成为重要的数学研究方向.这些理论包括方程解的存在性、唯一性、有界性、有限时刻爆破以及大时间渐近行为等.特别地,趋化性是细胞或生物体对化学刺激产生的定向运动,在胚胎发育、伤口愈合和肿瘤入侵等各种生物过程中发挥重要作用.本文研究如下三类趋化模型:(ⅰ)具有一般旋转灵敏性的吸引-排斥趋化模型(?)(ⅱ)具有一般旋转灵敏
现代微电子和集成电路技术飞速发展的同时也面临着集成度过高引起的局部热功耗过高等不可避免的问题。兴起的自旋电子学利用了电子的自旋自由度,成为了构造高密度、低功耗、非易失的新一代自旋电子学器件的基础,也推动了电子信息技术的全新变革。对自旋流的激发、传输和探测成为当今自旋电子学的研究主题。本文以液相外延法制备的超低损耗亚铁磁性绝缘体钇铁石榴石(YIG)薄膜为基础,研究了YIG/金属自旋异质结中的自旋泵浦
合作行为广泛存在于自然界和社会系统中。复杂网络上的演化博弈理论为解释自私群体中合作行为的涌现提供了有效的理论工具。复杂网络刻画了种群的拓扑结构,博弈模型描述了智能体之间的交互模式。网络结构和个体行为呈现多种多样的形式,这种多样性为选择提供了基础,因此对合作涌现至关重要。零行列式策略作为重复博弈新理论,可以单方面保证双方长期收益满足线性关系。剥削策略是零行列式策略的一个重要子集,保证了自身的收益不低
作为应用于电子通信领域的基础关键性材料,微波介质陶瓷具有广阔的应用前景。发展至今,陶瓷材料在微波频段介电性能的改善问题始终是研究的重点以及难点。然而目前该领域始终存在以“经验”作为指导思想展开实验,尽管存在可借鉴性,却无法从根本上利用理论指导实验,在研发新体系或对材料进行改性研究时没有“指挥棒”的作用发挥,因此寻找一种行之有效的理论成为研究学者的迫切需求。众所周知介质陶瓷的独特性体现在微波特性随晶
金属薄膜电爆炸驱动飞片技术是一种重要的动高压加载技术,以该技术为核心的冲击片雷管已应用于多个武器型号中。但由于现有的金属电爆炸薄膜电能转化为飞片动能效率较低,使得该技术难以适应未来武器系统小型化的发展需求,亟需开展能量转换效率高、易集成、可靠、安全的新型电爆炸薄膜设计、制备及应用技术研究。为此,本论文从电爆炸金属薄膜材料设计出发,将Al/Ni含能多层膜作为电爆炸材料,系统地开展了Al/Ni含能多层
半个多世纪以来,描述生物趋化现象的偏微分方程越来越受生物学家和数学家们的关注。考虑到趋化实验的条件设置和现实生活中的趋化背景,本文主要研究了趋化-流体耦合模型在带有边界的无界和有界区域上的初边值问题,具体研究内容如下:1.研究了三维带边无界区域上的趋化-Navier-Stokes耦合方程在Neumann-Neumann-Dirichlet边界条件下的初边值问题。首先利用各向异性的Lp插值不等式和常
电磁积分方程方法因其具有较高的计算精度和较少的未知量,已被广泛应用于电磁仿真设计领域中。使用矩量法离散积分方程会产生一稠密矩阵,为此诸多快速算法被提出以进一步提升了积分方程方法的计算能力。绝大多数快速算法通常用于加速积分方程离散形成的线性矩阵与右端向量之间的乘积。所以这些快速算法常结合迭代方法,求解积分方程离散形成的线性矩阵方程。然而,迭代方法仍然有严重的弊端。首先是对于复杂问题,迭代方法的收敛很
随着当代武器装备和电子器件的迅速增长,例如大功率真空电子器件、军舰和装甲导弹等系统,在微波器件设计、卫星通信及雷达等领域都各自发挥着重要的作用。实际上,这些设备本身表面可能设置有各类天线、传感器等细小装置,同时组成的介质材料往往是各不相同的,使得整个设备的物理特性变得非常复杂,因此具有几何及材料的多尺度特征。此外,在现代战场中,为了发挥不同的战场功效,辐射源的数量变得越来越大,而这导致电磁环境日趋
按疾病诊断相关分组付费(diagnosis-related groups, DRG)和按病种分值付费(diagnosis-intervention packet, DIP)两大医保费用支付方式已开始在我国多地试点实施,住院患者医保费用支付方式改革对医保基金管理和医疗机构运营将产生重要影响。药品费用作为医保费用的重要组成部分也成为医保管理部门及医疗机构的关注重点,但DRG和DIP对于药品费用的影响程