数值分析二.pptVIP

数值分析 几点要求: 记好课堂笔记 课程安排 第二章 数值计算的基本概念 （2学时） 误差与有效数字 计算机计算的几个问题 算法的稳定性问题 算法设计原则 课程安排 第三章线性方程组求解的数值方法（4学时） 高斯消元法 矩阵分解 向量范数与矩阵范数 迭代法求解 方程组的病态问题与误差分析 课程安排 第四章函数的数值逼近 （5学时） 代数多项式插值问题 线性插值 三次样条插值 曲线拟合的最小二乘法 课程安排 第五章数值积分 （4学时） 插值型求积公式 积分方程的数值求解 课程安排 第六章常微分方程初值问题 （3学时） 欧拉方法 稳定性与收敛性分析 课程安排 第七章非线性方程求根方法 （4学时） 方程求根问题 简单迭代法 牛顿迭代法 课程安排 课程总结和复习 （2学时） 上机实验 （8学时） 课程基本要求 掌握数值方法的基本原理 构造算法的基本思想和技巧 明确如何在计算机上使用 教材和参考书 教材： 《数值计算引论》，白峰杉，高等教育出版社， 参考资料： 《科学计算引论—基于MATLAB的数值分析》, Shoichiro Nakamura,电子工业出版社 《数值分析基础教程》, 李庆杨, 高等教育出版社 联系方式 Email: wangjg@uestc.edu.cn Call: Address: 科研楼B505房 教师社区 电子工程学院 王建国 数值分析有什么用？ 研究工作需要什么? 研究工作需要什么? 前 言 数值分析的任务 数值分析是研究求解各类数学问题的数值方法和有关理论的学科 例： 地球的表面积有多大？ 数学工具 前 言 例：线性方程组求解 第二章 数值计算的基本概念 本章内容 数的表示方法 误差概念和分析 数值计算的算法问题 数值计算应注意的问题 §2.1 数的表示方法 计数 1. 什么叫数？ 通用的记数法 定点形式（定点数) 通用的记数法 浮点形式（浮点数） 通用的记数法 在下面三数中： 二进制记数法 §2.2 误差的概念 绝对误差: 设x是精确值，p是近似值，则定义两者之差是绝对误差: 为近似数p的绝对误差,简称误差 绝对误差限 相对误差 相对误差定义为绝对误差与精确值之比 表示单位长度中含有多少误差 相对误差限 误差来源：舍入误差 截断近似法 舍去低位部分，即取 a= a0 a1 … am . am+1 … am+n ，那么 四舍五入法 四舍五入法 五入情况 当am+n+1 =5,6,7,8,9时，取 a= a0 a1 … am . am+1 … (am+n+1) ,那么 例：多项式值的数值计算 x=0.988:0.0001:1.012; y=x.^7-7*x.^6+21*x.^5-35*x.^4+35*x.^3-21*x.^2+7*x-1; plot(x,y) 有效数字概念 对于 a=a0 a1 … am . am+1 … am+n(a0≠0) 的近似数， 若|Δ|≤0.5x10-n， 则称a为具有m+n+1位有效数字的有效数，其中每一位数字都叫做a的有效数字。有效数和可靠数的最末位数字称为可疑数字 关于有效数字的推论 关于有效数字的推论 截断误差 用数值法求解数学模型时，往往用简单代替复杂,或者用有限过程代替无限过程所引起的误差 解：取Rn =ε/2=0.0005, ? =ε/2 =0.0005 设计算公式为 令Rn ≤1/4n4≤0.0005 则n应满足 n≥(500)4≈4.7 取n=5，设计算每项数值的舍入误差为Δ， 令4Δ≤?=0.0005 Δ≤0.0002,取Δ=0.00005=0.5*10-40.0002 用四位小数计算得 S5=1+0.0313+0.0041+0.0010+0.0003=1.0367 按?= 0.5*10-3，将S5舍入为1.037 §2.3 数值计算问题 数值计算问题的适定性 1. “良态”问题和“病态”问题 在适定的情况下，若对于原始数据很小的变化δX，对应的参数误差δy也很小，则称该数学问题是良态问题；若δy很大，则称为病态问题 病态问题中解对于数据的变化率都很大，因此数据微小变化必将导致参数模型精确解的很大变化 数学问题的性态完全取决于该数学问题本身的属性，在采用数值方法求解之前就存在，与数值方法无关 数值计算问题的适定性 2. 稳定算法和不稳定算法 如果用数值方法计算时,初始误差对结果影响小的算法称为稳定算法。否则称为不稳定算法 §