第一章数值分析2008(绪论).ppt
《第一章数值分析2008(绪论).ppt》由会员分享,可在线阅读,更多相关《第一章数值分析2008(绪论).ppt(66页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、海军航空工程学院二0五教研室数值分析数值分析 海军航空工程学院二0五教研室教材及参考书目&教材教材 (Text Book)数值分析与算法数值分析与算法 徐士良徐士良 编著编著(机械工业出版社)(机械工业出版社)(华中科技大学出版社)(华中科技大学出版社)&参考书目参考书目(Reference)Numerical Analysis (Seventh Edition)数值分析数值分析 (第七版(第七版 影印版)影印版)Richard L.Burden&J.Douglas Faires (高等教育出版社)高等教育出版社)数值方法数值方法 金一庆、陈越金一庆、陈越 编著编著(机械工业出版社)(机械工业
2、出版社)数值分析数值分析 李庆扬、易大义、王能超李庆扬、易大义、王能超 编著编著海军航空工程学院二0五教研室数值分析的作用数值数值分析分析输入复杂问题或运算输入复杂问题或运算 计算机计算机近似解近似解海军航空工程学院二0五教研室绪论1 1 绪论绪论 1.1 1.1 数值分析数值分析 研究对象与特点研究对象与特点 1.2 1.2 数值计算的误差数值计算的误差 1.2.1 1.2.1 误差来源与分类误差来源与分类 1.2.2 1.2.2 误差与有效数字误差与有效数字 1.2.3 1.2.3 数值计算中应注意的几个原则数值计算中应注意的几个原则 1.3 1.3 算法算法 1.3.1 1.3.1 算法
3、的基本概念算法的基本概念 1.3.2 1.3.2 数值型算法的特点数值型算法的特点 1.3.3 1.3.3 算法设计基本方法算法设计基本方法 1.3.4 1.3.4 算法的复杂度算法的复杂度海军航空工程学院二0五教研室1.11.1数值分析数值分析 研究对象与特点研究对象与特点 数值分析是研究适合于在计算机上使用的实际可数值分析是研究适合于在计算机上使用的实际可行、理论可靠、计算复杂性好的数值计算方法行、理论可靠、计算复杂性好的数值计算方法 第一,面向计算机,要根据计算机特点提供实际第一,面向计算机,要根据计算机特点提供实际可行的算法可行的算法 。第二,要有可靠的理论分析,数值分析中的算法第二,
4、要有可靠的理论分析,数值分析中的算法理论主要是连续系统的离散化及离散型方程数值理论主要是连续系统的离散化及离散型方程数值求解。求解。第三,要有良好的复杂性及数值试验第三,要有良好的复杂性及数值试验 海军航空工程学院二0五教研室要求要求 本课程简单介绍数值方法及其理论,大部分本课程简单介绍数值方法及其理论,大部分定理的证明需要大家自己去学习掌握,着重介绍定理的证明需要大家自己去学习掌握,着重介绍具体的算法设计及编程技巧,作为基本要求让大具体的算法设计及编程技巧,作为基本要求让大家做一些计算机上的数值试验,它对加深算法的家做一些计算机上的数值试验,它对加深算法的理解是很有好处的。理解是很有好处的。
5、海军航空工程学院二0五教研室1.21.2数值计算的误差数值计算的误差 1.2.11.2.1误差来源与分类误差来源与分类 本课程只研究数值求解数学问题时产生的误差,本课程只研究数值求解数学问题时产生的误差,它们主要有以下三类。它们主要有以下三类。1 1截断误差或方法误差截断误差或方法误差 2 2舍入误差舍入误差 3 3输入数据误差输入数据误差 海军航空工程学院二0五教研室1.2.11.2.1误差来源与分类误差来源与分类1 1 1 1截断误差或方法误差截断误差或方法误差它它是是指指将将数数学学问问题题转转化化为为数数值值计计算算问问题题时时产产生生的的误误差差,通常是用有限过程近似无限过程时产生的
6、误差。通常是用有限过程近似无限过程时产生的误差。海军航空工程学院二0五教研室1.2.11.2.1误差来源与分类误差来源与分类2 2 2舍入误差 数值计算时由于计算是有限位的,所以原始数据、中间结果和最后结果都要舍入,这就产生舍入误差。在十进制运算中一般采用四舍五入。例如例如1/31/3写成写成0.333 30.333 3,3.141 63.141 6等等,都有等等,都有舍入误差。舍入误差。海军航空工程学院二0五教研室1.2.11.2.1误差来源与分类误差来源与分类3 33 3输入数据误差输入数据误差 称为初始误差,这些误差对计算也将造成影响,称为初始误差,这些误差对计算也将造成影响,但分析初始
7、误差与对舍入误差分析相似,因此可将但分析初始误差与对舍入误差分析相似,因此可将它归入第二类。它归入第二类。由于对大规模数值计算问题舍入误差目前尚无由于对大规模数值计算问题舍入误差目前尚无有效方法进行定量估计,所以我们主要进行定性分有效方法进行定量估计,所以我们主要进行定性分析析.但对误差估计的基本概念及较简单的数值运算误但对误差估计的基本概念及较简单的数值运算误差估计还需作简单介绍。差估计还需作简单介绍。海军航空工程学院二0五教研室1.2.2 1.2.2 误差与有效数字误差与有效数字1 1 1 1误差的定义误差的定义 绝对误差绝对误差 /*absolute error*/*absolute e
8、rror*/其中其中x x为精确值,为精确值,x x*为为x x的近似值。的近似值。,例如:,例如:工程上常记为工程上常记为,称为绝对误差限,称为绝对误差限 /*accuracy*/*accuracy*/,的上限记为的上限记为注:注:注:注:e*e*e*e*理论上讲是唯一确定的,可能取正理论上讲是唯一确定的,可能取正理论上讲是唯一确定的,可能取正理论上讲是唯一确定的,可能取正,也可能取负。也可能取负。也可能取负。也可能取负。e*e*e*e*0 0 0 0 不唯一,当然不唯一,当然不唯一,当然不唯一,当然 e*e*e*e*越小越具有参考价值。越小越具有参考价值。越小越具有参考价值。越小越具有参考
9、价值。海军航空工程学院二0五教研室1.2.2 误差与有效数字2 1 1误差的定义误差的定义 相对误差相对误差 /*relative error*/x 的的相对误差上限相对误差上限/*relative accuracy*/定义为定义为注:从注:从注:从注:从 的定义可见,的定义可见,的定义可见,的定义可见,实际上被实际上被实际上被实际上被偷换偷换偷换偷换成了成了成了成了 ,而后才考察,而后才考察,而后才考察,而后才考察其上限。那么这样的偷换是否其上限。那么这样的偷换是否其上限。那么这样的偷换是否其上限。那么这样的偷换是否合法合法合法合法?严格的说法是,严格的说法是,严格的说法是,严格的说法是,与
10、与与与 是否反映了是否反映了是否反映了是否反映了同一数量级同一数量级同一数量级同一数量级的误差?的误差?的误差?的误差?关于此问题的详细讨论可见教材。关于此问题的详细讨论可见教材。关于此问题的详细讨论可见教材。关于此问题的详细讨论可见教材。海军航空工程学院二0五教研室1.2.2 误差与有效数字3 2 2有效数字有效数字 用科学计数法,记用科学计数法,记 (其中(其中 )。若)。若 (即(即 的截取按四舍五入规则),则称的截取按四舍五入规则),则称 为有为有n 位有效数字,精确到位有效数字,精确到 。例:例:问:问:有几位有效数字?有几位有效数字?证明:证明:有有 位有效数字,精确到小数点后第位
11、有效数字,精确到小数点后第 位。位。43海军航空工程学院二0五教研室1.2.2 误差与有效数字误差与有效数字4 4 3.有效数字与相对误差的关系有效数字与相对误差的关系 有效数字有效数字 相对误差限相对误差限已知已知 x*有有 n 位位有效数字有效数字,则其,则其相对误差限相对误差限为为海军航空工程学院二0五教研室1.2.2 误差与有效数字误差与有效数字5 5 相对误差限相对误差限 有效数字有效数字已知已知 x*的的相对误差限相对误差限可写为可写为则则可见可见 x*至少至少有有 n 位位有效数字有效数字。海军航空工程学院二0五教研室1.2.2 误差与有效数字64.函数的误差估计函数的误差估计
12、/*Error Estimation for Functions*/问题问题:对于:对于 y =f(x),若用若用 x*取代取代 x,将对将对y 产生什么影响?产生什么影响?分析分析:e*(y)=f(x*)f(x)e*(x)=x*x=f()(x*x)x*与与 x 非常接近时,可认为非常接近时,可认为 f()f(x*),则有:则有:|e*(y)|f(x*)|e*(x)|即:即:x*产生的误差经过产生的误差经过 f 作用后被放大作用后被放大/缩小了缩小了|f(x*)|倍。故称倍。故称|f(x*)|为为放大因子放大因子 或或 绝对条件数绝对条件数海军航空工程学院二0五教研室1.2.2 误差与有效数字
13、误差与有效数字7 7 相对误差条件数相对误差条件数 f 的条件数在某一点是的条件数在某一点是小小大大,则称,则称 f 在该点是在该点是好条件的好条件的 坏条件的坏条件的。注:关于多元函数注:关于多元函数 的讨论,请参阅教材。的讨论,请参阅教材。海军航空工程学院二0五教研室1.2.3 1.2.3 数值计算中应注意的几个原则数值计算中应注意的几个原则1 1 1.1.关于数值稳定性的算法。关于数值稳定性的算法。一一个个程程序序往往往往要要进进行行大大量量的的四四则则运运算算才才能能得得出出结结果果,每每一一步步的的运运算算均均会会产产生生舍舍入入误误差差。在在运运算算过过程程中中,舍舍入入误误差差能
14、能控控制制在在某某个个范范围围内内的的算算法法称称之之为为数数值稳定的算法值稳定的算法,否则就称之为不稳定的算法。否则就称之为不稳定的算法。海军航空工程学院二0五教研室1.2.3 1.2.3 数值计算中应注意的几个原则数值计算中应注意的几个原则2 2 例:计算例:计算 公式一:公式一:注意此公式注意此公式精确精确成成立立记为记为则初始误差则初始误差海军航空工程学院二0五教研室1.2.3 1.2.3 数值计算中应注意的几个原则数值计算中应注意的几个原则3 3考察第考察第n步的误差步的误差我们有责任改变。我们有责任改变。造成这种情况的是造成这种情况的是不稳定的算法不稳定的算法 /*unstable
15、 algorithm*/迅速积累,误差呈递增走势。迅速积累,误差呈递增走势。可见初始的小扰动可见初始的小扰动 公式二:公式二:注意此公式与公式一注意此公式与公式一在理论上在理论上等价等价。方法:先估计一个方法:先估计一个IN,再反推要求的再反推要求的In(n N)。可取可取海军航空工程学院二0五教研室1.2.3 1.2.3 数值计算中应注意的几个原则数值计算中应注意的几个原则4 4取取海军航空工程学院二0五教研室1.2.3 1.2.3 数值计算中应注意的几个原则数值计算中应注意的几个原则5 5 考察反推一步的误差:考察反推一步的误差:以此类推,对以此类推,对 n 1n1时)。时)。写成递归函数
16、有:写成递归函数有:intint fib(intfib(int n)n)if(n=0)return 0;if(n=0)return 0;if(n=1)return 1;if(n=1)return 1;if(n1)return fib(n-1)+fib(n-2);if(n1)return fib(n-1)+fib(n-2);海军航空工程学院二0五教研室1.3.3 1.3.3 算法算法-算法设计基本方法算法设计基本方法-递归法递归法3 3 递递归归算算法法的的执执行行过过程程分分递递推推和和回回归归两两个个阶阶段段。在在递递推推阶阶段段,把把较较复复杂杂的的问问题题(规规模模为为n n)的的求求解
17、解推推到到比比原原问问题题简简单单一一些些的的问问题题(规规模模小小于于n n)的的求求解解。例例如如本本题题中中,求求解解fib(n)fib(n),把把它它推推到到求求解解fib(n-1)fib(n-1)和和fib(n-2)fib(n-2)。也也就就是是说说,为为计计算算fib(n)fib(n),必必须须先先计计算算fib(n-1)fib(n-1)和和fib(n-fib(n-2)2),而而计计算算fib(n-1)fib(n-1)和和fib(n-2)fib(n-2),又又必必须须先先计计算算fib(n-3)fib(n-3)和和fib(n-4)fib(n-4)。依依次次类类推推,直直至至计计算算
18、fib(1)fib(1)和和fib(0)fib(0),分分别别能能立立即即得得到到结结果果1 1和和0 0。在在递递推推阶阶段段,必必须须要有终止递归的情况。要有终止递归的情况。海军航空工程学院二0五教研室1.3.3 1.3.3 算法算法-算法设计基本方法算法设计基本方法-递归法递归法4 4 在编写递归函数时要注意,函数中的局部变在编写递归函数时要注意,函数中的局部变量和参数只局限于当前调用层,当递推进入量和参数只局限于当前调用层,当递推进入“简简单问题单问题”层时,原来层次上的参数和局部变量便层时,原来层次上的参数和局部变量便被隐蔽起来。在一系列被隐蔽起来。在一系列“简单问题简单问题”层,它
19、们各层,它们各有自己的参数和局部变量。有自己的参数和局部变量。由于递归引起一系列的函数调用,并且可能由于递归引起一系列的函数调用,并且可能会有一系列的重复计算,递归算法的执行效率相会有一系列的重复计算,递归算法的执行效率相对较低。当某个递归算法能较方便地转换成递推对较低。当某个递归算法能较方便地转换成递推算法时,通常按递推算法编写程序。例如计算斐算法时,通常按递推算法编写程序。例如计算斐波那契数列的第波那契数列的第n n项的函数项的函数fib(n)fib(n)应采用递推算法,应采用递推算法,即从斐波那契数列的前两项出发,逐次由前两项即从斐波那契数列的前两项出发,逐次由前两项计算出下一项,直至计
20、算出要求的第计算出下一项,直至计算出要求的第n n项。项。海军航空工程学院二0五教研室1.3.3 1.3.3 算法算法-算法设计基本方法算法设计基本方法-回溯法回溯法1 1 回回溯溯法法也也称称为为试试探探法法,该该方方法法首首先先暂暂时时放放弃弃关关于于问问题题规规模模大大小小的的限限制制,并并将将问问题题的的候候选选解解按按某某种种顺顺序序逐逐一一枚枚举举和和检检验验。当当发发现现当当前前候候选选解解不不可可能能是是解解时时,就就选选择择下下一一个个候候选选解解;倘倘若若当当前前候候选选解解除除了了还还不不满满足足问问题题规规模模要要求求外外,满满足足所所有有其其他他要要求求时时,继继续续
21、扩扩大大当当前前候候选选解解的的规规模模,并并继继续续试试探探。如如果果当当前前候候选选解解满满足足包包括括问问题题规规模模在在内内的的所所有有要要求求时时,该该候候选选解解就就是是问问题题的的一一个个解解。在在回回溯溯法法中中,放放弃弃当当前前候候选选解解,寻寻找找下下一一个个候候选选解解的的过过程程称称为为回回溯溯。扩扩大大当当前前候候选选解解的的规规模模,以以继继续续试探的过程称为向前试探。试探的过程称为向前试探。海军航空工程学院二0五教研室1.3.3 1.3.3 算法算法-算法设计基本方法算法设计基本方法-回溯法回溯法2 2【问题问题】组合问题组合问题问问题题描描述述:找找出出从从自自
22、然然数数1 1,2 2,n n中中任任取取r r个数的所有组合。个数的所有组合。采采用用回回溯溯法法找找问问题题的的解解,将将找找到到的的组组合合以以从从小小到到大大顺顺序序存存于于a0a0,a1a1,ar-1ar-1中中,组组合合的元素满足以下性质:的元素满足以下性质:(1 1)ai+1aiai+1ai,后一个数字比前一个大;后一个数字比前一个大;(2 2)ai-i=n-r+1ai-i=n-r+1。海军航空工程学院二0五教研室1.3.3 1.3.3 算法算法-算法设计基本方法算法设计基本方法-回溯法回溯法3 3 首首先先放放弃弃组组合合数数个个数数为为r r的的条条件件,候候选选组组合合从从
23、只只有有一一个个数数字字1 1开开始始。因因该该候候选选解解满满足足除除问问题题规规模模之之外外的的全全部部条条件件,扩扩大大其其规规模模,并并使使其其满满足足上上述述条条件件(1 1),候候选选组组合合改改为为1 1,2 2。继继续续这这一一过过程程,得得到到候候选选组组合合1 1,2 2,3 3。该该候候选选解解满满足包括问题规模在内的全部条件,因而是一个解。足包括问题规模在内的全部条件,因而是一个解。在在该该解解的的基基础础上上,选选下下一一个个候候选选解解,因因a a 22上上的的3 3调调整整为为4 4,以以及及以以后后调调整整为为5 5都都满满足足问问题题的的全全部部要要求求,得得
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第一章 数值分析2008绪论 数值 分析 2008 绪论
限制150内