若干非线性问题特殊解的迭代解法

摘要

在现代科学工程计算领域中，会遇到一些具有许多解（或方案）的问题.对于这类问题，人们感兴趣的往往是它的一些特殊解.本文研究的问题便属于这类问题.对核子物理中提出的非对称代数Riccati方程，人们关心的是它的最小非负解.对多材料鼓面形状设计问题和非线性磁材料中的裂纹和不均匀性检测问题，人们想得到的是最优形状.如何快速有效地求得人们所关心的解，吸引了众多研究者的兴趣.
　　关于非线性方程的求解，已经存在很多成熟的迭代方法，如牛顿法，Eu-ler迭代族和Halley迭代族等.对于同一个问题，不同迭代法的优劣性也会因Jacobian矩阵的计算和函数求值的复杂性而不同.King-Werner迭代法自被King和Werner分别提出后，因其每步迭代中Jacobian矩阵和函数值的计算次数与牛顿法相同，而其收敛阶在非奇异时为1+√2的优势而引起了研究者的注意.对于最佳形状设计和形状重构问题，其用数学刻画是个约束优化问题，这可以有许多算法求解，如拉格朗日乘子法，增广拉格朗日法等.
　　本文将求解三个问题的特殊解.第一个问题是在中子迁移理论中提出的非对称代数Riccati方程的最小正解，第二个问题是鼓面设计中的符合特定设计要求的最佳形状问题，第三个问题是非线性磁材料中的裂纹和材料不均匀性检测问题.
　　在第二章中，我们首先介绍一些已有的求解非对称代数Riccati方程的最小正解的算法，然后用King-Werner方法求非对称代数Riccati方程的最小正解并给出King-Werner迭代法求解该非线性方程的最小正解的收敛性分析以及误差分析.数值例子表明King-Werner方法求奇异情况下的最小正解具有一定的优势.在求解奇异的非对称代数Riccati方程的最小正解时，许多已有的算法迭代速度会比较慢，借助Guo等提出的移技术，一些算法收敛速度慢的问题得以解决.在移技术中，会涉及到一个正实参数，并且该参数的不同取值对迭代算法的收敛速度的影响会有所不同.我们将算法得到的迭代点列看作是该参数的函数列，通过分析该参数对函数列的影响来分析该参数对简单迭代法，修改的简单迭代法，非线性块雅可比迭代法和非线性块高斯-赛德尔迭代法四种迭代法求最小正解时的收敛速率的影响.随后的数值例子验证了理论分析的正确性.
　　在第三章中，我们将研究鼓面的最佳材料分布问题.我们首先简单介绍目标函数关于密度函数的G(a)teaux导数，然后选取密度函数的一个能使目标函数的增量尽量大于零的更新方向，通过这一方式来设计一个算法，最后通过数值实验来验证算法求解该类问题的有效性和合理性.
　　在第四章里，我们解决的问题是非线性磁材料中裂纹和材料不均性的检测.这个问题属于形状重构问题.我们首先用分片常数水平集方法来表示材料的组成，然后将该问题转化为一个约束优化问题，并通过拉格朗日乘子法将该约束优化问题转化成无约束优化问题，最后通过梯度法来求解该无约束问题.数值实验显示我们的算法不但对水平集函数的初始值依赖不大，而且有效稳定.

目录

声明

摘要

第一章 绪论

1．1 引言

1．2 迁移理论中提出的非对称代数Riccati方程

1．3 鼓面材料设计中的优化问题

1．4 非线性磁材料的缺陷检测

第二章 迁移理论中的非对称代数Riccati方程的最小非负解

2．1 研究现状

2．2 King-Werner迭代法求解方程(1．8)的最小正解

2．2．1 King-Werner迭代法

2．2．2 收敛性分析

2．2．3 误差分析

2．2．4 数值实验

2．3 一些方法移位技术中的参数分析

2．3．1 迭代格式及理论分析

2．3．2 数值实验

第三章 最佳鼓面设计中的优化问题

3．1 迭代格式推导与介绍

3．2 数值实验

第四章 非线性磁材料无损检测中的优化问题

4．1 分片常数水平集方法

4．2 正问题和反问题

4．3 算法的推导与实现

4．4 数值实验

第五章 工作展望

5．1 工作展望

第六章 附录

参考文献

发表及投稿文章目录

简历

致谢