[研究生入学考试]工程数学—第一章.pptVIP

误差分析 误差分析 分部积分 误差分析 泰勒展开 误差分析 误差分析 n In n In 0 0.6321 5 0.1480 1 0.3679 6 0.1120 2 0.2642 7 0.2160 3 0.2074 8 -0.7280 4 0.1704 9 7.552 误差分析 误差分析 n In n In 0 0.6321 5 0.1480 1 0.3679 6 0.1120 2 0.2642 7 0.2160 3 0.2074 8 -0.7280 4 0.1704 9 7.552 积分结果： 误差分析 误差分析 误差分析 误差分析 误差随计算次数增加而增加 误差分析 误差随计算次数增加而减小 误差分析 误差随计算次数增加而减小 误差分析 误差随计算次数增加而减小 误差分析 n In n In 0 0.6321 5 0.1455 1 0.3679 6 0.1268 2 0.2644 7 0.1121 3 0.2073 8 0.1035 4 0.1708 9 0.0684 误差分析 n In n In 0 0.6321 5 0.1455 1 0.3679 6 0.1268 2 0.2644 7 0.1121 3 0.2073 8 0.1035 4 0.1708 9 0.0684 n In n In 0 0.6321 5 0.1480 1 0.3679 6 0.1120 2 0.2642 7 0.2160 3 0.2074 8 -0.7280 4 0.1704 9 7.552 误差分析 重要性：对于算法进行评估 算法的数值稳定性 初始数据的误差在计算中的传递特性 避免误差危害 要避免除数绝对值远远小于被除数绝对值的除法。 要避免两相近数相减 要防止“大数”吃掉小数 减少运算次数 避免误差危害 避免误差危害 避免误差危害 避免误差危害 小结 误差分类：模型、观测、截断、舍入 误差度量：绝对、相对、有效数字 误差传递：一元、多元 数值计算的稳定性：误差递增、递减 避免误差危害：加、减、乘、除、运算量 作业1-1 作业1-2 利用作业1-1的结果，计算要令 的近似值的相对误差限小于0.1%，所需要的 最少有效数字位数，并计算在此有效数字下 的相对误差。 作业1-3 1）F(x)是一元函数，利用泰勒展开式推导函数值F的绝对误差与自变量x的绝对误差之间的关系。 2）利用上述关系，估计 当x0, x*(0)是x的近似值，其相对误差限是ε时， 函数值ln(x*)作为ln(x)的近似值的绝对误差限 作业1-4 作业1-5 工程数学 厦门大学物理与机电工程学院 课程介绍 教材：《高等工程数学》 姚仰新、罗家洪、庄楚强，华南理工大学出版社，2007 参考书目：《数值方法》关治、陆金甫，清华大学出版社，2006 学时：54 学时 成绩：作业（40%）+考勤（10%）+期末（50%） 任课教师：张宇锋 电子邮箱：yufengzhang@xmu.edu.cn 数值分析绪论 数值分析研究对象、内容与特点 数值计算的误差 误差定性分析 避免误差危害 数值分析研究对象 工程和科学研究的一般过程 发现问题 寻找解法 问题获解 解析解 数值解 MatLab数值计算结果 矿道中移动梯子的问题 7 m 9 m 123o 两个矿道以123度角相交，直道宽7米，斜道宽9米，能够转过这个弯的梯子不能超过多长？ 矿道中移动梯子的问题 两个矿道以123度角相交，直道宽7米，斜道宽9米，能够转过这个弯的梯子不能超过多长？ w1 w2 C B A l1 l2 矿道中移动梯子的问题 两个矿道以123度角相交，直道宽7米，斜道宽9米，能够转过这个弯的梯子不能超过多长？ 数值分析研究对象 工程和科学研究中的数值方法 发现问题 问题获解 设计模型 算法设计 程序编写 执行运算 数值分析研究对象 工程和科学研究中的数值方法 发现问题 问题获解 设计模型 算法设计 程序编写 执行运算 求 方程求根 牛顿法 解 数值分析研究内容 函数数值求解 数值微分 数值积分 常微分方程数值解 偏微分方程数值解 线性方程组数值解 …… 应用软件：MatLab, Fortran, Origin … 数值分析研究特点 算筹 算盘 计算机 计算机的应用：加、减、乘、除；逻辑运算；储存量、运算能力有限 算法的时空复杂性，收敛性，数值稳定性 数值计算的误差 数值计算的误差 数值计算的误差 数值解 不等于 解析解 或 精确值 数值解 等于 精确值 + 误差 （合理性） 误差来源 模型误差：在建立数学模型时的 “简单化”导致的误差 误差来源 模型误差：在建立数学模型时的 “简单化”导致的误差 例如：运动中的时钟： 经典力学： 相对论： 误差来源