广义牛顿型算法求解两类离散非光滑问题

来源 :湖南大学 | 被引量 : 0次 | 上传用户:jiaoyang_204
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
障碍问题和Hamilton-Jacobi-Bellman方程(简称HJB方程)问题产生于机械、工程技术、物理、金融、最优控制等领域.它们的数值解,尤其是大规模问题数值解的研究是工程界和计算数学界一个非常热门课题.近几十年来,取得了许多成果.既然障碍问题和HJB方程是两类典型的非光滑问题,在本文中,我们将研究求解这两类问题的广义牛顿型算法.   在第2章,我们提出了广义牛顿Schwarz迭代法来求解离散的单边障碍问题.该算法的优点如下:(1)算法在每个牛顿迭代步,只采用有限步加性或乘性Schwarz迭代来求解一个低维线性方程组的近似解而不需要求解该线性方程组的精确解,从而可以大大减少计算工作量;(2)算法具有单调收敛性且在适当条件下超线性收敛到问题的解.此外,与其它具有单调收敛的Schwarz算法相比较,该算法的初始迭代很容易选取.   在第3章,我们提出了求解离散HJB方程的广义牛顿法并证明了算法的单调收敛及局部超线性收敛性.该算法的优点是每个牛顿步只求解一个线性方程组从而便于采用线性方程组的快速求解器进行求解.特别地,我们验证了Lions以及Mercier于1980年提出的迭代格式Ⅱ是一类特殊的广义牛顿法,所以该迭代格式具有局部超线性收敛性.进一步,我们研究了求解离散HJB方程的广义牛顿迭代法.该算法在每个牛顿步均采用迭代法来求解线性子问题的一个近似解从而大大地减小了计算工作量.在适当的条件下我们证明了算法具有局部超线性收敛性.数值实验表明了算法是非常有效的.   在第4章,我们研究了求解离散的双边障碍问题的阻尼广义牛顿法.与离散的单边障碍问题的情形相比,离散的双边障碍问题的求解难度更大.当采用古典的有效集策略或增广拉格朗日策略进行求解,常常得不到算法的单调收敛性.本章中,通过选取适当的初始迭代以及在每个迭代步选取一个适当的阻尼因子,我们证明了阻尼广义牛顿法是单调收敛的且具有有限步终止性.而且当问题退化为单边障碍问题时,阻尼广义牛顿法等价于古典的有效集策略算法(或增广拉格朗日策略算法).   在第5章,我们提出了求解离散的双边障碍问题的阻尼广义牛顿迭代法.注意到在前一章中得到的阻尼广义牛顿法在每个牛顿步均需要求解一个低维的线性方程组.当离散问题的规模很大时,精确求解子问题需要很大计算工作量.为了减少计算工作量,本章中我们将在每个牛顿步均采用迭代法来求解线性子问题的一个近似解.在适当条件下,我们证明了算法是超线性收敛的.数值结果表明算法是十分有效的.
其他文献
本文研究了拟复射影空间CQn+p中的全实伪脐子流形Mn的一些性质,采用活动标架场,通过估算子流形第二基本形式模长的平方的Laplacian,得到了一些Simons型积分不等式,从而推广了拟
本文研究Hamilton-Jacobi方程和对流扩散方程的一些新的数值解法,建立这些方法的稳定性和收敛性,并通过大量的数值实验对所提出的算法进行检验.   在第2章基于单调数值通
混沌科学是随着现代科学技术的迅猛发展,尤其是在计算机技术的出现和普遍应用的基础上发展起来的新兴交叉学科。混沌研究已是非线性科学领域的热点问题之一,而混沌的同步与控制
本文主要研究内容为一类退化的Kirchhoff方程(公式省略).其中b>0,2
通过等式((a060)0c0)0=(a0(b0c0)0)0定义了一类特殊的完全正则半群,称其为带相关群并半群.首先应用半群的基本半格刻画了带相关群并半群的结构,然后给出了几类特殊带相关群并半