相位恢复在多个不同领域均被提出，如量子力学、光学成像等.相位恢复即具有多种应用背景，亦具有丰富的数学内涵，因而近期该问题吸引了多个不同领域专家的关注，如计算数学、数据科学、最优化、代数几何等.本文将主要介绍相位恢复中的理论基础问题，特别是最少观测次数问题，并介绍求解相位恢复的模型性能，以及求解算法等.本文也介绍了一些当前相位恢复中研究的热点方向.

原子间相互作用建模是分子动力学模拟的核心问题之一.基于第一性原理的建模准而不快，经验势模型快而不准，因此人们长期面临精度和效率只得其一的两难困境.基于机器学习的原子间相互作用建模在达到第一性原理精度的同时，计算开销大大降低，因而有希望解决这一两难困境.本文将介绍构造基于机器学习的原子间相互作用模型的一般框架，归纳近年来的主要建模工作，并探讨这些工作的优势和劣势.

玻尔兹曼方程作为空气动理学中最基本的方程之一,是连接微观牛顿力学和宏观连续介质力学的重要桥梁.该方程描述了一个由大量粒子组成的复杂系统的非平衡态时间演化:除了基本的输运项,其最重要的特性是粒子间的相互碰撞由一个高维,非局部且非线性的积分算子来描述,从而给玻尔兹曼方程的数值求解带来非常大的挑战.在过去的二十年间,基于傅里叶级数的谱方法成为了数值求解玻尔兹曼方程的一种很受欢迎且有效的确定性算法.这主要归功于谱方法的高精度及它可以被快速傅里叶变换加速的特质.本文将回顾玻尔兹曼方程的傅里叶谱方法,具体包括方法的导出,稳定性和收敛性分析,快速算法,以及在一大类基于碰撞的空气动理学方程中的推广.

许多物理现象可以在数学上描述为受曲率驱动的自由界面运动,例如薄膜和泡沫的演变、晶体生长,等等.这些薄膜和界面的运动常依赖于其表面曲率,从而可以用相应的曲率流来描述,其相关自由界面问题的数值计算和误差分析一直是计算数学领域中的难点.参数化有限元法是曲率流的一类有效计算方法,已经能够成功模拟一些曲面在几类基本的曲率流下的演化过程.本文重点讨论曲率流的参数化有限元逼近,它的产生、发展和当前的一些挑战.

在电子结构计算领域, Kohn-Sham方程是最为广泛使用的数学模型之一. 然而, 由于现有的交换关联能近似仍存在缺陷, Kohn-Sham方程无法较好地描述强关联多电子体系. 近年来, 有学者从密度泛函理论的强相关极限出发, 提出了严格关联电子能量的优化模型. 该模型有望弥补Kohn-Sham方程的缺陷, 从而拓宽密度泛函理论的应用面. 由于在该模型中存在维数灾难, 近年来, 它的一些低维转化模型陆续被提出. 在本文中, 我们将介绍严格关联电子能量的优化模型、它的研究重点以及现有的一些低维转化模型. 我们也将介绍这些转化模型的数值求解方法, 并探讨未来的研究方向.

本文讨论双曲型守恒律方程的熵稳定格式.对于给定的熵对，格式所满足的熵条件中的数值熵通量是不唯一的.Tadmor的充分条件可以唯一地确定标量方程的熵守恒通量，但不能唯一确定方程组的熵守恒通量，却可以给出方程组的空间一阶精度的熵守恒格式.也讨论了在熵守恒通量上添加数值粘性得到的显式熵稳定格式需要满足的条件及常见的时间离散对熵守恒和熵稳定的影响.

很多交叉科学的实际问题在数学上都可以被归为求解具有多个变量的非线性函数或泛函的极小值问题, 如何有效地寻找其能量景观的全局极小和如何找到不同极小之间的关系是计算数学领域两个长久以来尚未解决的重要科学问题. 本文着重介绍近年来提出的“解景观”概念和方法. 我们将回顾解景观的概念、构建解景观的鞍点动力学方法、以及解景观在液晶和准晶方面的应用.

本文研究了数值求解非自治随机微分方程的正则Euler-Maruyama分裂（CEMS）方法，该方程的漂移项系数带有刚性且允许超线性增长，扩散项系数满足全局Lipschitz条件.首先，证明了CEMS方法的强收敛性及收敛速度.其次，证明了在适当条件下CEMS方法是均方稳定的.进一步，利用离散半鞅收敛定理，研究了CEMS方法的几乎必然指数稳定性.结果表明，CEMS方法在漂移系数的刚性部分满足单边Lipschitz条件下可保持几乎必然指数稳定性.最后通过数值实验，检验了CEMS方法的有效性并证实了我们的理论结果.

针对低秩稀疏矩阵恢复问题的一个非凸优化模型，本文提出了一种快速非单调交替极小化方法.主要思想是对低秩矩阵部分采用交替极小化方法，对稀疏矩阵部分采用非单调线搜索技术来分别进行迭代更新.非单调线搜索技术是将单步下降放宽为多步下降，从而提高了计算效率.文中还给出了新算法的收敛性分析.最后，通过数值实验的比较表明，矩阵恢复的非单调交替极小化方法比原单调类方法更有效.

平移对称幂法（SS-HOPM）在求解源自玻色-爱因斯坦凝聚态的非线性特征值问题时，不仅具有较高的计算效率，而且具有点列收敛性，但其收敛率尚未得到有效估计.本文通过将多项式Kurdyka-Łojasiewicz（K-Ł）指数界的相关结果应用到所涉及优化问题的Lagrange函数上，得到了平移对称幂法的次线性收敛率估计，从理论上解释了平移对称幂法的计算效率.

本文采用Modulus-based变换将非线性互补问题转化为非光滑方程组，并将一种多步自适应Levenberg-Marquardt方法推广应用于求解所得的非光滑方程组，从而得到原问题的解.在适当条件下，本文证明了算法的全局收敛性.与一种已有的参数自适应Levenberg-Marquardt方法（PSA-LMM）相比较，数值实验结果表明了本文所提出的算法具有更好的效率.

图像融合通常是指从多源信道采集同一目标图像,将互补的多焦点、多模态、多时相和/或多视点图像集成在一起,形成新图像的过程.在本文中,我们采用基于Huber正则化的红外与可见光图像的融合模型.该模型通过约束融合图像与红外图像相似的像素强度保持热辐射信息,以及约束融合图像与可见光图像相似的灰度梯度和像素强度保持图像的边缘和纹理等外观信息,同时能够改善图像灰度梯度相对较小区域的阶梯效应.为了最小化这种变分模型,我们结合增广拉格朗日方法(ALM)和量身定做有限点方法(TFPM)的思想设计数值算法,并给出了算法的收敛性分析.最后,我们将所提模型和算法与其他七种图像融合方法进行定性和定量的比较,分析了本文所提模型的特点和所提数值算法的有效性.

本文构建一类双参数拟Toeplitz分裂（TQTS）迭代方法求解变系数非定常空间分数阶扩散方程.TQTS迭代法是基于QTS迭代法引入双参技术建立而成，通过选取适当的参数使迭代矩阵谱半径变得更小，从而有效提升收敛的速度.然后对TQTS迭代法进行收敛性分析，获得相应的收敛区域，并对迭代法中涉及的参数进行讨论，获得使迭代矩阵谱半径上界达到最小的最优参数的表达式.最后通过数值仿真实验验证TQTS迭代法的有效性，实验结果表明TQTS迭代法改进效果十分突出，在迭代时间和步数上均有明显的减小.

对于求解大型稀疏连续Sylvester方程，Bai提出了非常有效的Hermitian和反Hermitian分裂（HSS）迭代法.为了进一步提高求解这类方程的效率，本文建立一种广义正定和反Hermitian分裂（GPSS）迭代法，并且提出不精确GPSS（IGPSS）迭代法从而可以降低计算成本.对GPSS迭代法及其不精确变体的收敛性作了详细分析.另外，建立一种超松弛加速GPSS（AGPSS）方法并且讨论了收敛性.数值结果表明了方法的高效性和鲁棒性.

本文针对带有随机杨氏模量和荷载的平面线弹性问题，提出了一类随机弱Galerkin有限元方法.先利用Karhunen-Loève展开把随机项参数化，将方程转化为一个确定性问题；再采用弱Galerkin有限元法和$k$-/$p$-型方法分别离散空间区域和随机场.在弱Galerkin离散中，用分片$s（s\geqslant 1$）和$s+1$次多项式逼近单元内部的应力和位移，用分片$s$次多项式逼近位移在单元边界上的迹.证明了该方法关于空间网格尺度最优且与Lamé常数$\lambda$一致无关的误差估计.最后通过数值算例验证了理论结果.

本文考虑了分布阶时间分数阶扩散波动方程，其中时间分数阶导数是在Caputo意义上定义的，其阶次$\alpha，\beta$分别属于（0，1）和（1，2）.文中提出了在计算上行之有效的数值方法来模拟分布阶时间分数阶扩散波动方程.在时间上，通过中点求积公式把分布阶项转换为多项的时间分数阶导数项，并且利用$L1$和$L2$公式来近似Caputo分数阶导数；空间上使用Galerkin有限元方法进行离散.给出了基于$H^1$范数的有限元解的稳定性和误差估计的详细证明，最后的数值算例结果说明了理论分析的正确性以及有效性.

针对高阶马尔科夫链的极限概率分布问题，本文提出了两种带动量项的高阶幂法，并在一定条件下建立了所提算法的收敛性定理.数值实验结果表明动量项能够有效改善原幂法的计算效率.

信赖域方法是求解非线性方程组的一种重要方法.本文研究了求解非线性方程组的信赖域半径趋于零的信赖域算法在Jacobi矩阵Hölderian连续条件下的全局收敛性质，以及其在Hölderian局部误差界和Jacobi矩阵Hölderian连续条件下的收敛速度.

本文研究跳适应向后Euler方法求解跳扩散随机微分方程在非全局Lipschitz条件下的强收敛性.通过克服方程非全局Lipschitz系数给收敛性分析带来的主要困难,我们成功地建立了跳适应后向Euler方法的强收敛性结果并得到相应的收敛率.最后,我们通过数值试验对前文所得理论结果做进一步的验证.

本文利用多边形网格上的间断有限元方法离散二阶椭圆方程,在曲边区域上,采用多条直短边逼近曲边的以直代曲的策略,实现了高阶元在能量范数下的最优收敛.本文还将这一方法用于带曲边界面问题的求解,同样得到高阶元的最优收敛.此外我们还设计并分析了这一方法的\linebreakW-cycle和Variable V-cycle多重网格预条件方法,证明当光滑次数足够多时,多重网格预条件算法一致收敛.最后给出了数值算例,证实该算法的可行性并验证了理论分析的结果.

本文提出求解数学物理方程大型离散特征值问题的几何网格预变换块因式分解算法(简称GPA算法).
通过长期研究我们发现:结构化网格矩阵$G$满足幂等方程$G^m=I_N,(m\ll N={\rm dim}(G))$,故可在实数域或复数范围内进行因式分解;且$G$与有限元刚度矩阵$A$之间乘法存在互易性:$A\cdot G=G\cdot A$,利用$G$的几何不变性可把$N$阶大型矩阵$A$正交分解为$m-$块对角块矩阵异步并行是我们算法的计算数学基础.
本文以正三角形、方形、平行六边形及正十七边形等结构化网格为例,特别是详细分析了六边形上的离散特征值异步并行算法及程序实现细节.文后附有若干2-3万阶量级离散矩阵特征值的桌面电脑数值计算例子(正三角形与方形网格,串行加速比分别为3-4倍),符合本文算法分析得出的"几何网格预处理的并行度与正多边形边数成正比"的结论.这类几何网格因式分解算法原则上可推广到三维乃至高维数学物理方程离散特征值计算问题，也可用于大型线性方程组的高效并行求解.

为了高效求解中小型线性互补问题，本文提出了改进的分块模方法，并证明了关于严格对角占优（对角元素均为正数）线性互补问题的收敛性.对于广义对角占优线性互补问题，先将其转化为严格对角占优线性互补问题，再采用改进的分块模方法求解.数值结果表明，改进的分块模方法在求解广义对角占优线性互补问题时在内迭代次数和计算时间上均明显优于分块模方法.

本文研究球面上的$\ell_1$正则优化问题，其目标函数由一般光滑函数项和非光滑$\ell_1$正则项构成，且假设光滑函数的随机梯度可由随机一阶oracle估计.这类优化问题被广泛应用在机器学习，图像、信号处理和统计等领域.根据流形临近梯度法和随机梯度估计技术，提出一种球面随机临近梯度算法.基于非光滑函数的全局隐函数定理，分析了子问题解关于参数的Lipschtiz连续性，进而证明了算法的全局收敛性.在基于随机数据集和实际数据集的球面$\ell_1$正则二次规划问题、有限和SPCA问题和球面$\ell_1$正则逻辑回归问题上数值实验结果显示所提出的算法与流形临近梯度法、黎曼随机临近梯度法相比CPU时间上具有一定的优越性.

本文研究带有五次项的非线性Schrödinger方程初边值问题的有限差分法，其中方程中二阶偏导数项的系数、五次项的系数及初值满足下面的条件（1.6）.针对此问题，我们研究了一个守恒差分格式，在条件（1.6）下，差分解的$L^{\infty}$模先验估计被得到.在此基础上，我们得到了差分解最优$L^2$模的误差估计.

2017年，李昭祥等提出了一种偏牛顿-校正法（Partial Newton-Correction Method，简记为PNC方法），并利用它成功地计算出了三类非线性偏微分方程的多重不稳定解.本文在PNC方法的基础上，提出并发展了一种改进的PNC方法.首先，利用Nehari流形$\mathcal{N}$与零平凡解的可分离性，建立并证明了$\mathcal{N}$的某特殊子流形$\mathcal{M}$上的全局分离定理及其推广（即局部分离定理）.全局分离定理只跟非线性偏微分算子或相应的非线性泛函本身有关，而与具体的计算方法无关.对一些典型的非线性偏微分方程多解问题（比如，Henon方程问题），该全局分离定理的分离条件，经验证是成立的.另一个方面，通过修改或补充原辅助变换的定义，去掉了原辅助变换的奇异性；接着建立并证明了某些非线性偏微分方程问题的新未知解与该非线性偏微分算子零核空间的密切关系；在证明中，去掉了在原奇异变换下所需的标准收敛（standard convergence）假设.最后，计算实例与数值结果验证了改进的PNC方法的可行性和有效性；同时表明子流形$\mathcal{M}$与已知解的可分离性是PNC方法和本文新方法能成功找到多解的关键.

基于已有文献的研究成果及前期工作，我们考察了非线性弱奇性Volterra积分方程（VIE）的谱配置法，并对该方法进行了收敛性分析.得到的结论是数值误差呈谱收敛.误差收敛阶与配置点个数及方程解的正则性相关.数值实验也证实了这一结论.本文的方法解决了已有文献中类似数值方法（Allaei（2016），Sohrabi（2017））存在的问题.

卷积位势广泛存在于科学和工程领域, 它的高效高精度计算往往是数值仿真的瓶颈.卷积位势是典型的非局部积分, 卷积核函数通常在原点或者无穷远处具有奇异性,密度函数是光滑速降函数并可能具有较强的各向异性. 无论是从卷积还是从傅里叶积分出发,我们首先将全空间截断到有界矩形区域并将其等距离散, 再应用傅里叶谱方法来高精度逼近密度函数.理想的求解器需要在保证高精度的同时, 尽可能提高计算效率, 并妥善处理各向异性密度函数的情形.本文详细回顾了目前流行的三类基于积分方程的高精度快速算法, 包括基于非均匀快速傅里叶变换的算法、基于高斯和的算法与核截断算法. 它们都能达到谱精度, 计算效率都类似于离散快速傅里叶变换(FFT), 并都能处理各向异性的密度函数. 这三类算法具有离散卷积结构;一旦生成了离散张量, 位势的计算将转化为两倍长度向量的傅里叶变换, 计算效率达到了近似最优,且与各向异性强度无关.最后我们介绍了误差估计的已有结果, 并用实例从精度、效率和各向异性等方面展示了算法能力.

在"求解加权线性最小二乘问题的一类预处理GAOR方法"一文中，作者提出了求解加权线性最小二乘问题等价$2\times 2$块线性系统的一类预处理GAOR方法，并给出了几个比较定理来说明新提出预处理GAOR方法的优越性.本文我们将指出该文中几个比较定理的不完善之处和证明的错误之处，并给出正确的证明.

应用求解算子方程的Ulm方法构造了求解一类矩阵特征值反问题（IEP）的新算法.所给算法避免了文献[Aishima K.，A quadratically convergent algorithm based on matrix equations for inverse eigenvalue problems，Linear Algebra and its Applications，2018，542：310-33]中算法在每次迭代中要求解一个线性方程组的不足，证明了在给定谱数据互不相同的条件下所给算法具有根收敛意义下的二次收敛性.数值实验表明本文所给算法在矩阵阶数较大时计算效果优于上文所给算法.

设$n+1$个$m\times n（m\geq n）$实矩阵$\{A_i\}_{i=0}^n$和给定的$n$个正数$\{\sigma_i^{*}\}_{i=1}^n$.本文研究如下的逆奇异值问题：求$n$个实数$\{c_i^{*}\}_{i=1}^n$，使得矩阵$A_0+c_1^{*}A_1+\cdots +c_n^{*}A_n$有奇异值$\{\sigma_i^*\}_{i=1}^n.$基于矩阵方程，我们给出了求解逆奇异值问题的一个新的算法，并证明了它的二阶收敛特性.该算法可以看成是Aishima[Linear Algebra and its Applications，2018，542：310-333]中逆对称特征值问题算法的推广.数值例子表明算法的有效性.

针对一类变延迟微分方程，应用全隐式方法—平衡方法，研究了其收敛性和稳定性.结果表明平衡方法以$\frac{1}{2}\gamma，\gamma\in（0，1]$阶收敛到精确解；并且强平衡方法和弱平衡方法都能保持解析解的均方稳定性；进一步数值实验验证了算法理论分析的正确性，并且表明全隐式的平衡方法比显式方法—Euler方法具有更好的稳定性.

本文基于一维椭圆型界面问题提出了一种有效的数值方法.首先,根据模型构建一个崭新的破裂再生核空间.其次,应用破裂再生核方法给出了此类界面问题的近似解,并讨论该方法的收敛性.最后,通过几个有效的数值算例来说明该方法的精确性和稳定性.

四阶不完全对称张量的M-特征值在非线性弹性材料分析中有着广泛的应用.本文的目的是给出四阶不完全对称张量M-特征值的新包含域，得到最大M-特征值上界更精确的估计，并将得到的上界估计值应用到计算最大M-特征值的WQZ算法中，数值例子验证了结果的有效性.最后，基于得到的包含域，给出了四阶不完全对称张量正定性判定的充分条件.

增广拉格朗日方法是求解带线性约束的凸优化问题的有效算法.线性化增广拉格朗日方法通过线性化增广拉格朗日函数的二次罚项并加上一个临近正则项,使得子问题容易求解,其中正则项系数的恰当选取对算法的收敛性和收敛速度至关重要.较大的系数可保证算法收敛性,但容易导致小步长.较小的系数允许迭代步长增大,但容易导致算法不收敛.本文考虑求解带线性等式或不等式约束的凸优化问题.我们利用自适应技术设计了一类不定线性化增广拉格朗日方法,即利用当前迭代点的信息自适应选取合适的正则项系数,在保证收敛性的前提下尽量使得子问题步长选择范围更大,从而提高算法收敛速度.我们从理论上证明了算法的全局收敛性,并利用数值实验说明了算法的有效性.

研究含参数$l$非方矩阵对广义特征值极小扰动问题所导出的一类复乘积流形约束矩阵最小二乘问题.与已有工作不同,本文直接针对复问题模型,结合复乘积流形的几何性质和欧式空间上的改进Fletcher-Reeves共轭梯度法,设计一类适用于问题模型的黎曼非线性共轭梯度求解算法,并给出全局收敛性分析.数值实验和数值比较表明该算法比参数$l=1$的已有算法收敛速度更快,与参数$l=n$的已有算法能得到相同精度的解.与部分其它流形优化相比与已有的黎曼Dai非线性共轭梯度法具有相当的迭代效率,与黎曼二阶算法相比单步迭代成本较低、总体迭代时间较少,与部分非流形优化算法相比在迭代效率上有明显优势.

本文针对广义绝对值方程,提出了基于牛顿法的矩阵多分裂方法.并在该方法的基础上进一步改进,得到了基于牛顿法的交替矩阵多分裂方法.给出两种算法在一定条件下的全局收敛性,并分析当分裂为H分裂时，基于牛顿法的矩阵多分裂方法的收敛条件.通过数值实验验证了所提出的算法的可行性和有效性.

基于指数罚函数,对最近提出的一种求解无约束优化问题的三项共轭梯度法进行了修正,并用它求解更复杂的大规模极大极小值问题.证明了该方法生成的搜索方向对每一个光滑子问题是充分下降方向,而且与所用的线搜索规则无关.以此为基础,设计了求解大规模极大极小值问题的算法,并在合理的假设下,证明了算法的全局收敛性.数值实验表明,该算法优于文献中已有的类似算法.

首先利用变分原理和最优化理论得到了原问题的等价最优性条件；其次构造了椭圆最优控制问题分裂正定混合有限元方法的逼近格式；再次通过引入一些重要的中间变量和投影算子，并利用投影算子的相关性质，结合分裂正定混合有限元本身的逼近结果，得到了椭圆最优控制问题分裂正定混合有限元方法的超收敛性；最后数值实验结果验证了所得理论结果的正确性.

基于HS共轭梯度法的结构，本文在弱假设条件下建立了一种求解凸约束伪单调方程组问题的迭代投影算法.该算法不需要利用方程组的任何梯度或Jacobian矩阵信息，因此它适合求解大规模问题.算法在每一次迭代中都能产生充分下降方向，且不依赖于任何线搜索条件.特别是，我们在不需要假设方程组满足Lipschitz条件下建立了算法的全局收敛性和R-线收敛速度.数值结果表明，该算法对于给定的大规模方程组问题是稳定和有效的.

本文将局部投影稳定化（LPS）方法和连续时空有限元方法相结合研究对流扩散反应方程，给出稳定化连续时空有限元离散格式.与传统的时空有限元研究思路不同，时间方向利用Lagrange插值多项式，解耦时间和空间变量，降低时空有限元解的维数，具有减少计算量和简化理论分析的优点.通过引入Legendre多项式给出了有限元解的稳定性分析，进一步引进Lobatto多项式证明了有限元解的全局L^∞（L²）和局部L²（J_n；LPS）范数误差估计.最后给出数值算例验证理论分析的正确性，以及稳定化格式的可行性和有效性.