值方法的收敛性和绝对稳.PDF

数值解中误差的积累、数 值方法的收敛性和绝对稳 定性 一阶微分方程组与高阶方 程的数值解法 计计计算算算方方方法法法 第27讲 • – 1 / 22 2018.12.21 数值解中误差的积 累、数值方法的收敛 ⊲ 性和绝对稳 定性 数值解中误差的积累和数 值方法的收敛性 绝对稳 定性 一阶微分方程组与高阶方 程的数值解法 9.2 数数数值值值解解解中中中误误误差差差的的的积积积累累累、、、数数数值值值方方方法法法的的的收收收敛敛敛性性性 和和和绝绝绝对对对稳稳稳定定定性性性 第27讲 • – 2 / 22 2018.12.21 9.2 数数数值值值解解解中中中误误误差差差的的的积积积累累累、、、数数数值值值方方方法法法的的的收收收敛敛敛性性性和和和绝绝绝对对对稳稳稳定定定性性性 数值解中误差的积累、数 9.2.1 数数数值值值解解解中中中误误误差差差的的的积积积累累累和和和数数数值值值方方方法法法的的的收收收敛敛敛性性性 值方法的收敛性和绝对稳 定性 数值解中误差的积累 显式单步法可以写为如 下的一般形式 ： ⊲ 和数值方法的收敛性 绝对稳 定性 一阶微分方程组与高阶方 y = y + hφ(x , y , h). 程的数值解法 i+1 i i i 假设该法是p 阶方法，即局部截断误差R [y] = O(hp+1) . 则 y (x ) = y (x ) + hφ(x , y(x ), h) + R [y]. i+1 i i i 并假设max |y (p+1) (x)| ≤ Mp+1 ，则|R [y] | ≤ cMp+1hp+1 . a≤x≤b 在 实际计算中由于存在舍入误差，实际计算得到的y˜ 不是理论上 i+1 的数值解yi+1 ，而是yi+1 的近似值，它应 满足 y˜ = y˜ + hφ(x , y˜ , h) + η . i+1 i i i i 其中η 表示第i 步计算所产生的舍入误差 成的误差. i 第27讲 • – 3 / 22