计算方法(B) 主讲:张明波 email: mbzhang@ustc.edu.cn Tel: 3600365 (O), 13721101665.

Slides:



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

一、微分的定义 二、微分的几何意义 三、微分公式及微分法则 四、微分在近似计算中的应用 五、小结 思考题.
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量. 再例如, 既容易计算又是较好的近似值 问题 : 这个线性函数 ( 改变量的主要部分 ) 是否 所有函数的改变量都有 ? 它是什么 ? 如何求 ?
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
一、问题提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、 微分的求解 六、 微分的应用 七、 小结.
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
常系数线性微分方程组 §5.3 常系数线性方程组. 常系数线性微分方程组 一阶常系数线性微分方程组 : 本节主要讨论 (5.33) 的基解矩阵的求法.
第 4 章 数值微积分. 4.1 内插求积 Newton-Cotes 公式 第 4 章 数值微积分 4.1 内插求积 Newton-Cotes 公式.
计算机数学基础(下) --数值分析 教师:孙继荣 电话: 028 -
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
2.3 函数的微分. 四川财经职业学院 课前复习 高阶导数的定义和计算方法。 作业解析:
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
软件学院11级实习前培训-论文和学位申请 任皖英 Tel: (办)
计算方法教程 数学与统计学院 马军 理科楼 338 QQ :
1.2 误差 误差的来源与分类 误差与有效数字 函数求值的误差估计
1.非线性振动和线性振动的根本区别 §4-2 一维非线性振动及其微分方程的近似解法 方程
第一讲 数值计算的误差.
数 值 分 析 (第 4 版) 李庆扬 王能超 易大义 编 清华大学出版社 施普林格出版社.
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
第三章 函数逼近 — 最佳平方逼近.
一、利用导数作近似计算 1. 近似计算 是用计算方法得到一定精度的计算结果. y 于是 o x.
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
9.1 数值积分基本方法 9.2 梯形积分 9.3 Simpson积分 9.4 Newton-Cotes积分 9.5 Romberg积分
§5.3 定积分的换元法 和分部积分法 一、 定积分的换元法 二、 定积分的分部积分法 三、 小结、作业.
第四章 数值积分与数值微分 — 基本概念 — Newton-Cotes 公式.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
陈研 Tel: 新学科综合楼 数 值 计 算 方 法 陈研 Tel: 新学科综合楼 中国农业大学资源和环境学院 2005年9月.
第五节 函数的微分 一、微分的概念 二、微分运算法则 三、微分在近似计算中的应用 四、微分在估计误差中的应用 第二章
不确定度的传递与合成 间接测量结果不确定度的评估
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
§4.3 常系数线性方程组.
Examples for transfer function
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
§3 微分及其运算 一、微分的定义 二、基本初等函数的微分公式与 微分运算法则.
第2章 Z变换 Z变换的定义与收敛域 Z反变换 系统的稳定性和H(z) 系统函数.
Computer Graphics 计算机图形学基础 张 赐 Mail: CSDN博客地址:
第四章 数值积分与数值微分 — 复合求积公式 — Romberg 算法.
Tel: : 授课: 68 学分:4.
计算机数学基础 主讲老师: 邓辉文.
第三节 泰勒 ( Taylor )公式 — 应用 一、泰勒公式的建立 二、几个初等函数的麦克劳林公式 三、泰勒公式的应用 第三章 理论分析
第三章 线性代数方程组的解法 在自然科学和工程技术中很多问题的解决常常归结为解线性代数方程组。例如: (蓝色)建筑工程中的结构力学问题;
第5章 线性代数 矩阵分析 矩阵分解 线性方程组的求解 符号矩阵.
第二章 函数 插值 — 分段低次插值.
课题:1.5 同底数幂的除法.
主讲:张瑞 Tel: (O) 计算方法(B) 主讲:张瑞 Tel: (O)
概 率 统 计 主讲教师 叶宏 山东大学数学院.
成绩是怎么算出来的? 16级第一学期半期考试成绩 班级 姓名 语文 数学 英语 政治 历史 地理 物理 化学 生物 总分 1 张三1 115
第16讲 相似矩阵与方阵的对角化 主要内容: 1.相似矩阵 2. 方阵的对角化.
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
第4课时 绝对值.
学习任务三 偏导数 结合一元函数的导数学习二元函数的偏导数是非常有用的. 要求了解二元函数的偏导数的定义, 掌握二元函数偏导数的计算.
第一部分:概率 产生随机样本:对分布采样 均匀分布 其他分布 伪随机数 很多统计软件包中都有此工具 如在Matlab中:rand
海报题目 简介: 介绍此项仿真工作的目标和需要解决的问题。 可以添加合适的图片。
建模常见问题MATLAB求解  .
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
2019/5/20 第三节 高阶导数 1.
§2 方阵的特征值与特征向量.
§7.3 离散时间系统的数学 模型—差分方程 线性时不变离散系统 由微分方程导出差分方程 由系统框图写差分方程 差分方程的特点.
第三节 函数的微分 3.1 微分的概念 3.2 微分的计算 3.3 微分的应用.
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
Volterra-Lotka方程 1925年, A. Lotka(美)和V. Volterra(意)给出了第一个两物种间的捕食模型。
   小数的产生和意义      小站三小        刘宝霞.
海报题目 简介: 介绍此项仿真工作的目标和需要解决的问题。 可以添加合适的图片。
國立臺灣師範大學附屬高級中學 102學年度第一學期 家長日 校務方針報告
Presentation transcript:

计算方法(B) 主讲:张明波 email: mbzhang@ustc.edu.cn Tel: 3600365 (O), 13721101665

第0章 绪论 计算方法干什么? 计算方法的内容 误差

计算方法是一种研究并解决数学问题的数值近似解方法 实际问题 重点为:设计算法,并分析算法的稳定性﹑收敛性 ﹑复杂度﹑误差等 物理模型 数学模型 计算方法是一种研究并解决数学问题的数值近似解方法 计算方法 计算机求结果

内容 1、数值逼近-插值﹑拟合;数值微分与数值积分. 2、数值线性代数-解线性方程组、求逆矩阵、特征值、 特征向量等. 2、数值线性代数-解线性方程组、求逆矩阵、特征值、 特征向量等. 100亿/秒,算3,000年,而Gauss消元法2660次 3、(非线性﹑微分)方程数值求解.

学习的目的、要求 会应用、修改、创建公式 编制程序完成计算 课程评分方法  总分 (100) = 平时作业(20)+上机作业(10)+期末 (70)

上机作业要求 1、编程可以用任何语言C,C++,Matlab,Mathematica,Maple,等 不允许使用内置函数完成主要功能 2、以email形式发给助教 内容:一次作业一个附件,并在内容中写出运行结果

误差 绝对误差 设 为准确值, 为近似值, 称为近似值 的绝对误差;若| | ,则称 为 的一个绝对误差限. ,则称 为 的一个绝对误差限. 例:3.14作为 的近似值,其绝对误差为: 而 0.002 为 3.14 的一个绝对误差限.

相对误差 为准确值, 为近似值 设 称为 的相对误差. 若| | ,则称 为 的一个相对误差限. 称为 的相对误差. 若| | ,则称 为 的一个相对误差限. 例:将10,000米的跑道建成10,010与 100米的跑道建成101, 两者的绝对误差分别为10和1米,优劣如何?

前者相对误差 (10010-10000)/10000=0.001, 后者相对误差(101-100)/100=0.01. 故虽前者绝对误差较大,但在某种意义上前者更精确。 绝对误差与相对误差是描述误差大小的两种方式,都很重要.

函数的误差(误差的传播) 为准确值, 为近似值 , 函数值f(x)的误差为? 绝对误差限 相对误差限

多变元函数的误差 为准确值, 为近似值 , 函数值 的误差为? 其中 为 的绝对误差, 在 处取值.

故绝对误差限

相对误差限为

有效位的多少直接影响到近似值的绝对误差和相对误差 有效位数 当x的绝对误差不超过某一位的半个单位,则这一位到左数第一个非零位的位数称位x的有效位数。 有效位的多少直接影响到近似值的绝对误差和相对误差 例:3.14作为 的近似值有3位有效数字; 0.031或0.030作为0.0305的近似值均有两位有效 数字

计算中注意事项 分别令 和 知: 当 时,相对误差会较大; 当 时,绝对误差会较大. 算法中应避免相近两数相减,避免小数作分母.

减少运算次数 例:1. 2. 秦九昭算法 避免大数吃小数(必要时提高精度) 例:计算 设计稳定收敛的算法.

例子 要计算 ,我们有 1. 构造方法如下: 2.

n 0.182 1 0.088 0.090 2 0.058 0.050 3 0.0431 0.083 4 0.0343 -0.165 5 0.0284 1.025 6 0.024 -4.958 7 0.021 24.933 8 0.019 -124.540

原因:对格式1,如果前一步有误差, 则被放大5倍加到这一步 称为不稳定格式 稳定格式,对舍入误差有抑制作用

2、 有时候,模型本身就是病态 (系数引入小变化,解产生大变化)