matlab_常微分方程数值解法.ppt

matlab_常微分方程数值解法 ppt课件 2021/3/26 matlab_常微分方程数值解法 ppt课件 例题： 2021/3/26 matlab_常微分方程数值解法 ppt课件 3、龙格-库塔(R-K)方法 1 欧拉方法 对于一般的常微分方程： 几种方法的比较 取前一点的斜率作为平均斜率，即： 整体截断误差： O(h) 2021/3/26 matlab_常微分方程数值解法 ppt课件 2 向后的欧拉方法 取后一点的斜率作为平均斜率，即： 整体截断误差： O(h) 2021/3/26 matlab_常微分方程数值解法 ppt课件 3 改进的欧拉方法 取两点斜率平均值，即： 整体截断误差： O(h2) 2021/3/26 matlab_常微分方程数值解法 ppt课件 欧拉方法（前、后）和改进的欧拉方法优点是算法简单，但计算精度较低，不能满足实际问题求解对精度的要求，所以使用较少。 龙格-库塔(R-K)方法就是一种精度较高的较为实用计算常微分方程的方法。 2021/3/26 matlab_常微分方程数值解法 ppt课件 从以上三种方法的比较，可以推测：若取多点处斜率的加权平均作为平均斜率，精度会更高，这就是龙格—库塔法的基本思想。 2021/3/26 matlab_常微分方程数值解法 ppt课件 四阶龙格—库塔法计算公式为： 整体截断误差： O(h4) 2021/3/26 matlab_常微分方程数值解法 ppt课件 2021/3/26 matlab_常微分方程数值解法 ppt课件 例题： time_Euler = 0.2500 time_EulerPro = 0.5000 time_RK4 = 1.0150 2021/3/26 matlab_常微分方程数值解法 ppt课件 4、 MATLAB的常微分方程函数 ode45 ode23 ode113 ode15s ode23s ode23t ode23tb 格式 [x,y]=ode45(′ fun′, [x0,xn], y0,option] 说明： 适用于求解一阶常微分方程组 fun定义微分方程组的函数文件名 [x0,xn]求解区域 y0初始条件向量 option可选参数，由ODESET函数设置，比较复杂 x输出自变量向量，y输出[y, y ′, y ″,..] 2021/3/26 matlab_常微分方程数值解法 ppt课件 没有一种算法可以有效地解决所有的 ODE 问题，因此MATLAB 提供了多种ODE函数。 函数 ODE类型 特点 说明 ode45 非刚性 单步法；4，5 阶 R-K 方法；累计截断误差为 (△x)3 大部分场合的首选方法 ode23 非刚性 单步法；2，3 阶 R-K 方法；累计截断误差为 (△x)3 使用于精度较低的情形 ode113 非刚性 多步法；Adams算法；高低精度均可到 10-3～10-6 计算时间比 ode45 短 ode23t 适度刚性 采用梯形算法 适度刚性情形 ode15s 刚性 多步法；Gear’s 反向数值微分；精度中等 若 ode45 失效时，可尝试使用 ode23s 刚性 单步法；2 阶Rosebrock 算法；低精度 当精度较低时，计算时间比 ode15s 短 ode23tb 刚性 梯形算法；低精度 当精度较低时，计算时间比ode15s短 2021/3/26 matlab_常微分方程数值解法 ppt课件 例题 二、三阶龙格库塔法 matlab_常微分方程数值解法 ppt课件 matlab_常微分方程数值解法 ppt课件 matlab_常微分方程数值解法 ppt课件 matlab_常微分方程数值解法 ppt课件 2021/3/26 matlab_常微分方程数值解法 ppt课件 科学计算与MATLAB 2021/3/26 matlab_常微分方程数值解法 ppt课件 　常微分方程数值解法 2021/3/26 matlab_常微分方程数值解法 ppt课件 精品资料 2021/3/26 matlab_常微分方程数值解法 ppt课件 你怎么称呼老师？ 如果老师最后没有总结一节课的重点的难点，你是否会认为老师的教学方法需要改进？ 你所经历的课堂，是讲座式还是讨论式？ 教师的教鞭 “不怕太阳晒，也不怕那风雨狂，只怕先生骂我笨，没有学问无颜见爹娘 ……” “太阳当空照，花儿对我笑，小鸟说早早早……” 2021/3/26 matlab_常微分方程数值解法 ppt课件 内容提要 引言 欧拉近似方法 龙格-库塔(R-K)方法 MATLAB的常微分方程函数 小结 2021/3/26 matlab_常微分方程数值解法 ppt课件 1、引言 物理学所研究的各