1 第八章 常微分方程数值解法. 2 1 .微分方程的数值解法 3 在这些节点上把常微分方程的初值问题离散化为差 分方程的相应问题,再求出这些点上的差分方程的解 作为相应的微分方程的近似值(满足精度要求)。

Slides:



Advertisements
Similar presentations
数值分析 第五节 数值微分 在实际问题中,往往会遇到某函数 f(x) 是用表格 表示的, 用通常的导数定义无法求导, 因此要寻求其他 方法近似求导。常用的数值微分方法有 : 一. 运用差商求数值微分 二.运用插值函数求数值微分 三. 运用样条插值函数求数值微分 四. 运用数值积分求数值微分.
Advertisements

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第 8 章 常微分方程 实际中,很多问题的数学模型都是微分方程。我们可以研究它们的一些 性质。但是,只有极少数特殊的方程有解析解。对于绝大部分的微分方程是.
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
新疆医科大学 主讲人:张利萍 计 算 方 法. zlp 第五章 常微分方程数值解 5.1 引言 ( 基本求解公式 ) 5.2 Runge-Kutta 法 5.3 微分方程组和高阶方程解法简介.
第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
第二章 导数与微分 习题课 主要内容 典型例题 测验题. 求 导 法 则求 导 法 则 求 导 法 则求 导 法 则 基本公式 导 数 导 数 微 分微 分 微 分微 分 高阶导数 高阶微分 一、主要内容.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
第七节 函数的微分 一 、微分 概念 二、微分的几何意义 三、 基本初等函数的微分公 式与 微分运算法则 四 、小结.
1 第八章常微分方程初值问题的数值解法. 2 第八章 常微分方程数值解法 8.1 引言 ( 基本求解公式 )8.1 引言 ( 基本求解公式 ) 8.2 Runge-Kutta 法8.2 Runge-Kutta 法 8.3 微分方程组和高阶方程解法简介8.3 微分方程组和高阶方程解法简介.
第 4 章 不定积分 4.1 不定积分的概念与基本积分公式 4.2 换元积分法 4.3 分部积分法.
第 4 章 数值微积分. 4.1 内插求积 Newton-Cotes 公式 第 4 章 数值微积分 4.1 内插求积 Newton-Cotes 公式.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
计算机数学基础(下) --数值分析 教师:孙继荣 电话: 028 -
5.4 微 分 一、微分概念 二、微分的运算法则与公式 三、微分在近似计算上的应用. 引例 一块正方形金属片受热后其边长 x 由 x 0 变到 x 0  x  考查此薄片的面积 A 的改变情况  因为 A  x 2  所以金属片面 积的改变量为  A  (x 0 
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二章 导数与微分 一. 内 容 要 点 二. 重 点 难 点 三. 主 要 内 容 四. 例 题与习题.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
2.3 函数的微分. 四川财经职业学院 课前复习 高阶导数的定义和计算方法。 作业解析:
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
第三章 函数逼近 — 最佳平方逼近.
第一章 行列式 第五节 Cramer定理 设含有n 个未知量的n个方程构成的线性方程组为 (Ⅰ) 由未知数的系数组成的n阶行列式
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
一、原函数与不定积分 二、不定积分的几何意义 三、基本积分公式及积分法则 四、牛顿—莱布尼兹公式 五、小结
第二节 微积分基本公式 1、问题的提出 2、积分上限函数及其导数 3、牛顿—莱布尼茨公式 4、小结.
第九章 常微分方程数值解  考虑一阶常微分方程的初值问题
第四章 定积分及其应用 4.3 定积分的概念与性质 微积分基本公式 定积分的换元积分法与分部积分法 4.5 广义积分
第四章 函数的积分学 第六节 微积分的基本公式 一、变上限定积分 二、微积分的基本公式.
微积分基本定理 2017/9/9.
9.1 数值积分基本方法 9.2 梯形积分 9.3 Simpson积分 9.4 Newton-Cotes积分 9.5 Romberg积分
§5.3 定积分的换元法 和分部积分法 一、 定积分的换元法 二、 定积分的分部积分法 三、 小结、作业.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
定积分习题课.
不确定度的传递与合成 间接测量结果不确定度的评估
Ordinary Differential Equations
计算机数学基础(下) 第5编 数值分析 第14章 常微分方程的数值解法.
第四节 一阶线性微分方程 线性微分方程 伯努利方程 小结、作业 1/17.
第三节 格林公式及其应用(2) 一、曲线积分与路径无关的定义 二、曲线积分与路径无关的条件 三、二元函数的全微分的求积 四、小结.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
第二章 导数与微分 第二节 函数的微分法 一、导数的四则运算 二、复合函数的微分法.
非线性物理——混沌 向前 胡冰清
全 微 分 欧阳顺湘 北京师范大学珠海分校
第三章 导数与微分 习 题 课 主要内容 典型例题.
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
§3 微分及其运算 一、微分的定义 二、基本初等函数的微分公式与 微分运算法则.
数值计算方法 第八章 常微分方程初值问题数值解法  重庆邮电大学.
/* Numerical Methods for Ordinary Differential Equations */
计算机数学基础 主讲老师: 邓辉文.
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
第三模块 函数的微分学 第五节 隐函数及参数方程的求 导方法、高阶导数 一、隐函数的微分法 二、由参数方程所确定的函数的微分法
第三节 泰勒 ( Taylor )公式 — 应用 一、泰勒公式的建立 二、几个初等函数的麦克劳林公式 三、泰勒公式的应用 第三章 理论分析
高等数学 西华大学应用数学系朱雯.
Partial Differential Equations §2 Separation of variables
概 率 统 计 主讲教师 叶宏 山东大学数学院.
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
3.1 变化率与导数   3.1.1 变化率问题 3.1.2 导数的概念.
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
学习任务三 偏导数 结合一元函数的导数学习二元函数的偏导数是非常有用的. 要求了解二元函数的偏导数的定义, 掌握二元函数偏导数的计算.
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
第三节 函数的微分 3.1 微分的概念 3.2 微分的计算 3.3 微分的应用.
选修1—1 导数的运算与几何意义 高碑店三中 张志华.
第四章 函数的 积分学 第七节 定积分的换元积分法     与分部积分法 一、定积分的换元积分法 二、定积分的分部积分法.
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
《偏微分方程》第一章 绪论 第一章 绪论 1.1.
第一节 不定积分的概念与性质 原函数与不定积分的概念 基本积分表 不定积分的性质 小结、作业 1/22.
Presentation transcript:

1 第八章 常微分方程数值解法

2 1 .微分方程的数值解法

3 在这些节点上把常微分方程的初值问题离散化为差 分方程的相应问题,再求出这些点上的差分方程的解 作为相应的微分方程的近似值(满足精度要求)。

4 2 .欧拉法与改进欧拉法

5

6

7

在 x 0 处:以( x 0,y 0 )为点, f(x 0,y 0 ) 为斜率作直线 在 x 1 处:以( x 1,y 1 )为点, f(x 1,y 1 ) 为斜率作直线

9

10

11

12

13

14

15

16

17

18

20 考虑用函数 f ( x, y ) 在若干点上的函数值的线性组合来构 造近似公式,构造时要求近似公式在 ( x i, y i ) 处的 Taylor 展 开式与解 y ( x ) 在 x i 处的 Taylor 展开式的前面几项重合,从 而使近似公式达到所需要的阶数。既避免求偏导, 又提高 了计算方法精度的阶数。或者说, 在 [x i,x i+1 ] 这一步内多预 报几个点的斜率值,然后将其加权平均作为平均斜率, 则可构造出更高精度的计算格式,这就是龙格 — 库塔 ( Runge-Kutta )法的基本思想。

21 1 二阶龙格 — 库塔法 在 上取两点 x i 和, 以该两点处的斜率值 k 1 和 k 2 的加权平均 ( 或称为线性组合 ) 来求取平均斜率 k * 的近似值 K ,即 式中 : k 1 为 x i 点处的切线斜率值, k 2 为 点处的切线斜率值, 对照改进的欧拉法, 将 视为 ,即可得

22 对常微分方程初值问题的解 y = y ( x ), 根据微分中值定理, 存在点 ,使得 式中 K 可看作是 y = y ( x ) 在区间 上的平均斜率。所以 可得计算公式为: 也即

23 将 y ( x i ) 在 x = x i 处进行二阶 Taylor 展开: 将 可得: 将 f(x,y) 在 ( x i, y i ) 处进行一阶 Taylor 展开:

24 成立, 格式 的局部截断误差就等于 进行比较系数后可知, 只要 有 2 阶 精度

25 上式中具有三个未知量, 但只有两个方程, 因而有无穷多 解。若取, 则 p =1 ,这是无穷多解中的一个解,改 写可得 不难发现,上面的格式就是改进的欧拉格式。凡满足 条件上式有一簇形如上式的计算格式,这些格式统称为 二阶龙格 — 库塔格式。因此改进的欧拉格式是众多的二 阶龙格 — 库塔法中的一种特殊格式。

26 2 三阶龙格 - 库塔法 为了进一步提高精度,设除 外再增加一点 并用三个点,, 的斜率 k 1, k 2, k 3 加权平均 得出平均斜率 k * 的近似值,这时计算格式具有形式 :

27 为了预报点 的斜率值 k 3, 在区间 内有两 个斜率值 k 1 和 k 2 可以用, 可将 k 1, k 2 加权平均得出 上的平均斜率, 从而得到 的预报值 于是可得 运用 Taylor 展开方法选择参数, 可以使上述格 式的局部截断误差为, 即具有三阶精度,这类格式 统称为三阶龙格 — 库塔方法。

28 三阶龙格 — 库塔方法

29 3 四阶龙格 — 库塔法 如果需要再提高精度,用类似上述的处理方法, 只需在区间 上用四个点处的斜率加权平均作为 平均斜率 k * 的近似值,构成一系列四阶龙格 — 库塔公 式。具有四阶精度,即局部截断误差是 。 由于推导复杂,这里从略,只介绍最常用的一种四 阶经典龙格 — 库塔公式。

30 四阶龙格 — 库塔法

31 或

32 四阶龙格 — 库塔法的优点是 : ⑴ 它是一步法,即已知 y k 就可以求出 y k+1 ,不需 要知道其它的数据。 ⑵ 精确度高。 缺点是 : 计算量大。在一个步长的计算中,要四次计算 f(x,y) 的值。

33 例 4 用四阶龙格 — 库塔法解 在 处的近似值。 解

34 先用公式计算出 ,再计算出 ,然后步长增加, 计算下一步长。 计算结果为:

35

36

37

38

39

40

41

42

43

44

45

46