计算差分方程的收敛点_数值计算中的有限差分法

早期,在科学技术的发展过程中,理论分析和实验研究一直是两种主要的科学方法,但是随着计算机的出现与发展,情况就大大不同了,研究者可以采用计算机去计算过去那些根本不能求解的问题,模拟一些难以观测的现象,第三种科学方法——数值计算诞生了。

有限差分法作为一种最为朴素的计算方法,发展时间也最为悠久,至今依然活跃。有限差分法的朴素思想就是,在求解连续性问题时,将时空域离散后,利用差商代替微商求解问题。

从应用上看,有限差分法常使用在求解时间依赖问题(双曲型和抛物型方程),而有限元主要求解稳态问题(椭圆型方程)。从建立差分格式的方法上看,有限差分法既可以根据Taylor级数展开建立差分格式,也可以根据积分方法建立差分格式,其中,积分方法由于其本身的特殊性,也被单独称为有限体积法。

从微分方程建立差分方程,方法格式层出不穷,但三个基本性质确是核心所在,即差分格式的相容性、收敛性和稳定性。

  • 相容性:当时间和空间离散步长足够小时,差分方程将与微分方程充分的接近。
  • 收敛性:当时间和空间离散步长足够小时,差分方程的解将逼近微分方程的解。
  • 稳定性:有限差分格式是按照时间层逐层推进,如果初始计算引入的误差,随着逐层计算被控制,那么这个差分格式就是稳定的。

上述三个性质,通常是互相独立的,也就是相容的差分格式并不一定收敛。

通常,判断一个差分格式

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
有限分初学者必备】如何根据问题的特点将定解区域作网格剖分;如何把原微分方程离散化为方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确,还需从理论上分析方程组的性态,包括解的唯一性、存在性和分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种分格式,为了有实用意义,一个基本要求是它们能够任意逼近微分方程,这就是相容性要求。另外,一个分格式是否有用,最终要看方程的精确解能否任意逼近微分方程的解,这就是收敛性的概念。此外,还有一个重要的概念必须考虑,即分格式的稳定性。因为分格式的计算过程是逐层推进的,在计算第n+1层的近似值时要用到第n层的近似值 ,直到与初始值有关。前面各层若有舍入误,必然影响到后面各层的值,如果误的影响越来越大,以致分格式的精确解的面貌完全被掩盖,这种格式是不稳定的,相反如果误的传播是可以控制的,就认为格式是稳定的。只有在这种情形,分格式在实际计算的近似解才可能任意逼近方程的精确解。关于分格式的构造一般有以下3种方。最常用的方数值微分,比如用商代替微商等。另一方叫积分插值,因为在实际问题得出的微分方程常常反映物理上的某种守恒原理,一般可以通过积分形式来表示。此外还可以用待定系数构造一些精度较高的分格式。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值