微分方程数值解第一章答案课件.ppt

《微分方程数值解法》 参考书: 微分方程数值解法 李荣华等编, 高教出版社 课堂授课+计算实验 考核方式: 平时作业+课堂+期末考试 任课教师 教学内容 第一章、常微分方程的数值解法 第二章、椭圆型方程的差分方法 第七章、椭圆型方程的有限元方法 第四章、抛物型方程的差分方法 第五章、双曲型方程的差分格式 第一章 常微分方程初值问题的数值解法 教学目标 教学重点 教学过程 教学目标 了解ODE数值解法的基本内容, 掌握Euler法和线性多步方法, 会判断常用方法的优劣之处. 教学重点 基本概念和Euler法 线性多步方法 稳定性 教学过程 常微分方程基本概念 常微分方程初值问题 Euler法及其基本问题 线性多步方法 数值稳定性 Runge-Kutta方法 1: 常微分方程的基本概念 2 初值问题：标量形式 常微分方程来源举例1 2 Euler方法 举例1 利用Euler方法计算初值问题 举例2 P55 习题1 利用Euler方法求数值解 Euler方法的三种解释 数值微分：用差商来代替导数 数值积分：把微分方程变成积分方程 幂级数展开：将u(t+h) 在t 做Taylor展开 数值方法的基本问题 截断误差(局部、整体) 相容性 收敛性 稳定性 局部截断误差 设u(t)是初值问题(1)的解, 在[t,t+h]上定义算子 整体截断误差 设u(t)是初值问题(1)的解, un是(2)的解,定义算子 截断误差 局部截断误差Rn：假设第n步精确计算的前提下，计算解un+1和精确解u(tn+1)的误差 整体截断误差?n：在考虑误差累积的效应下，计算解un和精确解u(tn)的误差 相容性和相容的阶 相容性针对的是建立差分格式时由差商代替微商所引起的局部截断误差. q阶相容: 若一个离散变量方法的局部截断误差对任意n满足： 稳定性 在利用公式(2)计算数值解的过程中,难免有舍入误差.稳定性就是讨论舍入误差是否会随着计算无限扩大地传递下去. 数值方法稳定性指对初始误差的连续依赖性,以线性k步方法为例，即为存在常数C和h00,使得当h∈(0,h0] 时 这里常数C不依赖于h。通常这里定义的稳定性指 h→0 情况下的稳定性。 Eular方法的性质 相容性 √ (1阶) 收敛性 √ (1阶) 稳定性 √ 绝对稳定区域 数值求解微分方程过程示意 作业 1 利用Euler方法求数值解 收敛性与收敛的阶 收敛性研究的是误差累积产生的整体截断误差. 收敛：对任意的t∈(t0,T] ，成立 若此时，整体截断误差满足 则称方法的收敛为p阶的 总结：基本步骤 ③ 解差分方程，求出格点函数 ① 对区间作分割： 求 y(x) 在xi 上的近似值yi。 ② 由微分方程出发，建立求格点函数的差分方程。 这个方程应该满足： A、解存在唯一；B、稳定，收敛；C、相容 目的 关键 为了考察数值方法提供的数值解，是否有实用价值， 需要知道如下几个结论： ① 步长充分小时，所得到的数值解能否逼近 问题得真解；即收敛性问题 ② 误差估计 产生得舍入误差，在以后得各步计算中，是否会 无限制扩大；稳定性问题 微分方程 区域剖分 离散系统的 性态研究 递推计算或解线性代数方程组 微分方程离散 初始和边界条件处理 解的存在性、唯一性 解的收敛性和收敛速度 解的稳定性 得到数值解 步长h=0.1, 解区间[0,0.5] 绘制折线，与真解比较 2 证明梯形法的收敛性, 并估计整体截断误差. * * * * 北京·中国地质大学China University of Geosciences，Beijing 教材： 微分方程数值方法 (第二版)， 胡健伟,汤怀民著, 科学出版社, 2007,2 参考书: 微分方程数值解法 李荣华等编, 高教出版社 第一章 基本概念 第一章 基本概念 第一章 基本概念 微分方程: 常微分方程和偏微分方程 阶 解，通解和特解 定解问题: 初值问题和边值问题 常微分方程 偏微分方程 联系着自变量, 未知函数及其导数(微分)的方程, 称为微分方程 . :未知函数是一元函数 分类 微分方程: 常微分方程和偏微分方程 :未知函数是多元函数 方程中未知函数导数的最高阶数叫做微分方程的阶. 一阶微分方程 三阶微分方程 一阶微分方程 例如： 微分方程的阶 — 是使方程成为恒等式的函数. 通解 — 解中所含独立的任意常数的个数与方程 的阶数相同. 特解 微分方程的解 — 不含任意常数的解. （微分方程的绝大部分解） 解, 通解, 特解 — 确定通解中任意常数的条件. 1