新疆医科大学 主讲人:张利萍 计 算 方 法. zlp 第五章 常微分方程数值解 5.1 引言 ( 基本求解公式 ) 5.2 Runge-Kutta 法 5.3 微分方程组和高阶方程解法简介.

Slides:



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

第六章 数值微分 6.1 插值型数值微分公式 6.2 插值型数值积分. 6.1 插值型数值微分公式 当 x 为插值节点 时,上式简化为 故一般限于对节点上的导数值采用插值多项式的相应导数 值进行近似计算,以便估计误差。 一般地 这类公式称为插值型数值微分公式。
数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第 8 章 常微分方程 实际中,很多问题的数学模型都是微分方程。我们可以研究它们的一些 性质。但是,只有极少数特殊的方程有解析解。对于绝大部分的微分方程是.
高等数学( XJD ) 第二章 导数与微分 返回 高等数学( XAUAT ) 高等数学( XJD ) 求导法则 基本公式 导 数 导 数 微 分微 分 微 分微 分 求导方法 高阶导数 微分法则 导数与微分关系图导数与微分关系图.
一、微分的定义 二、微分的几何意义 三、微分公式及微分法则 四、微分在近似计算中的应用 五、小结 思考题.
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
第五节 全微分方程 一、全微分方程及其求法 二、积分因子法 三、一阶微分方程小结. 例如 所以是全微分方程. 定义 : 则 若有全微分形式 一、全微分方程及其求法.
第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
第二章 导数与微分 习题课 主要内容 典型例题 测验题. 求 导 法 则求 导 法 则 求 导 法 则求 导 法 则 基本公式 导 数 导 数 微 分微 分 微 分微 分 高阶导数 高阶微分 一、主要内容.
目录 上页 下页 返回 结束 习题课 一、导数和微分的概念及应用 二、导数和微分的求法 导数与微分 第二章.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
1 第八章 常微分方程数值解法. 2 1 .微分方程的数值解法 3 在这些节点上把常微分方程的初值问题离散化为差 分方程的相应问题,再求出这些点上的差分方程的解 作为相应的微分方程的近似值(满足精度要求)。
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
常系数线性微分方程组 §5.3 常系数线性方程组. 常系数线性微分方程组 一阶常系数线性微分方程组 : 本节主要讨论 (5.33) 的基解矩阵的求法.
第七节 函数的微分 一 、微分 概念 二、微分的几何意义 三、 基本初等函数的微分公 式与 微分运算法则 四 、小结.
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 -
1 热烈欢迎各位朋友使用该课件! 广州大学数学与信息科学学院. 2 工科高等数学 广州大学袁文俊、邓小成、尚亚东.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
理学院 张立杰 《数值分析》第四讲 数值积分与微分. §4.1 引言 第四章:数值积分与数值微分 1 、积分的概念 设 任取 做 如果 存在, 则称 可积,极限值称为函数 在区间 [a,b] 上的 定积分,记为 : Riemann 积分.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
2.3 函数的微分. 四川财经职业学院 课前复习 高阶导数的定义和计算方法。 作业解析:
圆的一般方程 (x-a)2 +(y-b)2=r2 x2+y2+Dx+Ey+F=0 Ax2+Bxy+Cy2+Dx+Ey+ F=0.
18.2一元二次方程的解法 (公式法).
6.9二元一次方程组的解法(2) 加减消元法 上虹中学 陶家骏.
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
第三章 函数逼近 — 最佳平方逼近.
《高等数学》(理学) 常数项级数的概念 袁安锋
第二章 数值微分和数值积分.
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
高等数学电子教案 第五章 定积分 第三节 微积分基本定理.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
一、原函数与不定积分 二、不定积分的几何意义 三、基本积分公式及积分法则 四、牛顿—莱布尼兹公式 五、小结
第二节 微积分基本公式 1、问题的提出 2、积分上限函数及其导数 3、牛顿—莱布尼茨公式 4、小结.
第九章 常微分方程数值解  考虑一阶常微分方程的初值问题
9.1 数值积分基本方法 9.2 梯形积分 9.3 Simpson积分 9.4 Newton-Cotes积分 9.5 Romberg积分
定积分的换元法 和分部积分法 换元公式 分部积分公式 小结 1/24.
§5.3 定积分的换元法 和分部积分法 一、 定积分的换元法 二、 定积分的分部积分法 三、 小结、作业.
第四章 数值积分与数值微分 — 基本概念 — Newton-Cotes 公式.
计算方法 第2章 数值微分与数值积分 2.1 数值微分.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
第三节 函数的求导法则 一 函数的四则运算的微分法则 二 反函数的微分法则 三 复合函数的微分法则及微分 形式不变性 四 微分法小结.
计算机数学基础(下) 第5编 数值分析 第14章 常微分方程的数值解法.
第四节 一阶线性微分方程 线性微分方程 伯努利方程 小结、作业 1/17.
第三节 格林公式及其应用(2) 一、曲线积分与路径无关的定义 二、曲线积分与路径无关的条件 三、二元函数的全微分的求积 四、小结.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
第三章 导数与微分 习 题 课 主要内容 典型例题.
Examples for transfer function
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
穩定是指偏離平衡時能夠回復平衡的特性,控制則是改變飛行狀態的機制。
数值计算方法 第八章 常微分方程初值问题数值解法  重庆邮电大学.
第四章 数值积分与数值微分 — 复合求积公式 — Romberg 算法.
计算机数学基础 主讲老师: 邓辉文.
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
高等数学 西华大学应用数学系朱雯.
Partial Differential Equations §2 Separation of variables
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
一元二次不等式解法(1).
第四节 第七章 一阶线性微分方程 一、一阶线性微分方程 *二、伯努利方程.
第六章 数值积分与数值微分.
§2 方阵的特征值与特征向量.
第三节 函数的微分 3.1 微分的概念 3.2 微分的计算 3.3 微分的应用.
我们能够了解数学在现实生活中的用途非常广泛
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
第一节 不定积分的概念与性质 原函数与不定积分的概念 基本积分表 不定积分的性质 小结、作业 1/22.
Presentation transcript:

新疆医科大学 主讲人:张利萍 计 算 方 法

zlp 第五章 常微分方程数值解

5.1 引言 ( 基本求解公式 ) 5.2 Runge-Kutta 法 5.3 微分方程组和高阶方程解法简介

本章要点: 本章作业 本章主要研究基于微积分数值解法的常 微分方程数值解,主要方法有 线性单步法中的 Euler 方法、 Simpson 方法、 Runge-Kutta 方法 高阶微分方程和微分方程组的数值解法 P

本章应用题 : 驱逐舰在浓雾中搜索潜艇,其时发现潜艇在 3 英里 的海面上,但潜艇立即下潜,驱逐舰速度两倍于潜 艇,且已知潜艇下潜后即以全速朝某一未知方向直 线前进,问驱逐舰应采取什么路线才能保证它会开 过潜艇的上方以投放深水炸弹? 提示 取极坐标,并以发现潜艇时潜艇的位置为原点 ——— 反潜

5.1 引言 ( 基本求解公式 ) 在工程和科学技术的实际问题中,常需要求解微分方程 只有简单的和典型的微分方程可以求出解析解 而在实际问题中的微分方程往往无法求出解析解 在高等数学中我们见过以下常微分方程: (1) (2)

(3) (1),(2) 式称为初值问题,(3) 式称为边值问题 (4) 另外, 在实际应用中还经常需要求解常微分方程组 : 本课程主要研究问题 (1) 的数值解法, 对 (2)~(4) 只作简单介绍 我们首先介绍初值问题 (1) 的解存在的条件

定理 1. 对于问题 (1), 要求它的数值解

(1) 从 (1) 的表达式 可以看出, 求它的数值解的关键在于 而数值微分或数值积分问题我们都已经学习过

一、基于数值微分的常微分方程数值解法 (1) 对于初值问题 (1) 在下列子区间上分别应用两点数值微分公式 为了讨论方便, 假设以下节点为等距节点

(5) ( 一 ) Euler 公式

由 (5) 式每组的前一半可得 (6) (7) 记 其中 (6) 和 (7) 式称为求解初值问题 (1) 的 ( 前进 )Euler 公式和误差项

由 (5) 式每组的后一半可得 记 其中 (8) (9) (8) 和 (9) 式称为求解初值问题 (1) 的后退 Euler 公式和误差项

从 (6) 或 (8) 式不难看出, 这种类型的方法称为单步格式或单步法 Euler 方法的几何体现 : 前进 Euler 公式 后退 Euler 公式

Euler1.m 例 1. 解:解: 由前进 Euler 公式

得 依此类推, 有

由于后退 Euler 公式是隐形公式, 计算例 1 将很麻烦 事实上大多数情况下用后退 Euler 公式都较困难 就可得到新的 Euler 公式 (10) 此方法称为预测 — 校正系统

用 Euler 公式的预测 —— 校正系统求解例 1.例 1 例 2. 解:解: 由 (10) 式, 有 Euler1.m

依此类推, 得 比较不同的结果

( 二 ) 常微分方程数值解的截断误差 评价一个微分方程求解公式的标准当然是其精度 而在求解公式 中 误差项

定义 1. 因为一般情况下, 求解公式的每一步都存在误差, 因此有 定义 2. 定义 3.

Euler 公式的局部截断误差为 具有 1 阶精度 后退 Euler 公式的局部截断误差为 也具有 1 阶精度 显然一个求解公式的精度越高, 计算解的精确性也就越好 从前面的分析可知,Euler 法的精度并不算高 因此有必要找寻精度更高的求解公式

二、基于数值积分的常微分方程数值解法 (1) 对于初值问题 (11)

矩形求积公式 梯形求积公式, 误差为 Simpson 求积公式, 误差为 将以上求积公式代入 (11) 式, 并加以处理就 可得到相对应的求解公式

( 一 ) 矩形求解公式 由 可得 令 (12) (12) 式称为矩形公式 ( 矩形法 ) 实际上就是 Euler 求解公式

( 二 ) 梯形求解公式 由 可得 令 (13) 称 (13) 式为梯形求解公式 ( 梯形法 ) 注意 :(13) 式是隐形公式

则梯形公式第 k 步的截断误差为 显然梯形法具有二阶精度 由于梯形公式为隐形公式, 一般情况下不易显化

------(14) 以上公式称为改进的 Euler 求解公式 ( 改进 Euler 法 ), 即 (15)

例 3. 用 Euler 公式、梯形公式和改进 Euler 公式求 解初值问题, 并比较结果的精度 解:解: (1)Euler 公式

(2) 梯形公式

(3) 改进 Euler 公式 改进 Euler 公式 x y 使用 MATLAB 软件 Euler2.m 结果为

Euler 公式 梯形公式 改进 Euler 公式 结果比较 Euler 法的精度不如梯形公式

( 三 ) Simpson 求解公式 将 Simpson 求积公式 代入( 11 ) 简化后, 得 (16)

由 Simpson 求积公式的误差 可以近似得到 (16) 式的截断误差为 仔细分析 (16) 式 如何求 ? 求积公式 (16) 的精度为 4 阶

考虑将 (16) 式改为下面的形式 (17) (18) (17) 式称为 Simpson 求解公式,(18) 为相应的截断误差项 (17) 式是一个隐式求解公式 这种形式称为多步法, 在本课程中将不予介绍

张利萍制作 39 小结 这节课我们主要学习了: 1 、基于数值微分的常微分方程数值解法 2 、基于数值积分的常微分方程数值解法 要求熟练掌握欧拉公式、改进的欧拉公 式、校正欧拉公式 作业:课后相应习题