最优控制和图像处理问题中的区域分解算法与并行计算

来源 :华东师范大学 | 被引量 : 2次 | 上传用户:q1q3q5
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着科技的进步,越来越多的领域提出了对大规模并行计算的要求.比如在航空航天,地球物理,药物设计,纳米材料,卫星遥感等等复杂系统在各种现实条件下的形态和行为做出精确模拟和预测.这些需求的出现和增长,催生了计算机硬件的快速发展.而如何使用超级计算机来解决现实中的大规模计算问题成为当前日益关注的问题.在数值计算领域,适用于大规模计算的并行算法的研究就变得非常必要.一个好的并行算法,必须具备稳健的性质,能够实现长时间大规模的模拟计算.同时,该算法需要具有很好的加速比和加速效率,具备良好的扩展性.而区域分解算法,恰恰具备以上的特性.在本文中,我们主要研究了在图像处理中的非线性变分问题和最优控制领域的区域分解算法和并行实现.这两类问题在图像处理领域和实际工程计算问题中有大量应用.因此如何高效并行求解,是迫切需要解决的问题.对这两类问题,我们主要做了如下的工作.首先是关于偏微分方程约束的最优控制问题的区域分解算法,我们考虑了如下几类算法:1.1椭圆方程约束的最优控制的施瓦兹交替梯度投影算法.实现重叠区域算法的可控性,即在内迭代有限步停止迭代.我们给出了施瓦兹交替梯度投影算法.该算法包括两层迭代,内迭代采用施瓦兹交替法来求解状态变量的偏微分方程,外迭代采用梯度投影方法求解控制变量.我们得到了算法的收敛性和收敛速度,同时给出了合适的内迭代次数的上界.采用给出的内迭代次数,有效地节省了计算量,同时保证了算法的整体收敛性.1.2对椭圆方程约束的最优控制时分布受限非重叠区域分解算法.我们考虑了Robin型的非重叠区域分解算法,同时把求解等价的变分不等式系统转化到求解子区域上的变分不等式系统.更进一步,我们可以通过选取合适的松弛参数来获取最优的收敛速度.这是在理论上首次得到区域分解算法求解非线性系统的收敛速度.1.3对椭圆方程约束的最优控制时积分受限非重叠区域分解算法.不同的限制集对区域分解算法有较大的影响,有些可以直接并行,有些不可以.1.2的工作是建立在满足相容性条件的控制受限的约束集上的.我们提出了更一般约束集下面的Robin型区域分解算法,同时得到了算法的收敛性.针对图像处理中的非线性变分问题,主要是基于ROF模型的相关能量极小化问题,我们有如下的工作.2.1针对图像处理中的Rudin, Osher和Fatemi (ROF)模型,我们研究了其对偶模型的区域分解算法.首先我们采用稳定的单位分解来分解不等式约束,在此基础上构造了连续子空间和并行子空间校正算法.更进一步,我们分析了算法的收敛性,得到了O(n-1)的收敛速度,其中n为迭代次数.据我们所知,这个结果是区域分解在图像处理领域应用中首次得到的.2.2非局部总变差模型利用了更多的平均信息,所以也就需要比ROF模型更多的计算量,有效地实现并行计算非常必要.基于Xu, Tai和(?)Wang [96]的工作,我们采用连续子空间校正算法将问题转化成子区域上的问题.对每个子问题,我们采用Bregmanized算子分裂算法来求解.数值试验显示,该算法成功实现了图像去噪,图像去模糊以及图像修复.最后的并行算例展示了算法良好的加速比和加速效率.
其他文献
近几年来,对于非单的C*-代数的分类研究取得了许多重要进展。H.Lin和H.Su对AT-代数进行了分类。这一工作的重要性在于,AT-代数通常可以用AT-代数通过AF-代数的扩张所得到。在本文中,我们将指出相反的结果是不成立的,存在着大量的例子说明A??-代数通过AF-代数的扩张不是AT-代数。并利用K-理论给出它们是AT-代数的一个必要条件。我们给出AT-代数一个等价定义并讨论一类更广泛的C*-代
应变是研究材料物性和结构的重要手段,同时也是调制和改良材料物性的重要方法,本文基于第一性原理计算方法,研究了稀磁半导体(Ga,Mn)As体系中的各种Mn掺杂的形成能和扩散势垒与应变之间的关系,我们得出Mn掺杂的形成能和扩散势垒是否存在线性变化,及应变对氧化物(LaAlO3/SrTiO3)界面处氧空位的形成能和电子结构的影响。本文的主要工作有以下三部分。一、应变对半导体的影响往往是通过连续弹性模型来
小鲵科为亚洲特有的有尾两栖动物。中国西部小鲵为小鲵科动物中的重要类群,阐明其各属及种间的系统发育关系,澄清其属、种分类问题,并探讨其与地理分布格局的关系,对整个小鲵科的系统演化与分布格局关系的研究具有关键性的意义。其次,由于中国西部小鲵的主要分布于青藏高原东缘和邻近周边地区,其独特的地理分布格局可为探讨青藏高原的隆升、地史变迁与相关生物类群演化历史间的相互关系提供有价值的线索。此外,中国西部小鲵中
塞加羚羊(Saiga tatarica)因拥有一对通透竖直的双角,一直是我国传统中医药中著名的物种。它曾广泛地分布在中亚半干旱草原上,数量超过100万头。但是,在过去短短的十几年时间中,由于没有节制的猎杀,这种动物的数量下降了95%以上,并且破坏了种群的自然性比。而在我国,由于捕杀,早在20世纪60年代之后就没有了该种动物的野外纪录。 考虑到物种多样性保护和国药发展的需求,从1987年开始,
学位
在本文中,我们主要研究复数域及素特征代数闭域上有限W-代数和有限W-超代数结构及其表示理论的相关问题.研究成果主要有:1.有限W-代数方面的工作第一部分是对复数域上B2型有限W-代数的研究.我们通过计算得到了其具体生成元和关系式,从而完全刻画了此代数的结构.第二部分是对素特征域上限制李代数sl2的约化包络代数的中心结构的研究.本部分利用素特征域上有限W-代数的相关理论,得到了限制李代数sl2关于p
健美操是一项集体操、音乐、舞蹈为一体,结合有氧运动和无氧运动,充分展示健、力、美的体育项目。健美操可分为健身类、表演类、竞技类3种类型。健身类健美操是以锻炼身体为目的,其动作简单、音乐缓慢。表演类健美操是专门为表演设计的,其动作较快,音乐可快可慢。竞技健美操属于健美体育的范畴,是一种要求运动员2分钟内在音乐伴奏下完成高强度、连续、成套动作的体育项目,其对运动员身体素质、技术能力以及艺术表现力
期刊
文章的主要研究对象是量子拟shuffle代数.本博士论文分为三个部分:第一部分给出多参数量子群的量子拟对称代数实现,这是一种公理化的构造方式.第二部分考虑一类特殊的量子拟shuffle代数,q-拟对称函数代数.特别地,我们研究了奇拟对称函数代数的组合性质.第三部分详尽刻画了一个tame表示型阶化Hopf代数(?)-1(2)的不可分解模结构以及其张量积分解结构,从而确定其表示环(Green环),以及
学位
基于超冷中性原子和单个囚禁离子的光频标频率不确定度远远低于当前最好的铯喷泉原子钟,是下一代时间频率基准的有力竞争者。光钟不仅可以用来测量基本物理常数的变化,也是全球定位系统和高速通信网络中的重要技术。目前最好的光钟是铝离子光钟,其不确定度已经达到了10-18量级,但是离子光钟的稳定度受到量子投影噪声的限制。而中性原子光晶格钟能同时探测大量囚禁在光晶格中的原子,其稳定度远远高于离子光钟。镱(Yb)原
棉藓属(Plagiothecium Bruch&Schimp.)是棉藓科的核心属,为温带山地侧蒴藓类中重要的组成部分。尽管目前基于数个形态特征性状的综合可以将其与近缘属区分,但由于该属缺乏独征,对于其是否为一个自然类群仍有争议。东亚地区有大量棉藓属及其近缘属物种分布,然而由于小范围修订造成的片面理解,属内存在诸多物种界定的模糊、不一致甚至相互矛盾之处,也遗留了一定数量仅有模式标本记载的“幽灵种”。
复流形是复几何所考虑的基本对象。我们设M是一紧复流形,E为其上的可微分向量丛。借助于它们上面的Hermitian度量我们可以定义两个模空间:一个联系着代数几何;一个联系着微分几何。Kobayashi-Hitchin correspondence说:这两大数学的分支有一座桥梁,由上述两个模空间的一一对应关系而产生。 上面的这一设想在大约二十年前被证明是对的,在基流形M是Kahler流形的情形:
学位