抛物型方程快速计算的新算法

来源 :湖南师范大学 | 被引量 : 1次 | 上传用户:lala_
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
目前在统计物理、概率论、量子力学、生物化学、声热同时传播、弹性振动、水文地质、石油开采等问题的研究时都会涉及到抛物型方程。因此,对于这类方程的快速计算研究有重要的意义,但是这类方程求解的传统思想都是按时间方向分层计算,逐层推进的,所以本文研究的方法就是提高每层的计算速度,从而得到一种对这类方程快速推进求解的新算法。本文采用的是精度高,稳定性好的crank-Nicolson格式进行离散,有限差分方法离散矩阵所需矩阵带宽小,有利于今后高维问题的研究。   本文采用的是CG迭代法求解,对于一般的迭代法来说,迭代误差都是由初始误差和压缩比两个因素来决定的。但是在网格比r很大的情况下,条件数也变大,从而使压缩比也随之接近1,那么此时很难依靠压缩比来加快计算速度,所以本文所研究的重点是如何为时间方向下一层的计算提供一个好的初值,以降低初始误差来加快计算速度。主要讨论如下:   1通过对常微初值问题和抛物初边值问题的研究,发现经典外推公式能够很准确地去逼近原方程的解,因此采用经典外推公式和共轭梯度法以固定迭代次数的模式来对抛物型方程求解,发现其比单纯CG法的计算结果要好,但是发现高次外推公式存在着误差放大和积累,不稳定性等问题;   2.针对前一阶段出现的问题,我们对经典公式做出改变,提出了一些新的外推公式,采用以真实误差控制的模式进行研究,发现外推公式结合CG迭代法求解的平均迭代次数只有单纯CG迭代法的1/3~1/6,而且稳定性高;   3.为了更好地考察CG迭代的效果,我们改用迭代余量控制模式进行计算,这样更符合解决实际问题的计算,计算的结果比较好。
其他文献
图的对称性研究一直是群与图研究中的热门课题.随着计算机及网络的发展,图的对称性研究的应用价值日益显露,这也引起了学者们极大的研究兴趣.图的对称性是通过图的全自同构群在
类似于遗传算法,群智能算法也属于启发式算法的一种.早在20世纪90年代初,已存在通过模拟自然界生物的社会行为来构造随机算法的思想.研究者对生物的群体行为进行模拟,提出群智能算
一、项目背景随着自动化技术的不断发展,自动化检测是未来计量发展的方向。近几年,国外主要比较仪生产厂家纷纷致力于研制毫克组、克组和公斤组砝码的自动检测机器人,以提高
数字图像处理的很多应用,需要对图像进行多尺度操作。尤其是在计算摄影学中,经常需要把一幅图像分解为一个分片光滑的基本层和一个或者多个细节层。图像的这种多尺度分解的关
带有causal算子的函数方程能够把常微分方程,积分微分方程,有限或无穷时滞微分方程,Volterra积分方程和中立型微分方程诸如此类微分方程进行有机的统一,在工程等技术领域具有广泛
在KNT-PGK3型高级过程控制实验装置的基础上,研究了基于超声波传感器的水箱液位模糊PID控制。具体是在现有实验装置与i FIX组态软件环境下采用模糊PID控制方法对实验台的水箱
本文研究了三维带粘性的液体-气体两相流模型古典解的全局存在性和唯一性.在这篇文章中我们假设了初始数据的能量模足够小.我们的主要想法是来源于参考文献[16],用连续性技巧
随着“互联网+”的到来,物联网、云计算、大数据等这些新一代的信息技术逐渐出现人们的面前.虽然这给人们的生活带来了极大的便利,但是面对规模如此庞大的数据信息,特别针对求解
现阶段我国2G移动通信管制政策的问题rn1、移动市场的不均衡格局rn从市场结构看,目前我国移动市场只有两家移动运营商.其中,中国移动占据了近75%的市场价值份额和66%的用户市
删失回归模型,又称为Tobit模型,是一种响应变量受非负限制的模型,删失回归模型是一种非常重要的模型,在计量经济学中有许多的应用,本文首先用经验过程方法,基于p分位数获得了删失回