《数值分析第七章常微分方程的数值解法二.优秀PPT.ppt》由会员分享，可在线阅读，更多相关《数值分析第七章常微分方程的数值解法二.优秀PPT.ppt（18页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第七章第七章 常微分方程的数值解法（二）常微分方程的数值解法（二）第四节第四节 线性多步法线性多步法第五节第五节 单步法的收敛性与稳定单步法的收敛性与稳定性性第六节第六节 一阶方程组和高阶方一阶方程组和高阶方程程第七节第七节 边值问题的数值解法边值问题的数值解法上页上页上页上页下页下页下页下页返回返回返回返回4 线性多步法线性多步法用用若干若干节点处的节点处的 y 及及 y 值的值的线性组合线性组合来近似来近似y(xi+1).).(.110111101kikiiikikiiiffffhyyyy +=a aa aa a其通式可写为：其通式可写为：当当 1 0 时，为时，为隐式公式隐式公式;1=0
2、 则为则为显式公式显式公式.基于数值积分的构造法基于数值积分的构造法将将 在在 上积分，得到上积分，得到只要只要近似地算出右边的积分近似地算出右边的积分 ，则可通过，则可通过 近似近似y(xi+1).而而选用不同近似式选用不同近似式 Ik，可得到不同的计算公式，可得到不同的计算公式.上页上页上页上页下页下页下页下页返回返回返回返回 亚当姆斯显式公式亚当姆斯显式公式利用利用k+1 个节点上的被积函数值个节点上的被积函数值 构造构造 k 阶牛顿阶牛顿后后插插多项式多项式 ,有有Newton插值余项插值余项/*显式计算公式显式计算公式*/局部截断误差为：局部截断误差为：例例1 k=1 时有时有上页上
3、页上页上页下页下页下页下页返回返回返回返回注：注：一般有一般有 ，其中，其中Bk 与与yi+1 计算公式中计算公式中 fi,fi k 各项的各项的系数系数均可查表得到均可查表得到.10123kfifi 1fi 2fi 3BkMisprint on p.106常用的是常用的是 k=3 的的4阶亚当姆斯显式公式阶亚当姆斯显式公式上页上页上页上页下页下页下页下页返回返回返回返回 亚当姆斯隐式公式亚当姆斯隐式公式利用利用k+1 个节点上的被积函数值个节点上的被积函数值 fi+1,fi,fi k+1 构造构造 k 阶阶牛顿牛顿前插前插多项式。与显式多项式完全类似地可得到一系列多项式。与显式多项式完全类似
4、地可得到一系列隐式公式隐式公式，并有，并有 ，其中，其中 与与 fi+1,fi,fi k+1 的系数亦可查表得到的系数亦可查表得到.10123kfi+1fifi 1fi 2Bk常用的是常用的是 k=3 的的4阶亚当姆斯隐式公式阶亚当姆斯隐式公式小于小于Bk较同阶显较同阶显式式稳定稳定上页上页上页上页下页下页下页下页返回返回返回返回 亚当姆斯预料亚当姆斯预料-校正系统校正系统 Step 1:用用Runge-Kutta 法法计算前计算前 k 个初值；个初值；Step 2:用用Adams 显式计算预料值；显式计算预料值；Step 3:用同阶用同阶Adams 隐式隐式计算计算校正校正值值.留意：三步所
5、用公留意：三步所用公式的精度必需相同。式的精度必需相同。通常用经典通常用经典Runge-Kutta 法协作法协作4阶阶Adams 公式公式.4阶阶Adams隐式公式的截断误差为隐式公式的截断误差为4阶阶Adams显式公式的截断误差为显式公式的截断误差为当当 h 充分小充分小时，可近似认为时，可近似认为 i i，则：，则：Predicted value pi+1Modified value mi+1Corrected value ci+1Modified final value yi+1外推技术外推技术/*extrapolation*/上页上页上页上页下页下页下页下页返回返回返回返回上页上页上页
6、上页下页下页下页下页返回返回返回返回5 收敛性和稳定性收敛性和稳定性一、一、收敛性收敛性 若若某某算算法法对对于于任任意意固固定定的的 x=xi=x0+i h，当当 h0(同时同时 i )时有时有 yi y(xi)，则称该算法是，则称该算法是收敛收敛的的.例例3 就初值问题就初值问题 考察欧拉显式格式的收敛性考察欧拉显式格式的收敛性.解：解：该问题的精确解为该问题的精确解为 尤拉公式为尤拉公式为对任意固定的对任意固定的 x=xi=i h，有，有 上页上页上页上页下页下页下页下页返回返回返回返回上页上页上页上页下页下页下页下页返回返回返回返回0.00.10.20.30.40.5精确解精确解改进尤
7、拉法改进尤拉法 尤拉隐式尤拉隐式尤拉显式尤拉显式 节点节点 xi二、二、稳定性稳定性例例4 考察初值问题考察初值问题 在区间在区间0,0.5上的解上的解.分别用欧拉显、隐式格式和改进的欧拉格式计算数值解分别用欧拉显、隐式格式和改进的欧拉格式计算数值解.1.0000 2.0000 4.0000 8.0000 1.6000 101 3.2000 101 1.00002.5000 10 1 6.2500 10 21.5625 10 23.9063 10 39.7656 10 41.00002.50006.25001.5626 1013.9063 1019.7656 1011.00004.9787 1
8、0 22.4788 10 31.2341 10 46.1442 10 63.0590 10 7上页上页上页上页下页下页下页下页返回返回返回返回若若某某算算法法在在计计算算过过程程中中任任一一步步产产生生的的误误差差在在以以后后的的计计算中都算中都逐步衰减逐步衰减，则称该算法是，则称该算法是绝对稳定的绝对稳定的.一般分析时为简洁起见，只考虑试验方程一般分析时为简洁起见，只考虑试验方程常数，可以常数，可以是复数是复数当步长取为当步长取为 h 时，将某算法应用于上式，并假设只在初值时，将某算法应用于上式，并假设只在初值产生误差产生误差 ，则若此误差以后逐步衰减，就称该算法，则若此误差以后逐步衰减，就
9、称该算法相对于相对于 绝对稳定绝对稳定，的全体构成的全体构成绝对稳定区域绝对稳定区域.我们称我们称算法算法A 比算法比算法B 稳定稳定，就是指，就是指 A 的绝对稳定区域比的绝对稳定区域比 B 的的大大.hl l h=h上页上页上页上页下页下页下页下页返回返回返回返回例例5 考察显式欧拉法考察显式欧拉法由此可见，要保证初始误差由此可见，要保证初始误差 0 以后逐步衰减，以后逐步衰减，必须满足：必须满足：0-1-2ReImg例例6 考察隐式欧拉法考察隐式欧拉法可见确定稳定区域为：可见确定稳定区域为：210ReImg注：一般来说，隐式尤拉法的确定稳定性比同阶的显式法注：一般来说，隐式尤拉法的确定稳
10、定性比同阶的显式法的好的好.上页上页上页上页下页下页下页下页返回返回返回返回例例7 隐式龙格隐式龙格-库塔法库塔法而显式而显式 1 4 阶方法的确定稳定阶方法的确定稳定区域为区域为其中其中2阶方法阶方法 的绝对稳定区域为的绝对稳定区域为0ReImgk=1k=2k=3k=4-1-2-3-123ReImg无条件稳定无条件稳定上页上页上页上页下页下页下页下页返回返回返回返回6 微分方程组与高阶方程微分方程组与高阶方程一、一、一阶微分方程组一阶微分方程组IVP的一般形式为：的一般形式为：=)(,.),(,()(.)(,.),(,()(1111xyxyxfxyxyxyxfxymmmm初值初值000202
11、0101)(,.,)(,)(mmyxyyxyyxy=将问题记作将问题记作向量形式向量形式，令：，令：前述全部公式皆前述全部公式皆适用于向量形式适用于向量形式.上页上页上页上页下页下页下页下页返回返回返回返回二、二、高阶微分方程高阶微分方程=10)1(1000)1()()(,.,)(,)(),.,(nnnnaxyaxyaxyyyyxfy化作化作一阶微分方程组一阶微分方程组求解求解.引入新变量引入新变量初值条件为：初值条件为：上页上页上页上页下页下页下页下页返回返回返回返回上页上页上页上页下页下页下页下页返回返回返回返回7 边值问题的数值解法边值问题的数值解法2 阶常微分方程边值问题阶常微分方程边值问题 打靶法打靶法先猜测一个初始斜率先猜测一个初始斜率 y (a)=s，通过解初值，通过解初值问题问题y(b)=(s)找出找出s*使得使得(s*)=，即把问，即把问题转化为求方程题转化为求方程 (s)=0 的根的根.yx0aby x()斜率斜率=s0()s0斜率斜率=s1()s1每计算一个每计算一个(s)都必需解一个都必需解一个ODE.上页上页上页上页下页下页下页下页返回返回返回返回 有限差分法有限差分法将求解区间将求解区间a,b 等分为等分为N 份，取节点份，取节点 xi=a+ih (i=0,N)，在每一个节点处将，在每一个节点处将 y 和和 y 离散化离散化.泰勒展开泰勒展开