几类反问题的有效数值解法

来源 :大连理工大学 | 被引量 : 0次 | 上传用户:dwwn123456
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
反问题在许多领域有着非常广泛的重要应用,例如医学成像、现代工业、气象预报等等.也因此,数值求解反问题引起了广泛关注.在本论文中,我们考虑以下几类反问题:带箱型状态约束的椭圆型最优控制问题;带积分型状态梯度约束和箱型控制约束的椭圆型最优控制问题;图像恢复中的病态的线性反问题.受到一阶算法成功用于求解大规模有限维优化问题的启发,我们试图研究如何利用一些有效的一阶算法并结合这些问题的结构来更有效地求解它们,取得的主要成果如下:1.带箱型状态约束的椭圆型最优控制问题.对于逐点状态约束,相应的拉格朗日乘子是Borel测度函数,这使得最优性条件中的互补性条件无法表示成逐点的形式.Lavrentiev正则化被用来解决这一困难,然后我们用分片线性有限元全离散方法来离散正则化的问题.我们关于正则化和离散化带来的误差给出了误差分析,由此可以看到全离散的误差阶并不比变分离散的误差阶差.为了数值求解离散问题,我们提出了一种异构交替方向乘子法(hADMM)以及一种两阶段策略,它们分别用来得到中等精度和高精度的解.2.带积分型状态梯度约束的椭圆型最优控制问题.我们首先证明了该问题的最优性条件并提出了一种有限元对偶的不精确majorized加速块坐标下降(FE-dABCD)算法用于求解该问题.具体地,我们用标准的分片线性有限元方法来离散该问题并给出了离散问题的对偶问题,该对偶问题是一个多块无约束凸优化问题.然后一种不精确majorized ABCD算法被用于求解该对偶问题.在求解子问题时,我们使用了一些有效的策略,这使得FE-dABCD算法享有O(1/k2)的迭代复杂度.此外,基于ABCD算法的网格独立性,我们还提出了一种多层FE-dABCD(mFE-dABCD)算法.3.图像恢复中的线性反问题.我们利用问题中隐藏的两个结构来改进快速迭代收缩阈值算法(FISTA)并称所提出的改进算法为结构化FISTA(sFISTA).我们探索的第一个结构是系数矩阵的Kronecker积近似,这一结构不仅可以引入一层并行,还使得在后续的优化进程中可以使用计算密集的矩阵-矩阵乘法.第二个结构则是在Kronecker积近似中出现的所有矩阵均为具有特殊结构的矩阵(Toeplitz,Hankel等),这使得我们在每次迭代时可以探索它们的快速矩阵-向量乘法.此外,我们证明了 Kronecker积近似带来的误差是可控的并且sFISTA可以达到与FISTA相同的精度水平.
其他文献
毫米波无线通信技术是未来无线通信技术发展的重要领域之一。传统的sub-6GHz通信频段和无线通信技术,已经无法满足人类快速增长的数据需求,而毫米波波段有着更为丰富的带宽资源,以及更大的信道容量,使得更加快速、可靠的通信成为可能,因此得到越来越多的人的关注,同时也有着极大的发展潜力。毫米波通常指频率处于30GHz到300GHz之间的电磁波,虽然有着丰富的带宽资源,但缺点是需要承受严重的路径损耗,导致
锚泊定位系统作为浮式结构物最常用的被动式定位方式,其工作原理是利用锚泊线把水面浮式结构物与固定在海底的锚固基础连接起来,约束浮式结构物的运动响应在允许的范围内,同时还需要保证锚泊系统自身的强度设计使其不产生过大的锚泊张力。按照锚泊线的线形形状可以分为两类:悬链式和张紧式。悬链式锚泊系统是浮式结构物传统且最常用的定位方式,通过自身重力的变化为浮式结构物提供回复力,但是对于深水或超深水浮式结构物,随着
第五代(5G)移动通信系统,由于移动数据流量的爆炸式增长和设备的超量连接,将在超大范围、海量容量、海量连接、超低延迟等方面面临新的要求和挑战。满足这些需求的主要限制来自可用频率资源的有限性。利用认知无线电(CR)技术提高现有频段的利用率,是解决频谱稀缺问题的一个有前途的解决方案。频谱感知作为认知无线电网络的主要环节,允许次级用户利用授权的主用户的空闲频谱进行机会通信。在实际应用中,低信噪比、无线信
车载自组织网络(VANET)作为智能交通系统的关键技术之一,可以为智能交通提供理想而有用的通信网络技术。在确保交通安全,效率,商业和娱乐方面,它将发挥重要作用。由于VANET广泛的应用前景,研究VANET的网络设计和实现对于交通信息化的发展具有深远的意义。网络可靠性和连通性是VANET应用的基础和前提,它包括物理连通性和逻辑连通性,涉及物理链路连通性,MAC协议和路由协议的有效性等。本文主要研究V
室内土工试验是探究岩土材料力学特性的必要手段,其中三轴试验是测定岩土材料应力应变性质、本构关系模型与强度参数的最重要、最常用的试验。但三轴试验系统的加载跟踪精度易受土样非线性力学特性及系统中非线性摩擦等因素的影响,同时,常规三轴试验中试样变形的测量一般采用外置式传感器,只能测量试样整体的变形,无法获得土样的局部变形特征。针对三轴试验系统存在问题,本文实现了三轴试验系统轴向加载的非线性补偿控制和试样
对于成分复杂且宏观离散的岩土材料,不能简单借鉴单相连续介质力学中的应力定义,因而Terzaghi提出了饱和土有效应力原理,被誉为现代土力学的基石。但拓展至非饱和土时,碍于复杂的相间作用,有效应力原理面临诸多挑战,迄今仍是土力学研究的热点问题。为探寻非饱和土中具有明确物理本质的应力体系,首先从合理建模入手,将非饱和土细分为土颗粒、胶结物质、固态孔隙水、液态孔隙水、收缩膜和孔隙气,提出了非饱和土6相模
本文主要研究了几类环上的图的自同构的问题。对于一个图G,如果点集(记作V(G))上的一个双射σ保持点与点的连接关系,那么σ被称作是图G的一个图自同构。如果把映射的合成看作群的乘法,那么图G的所有自同构在这个乘法下可以形成一个群,这个群被称为图G的自同构群,记作Aut(G)。给定一个环S,这个环的零因子图Γz(S)是一个有向图,图的点集是环中所有的非零零因子,存在一条边使得不同的点A和点B相连当且仅
在本文中,我们引入变分迭代法、同伦扰动法和拉普拉斯变换来求解各种类型的偏微分方程,并且通过这几种解析方法的比较研究,我们发现变分迭代法就是在校正函数中补充迭代步骤,从而提供对解的连续逼近;同伦扰动法就是克服方程中由非线性现象引起的非线性项所带来的困难,从而使逼近解快速收敛到精确解。本文以一些偏微分方程为例,如Klein-Gorden方程,Boussinesq方程和电报方程,来说明这些方法的有效性。
本文筛选出公众、科学、政府和法律为湿地保护制度的四个要素,以此构建了一个使不同国家湿地保护制度研究能在统一语境下进行比较研究的理论框架。通过研究美国湿地保护政策百年发展的特征、逻辑、执行及当前挑战,建立了“竹节模式”经验性模型。该模型揭示了美国湿地保护制度政策变迁的动因及各要素间相互作用从而保护湿地的机理。最后讨论了我国湿地保护制度变迁及当前挑战,提出了我国湿地保护制度发展的建议。第一章是绪论。规
在大型强子对撞机(LHC)能量(核子对质心能量在TeV量级)下,重离子碰撞产生的π介子多重数可数以千计。在这种情况下,粒子发射源内的全同π介子很可能发生玻色-爱因斯坦凝聚,从而对全同π介子的Hanbury Brown-Twiss(HBT)关联造成抑制。最近,欧洲核研究机构(CERN)的ALICE合作组在对LHC上运行的(?)=2.76 TeV Pb-Pb碰撞的实验数据分析中,发现实验测量到的全同π