一类非稳态分数阶扩散方程最优控制问题的快速算法

来源 :山东大学 | 被引量 : 0次 | 上传用户:sangsang126
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
最优控制问题的有限元逼近是工程设计中的重要课题,而分数阶扩散方程在数学物理领域中的应用也非常广泛,相比于整数阶方程,分数阶扩散方程更能准确恰当的描述反常扩散过程,比如模拟溶质的运动过程,湍流,地下水污染物运移以及古典保守系统的混沌动力学等。因此,对分数阶扩散方程最优控制问题的算法研究有着重要意义。而由于分数阶差分算子的非局部性质,有限元数值方法会产生稠密系数矩阵,针对该线性系统的直接算法往往需要O(N2)的存储需求和O(N3)的计算量。所以,寻找解决该线性系统的快速算法就意义重大。  本文主要研究一类一维非稳态分数阶扩散方程最优控制问题的快速算法,共分六章。  第一章,给出本文要研究的问题,minu∈L2[a,b]J(u)={1/2∫T0(∫ba(y-(y))2dx+∫bau2dx)dt},满足如下方程{(e)y/(e)t-D(k1aD-βx+k2xD-βb)Dy=f+u,a≤x≤b,0≤t≤T,y(a,t)=y(b,t)=0,0≤t≤T,y(x,t=0)=y0,a≤x≤b.其中u(x,t)为控制函数,y(x,t)为状态函数,(y)(x,t)为已知的某一状态函数。说明该问题解的存在唯一性,并给出该最优控制问题的等价形式。  第二章,给出该问题的向后差分-有限元逼近格式,并给出一种相应的先验误差估计,得到如下误差结果:|||u-uh|||L∞([0,T],L2(Ω))≤C(△t+h2),|||y-yh|||L∞([0,T],L2(Ω))≤C(△t+h2),|||p-ph|||L∞([0,T],L2(Ω))≤C(△t+h2).  第三章,给出该问题的中心差分-有限元逼近格式,同样给出对应的先验误差估计,得到如下误差结果:|||u-uh|||L∞([0,T],L2(Ω))≤C(△t2+h2),|||y-yh|||L∞([0,T],L2(Ω))≤C(△t2+h2),|||p-ph|||L∞([0,T],L2(Ω))≤C(△t2+h2).  第四章,根据两种不同的有限元逼近格式,分析得到的线性方程组矩阵特性,根据系数矩阵的Toeplitz性质,在对称情况下,应用PCG算法和Superfast算法将计算量从一般高斯消去法的O(N3)减少到O(Nlog2N)。在非对称情况下,应用预处理共轭梯度平方法(PCGS)将计算量控制在O(Nlog2N)内。  第五章,分别给出对称和不对称两种数值算例,并分别用Superfast算法,PCG(PCGS)算法以及高斯算法解两种有限元逼近格式,比较最后的收敛情况及CPU时间,结果表明,与传统的Guass算法相比,在保持了同样精度的收敛效果下,两种快速算法大大提高了计算效率。  第六章,给出全文的总结。
其他文献
本文主要讨论了广义Cholesky分解,Cholesky-like分解和SR分解的扰动分析的问题。这些矩阵分解是线性数值计算的重要的基础工具。本文主要讨论了这些矩阵分解受到加法或者乘法
本文是一篇关于复数域C上的光滑射影曲面的双有理分类的综述性文章。本文按照文献[10]的思路,采用Mori纲领的方法来重新阐述光滑复射影曲面的双有理分类结果。对于任意的射影
上世纪五、六十年代由三位著名数学家Kolmogorov,Arnold和Moser建立起来的经典KAM理论是哈密顿系统理论发展的里程碑,具有划时代意义,它使人们能够以一种新的方法来研究哈密顿系
随着科学技术的发展,电子设备在日常生活中越来越普及,我们能很容易地用各种电子设备来获取视频。然而,运动载体的拍摄装置在运动过程中会出现由于载体自身机械振动或者自然环境
我国的木构古建筑是历史上珍贵的建筑文化遗产,具有很高的文化价值和艺术价值。但是由于外界环境和自身材料结构的影响,可能导致结构突发性失效,从而严重影响到建筑结构的安全。
时间序列分析是按照处理时间顺序记录的数据的学科。由于我们周围的世界,很多数据是以时间顺序的形式存在和记录的,所以这种方法的应用非常广泛。在理论发展和实践运用的过程中
数字水印技术作为信息隐藏技术的一个重要分支,是信息隐藏领域的重要研究课题,本文主要根据水印透明性,鲁棒性的需求,结合数字图像的特点,运用离散小波变换原理、灰色系统理
图像去噪是图像处理中非常重要的一个领域,原因在于,首先经图像去噪后可以为后续图像处理提供更加准确的图像信息,为图像识别、目标信息分析等打下坚实的基础;其次,探索图像去噪方
粗糙集理论是波兰数学家Pawlak在1982年提出来的,它是一种处理模糊、不确定、不精确数据的数学工具,基于论域上多个二元关系的多粒度粗糙集是粗糙集理论研究的一个新方向.本