ODEs系统的快速求解及预处理技术

来源 :电子科技大学 | 被引量 : 0次 | 上传用户:chao1028
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
300多年来,常微分方程逐渐发展成为拥有自己独立研究对象的一门理论学科。自然界与工程技术中的很多现象,其数学表述可归结为常微分方程定解问题。像万有引力定律、人口发展规律、市场均衡价格的变化等,都可以建立一定的常微分方程模型来求解。然而在生产实际和科学研究中所遇到的微分方程往往很复杂,在很多情况下都不可能求出解的解析表达式,而且高次代数方程求根也并不容易,所以想要求出常微分方程的解析解的想法往往是无法实现的。在实际计算中,主要采用数值解法。  最近,白中治和任志茹等已经研究了三阶线性常微分方程由Sinc方法离散所得到的线性方程组的预处理方法。本文的研究就是建立在此研究成果的基础之上,同时又是他们的成果的延伸。首先介绍了Sinc方法,然后通过引入变量把三阶的线性常微分方程降阶为与之等价的由两个线性常微分方程组成的常微分方程组。用Sinc方法离散二阶线性常微分方程组后得到的是一个系数矩阵是2?2的分块矩阵的线性方程组,并且每一块都是由对角矩阵和Toeplitz矩阵组合而成。结合这些已有的理论知识和研究成果,本文做了以下的研究工作。  根据Sinc方法所得的线性方程组系数矩阵的结构性质,为了能够更有效的利用Krylov子空间方法求解线性方程组,本文给出了一种新的预处理子PN。我们同时选择使用了ILU预条件子,并且实验证明了这两种预条件子的预处理矩阵的谱分布都非常的密集。  通过数值实验结果证明Sinc方法能够很好的求解三阶线性常微分方程。对离散后的线性方程组用Krylov子空间方法进行求解。与前人的预条件子P做了比较,从实验图表可以看出用本文给出的预处理子PN时所需的迭代步数较少,且随着N的取值的增加用PN所需的迭代步数的将近减少了一半。因此可以看出本文提出的预条件子效果更好。
其他文献
本文首先用延拓结构理论分析(2+1)维海森堡铁磁链(HF)模型.之后进一步讨论可积的(2+1)维M(修正的)HF模型,首先用延拓结构理论对该模型进行分析,并通过赋予在闵氏空间中运动的空
本文的目的是利用Jacobi(雅克比)级数的Poisson积分在边界处的渐近性态来刻划函数的光滑性。 与Jacobi级数有关的函数理论是数学中的一个重要领域,相关问题的研究已取得了
进入新世纪以后,随着社会主义市场经济体制的确立和完善,我国房地产业迅猛发展起来,与此同时,施工管理在整个建筑工程中的地位尤为重要,同时,为保证工程施工的质量,就要及时解决掉施
期刊
在本文中,我们主要应用非线性泛函分析中的半序理论,锥拉伸与锥压缩不动点理论,对一些非线性边值问题进行讨论,全文共分为五章。 第一章是本文的绪论部分.主要介绍了本文的研