清华大学版数值分析绪论第1章
- 格式:doc
- 大小:833.50 KB
- 文档页数:25
教材:关治、陆金甫数值方法清华大学出版社,2006第1章引论§1 数值分析研究对象与特点(I)数学与科学、工程技术有非常密切关系,并相互影响。
科学与工程技术领域中的问题通过简化、抽象建立数学模型。
对数学模型的研究和求解,应用于科学与工程实践。
数学模型的建立与研究是数学研究的重要任务。
例如,设有一个质量为m 的质点作直线运动(设其在x 轴上运动),其坐标用x 表示。
在质点运动过程中,其坐标x 随时间t 而变动,要知道质点如何运动,就要知道x 对时间t 的依赖关系。
假定运动是在力F 的作用下进行的,而力F 又与时间t ,质量的位置x ,速度v dx dt =有关,即(,,)dx F f t x dt=。
由Newton 第二定律,22v ,d d x F ma a dt dt===,即在时刻t 有 22(,,)d x dx m F t x dt dt= 为使问题定解,还需加上初始条件 00()x t x =00v t t dxdt ==这样的初值问题,其解存在,唯一,连续依赖于初始数据…是数学工作者要研究和解决的问题。
作为科学与工程技术工作者,仅知道其解的存在,唯一…是不够的,还必须知道其数量是多少。
很多线性与非线性问题很难用解析方法来求得,甚至于看来简单的问题,也难于解析求解,如21sin xdxx不能用解析方法求得。
很多需要数字结果的科学与工程问题,得到数学模型后,必须采用数值方法来求解。
数值分析则是研究数值方法求解科学与工程问题的一门学科,主要是构造适合于不同类型问题的数值方法,分析方法的精度、稳定性、收敛性等一系列理论与实际问题。
改进计算方法,创造新的数值方法是数值分析很重要的任务。
数值分析具有两个显明特点:1.实践性数值求解的总是来源于科学与工程实践,解决之后又为科学与工程技术服务。
2.与计算机的密切相关虽然数值分析早于计算机的出现,但是数值分析的飞速发展是计算机出现之后,并随计算机发展而迅速发展。
第一章绪论e In X* =In X * -Inx :丄e*X*进而有;(In X *):2. 设X 的相对误差为2% ,求X n 的相对误差。
解:设f(χZ ,则函数的条件数为Cp=l fX+n _1X nχ I Xn n又;r ((X*) n) C P 7(X *)且 e r (χ*)为 2.7((χ*)n) 0.02 n3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指* * * * *出它们是几位有效数字: X 1 =1.1021, χ2 =0.031, χ3 =385.6, χ4 = 56.430,x 5 = 7".0.. *解:X I -1.1021是五位有效数字;X 2 = 0.031是二位有效数字;X 3 =385.6是四位有效数字;X 4 =56.430是五位有效数字;X 5 =7 1.0.是二位有效数字。
4. 利用公式(2.3)求下列各近似值的误差限: (1) X 1 X 2 X 4,(2) X 1 X 2X 3 ,(3) X 2 /X 4 .其中χl ,x 2,x 3,X 4均为第3题所给的数。
1设X 0, x 的相对误差为 解:近似值X*的相对误差为 、:,求InX 的误差。
e* X* -X而InX 的误差为 又 f '(χ) =nx n 」 C P解:* 1 4;(x 1) 102* 1 3 ;(x 2) 10 2* 1 1;(x 3) 10* 1 3;(x 4) 102* 1 1;(x 5) 102(1) ;(x ; x ; x *)* * *=;(%) ;(x 2) *x 4)1 A 12 1 j310 10 102 2 2 -1.05 10J 3* * *(2) S(X I X 2X 3)* * * * * * ** * =X1X 2 £(X 3)+ X 2X 3 ^(X J + X 1X 3 E (X 2):0.215 ⑶;(x 2/x ;)* Il * * I * X 2 E(X 4) + X 4 &(X 2)全 Γ"2X 41-3 1 30.031 10 56.430 10= ______________________ 256.430X56.430-10 54 3解:球体体积为V R3则何种函数的条件数为1.1021 0.031 11θ' 2 + 0.031X385.6 x 1><10* 2 +∣ 1.1021 X 385.6卜-×1^35计算球体积要使相对误差限为 1 ,问度量半径R 时允许的相对误差限是多少?C P 愕': C P “(R*) 9(R*)又γ(V*) -11故度量半径R 时允许的相对误差限为 ;r (R*) 1 : 0.3331 ____6.设 Y 0 =28,按递推公式 Yn =Ynd- ------- : 783 (n=1,2,…)100计算到Y oo 。
第一章绪论上世纪中叶诞生的计算机给科学、工程技术和人类的社会生活带来一场新的革命。
它使科学计算平行于理论分析和实验研究,成为人类探索未知科学领域和进行大型工程设计的第三种方法和手段。
在独创性工作的先行性研究中,科学计算更有突出的作用。
在今天,熟练地运用电子计算机进行科学计算,已成为科学工作者的一项基本技能。
然而,科学计算并不是计算机本身的自然产物,而是数学与计算机结合的结果,它的核心内容是以现代化的计算机及数学软件为工具,以数学模型为基础进行模拟研究。
近年来,它同时也成为数学科学本身发展的源泉和途径之一。
1. 数值分析的研究对象与特点数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。
一般地说,用计算机解决科学计算问题,首先需要针对实际问题提炼出相应的数学模型,然后为解决数学模型设计出数值计算方法,经过程序设计之后上机计算,求出数值结果,再由实验来检验。
概括为由实际问题的提出到上机求得问题的解答的整个过程都可看作是应用数学的任务。
如果细分的话,由实际问题应用有关科学知识和数学理论建立数学模型这一过程,通常作为应用数学的任务,而根据数学模型提出求解的数值计算方法直到编出程序上机计算出结果,这一过程则是计算数学的任务,即数值分析研究的对象。
因此,数值分析是寻求数学问题近似解的方法、过程及其理论的一个数学分支。
它以纯数学作为基础,但却不完全像纯数学那样只研究数学本身的理论,而是着重研究数学问题求解的数值方法及与此有关的理论,包括方法的收敛性,稳定性及误差分析;还要根据计算机的特点研究计算时间最省(或计算费用最省)的计算方法。
有的方法在理论上虽然还不够完善与严密,但通过对比分析,实际计算和实践检验等手段,被证明是行之有效的方法也可采用。
因此数值分析既有纯数学高度抽象性与严密科学性的特点,又有应用的广泛性与实际试验的高度技术性的特点,是一门与使用计算机密切结合的实用性很强的数学课程。
在电子计算机成为数值计算机的主要工具以后,则要求研究适合计算机使用的,满足精确要求,计算时间省的有效算法及其相关的理论。
在实现这些算法时往往还要根据计算机的容量、字长、速度等指标,研究具体的求解步骤和程序设计技巧。
有的方法在理论上虽还不够严格,但通过实际计算、对比分析等手段,证明是行之有效的方法,也应采用。
这些就是数值分析具有的特点,概括起来有四点:第一,面向计算机,要根据计算机特点提供切实可行的有效算法。
即算法只能包括加、减、乘、除运算和逻辑运算,这些运算是计算机能直接处理的运算。
.第二,有可靠的理论分析,能任意逼近并达到精确要求,对近似算法要保证收敛性和数值稳定性,还要对误差进行分析。
这些都建立在相应数学理论的基础上。
第三,要有好的计算复杂性,时间复杂性好是指节省时间,空间复杂性好是指节省存储量,这也是建立算法要研究的问题,他关系到算法能否在计算机上实现。
第四,要有数值试验,即任何一个算法除了从理论上要满足上述三点外,还要通过数值试验证明是行之有效的。
根据“数值分析”课程的特点,学习是我们首先要注意掌握方法的基本原理和思想,要注意方法处理的技巧及其与计算机的结合,要重视误差分析、收敛性及稳定性的基本理论;其次,要通过例子,学习使用各种数值方法解决实际计算问题;最后,为了掌握本课的内容,还应作一定数量的理论分析与计算练习。
由于本课内容包括了微积分、代数、常微分方程的数值方法,读者必须掌握这几门课的基本内容才能学好这门课。
2 数值计算的误差2.1 误差来源与分类我们算出数学模型的近似解和一个物理量的真的值往往不相等,它们之差称为误差。
用电子计算机进行解决实际问题的数值计算,误差是不可避免的。
引起数值结果中的误差的原因是多方面,通常来自固有误差和计算误差,如下面所示: 模型误差:用数学方法解决实际问题,首先必须建立该问题的数学模型。
即把实际问题经过抽象,忽略一些次要的因素,简化成一个确定的数学问题。
数学模型只是对实际问题的一种近似、一种粗糙的描述,因而它与实际问题或客观现象之间必然存在误差,这种误差称为“模型误差”。
这样的误差常常是可以忽略不计的。
例如在经典力学问题中,我们常常忽略相对论效应。
但是,如果这种误差不可忽略,说明数学模型选择得不好。
那么不论数值计算多么精确,其结果都将存在不可忽略的误差。
观测误差:数学问题中总包含一些参量(或物理量,如电压、电流、温度、长度等),它们的值(输入数据)往往是由观测得到的。
而观测的误差是难以避免的,由此产生的误差称为“观测误差”。
由于观测误差通常具有随机的性质,所以想用分析的方法来估计它们的影响常常是一件非常困难的事。
由于固有误差的产生往往涉及各专业知识及实验手段,所以它不是数值分析所研究的模型误差观测误差误差来源计算误差固有误差截断误差舍入误差内容。
截断误差:数值计算的本质就是用有限的过程、离散的数据来近似(刻画)无限的过程、连续的量。
这样对一些具有连续量、无限计算量的数学问题的求解的过程中,我们不得不对一些连续的量进行有限的离散化近似,不得截去无限的计算过程,只进行有限次的计算。
如求一个收敛的无穷级数之和,必须截去该级数后面的无穷多项,而用前面有限项的部分来近似代替,于是产生了有限过程代替无限过程的误差,称为“截断误差”,这是计算方法本身所出现的误差,所以也称为“方法误差”。
例如2462(1)cos 124!6!(2)!n nx x x x x n -=-+-+++ , 当x 很小时,可以用212x -作为cos x 近似值。
由交错级数判敛的莱布尼兹(Leibniz )准则,它的截断误差的绝对值不超过424x 。
有限过程代替无限过程的误差和计算量取决于其方法的收敛性及收敛速度。
舍入误差:计算中遇到的数据可能位数很多或是无穷小数,如1.41421356,=受机器字长的限制,无穷小数和位数很多的数必须舍入成一定的位数(机器字长)。
舍入方法:a) 舍入机。
采用四舍五入的办法,如将1.41421356 舍为1.4142136;b) 截断机。
在八位字长的截断机里取成1.4142135。
这样产生的误差称为“舍入误差”。
少量的舍入误差是微不足道的,但在计算机是作了成千上万次运算后,舍入误差的累积有时可能是十分惊人的。
它取决于方法的稳定性。
如果方法能够累积大量的误差,此算法是不稳定的,反之为稳定算法。
如算法1:1n n y a ny -=+。
显然是不稳定的;算法2:11n n y a y n-=+, 研究计算结果的误差是否满足精度要求就是误差估计问题,本书主要讨论算法的截断误差与舍入误差,而截断误差将结合具体算法讨论。
为分析数值运算的舍入误差,先要对误差基本概念做简单介绍。
2.2 绝对误差与相对误差定义1 设x 为准确值,*x 为x 的一个近似值,称**e x x =-为近似值的绝对误差,简称误差。
注意这样定义的误差*e 可正可负,当绝对误差为正时近似值偏大,叫强(赢)近似值;当绝对误差为负时近似值偏小,叫弱(亏)近似值。
通常我们不能算出准确值x ,当然也不能算出误差*e 的准确值,只能根据测量工具或计算情况估计出误差的绝对值不超过某正数*ε,也就是误差绝对值的一个上界。
*ε叫做近似值的误差限,它总是正数。
一般情形**x x ε-≤,即****x x x εε-≤≤+。
这个不等式有时也表示为**x x ε=±。
我们把近似值的误差*e 与准确值x 的比值**e x x x x -= 称为近似值*x 的相对误差,记作*r e 。
在实际计算中,由于真值x 总是不知道的,通常取*****r e x x e x x -==作为*x 的相对误差,条件是***r e e x =较小,此时*****2**2*******()()(/)()1(/)e e e x x e e x x x x x x x e e x --===-- 是*r ε的平方项级,故可忽略不计。
相对误差也可正可负,它的绝对值上界叫做相对误差限,记作*r ε,即***r x εε=。
注1:1)绝对误差是有量纲的量,它与所研究问题的背景有关,因此我们不能单单从绝对误差值的大小来判断计算结果的精度,还必须考虑到实际问题的应用背景;2)通过引入相对误差,我们可比较不同算法,不同应用问题的计算精度,这一点是无法通过比较绝对误差做到的,因为不同值的东西,在量上是无法比较的。
正如我们不能通过比较两个商品的使用价值来判定商品的真贵,只能通过抽象的价值。
2.3 有效数字有效数字是近似值的一种表示法,它既能表示近似值的大小,又能表示其精确程度。
定义2 若近似值*x 的误差限是某一位的半个单位,该位到*x 的第一位非零数字共有n 位,则称近似值有n 位有效数字。
在科学记数法中,将近似值*x 写成规格化形式为120.10m n x a a a =±⨯ (1-1)其中m 为整数,10,(2,)i a a i ≠= 为0到9之间的整数。
按照定义2,近似值*x 有n 位有效数字当且仅当*1102m n x x --≤⨯ (1-2) 因此在m 相同的情形下,n 越大则误差越小,亦即一个近似值的有效位数越多其误差限越小。
例1 按四舍五入原则写出下列各数具有5位有效数字的近似数:187.9325,0.03785551, 8.000033, 2.7182818。
按定义,上述各数具有5位有效数字的近似数分别是187.93, 0.037856, 8.0000, 2.7183 注意8.000033x =的5位有效数字的近似数是8.0000而不是8,因为8只有1位有效数字。
例2 重力常数g ,如果以2/m s 为单位,120.98010/g m s ≈⨯;若以2/k m s 为单位,220.98010/g km s -≈⨯,它们都具有3位有效数字,因为按第一种写法213119.80101022g ---≤⨯=⨯, 按第二种写法523110.00980101022g ----≤⨯=⨯, 他们虽然写法不同,但都具有3位有效数字。
至于绝对误差限,由于单位不同结果也不同,*22*52121110/,10/22m s km s εε--=⨯=⨯,而相对误差都是*0.005/9.800.000005/0.00980r ε==。
例2说明有效位数与小数点后有多少位数有关。
然而,从(2-2)可以得到具有n 位有效数字的近似数*x ,其绝对误差限为*1102m n ε-=⨯,在m 相同的情况下,n 越大则10m n -越小,故有效位数越多,绝对误差限越小。
关于一个近似数的有效位数与其相对误差的关系,有下面的定理定理1.1 设近似数*x 具有规格化形式(1-1), (1)若*x 具有n 位有效数字,则其相对误差限为*111102n r a ε-+≤⨯ (1-3) (2) 如果*111102(1)n ra ε-+≤⨯+ (1-4) 则*x 至少具有n 位有效数字。