常微分方程数值解法yjs.pptVIP

第6章 常微分方程的数值解法 例6.1：考虑初值问题： 取步长h = 0.1，并与准确解 比较。 解：因为xi = 0.1i，而f(x, y) = -y + x，故Euler格式 yi+1=yi+hf(xi, yi) = (1-h)yi + hxi (i=0,1,…) 隐式Euler格式为 yi+1=yi+hf(xi+1, yi+1) = yi – hyi+1+hxi+1 (i=0,1,…) 整理得 yi+1= (yi+hxi+1 ) / (1+h) (i=0,1,…) 用Matlab计算：clear;x=0,y=0; for i=1:5 y=(1-h)*y+h*x x=x+h; end 当?-1?0时,公式为隐式公式,反之为显式公式.参数??i,?i?的选择原则是使方法的局部截断误差为 y(xn+1)-yn+1=O(hr+2) 例 选取参数?,?0,?1,?2,使三步方法 yn+1=?yn+h(?0?n+?1?n-1+?2?n-2) 这里,局部截断误差是指,在yn-i=y(xn-i), i=0,1,…,r的前提下,误差y(xn+1)-yn+1. 为三阶方法. 解 设yn=y(xn),yn-1=y(xn-1),yn-2=y(xn-2),则有 ?n=?(xn , y(xn))=y?(xn) y(xn+1)=y(xn)+hy?(xn)+1/2h2y??(xn)+1/6h3y???(xn) 于是有 若使: y(xn+1)-yn+1=O(h4) ,只要?,?0,?1,?2满足: ?n-1=?(xn-1,y(xn-1))=y?(xn-1)=y?(xn-h) =y?(xn)-hy??(xn)+1/2h2y???(xn)-1/6h3y(4)(xn)+O(h4) ?n-2=y?(xn)-2hy??(xn)+2h2y???(xn)-4/3h3y(4)(xn)+O(h4) yn+1=?y(xn)+h(?0+?1+?2)y?(xn)-h2(?1+2?2)y??(xn) +h3(1/2?1+2?2)y???(xn)-h4/6(?1+8?2)y(4)(xn)+O(h5) +1/24h4y(4)(xn)+O(h5) ?=1, ?0+?1+?2=1, ?1+2?2=-1/2 , ?1+4?2=1/3 解之得: 于是有三步三阶显式差分公式 yn+1=yn+h/12(23?n-16?n-1+5?n-2) 科大研究生学位课程 数值分析 * 数值分析 数值分析 实际中，很多问题的数学模型都是微分方程。我们可以研究它们的一些性质。但是，只有极少数特殊的方程有解析解。对于绝大部分的微分方程是没有解析解的。 因此，需要研究和掌握微分方程的数值解法，即计算解域内离散点上的近似值的方法。 ? 考虑一阶常微分方程的初值问题 首先考虑微分方程解的存在性。 只要 f (x, y) 在[a, b] ? R1 上连续，且关于 y 满足 Lipschitz 条件，即存在与 x, y 无关的常数 L 使 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则上述问题存在唯一解。 称(6.1)的解y(x)在节点xi处的近似值 yi ? y(xi) a x1 x2 ... xn = b. 为其数值解，方法称为数值方法。 注：① 考虑等距节点: xi = a + ih，h = (b – a)/n. ② 从初始条件y(a) = y0出发，依次逐个计算y1，y2，…，yn的值，称为步进法。 两种：单步法、多步法。 微分方程的数值方法。求的是在一系列离散点列上，未知函数y(x)在这些点上的值的近似。 为了考察数值方法提供的数值解，是否有实用价值，需要知道如下问题： ① 步长充分小时，所得到的数值解能否逼近问题得真解； 即收敛性问题 ② 误差估计，主要考虑局部截断误差。 产生的舍入误差，在以后得各步计算中，是否会无限制扩大；稳定性问题 6.1 欧拉Euler 法与改进欧拉法 1.欧拉法： x0 x1 向前差商近似导数 记为 定义 　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差