单步法的收敛性和稳定性.pptVIP

第1页，共11页，星期日，2025年，2月5日8.3.1单步法的收敛性数值解法的基本思想就是要通过某种离散化方法，将微分方程转化为某种差分方程（例如，（8.1.8）式）来求解。这种转化是否合理，还要看差分方程的解，是否收敛到微分方程的准确解。定义8.3对于任意固定的，若对于初值问题（8.1.1）的显式单步法（8.1.8）产生的近似解，均有，则称该方法是收敛的。在定义中，是固定的点，当时有，n不是固定的。显然，若方法是收敛的，则在固定点处的整体截断误差趋于零。下面给出方法收敛的条件。定理8.1设初值问题（8.1.11）的单步法（8.1.8）是p阶的（），且函数满足对y的Lipschitz条件即存在常数，使第2页，共11页，星期日，2025年，2月5日对一切成立，则方法（8.1.8）收敛，且。因为（8.1.8）是p阶的，所以存在，当时有。再用的Lipschitz条件有为了方便，记，即有。由此可推得证仍记，根据局部截断误差的定义将此式与（8.1.8）相减得第3页，共11页，星期日，2025年，2月5日利用关系式可以得到现在取，有，于是有。定理得证。容易证明，如果（8.1.1）的满足Lipschitz条件是,且初值是正确的，则显示Euler法、改进的Euler法和R-K方法是收敛的。由定理8.1说明,f关于y满足Lipschitz条件是使单步收敛的充分条件，而且，还说明一个方法的整体截断误差比局部截断误差低一阶。所以，常常通过求出局部截断误差去了解整体截断误差的大小。单步法的显式形式（8.1.8）可写成（8.3.1）第4页，共11页，星期日，2025年，2月5日称为增量函数。对于收敛的方法，固定，有从而。对于（8.3.1），我们自然要考虑是否成立。这就是相容性问题。则称方法（8.1.8）与初值问题（8.1.1）是相容的。相容性说明数值计算的差分方程（8.3.1）趋于（8.1.1）中微分方程。我们本章讨论的数值方法都是与原初值问题相容的。定义8.4若方法（8.1.8）的增量函数满足第5页，共11页，星期日，2025年，2月5日8.3.2单步法的稳定性对于一种收敛的相容的差分方程，由于计算过程中舍入误差总会存在，我们需要讨论其数值稳定性。一个不稳定的差分方程会使计算解失真或计算失败。为了讨论方便起见。将（8.1.1）中的在解域内某一点作Taylor展开并局部线性化，即令利用线性化的关系，可得。因此，我们通过如下的试验方程第6页，共11页，星期日，2025年，2月5日（8.3.2）讨论数值方法的稳定性。当某一步有舍入误差时，若以后的计算中不会逐步扩大，则称这种稳定性为绝对稳定性。现在讨论显式Euler法的稳定性。将显式Euler法用于试验方程（8.3.2），有。当有舍入误差时，其近似值为