数值分析全册完整课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《数值分析全册完整课件.ppt》由用户(金钥匙文档)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 完整 课件
- 资源描述:
-
1、 计计 算算 方方 法法 数值分析全册完整课件 教材和参考书教材和参考书 教材: 数值分析,电子科技大学应用数学学院,钟尔杰, 黄廷祝主编,高等教育出版社 参考书:参考书: 数值方法(数值方法(MATLAB版)(第三版),版)(第三版),John H. Mathews,Kurtis D. Fink 著,电子工业出版社;著,电子工业出版社; 数值分析(第四版),李庆扬,王能超,易大义编,清华数值分析(第四版),李庆扬,王能超,易大义编,清华 大学出版社;大学出版社; 计算方法(第二版),邓建中,刘之行,西安交通大学出计算方法(第二版),邓建中,刘之行,西安交通大学出 版社;版社; 数值分析与实验
2、学习指导,蔡大用编,清华大学出版社。数值分析与实验学习指导,蔡大用编,清华大学出版社。 教学要求教学要求 了解计算方法研究的主要内容; 掌握计算方法的基本概念和基本原理,进一步提掌握计算方法的基本概念和基本原理,进一步提 高抽象思维和逻辑推理的能力;高抽象思维和逻辑推理的能力; 掌握数值计算的各种方法(或算法)的基本思掌握数值计算的各种方法(或算法)的基本思 想,进一步提高数值计算能力想,进一步提高数值计算能力 ; 能够与实际问题相结合,利用所学算法解决一些能够与实际问题相结合,利用所学算法解决一些 实际的数学模型问题实际的数学模型问题 ; 能够利用数学软件编程实现所学算法(可用能够利用数学软
3、件编程实现所学算法(可用 MATLAB,MATHEMATICA等)。等)。 课程内容 引论(误差、有效数字等) 非线性方程求根 线性方程组的直接解法 线性方程组的迭代解法 数值揑值方法 数据拟吅方法 数值积分不数值微分 常微分方程的数值解法 说明:说明: 学时: 授课48上机16 (课程设计) 成绩评定: 期末考试(60%) 半期考试(20%) 平时成绩(20%) 第一章第一章 引论引论 第一章第一章 引论引论 1.1 计算方法研究的对象不特点 1.2 数值计算的误差不有效数字 1.3 数值运算的误差估计 1.4 数值计算中的一些基本原则 计算方法是做什么用的? 1.1 计算方法研究的对象不特
4、点 计算计算 方法方法 输入复杂问题或运算输入复杂问题或运算 .),(,)( ,ln, xf dx d dxxf bxAxax b a x 计算机计算机 近似解近似解 研究对象:用计算机求解各种数学问题的数值计算 方法及其理论不软件实现。 计算机解决科学计算的过程:计算机解决科学计算的过程: 实际问题实际问题数学模型数学模型数值计算方法数值计算方法程序设计程序设计结果结果 实际上述过程可分两部分:实际上述过程可分两部分: 1.由实际问题应用有关知识和数学理论建立模型,由实际问题应用有关知识和数学理论建立模型, -应用数学任务应用数学任务 2.由数学模型提出求解的数值计算方法直到编程出结果,由数
5、学模型提出求解的数值计算方法直到编程出结果, -计算数学任务计算数学任务 计算方法是计算数学的一个主要部分,研究的即是后半 部分,将理论与计算相结合。 特点: 面向计算机,提供切实可行的算法; 有可靠的理论分析,能达到精度要求,保证近 似算法的收敛性和数值稳定性; 要有好的计算复杂性,节省时间及存储量; 有数值实验,证明算法有效。 算法基本结构:算法基本结构:顺序,分支,循环顺序,分支,循环 算法描述:算法描述:程序或流程图程序或流程图 常采用的处理方法: 构造性方法 离散化方法 递推化方法 迭代法 近似替代方法 以直代曲法 化整为零的处理方法 外推法 微积分的若干定理: 罗尔定理和微分中值定
6、理; 介值定理及推论; 泰勒公式(一元、二元); 积分中值定理; 高等代数的若干概念和结论:高等代数的若干概念和结论: 多项式;多项式; 行列式;行列式; 初等矩阵;初等矩阵; 特殊三角阵。特殊三角阵。 数学基础:数学基础: 1.2 数值计算的误差与有效数字数值计算的误差与有效数字 1.2.1 误差来源与分类误差来源与分类: 按来源分,分为固有误差和计算误差。按来源分,分为固有误差和计算误差。 固有误差固有误差:建立模型时已存在。:建立模型时已存在。 模型误差模型误差:建立数学模型时所引起的误差;建立数学模型时所引起的误差; 观测误差观测误差:测量工具的限制或在数据的获取时测量工具的限制或在数
7、据的获取时 随机因素所引起的物理量的误差。随机因素所引起的物理量的误差。 计算误差计算误差:计算过程中出现的误差。:计算过程中出现的误差。 截断误差截断误差:用数值方法求解数学模型时,用简单用数值方法求解数学模型时,用简单 代替复杂代替复杂, ,或者用有限过程代替无限过程所引起的或者用有限过程代替无限过程所引起的 误差;误差; . ! . ! 2 1 2 n xx xe n x ! . ! 2 1)( 2 n xx xxS n n 10 , )!1( )( 1 x n n x e n x xSe 例如:例如: 由泰勒余项定理得其截断误差为:由泰勒余项定理得其截断误差为: 舍入误差舍入误差:计算
8、机表示的数的位数有限,通常用计算机表示的数的位数有限,通常用 四舍五入的办法取近似值,由此引起的误差。四舍五入的办法取近似值,由此引起的误差。 14159265. 3 414213562. 12 166666666. 0 6 1 ! 3 1 1415927. 3 4142136. 12 16666667. 0 ! 3 1 将将 作作Taylor展开后再积分展开后再积分 2 x e . 9 1 !4 1 7 1 !3 1 5 1 !2 1 3 1 1 ). !4!3!2 1( 1 0 864 2 1 0 dx xxx xdxe 2 x S4 R4 743002401033301 42 1 10
9、1 3 1 1 4 .S 0010200050. 舍入误差舍入误差 21 0 0 0050 0010 006 -x edx. 总体误差的的 (= 0.747) , 1 0 4 Sdxe 2 x 取取 则则 . 11 1 !5 1 9 1 !4 1 4 R 称为称为截断误差截断误差 0050 9 1 !4 1 4 .R 这里这里 dxe x 1 0 2 例:例:近似计算近似计算 解:解: 例:蝴蝶效应 纽约的一只蝴蝶翅膀一拍, 风和日丽的北京就刮起台风来了?! NY BJ 误差的传播与积累:误差的传播与积累: 该问题是一病态问题。该问题是一病态问题。 营长对值班军官营长对值班军官: : 明晚大约
10、明晚大约 8 8点钟左右,哈雷彗星将可点钟左右,哈雷彗星将可 能在这个地区看到,这种彗星每隔能在这个地区看到,这种彗星每隔 7676年才能看见一次。年才能看见一次。 命令所有士兵着野战服在操场上集合,我将向他们解释命令所有士兵着野战服在操场上集合,我将向他们解释 这一罕见的现象。如果下雨的话,就在礼堂集合,我为这一罕见的现象。如果下雨的话,就在礼堂集合,我为 他们放一部有关彗星的影片。他们放一部有关彗星的影片。 值班军官对连长值班军官对连长: : 根据营长的命令,明晚根据营长的命令,明晚8 8点哈雷彗星点哈雷彗星 将在操场上空出现。如果下雨的话,就让士兵穿着野战服将在操场上空出现。如果下雨的话
11、,就让士兵穿着野战服 列队前往礼堂,这一罕见的现象将在那里出现。列队前往礼堂,这一罕见的现象将在那里出现。 据说,美军据说,美军 1910 1910 年的一次部队的命令传递是这样的年的一次部队的命令传递是这样的: : 班长对士兵班长对士兵: : 在明晚在明晚8 8点下雨的时候,著名的点下雨的时候,著名的7676岁哈雷岁哈雷 将军将在营长的陪同下身着野战服,开着他那将军将在营长的陪同下身着野战服,开着他那“彗星彗星”牌牌 汽车,经过操场前往礼堂。汽车,经过操场前往礼堂。 连长对排长连长对排长: : 根据营长的命令,明晚根据营长的命令,明晚8 8点,非凡的哈雷点,非凡的哈雷 彗星将身穿野战服在礼堂
12、中出现。如果操场上下雨,营长彗星将身穿野战服在礼堂中出现。如果操场上下雨,营长 将下达另一个命令,这种命令每隔将下达另一个命令,这种命令每隔7676年才会出现一次。年才会出现一次。 排长对班长排长对班长: : 明晚明晚8 8点,营长将带着哈雷彗星在礼堂中点,营长将带着哈雷彗星在礼堂中 出现,这是每隔出现,这是每隔 7676年才有的事。如果下雨的话,营长将年才有的事。如果下雨的话,营长将 命令彗星穿上野战服到操场上去。命令彗星穿上野战服到操场上去。 1.2.2 误差与有效数字误差与有效数字: 定义定义1.1 设设x*为某一数据的准确值,为某一数据的准确值,x为为x*的一个近的一个近 似值,称似值
13、,称e(x)=x-x*(近似值准确值近似值准确值)为近似值为近似值x的的绝对绝对 误差误差,简称误差。,简称误差。 e(x) 可正可负,当可正可负,当e(x) 0时近似值偏大,叫强近似值;当时近似值偏大,叫强近似值;当e(x) 0时近似值偏小时近似值偏小,叫弱近似值。叫弱近似值。 由于由于x*通常无法确定,只能估计其绝对误差值通常无法确定,只能估计其绝对误差值 不超过某整数不超过某整数 (x),即即 )(| )(|xxxxe 则称则称 (x) 为为绝对误差限绝对误差限。 1. 绝对误差与相对误差绝对误差与相对误差: )()(xxxxx )( * xxx 由上式得由上式得 可知可知x*的范围。或
14、记为的范围。或记为 但误差但误差e(x)并不足以刻划并不足以刻划x的精度。的精度。 如:如: x*=152, x=15, (x) =2; y*=10005, y=1000, (y)=5 因此考虑精度时除看误差大小外,还应考虑精确值本因此考虑精度时除看误差大小外,还应考虑精确值本 身的大小,故引入相对误差概念。身的大小,故引入相对误差概念。 定义定义1.2 设设x*为某一数据的准确值,为某一数据的准确值,x为为x*的一个近似的一个近似 值,称值,称 )0(, )( )( * x x xx x xe xer 为近似值为近似值x的的相对误差相对误差。 实际计算时,由于实际计算时,由于x*不知,通常取
15、不知,通常取 x xx x xe xer * )( )( 如果存在一适当小的如果存在一适当小的正数正数 r , ,使得使得 rr x xx x xe xe )( )( 则称则称 r为 为相对误差限相对误差限。 例:例:x=15, (x) =2, r( (x)=2/15=13.33% )=2/15=13.33%; y=1000, (y)=5 , r( (y)=5/1000=0.5% )=5/1000=0.5% 2. 有效数字:有效数字: 定义定义1.3 若近似值若近似值x的误差限是某一位的半个单位,的误差限是某一位的半个单位, 该位到该位到x的第一位非零数字共有的第一位非零数字共有n位,就说位,
16、就说x有有n位有位有 效数字,效数字,x可表示为可表示为 m n aaax10. 0 21 其中,其中, a1, a2, , an 都是都是 0 9 中的任一整数,中的任一整数, 但但 a10。 其绝对误差限满足:其绝对误差限满足: nm xxxe 10 2 1 )( 如:如:x*= =3.14159265 (1)取)取x3.14,则则m=1, 31 10 2 1 005. 0002. 0 xx 即即n=3,有三位有效数字;有三位有效数字; (2)取)取x=3.1416,则则m=1, 51 10 2 1 00005. 0000008. 0 xx 即即n=5,有有5位有效数字。位有效数字。 定理
17、定理1:设近似数设近似数x表示为表示为 m n aaax10. 0 21 若若x具有具有n位有效数字,则其相对误差限为位有效数字,则其相对误差限为 )1( 1 10 2 1 )( n r a x 反之,若反之,若x的相对误差限为的相对误差限为 )1( 1 10 )1(2 1 )( n r a x 则则x至少具有至少具有n位有效数字。位有效数字。 例例1. 某零件质量取决于零件某参数某零件质量取决于零件某参数x,设参数标定值设参数标定值 为为1个单位,在生产过程中允许参数与标定值间有一定个单位,在生产过程中允许参数与标定值间有一定 误差,据此将零件分成误差,据此将零件分成A、B、C三等,等级由相
18、对误三等,等级由相对误 差限决定,差限决定,A:1%,B:5 %,C:10 %,试确定三个试确定三个 等级的零件参数允许变化的范围。等级的零件参数允许变化的范围。 解解 x*=1, 由由 r x xx 得得 故故 )1 (*)1 (* rr xxx 将三个相对误差限分别带入,得范围如下:将三个相对误差限分别带入,得范围如下: A:x 0.99,1.01 B:x 0.95,0.99)U( 1.01,1.05 C:x 0.9,0.95)U(1.05,1.1 rr x x 1 * 例例2. 测量一物体的长度为测量一物体的长度为954cm,问测量数据的相对问测量数据的相对 误差限多大?误差限多大? 解
19、解 因因实际问题所截取的近似数,其绝对误差限一般实际问题所截取的近似数,其绝对误差限一般 不超过最小刻度的半个单位不超过最小刻度的半个单位, 故当故当x=954cm时,有时,有 (x)0.5cm, 而而x的相对误差的相对误差 er(x) 0.5/954=0.0005241.0.00053=0.053 % 故故 r( (x)= )=0.053 %. 例例3. 200.1%要使的相对误差不超过,应取 几位有效数字? 解解: . 420 1 a的首位数是 .20位有效数字有的近似值设nx * 1 * 1 |1 |( )|10 |2 n r xx ex xa 则由定理则由定理1,相对误差满足,相对误差
20、满足 001. 010 42 1 1 n 097. 3n 即应取即应取4位有效数字,近似值的误差不超过位有效数字,近似值的误差不超过0.1%. 1.3 数值运算的误差估计 1. 函数运算的误差估计:函数运算的误差估计: 设设y=f(x)为一元函数,自变量准确值为一元函数,自变量准确值x*,对应函数准确对应函数准确 值值y*=f(x*),x误差为误差为e(x),误差限为误差限为 (x),函数近似值函数近似值 误差误差e(y),误差限为误差限为 (y)。则(可由则(可由Taylor公式推得)公式推得) )(| )( |)(xxfy )( | )(| | )( | )(x xf xxf y rr 对
21、于多元函数对于多元函数 ),( 21n xxxfz 设准确值设准确值 ),( * 2 * 1 * n xxxfz 由多元函数由多元函数Taylor公式,可得误差估计:公式,可得误差估计: 1 ( )() n k k k f zx x 相对误差限为:相对误差限为: 1 () ( ) n rk rk k k xf zx xz 2. 算术运算的误差估计:算术运算的误差估计: 两个近似数两个近似数x1, x2,其误差限分别为其误差限分别为 (x1), (x2),它它 们进行加、减、乘、除运算得到的误差限分别为:们进行加、减、乘、除运算得到的误差限分别为: )()()( 2121 xxxx )()()(
22、 122121 xxxxxx 2 2 1221 2 1 )()( )( x xxxx x x 例例4: 设设a=2.31,b=1.93,c=2.24都是三位有效数字的近似数,令都是三位有效数字的近似数,令 p=a+bc,求求 (p)和和 r( (p p) ), ,并判断并判断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
23、( (p p)= )= (p)/|p| 0.02585/6.6332 0.0039=0.39% 因为因为 (p) 0.025850.05=1/2 101-2 所以所以 p=6.6332中只有两位有效数字。中只有两位有效数字。 1. 避免除数的绝对值进小亍被除数的绝对值; 2. 避免两个相近的数相减; 3. 防止大数“吃掉”小数现象; 4. 简化计算步骤,减少运算次数; 5. 选用数值稳定性好的算法。 1.4 数值运算中的一些基本原则 1. 避免除数的绝对值远小于被除数的绝对值;避免除数的绝对值远小于被除数的绝对值; 设设 )0( x x y z若若|x|x2|xn| |,应按绝对值由小到大的顺
展开阅读全文