第一章引论素材.pptVIP

数 值 分 析 欧志英 兰州理工大学 教材和参考书 教材： 数值分析基础，欧志英等编，甘肃人民出版社 参考书： 数值分析（第四版），李庆扬，王能超，易大义编，清华大学出版社； 计算方法（第二版），邓建中，刘之行，西安交通大学出版社； 数值分析与实验学习指导，蔡大用编，清华大学出版社。 课程内容 引论（误差、有效数字等） 非线性方程求根 线性方程组的直接解法 线性方程组的迭代解法 数值插值方法 数据拟合方法 数值积分与数值微分 常微分方程的数值解法 说明： 学时： 授课48 成绩评定： 期末考试（80%） 平时成绩（20%） 第一章 引论 第一章 引论 1.1 数值分析研究的对象与特点 1.2 数值计算的误差与有效数字 1.3 数值运算的误差估计 1.4 数值计算中的一些基本原则 1.1 数值分析研究的对象与特点 数值分析是做什么用的？ 研究对象：用计算机求解各种数学问题的数值计算方法及其理论与软件实现。 特点: 面向计算机，提供切实可行的算法； 有可靠的理论分析，能达到精度要求，保证近似算法的收敛性和数值稳定性； 要有好的计算复杂性，节省时间及存储量； 有数值实验，证明算法有效。 常采用的处理方法： 构造性方法； 离散化方法； 递推化方法； 迭代法； 近似替代方法； 以直代曲法； 化整为零的处理方法； 外推法。 微积分的若干定理： 罗尔定理和微分中值定理； 介值定理及推论； 泰勒公式（一元、二元）； 积分中值定理； 1.2 数值计算的误差与有效数字 例：蝴蝶效应 —— 纽约的一只蝴蝶翅膀一拍，风和日丽的北京就刮起台风来了？！ 营长对值班军官: 明晚大约 8点钟左右，哈雷彗星将可能在这个地区看到，这种彗星每隔 76年才能看见一次。命令所有士兵着野战服在操场上集合，我将向他们解释这一罕见的现象。如果下雨的话，就在礼堂集合，我为他们放一部有关彗星的影片。 1.4 数值运算中的一些基本原则 1. 避免除数的绝对值远小于被除数的绝对值； 2. 避免两个相近的数相减； 3. 防止大数“吃掉”小数现象； 4. 简化计算步骤，减少运算次数； 5. 选用数值稳定性好的算法。 1.3 数值运算的误差估计 1. 函数运算的误差估计： 设y=f(x)为一元函数，自变量准确值x*，对应函数准确值y*=f(x*)，x误差为e(x)，误差限为ε(x)，函数近似值误差e(y)，误差限为ε(y)。则（可由Taylor公式推得） 对于多元函数 设准确值 由多元函数Taylor公式，可得误差估计： 相对误差限为： 2. 算术运算的误差估计： 两个近似数x1, x2，其误差限分别为ε(x1)， ε(x2)，它们进行加、减、乘、除运算得到的误差限分别为： 例4: 设a=2.31,b=1.93,c=2.24都是三位有效数字的近似数，令p=a+bc，求ε(p)和 εr(p)，并判断p有几位有效数字。 解 由题知， ε(a) =ε(b) =ε(c)=0.005 ε(p)= ε(a)+ ε(bc)? ε(a)+|b| ε(c)+|c| ε(b) =0.005+1.93?0.005+2.24 ?0.005=0.02585 p=a+bc=2.31+1.93 ?2.24=6.6332 又 故 εr(p)=ε(p)/|p| ? 0.02585/6.6332 ? 0.0039=0.39% 因为 ε(p) ? 0.025850.05=1/2 ?101-2 所以p=6.6332中只有两位有效数字。 1. 避免除数的绝对值远小于被除数的绝对值； 设 若|x||y|， 因 表明当|x|相对太小时，商的绝对误差可能很大。 2. 避免两个相近的数相减； 若y ? x，设z=y-x，则ε(z)= ε(y)+ ε(x)， 当y ? x时，z ? 0，相对误差限会很大，导致结果有效数字位数减少。 避免两个相近的数相减方法：常用一些恒等变形实现。 如当x充分大时， 当x绝对值很小时， 3. 防止大数“吃掉”小数现象； 原因：计算机表示的数位数有限，很大的数和很小的数相加减时，很小的数会被“吃掉”（舍去）。 方法：当绝对值悬殊的一系列数相加时，若有 |x1||x2|…|xn|，应按绝对值由小到大的顺序累加。 例 计算方程 的根。 算法一： 利用求根公式 在计算机内， 109存为0.1? 1010 ，1存为0.1? 101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。 * * E-mail：zhiyingou@163.com 数值 分析 输入复杂问题或运算 ? ? ? ? 计算机 近似解 计算机解决科学计算的过程：