数值分析第五版章课件.ppt
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《数值分析第五版章课件.ppt》由会员分享,可在线阅读,更多相关《数值分析第五版章课件.ppt(91页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、数值分析第第5 5版版 编编2023/2/21第1章 数值分析与科学计算引论数值分析与科学计算引论数值分析研究对象、作用与特点数值分析研究对象、作用与特点数值计算的误差数值计算的误差误差定性分析与避免误差危害误差定性分析与避免误差危害数值计算中算法设计的技术数值计算中算法设计的技术数学软件数学软件2023/2/22第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件1 1 研究对象研究对象 用计算机求解数学问题的数值计算方法、理论及软件实现用计算机求解数学问题的数值计算方法、理论及
2、软件实现 实际问题实际问题 数学模型数学模型 数值计算方法数值计算方法 程序设计(数学软件)程序设计(数学软件)上机计算求出结果上机计算求出结果应用数学应用数学计算数学即数值分析计算数学即数值分析数值分析(计算方法)数值分析(计算方法)的研究对象的研究对象(理论与计算结合理论与计算结合)插值与函数逼近(插值与函数逼近(2、3)数值微分与数值积分()数值微分与数值积分(4)非线性方程数值解(非线性方程数值解(7)数值线性代数()数值线性代数(5、6、8)常微与常微与偏微分方程的数值解偏微分方程的数值解(9)等)等1.1 1.1 数值分析的对象、作用与特点数值分析的对象、作用与特点2023/2/2
3、3第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件例例 考虑线性方程组数值解问题考虑线性方程组数值解问题相关理论与精确解法相关理论与精确解法根据方程的特点研究算法及相关理论根据方程的特点研究算法及相关理论纯数学纯数学计算数学计算数学3 数值分析特点数值分析特点 面向计算机提供有效算法;面向计算机提供有效算法;TSP问题问题 可靠的理论分析(精度、收敛、稳定)可靠的理论分析(精度、收敛、稳定);好的计算复杂性(时、空)好的计算复杂性(时、空);阿基米德、刘徽、祖冲之阿基米德、刘徽
4、、祖冲之 数值实验数值实验 B.C.3 A.D.3 A.D.52 数值分析作用数值分析作用 科学研究的手段:科学理论、科学实验与科学研究的手段:科学理论、科学实验与科学计算科学计算2023/2/24第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件1.2 1.2 数值计算的误差数值计算的误差1 误差分类误差分类模型误差模型误差:数学模型数学模型 实际问题实际问题观测误差观测误差:由观测产生由观测产生截断误差截断误差/方法误差方法误差:近似解近似解 精确解精确解舍入误差舍入误差:计
5、算机字长的限制计算机字长的限制(不讨论不讨论)2 误差与有效数字误差与有效数字记记x为准确值,为准确值,x*为为 x的一个近似值的一个近似值 定义定义1 称称 为近似值为近似值 的绝对误差,简称的绝对误差,简称误差误差。强近似值强近似值:当当弱近似值弱近似值:当当2023/2/25第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件误差限误差限:误差绝对值的上界误差绝对值的上界(不能完全反映近似值的好坏不能完全反映近似值的好坏)相对误差相对误差:相对误差限相对误差限:(可反映出近似
6、程度的好坏可反映出近似程度的好坏)定义定义2 有效数字有效数字 若近似值若近似值 的误差限是某一位的半个单位,该位到的误差限是某一位的半个单位,该位到 的的 第一位非零数第一位非零数 字共有字共有n位,就说位,就说 有有n位有效数字。位有效数字。2023/2/26第1章 数值分析与科学计算引论注注:有效位数与小数点后有多少位无关;有效位数与小数点后有多少位无关;m相同情况下,有效位数越多,误差限越小;相同情况下,有效位数越多,误差限越小;相对误差及相对误差限是无量纲的,绝对误差及误差限是有量纲的。相对误差及相对误差限是无量纲的,绝对误差及误差限是有量纲的。研究对研究对象象作用特点作用特点数值计
7、算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件2023/2/27第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件3 数值运算的误差估计数值运算的误差估计2023/2/28第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件1.3 1.3 误差定性分析及避免误差危害误差定性分析及避免误差危害概率分析法概率分析法 向后
8、误差分析法向后误差分析法 区间分析法区间分析法1.病态问题与条件数病态问题与条件数病态问题病态问题输入输入(微小的扰动微小的扰动)输出输出(相对误差很大相对误差很大)条件数条件数2023/2/29第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件2.算法的数值稳定性算法的数值稳定性定定义义3 一一个个算算法法如如果果输输入入数数据据有有误误差差,而而在在计计算算过过程程中中舍舍入入误误差差不不增增长长,则则称称此此算算法法是是数数值值稳稳定定的的,否否则则称称此此算算法法为为不不
9、稳稳定定 的。的。例例1.1:P.2023/2/210第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件2023/2/211第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件3.误差危害的避免误差危害的避免(1)避免除数的绝对值远远小于被除数绝对值的除法;)避免除数的绝对值远远小于被除数绝对值的除法;(2)避免两相近数相减,引起有效数字严重损失;)避免两相近数相减,
10、引起有效数字严重损失;(3)防止大数吃小数;防止大数吃小数;(4)简化计算简化计算 步骤,减少运算次数;步骤,减少运算次数;秦九韶、秦九韶、Horner(5)数值稳定性。数值稳定性。2023/2/212第1章 数值分析与科学计算引论研究对研究对象象作用特点作用特点数值计算数值计算误差误差误差分析误差分析避免危害避免危害数值计算数值计算算法设计算法设计数学软件数学软件算法设计的三种基本技术算法设计的三种基本技术(1)化大为小的缩减技术)化大为小的缩减技术 Zeno悖论悖论 古希腊哲学家古希腊哲学家 结绳记数(结绳记数(13)10=(1101)2(2)化难为易的校正技术)化难为易的校正技术 用四则
11、运算计算开方用四则运算计算开方(3)化粗为精的松弛技术)化粗为精的松弛技术1.1.数值计算中算法设计的技术数值计算中算法设计的技术周易.系辞下说:“上古结绳而治,后世圣人(伏羲)易之以书契“。2023/2/213第1章 数值分析与科学计算引论谢谢谢谢!第一章习题第一章习题1,5,12,12023/2/214第1章 数值分析与科学计算引论第2章 插值法插值法引言引言拉格朗日(拉格朗日(Lagrange)插值)插值均差与牛顿(均差与牛顿(Newton)插值)插值埃尔米特(埃尔米特(Hermite)插值)插值分段低次插值分段低次插值三次样条插值三次样条插值2023/2/215第2章 插值法引言引言拉
12、格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2.1 2.1 引言引言 在科学计算中经常要用较简单的函数来逼近较复杂的函数。按函数逼近问题提法的不同,通常有插值、函数逼近及曲线拟合等三种问题。(1)插值 设函数在点上的函数值分别为,求一简单函数 ,使 (2.1)0 0 2023/2/216第2章 插值法称点为插值节点,区间为插值区间,为被插值函数,为的插值函数。引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值 由于多项式具有结构简单,数值计算和理论分析都很方便
13、的优点,因此通常取 为多项式。相应的插值法称为多项式插值。2023/2/217第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值(2)函数逼近 设 为定义在区间 上的某类函数构成的线性空间,为 的子集,对于 ,求 ,使 与 的差在某种度量意义下最小,这就是函数逼近问题。通常取为连续函数空间,而 通常是多项式、有理函数或三角多项式函数等。另一种是2-范数对应的函数逼近分别称为最佳一致逼近和最佳平方逼近。常用的度量标准有两种:一种是-范数2023/2/218第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿
14、插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值(3)曲线拟合 设已知函数 的数据点为,在某一函数类中求函数其中为待定系数,使在节点 处的误差的平方和最小。2023/2/219第2章 插值法引言引言拉格朗拉格朗日插值日插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2.2 2.2 拉格朗日插值拉格朗日插值n次插值多顶式:次插值多顶式:1 线性插值(一次)线性插值(一次)y O x xk xk+1 yk yk+12023/2/220第2章 插值法引言引言拉格朗拉格朗日插值日插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次
15、分段低次插值插值三次样条三次样条插值插值 x y O xk xk+1 1 利用函数的零点及一次多项式性质,由待定系数法可求得2023/2/221第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2 抛物插值(二次)抛物插值(二次)由待定系数法可求得:由待定系数法可求得:1 y O x2023/2/222第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2023/2/223第2章 插值法3 拉格朗日插值多项式拉格朗日插值多项式将线
16、性与抛物型插值推广到一般的情形将线性与抛物型插值推广到一般的情形引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2023/2/224第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2023/2/225第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值证明 设所求的插值多项式为 则由插值条件得关于 的线性方程组从而线性方程组有唯一解。这就证明了 的存在唯一性
17、。2023/2/226第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值4 插值余项与误差估计插值余项与误差估计2023/2/227第2章 插值法例2.1 已知 (1)用线性插值及抛物插值计算 的近似值;(2)并问它们各有几位有效数字;(3)求抛物插值的误差。解 先用线性插值计算,取 ,由线性插值公式得 所以而因此用线性插值所得 的近似值具有3位有效数字。引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2023/2/228第2章 插值法引言
18、引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值再用抛物插值计算,取由抛物插值公式得:所以因此用抛物插值所得 的近似值具有4位有效数字。线性插值仅用两个节点上的信息,精度自然较低,而抛物插值用了三个节点上的信息,精度通常会有所提高。2023/2/229第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值这又一次说明了用抛物插值所得的近似值具有4位有效数字。2023/2/230第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃
19、尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2.3 2.3 均差与牛顿插值均差与牛顿插值一阶均差一阶均差:二阶均差二阶均差:K阶均差阶均差:承袭性承袭性1 均差及其性质均差及其性质2023/2/231第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值2 牛顿插值公式牛顿插值公式可克服可克服Lagrange插值法插值法无承袭性的缺点。无承袭性的缺点。2023/2/232第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插
20、值插值差商表一阶差商二阶差商三阶差商2023/2/233第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值一阶差商二阶差商例2.2解 的表达式中前两项为线性插值,加上第三项后为二次插值,与前例比较结果是相同的。2023/2/234第2章 插值法引言引言拉格朗日拉格朗日插值插值牛顿插值牛顿插值埃尔米特埃尔米特插值插值分段低次分段低次插值插值三次样条三次样条插值插值 前面讨论插值问题只提函数值条件,没有导数条件。有些实际问题不但要求插值函数与被插值函数在节点上函数值相同,即“过点”,而且要求导数值也相同,即“相切”
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 第五 课件
![提示](https://www.taowenge.com/images/bang_tan.gif)
限制150内