数值计算方法-复习-第六章
- 格式:pdf
- 大小:335.68 KB
- 文档页数:40
《数值计算方法》复习资料课程的性质与任务数值计算方法是一门应用性很强的基础课,在学习高等数学,线性代数和算法语言的基础上,通过本课程的学习及上机实习、使学生正确理解有关的基本概念和理论,掌握常用的基本数值方法,培养应用计算机从事科学与工程计算的能力,为以后的学习及应用打下良好基础。
第一章数值计算方法与误差分析一考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;绝对误差的传播。
二复习要求1. 知道产生误差的主要来源。
2. 了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系。
3. 知道四则运算中的误差传播公式。
三例题例1设x*= =3.1415926…近似值x=3.14=0.314×101,即m=1,它的绝对误差是-0.001 592 6…,有即n=3,故x=3.14有3位有效数字.x=3.14准确到小数点后第2位.又近似值x=3.1416,它的绝对误差是0.0000074…,有即m=1,n=5,x=3.1416有5位有效数字.而近似值x=3.1415,它的绝对误差是0.0000926…,有即m=1,n=4,x=3.1415有4位有效数字.这就是说某数有s位数,若末位数字是四舍五入得到的,那么该数有s位有效数字;例2 指出下列各数具有几位有效数字,及其绝对误差限和相对误差限:2.000 4 -0.002 00 9 000 9 000.00=2.000 4=0.200 04×101, 它的绝对误差限0.000 05=0.5×10 1―5,即解因为x1m=1,n=5,故x=2.000 4有5位有效数字. a=2,相对误差限1x 2=-0.002 00,绝对误差限0.000 005,因为m =-2,n=3,x 2=-0.002 00有3位有效数字. a 1=2,相对误差限εr ==0.002 5x 3=9 000,绝对误差限为0.5×100,因为m =4, n=4, x 3=9 000有4位有效数字,a =9,相对误差限εr ==0.000 056x 4=9 000.00,绝对误差限0.005,因为m =4,n=6,x 4=9 000.00有6位有效数字,相对误差限为εr ==0.000 000 56由x 3与x 4可以看到小数点之后的0,不是可有可无的,它是有实际意义的. 例3 ln2=0.69314718…,精确到10-3的近似值是多少?解 精确到10-3=0.001,意旨两个近似值x 1,x 2满足,由于近似值都是四舍五入得到的,要求满足,近似值的绝对误差限应是ε=0.0005,故至少要保留小数点后三位才可以。
第一章引论计算方法解决问题的主要思想计算方法的精髓:以直代曲、化繁为简1、采用“构造性”方法构造性方法是指具体地把问题的计算公式构造出来。
这种方法不但证明了问题的存在性,而且有了具体的计算公式,就便于编制程序上机计算。
2、采用“离散化”方法把连续变量问题转为求离散变量问题。
例:把定积分离散成求和,把微分方程离散成差分方程。
3、采用“递推化”方法将复杂的计算过程归结为简单过程的多次重复。
由于递推算法便于编写程序,所以数值计算中常采用“递推化”方法。
4、采用“近似代替”方法计算机运算必须在有限次停止,所以数值方法常表现为一个无穷过程的截断,把一个无限过程的数学问题,转化为满足一定误差要求的有限步来近似替代。
算法的可行性分析时间复杂度、空间复杂度分析算法的复杂性(包含时间复杂性和空间复杂性)。
时间复杂度是算法耗费时间的度量。
算法的空间复杂度是指算法需占用存储空间的量度算法的可靠性分析良态算法、病态算法一个算法若运算过程中舍入误差的积累对最后计算结果影响很大,则称该算法是不稳定的或病态算法,反之称为稳定算法或良态算法。
误差的来源1、模型误差我们所建立的数学模型是对实际问题进行抽象简化而得到的。
因而总是近似的,这就产生了误差。
这种数学模型解与实际问题的解之间出现的误差,称为模型误差。
2、观测误差观测到的数据与实际数据之差。
3、截断误差数学模型的准确解与计算方法的准确解之间的误差。
4、舍入误差由于计算机字长有限,原始数据在计算机上表示会产生误差,每次计算又会产生新的误差,这种误差称为舍入误差。
绝对误差、相对误差定义2 记x*为x的近似数,称E(x)=x-x*为近似数x*的绝对误差,|E(x)|为绝对误差限。
定义3 称Er(x)=(x-x*)/x为近似数x*的相对误差。
实际运算时也将Er*(x)=(x-x*)/x*称为近似数x*的相对误差。
“四舍五入”:即尾数是4或以下则舍去,尾数是6或以上则进1,如果尾数是5,则规定:前面一位数字是偶数则舍去,奇数则进1。
2016计算方法复习务必通过本提纲例子和书上例子掌握如下书本内容:1. 会高斯消去法;会矩阵三角分解法;会Cholesky 分解的平方根法求解方程组2. 会用插值基函数;会求Lagrange, 会计算差商和Newton 插值多项式和余项3. 会Jacobi 迭代、Gauss —Seidel 迭代的分量形式,迭代矩阵,谱半径,收敛性4. 会写非线性方程根的Newton 迭代格式;斯蒂芬森加速5. 会用欧拉预报-校正法和经典四阶龙格—库塔法求解初值问题6. 会最小二乘法多项式拟合7. 会计算求积公式的代数精度;(复化)梯形公式和(复化)辛普生公式求积分;高斯-勒让德求积公式第1章、数值计算引论(一)考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;误差的传播。
(二) 复习要求1。
了解数值分析的研究对象与特点。
2。
了解误差来源与分类,会求有效数字; 会简单误差估计. 3.了解误差的定性分析及避免误差危害。
(三)例题例1. 设x =0.231是精确值x *=0。
229的近似值,则x 有2位有效数字。
例2. 为了提高数值计算精度, 当正数x 充分大时, 应将)1ln(2--x x 改写为)1ln(2++-x x .例3. 3*x 的相对误差约是*x 的相对误差的1/3 倍.第2章、非线性方程的数值解法(一)考核知识点对分法;不动点迭代法及其收敛性;收敛速度; 迭代收敛的加速方法;埃特金加速收敛方法;Steffensen 斯特芬森迭代法;牛顿法;弦截法. (二) 复习要求1.了解求根问题和二分法.2。
了解不动点迭代法和迭代收敛性;了解收敛阶的概念和有关结论。
3。
理解掌握加速迭代收敛的埃特金方法和斯蒂芬森方法。
4。
掌握牛顿法及其收敛性、下山法, 了解重根情形. 5.了解弦截法. (三)例题1。
为求方程x 3―x 2―1=0在区间[1.3,1.6]内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( )(A )11,1112-=-=+k k x x x x 迭代公式 (B )21211,11kk x x x x +=+=+迭代公式(C ) 3/12123)1(,1k k x x x x +=+=+迭代公式 (D )231x x =-迭代公式11221+++=+k k kk x x x x 解:在(A)中,2/32)1(21)(,11)(,11--='-=-=x x x x x x ϕϕ2/3)16.1(21->=1.076故迭代发散。
《数值计算方法》复习资料课程的性质与任务数值计算方法是一门应用性很强的基础课,在学习高等数学,线性代数和算法语言的基础上,通过本课程的学习及上机实习、使学生正确理解有关的基本概念和理论,掌握常用的基本数值方法,培养应用计算机从事科学与工程计算的能力,为以后的学习及应用打下良好基础。
第一章数值计算方法与误差分析一考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;绝对误差的传播。
二复习要求1. 知道产生误差的主要来源。
2. 了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系。
3. 知道四则运算中的误差传播公式。
三例题例1设x*= π=3.1415926…近似值x=3.14=0.314×101,即m=1,它的绝对误差是-0.001 592 6…,有即n=3,故x=3.14有3位有效数字.x=3.14准确到小数点后第2位.又近似值x=3.1416,它的绝对误差是0.0000074…,有即m=1,n=5,x=3.1416有5位有效数字.而近似值x=3.1415,它的绝对误差是0.0000926…,有即m=1,n=4,x=3.1415有4位有效数字.这就是说某数有s位数,若末位数字是四舍五入得到的,那么该数有s位有效数字;例2 指出下列各数具有几位有效数字,及其绝对误差限和相对误差限:2.000 4 -0.002 00 9 000 9 000.00解因为x1=2.000 4=0.200 04×101, 它的绝对误差限0.000 05=0.5×10 1―5,即m=1,n=5,故x=2.000 4有5位有效数字. a1=2,相对误差限x2=-0.002 00,绝对误差限0.000 005,因为m=-2,n=3,x2=-0.002 00有3位有效数字. a1=2,相对误差限εr==0.002 5x3=9 000,绝对误差限为0.5×100,因为m=4, n=4, x3=9 000有4位有效数字,a=9,相对误差限εr==0.000 056x4=9 000.00,绝对误差限0.005,因为m=4,n=6,x4=9 000.00有6位有效数字,相对误差限为εr==0.000 000 56由x3与x4可以看到小数点之后的0,不是可有可无的,它是有实际意义的.例3ln2=0.69314718…,精确到10-3的近似值是多少?解精确到10-3=0.001,意旨两个近似值x1,x2满足,由于近似值都是四舍五入得到的,要求满足,近似值的绝对误差限应是ε=0.0005,故至少要保留小数点后三位才可以。
第6章 逐次逼近法一、考核知识点:向量范数与矩阵范数及其性质,谱半径,对角占优矩阵,迭代法的收敛性,雅可比迭代法与高斯-塞德尔迭代法及其收敛性。
简单迭代法,牛顿迭代法,割线法,收敛性。
二、考核要求:1.了解向量范数的定义、性质;了解矩阵范数的定义、性质,知道谱半径的定义。
2.了解对角占优矩阵;了解迭代法的收敛性。
3.熟练掌握雅可比迭代法其收敛性的判断。
4.熟练掌握高斯-塞德尔迭代法其收敛性的判断。
5.熟练掌握用牛顿迭代法、割线法求方程近似根的方法。
了解其收敛性。
6.掌握用简单迭代法求方程的方法近似根的方法。
了解其收敛性。
三、重、难点分析例1 已知向量X=(1,-2,3),求向量X 的三种常用范数。
解 3max ==∞i ix X,14,612211====∑∑==ni ini i xXx X例2证明 ,1∞∞≤≤Xn X X证明 因为 11max X x x x Xn i i p i i=<==∑=∞∞=≤=≤∑Xn x n x n x i ip ni i max 1所以,1∞∞≤≤Xn X X例3 已知矩阵⎥⎦⎤⎢⎣⎡-=2212A ,求矩阵A 的三种常用范数。
解 4max 31==∑=∞j ij i a A,∑===ni ij ja A 114max , 39)9)(4(36135228522822122122122===--=+-=--=-⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡-=λλλλλλλλA I A A A A T T例4 已知方程组⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛121212212321x x x a a a (1)写出解此方程组的雅可比法迭代公式 (2)证明当4>a 时,雅可比迭代法收敛 (3)取5=a ,TX )101,51,101()0(=,求出)2(X 。
解 (1)对3,2,1=i ,从第i 个方程解出i x ,得雅可比法迭代公式为:⎪⎪⎪⎩⎪⎪⎪⎨⎧=--=--=--=+++ ,1,0,)21(1)222(1)21(1)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1m x x a x x x a x x x a x m m m m m m m m m (2)当4>a 时,A 为严格对角占优矩阵,所以雅可比迭代法收敛。
第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。
[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%,求n x 的相对误差。
[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。
4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。
(1)*4*2*1x x x ++; [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。
数值计算方法思考题数值计算方法思考题第一章预篇1.什么是数值分析?它与数学科学和计算机的关系如何? 2.何谓算法?如何判断数值算法的优劣?3.列出科学计算中误差的三个来源,并说出截断误差与舍入误差的区别。
4.什么是绝对误差与相对误差?什么是近似数的有效数字?它与绝对误差和相对误差有何关系?5.什么是算法的稳定性?如何判断算法稳定?为什么不稳定算法不能使用? 6.判断如下命题是否正确:一个问题的病态性如何,与求解它的算法有关系。
无论问题是否病态,好的算法都会得到好的近似解。
解对数据的微小变化高度敏感是病态的。
高精度运算可以改善问题的病态性。
用一个稳定的算法计算良态问题一定会得到好的近似值。
用一个收敛的迭代法计算良态问题一定会得到好的近似值。
两个相近数相减必然会使有效数字损失。
计算机上将1000个数量级不同的数相加,不管次序如何结果都是一样的。
7.考虑二次代数方程的求解问题ax2 + bx + c = 0.下面的公式是熟知的bb24acx.2a与之等价地有x?对于2c?b?b?4ac2.a = 1,b = -100 000 000 ,c = 1应当如何选择算法?8.指数函数有著名的级数展开x2x3e?1?x2!3!x 如果对x 9.考虑数列xi, i = 1,…, n, 它的统计平均值定义为x?1n?xi xi?1 它的标准差2?1n2(xi?x)? n?1i?1??1 数学上它等价于1n222xinx n1i11 作为标准差的两种算法,你如何评价它们的得与失?第二章非线性方程求根1.判断如下命题是否正确:(a) 非线性方程的解通常不是唯一的;(b) Newton法的收敛阶高于割线法;(c) 任何方法的收敛阶都不可能高于Newton法; (d)Newton法总是比割线法更节省计算时间;(e) 如果函数的导数难于计算,则应当考虑选择割线法;(f) Newton法是有可能不收敛;(g) 考虑简单迭代法xk+1 = g(xk),其中x* = g(x*)。
第一章 绪论(12)1、设0>x ,x 的相对误差为δ,求x ln 的误差。
[解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x , 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%,求n x 的相对误差。
[解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε,相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字:1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字;0031.0*2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56*4=x 有5位有效数字;0.17*5⨯=x 有2位有效数字。
4、利用公式(3.3)求下列各近似值的误差限,其中*4*3*2*1,,,x x x x 均为第3题所给的数。
(1)*4*2*1x x x ++; [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε;(2)*3*2*1x x x ;[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε;(3)*4*2/x x 。
数值计算方法总复习第一章算法与误差第二章非线性方程求解第三章线性代数方程求解第四章函数插值与曲线拟合第五章数值积分与数值微分第六章常微分方程的数值解法Chap.1 (1)关于数值计算方法,What,特点一、关于《数值计算方法》数值计算方法是应用数学的一个分支,又称数值分析或计算方法,它是研究数字计算机求解各种数学问题的数值方法及其理论的一门科学,是程序设计和对数值结果进行分析的依据和基础。
应用计算机解决科学计算问题包括以下几个过程:提出实际问题;建立数学模型;选用数值计算方法;程序设计和上机计算。
可见数值计算方法是进行科学计算全过程的一个重要环节。
计算机计算的特点:(1)运算速度快;(2)但只能完成加、减、乘、除和一些逻辑运算。
所以,各种复杂的数学问题------→归结为四则运算------→编程指令。
把对数学问题的解法归纳为有加、减、乘、除等基本运算,并对运算顺序有完整而准确的描述的算法称为数值计算方法或简称数值算法。
研究各种算法和相关理论的一门课程。
§1.2 误差一、误差的来源数分为两类:精确数(准确数、真值);近似数/近似值。
1)模型误差或描述误差2)测量误差(观测误差)3)截断误差(方法误差)4)舍入误差(计算误差):数值计算关心的是截断误差(方法误差)和舍入误差(计算误差)二、误差限和有效数字1. 误差限的定义设Z 是准确值Z *的某个近似值,如果根据具体测量或计算的情况,可以事先估计出误差的绝对值不超过某个正数ε:即: |Z * - Z |≤ε则称ε为近似值的误差限。
或称在允许误差ε的情况下,结果Z 是“准确的”.2. 误差限和有效数字在表示一个近似数时,常常用到“有效数字”,有效数字和误差限都是用来定量表示误差的大小,且它们之间有对应关系。
有效数字的定义:设数x 的近似值m n x x x x 10.021*⨯= , 其中 xi 是0到9之间的任一个数,但x 1≠0,i=1,2,3…,n 正整数,m 整数,若nm *|x x |-⨯≤-1021 则称x *为x 的具有n 位有效数字的近似值,x *准确到第n 位,x1x2…xn 是x *的有效数字。