论文部分内容阅读
最优控制问题的数值近似是工程设计中的重要课题,而分数阶扩散方程在数学物理领域中的应用也非常广泛。相比于整数阶方程,分数阶扩散方程更能准确恰当地描述反常扩散过程。比如模拟溶质的运动过程,湍流,地下水污染物运移以及古典保守系统的混沌动力学等。因此,对分数阶扩散方程最优控制问题的算法研究有着重要的意义。由于分数阶差分算子的非局部性质,有限差分数值方法会产生稠密的系数矩阵,针对该线性系统的直接算法往往需要O(N3)的计算量和O(N2)的存储量。所以,寻找解决该系统的快速算法就意义重大。 本文主要研究基于非稳态空间分数阶扩散方程及点态受限约束的三维最优控制问题的快速差分算法,共分五章。 第一章,给出分数阶最优控制问题的研究背景及现状,并给出所要研究的问题模型:寻找u∈K={u∈L2(0,T;Ω),u≥ u0}使得如下目标泛函达到最小J(u)=1/2∫T0∫Ω(p-(p))2dxdydzdt+1/2∫T0∫Ωu2 dxdydzdt,其中状态方程满足(e)p/(e)t-d(k1a1D2-αx+k2xD2-αb1-e(l1a2D2-βy+l2yD2-βb2)p-s(m1a3D2-γz+ m2zD2-γb3)p=f+u,P|(e)Ω=0, P|t=0=p0(x,y,z),其中u(x,y,z,t)是控制,p(x,y,z,t)是状态,(p)(x,y,z,t)表示对状态的观测函数。 第二章,通过引入伴随状态方程,给出了梯度投影算法并介绍了CN-WSGD差分格式。 第三章,首先给出状态和伴随状态方程的交替方向法,然后给出求解线性方程组的PCG/PCGS算法,最后根据ADI-WSGD格式,分析得到的线性方程组的矩阵特性。根据系数矩阵的Toeplitz性质,在对称情况下,应用PCG算法将计算量从一般高斯消去法的O(N3)减少到O(NlogN)。在非对称情况下,应用PCGS算法将计算量控制在O(NlogN)内。 第四章,分别给出对称和不对称两种数值算例,采用PCG/PCGS方法和GAUSS消去法分别进行求解,并比较最后的收敛情况及CPU时间。结果表明,与传统的GAUSS消去法相比,在保持同样的收敛效果下,快速算法大大提高了计算效率。 第五章,给出全文的总结。