2016数值分析课件详解

数值分析电子课件工科研究生公共课程数学系列辽宁科技大学理学院2016年9月机动上页下页首页结束工科研究生公共课程数学系列第1章数值分析与科学计算引论内容提要：1.1数值分析研究对象与特点1.2数值计算的误差1.3误差定性分析与避免误差危害机动上页下页首页结束工科研究生公共课程数学系列1.1数值分析研究对象与特点一、数值分析研究对象计算机解决科学计算问题时经历的过程2实际问题模型设计算法设计问题的解上机计算程序设计求方程求根22x牛顿法)2(211KKkxxx程序设计解上机计算,417.1,5.1,1310xxx实例机动上页下页首页结束工科研究生公共课程数学系列数值分析的内容包括函数的数值逼近、数值微分与数值积分、非线性方程数值解、数值线性代数、常微和偏微数值解等。数值分析研究对象以及解决问题方法的广泛适用性，著名流行软件如Maple、Matlab、Mathematica等已将其绝大多数内容设计成函数，简单调用之后便可以得到运行结果。但由于实际问题的具体特征、复杂性,以及算法自身的适用范围决定了应用中必须选择、设计适合于自己特定问题的算法，因而掌握数值方法的思想和内容是至关重要的。本课程内容包括了微积分、代数、常微分方程的数值方法，必须掌握这几门课程的基础内容才能学好本门课程。机动上页下页首页结束工科研究生公共课程数学系列二、数值分析的特点•面向计算机，要根据计算机的特点提供切实可行的有效算法。•有可靠的理论分析，能任意逼近并达到精度要求，对近似算法要保证收敛性和数值稳定性，还要对误差进行分析。这些都是建立在数学理论的基础上，因此不应片面的将数值分析理解为各种数值方法的简单罗列和堆积。•要有好的计算复杂性，时间复杂性好是指节省时间，空间复杂性好是指节省存储量，这也是建立算法要研究的问题，它关系到算法能否在计算机上实现。•要有数值实验，即任何一个算法除了从理论上要满足上述三点外，还要通过数值实验证明是行之有效的。机动上页下页首页结束工科研究生公共课程数学系列三、数值分析的学习方法初学可能会觉得公式多，理论分析复杂。给出如下的几点学习方法。•认识建立算法和对每个算法进行理论分析是基本任务，主动适应公式多和讲究理论分析的特点。•注重各章节所研究算法的提出，掌握方法的基本原理和基本思想，要注意方法处理的技巧及其与计算机的结合。•理解每个算法建立的数学背景、数学原理和基本线索，而且对一些最基本的算法要非常熟悉。•要通过算例学习使用各种数值方法解决实际计算问题。•为掌握本课的内容，还应做一些理论分析和计算练习。机动上页下页首页结束工科研究生公共课程数学系列1.2数值计算的误差一、误差的来源与分类在运用数学方法解决实际问题的过程中，每一步都可能带来误差。1、模型误差在建立数学模型时，往往要忽视很多次要因素，把模型“简单化”，“理想化”，这时模型就与真实背景有了差距，即带入了误差。数学模型与实际问题之间出现的误差称为模型误差。2、观测误差（测量误差）数学模型中的已知参数，多数是通过测量得到。而测量过程受工具、方法、观察者的主观因素、不可预料的随机干扰等影响必然带入误差。机动上页下页首页结束工科研究生公共课程数学系列3、截断误差（方法误差）数学模型常难于直接求解，往往要用数值方法求近似解替代，这种简化带入误差称为方法误差或截断误差。()2()(Taylor)(0)(0)(0)()(0)1!2!!()nnnnfxfffPxfxxxnRx：可微函数用泰勒多项式近似代替，则数值方法的截断误差是泰勒余项例如。4、舍入误差计算机只能处理有限数位的小数运算，初始参数或中间结果都必须进行四舍五入运算，这必然产生舍入误差。3.141593.141590.0000026R：用近似代替，产生的误差例如机动上页下页首页结束工科研究生公共课程数学系列1-1xxxexxxe设为准定义绝对误差确值，为的一个近似值，称为近似值的，简称误差,记为。误差分析是一门比较艰深的专门学科。在数值分析中主要讨论截断误差及舍入误差。但一个训练有素的计算工作者，当发现计算结果与实际不符时，应当能诊断出误差的来源，并采取相应的措施加以改进，直至建议对模型进行修改。二、绝对误差、相对误差与有效数字1、绝对误差与绝对误差限误差是有量纲的量，量纲同x，它可正可负。误差一般无无法准确计算，只能根据测量或计算情况估计出它的误差绝对值的一个上界，这个上界称为近似值x*的误差限，记为ε*。它是正数，有量纲的。如用毫米刻度尺测量长度。误差限是0.5mm。机动上页下页首页结束工科研究生公共课程数学系列101,10001,10,1;1000,1xyxxxxxxxxyxy误差限的三种表示：（1）对于一般情形（2）（3）工程中常表示为：，有两个量则例如。2、相对误差与相对误差限2rexexxxxxe近似值的误差与准确值的比值称为近似值的，记作。相对误差无量定义纲，1-相对误差可正可负。机动上页下页首页结束工科研究生公共课程数学系列**,||||||||10%0.1%rrrryxxexxexxxxexxxyxyyyxx在计算中，由于真值总是不知道的，通常取。相对误差的绝对值上界叫相对误差限，记作即上例中与分别为与的对误差限，可见近似的程度比近似相的对误差限程度好。3、有效数字定义1-3如果近似值x*的误差限是某一位的半个单位，该位到x*的第一位非零数字共有n位，就说x*有n位有效数字.机动上页下页首页结束工科研究生公共课程数学系列5187.93250.037855518.0000332.7182818下列数按四舍五入原则写出下列各数的具有位有效数字的近似数？例如187.930.0378568.00002.7183解：1231.102156.4300.031xxx下列各数都是经过四舍五入原则得到的近似数指出它们有几个有效数字？再例123xxx解：五位五位二位机动上页下页首页结束工科研究生公共课程数学系列4、绝对误差，相对误差与有效数字的关系绝对误差与相对误差：由两者定义可知。*||rreexx，1(1)12110(1010)1102mnnmnxxaaanxx即，近似数表示为具有位有效数字，则绝对误差与有效数字：绝对误差不超过末位有效数字的半个单位。机动上页下页首页结束工科研究生公共课程数学系列1(1)121(1)1(1)110(1010)(1,2,,)090,110112110211)-(mllinrnrxxaaaailamxnaxxan设近似数表示为其中是到中的一个数字，为整数。若具有位有效数字，则其相对误差限为（）反之，若的相对误差限，定理则至少具有位有效数字。有效数字与相对误差限定理说明有效数位越多，相对误差限越小。定理也给出了相对误差限的求法。机动上页下页首页结束工科研究生公共课程数学系列1111111111110(1)10,0.510110;1021(1)10102(1)1102mmmnnrmmnrmnaxaxnxxaaxxxxaaxn证明：由（）可得当有位有效数字时反之，由故至少有位有效数字。定理证毕。1(1)1210(1010)mllxxaaa注：近似数表示为机动上页下页首页结束工科研究生公共课程数学系列1212()()xxxx两个近似数与，其误差限分别为及，则它们进行加、减、乘、除运算得到的误差限分别满足不等式2010.1%-1要使的近似值的相对误差限小于，至少要取几例位有效数字？三、数值运算的误差估计1、四则运算*1113311102204.4,4,4,0.12510100.1%2040.1%nrrnaan解：设取位有效数字，由定理，。由于知故只要取就有即只要对的近似值取位有效数字，其相对误差限就小于。机动上页下页首页结束工科研究生公共课程数学系列1212121221122112222()()();()()();()()()0xxxxxxxxxxxxxxxxxx，2201-2530010,VVRI若电压，电阻求电流并计算其误差限及例相对误差限。**22200.7333()300()=5()()=10()()()()2201030050.0411()900000.73330.0411()0.0411()6%0.7333rVIARVVRVRRVIRAIAI解：由题意知，，于是，（）所以机动上页下页首页结束工科研究生公共课程数学系列11080,-0.2-0.1llmddmllmddmsld已测得某场地长的值为，宽的值为已知，。试求面积的绝对误差限与相对误差限。例1-32()()()1100.1800.227()()()()270.31%8800rslddlmssslds解：由误差限乘法公式知，故2、函数误差当自变量有误差时计算函数值也产生误差，可以利用函数的泰勒展开式估计其误差界。机动上页下页首页结束工科研究生公共课程数学系列22()()()(())()()()=()()(),,2()()()()()()2()()()fxxxfxfxfxffxfxfxxxxxxxffxfxfxxxfxfxx设是一元可微函数，的近似值为，以近似，其误差界记作，由泰勒展开介于之间，取绝对值得假定与的比值不太大，可忽略的高阶项，于是可得计算函数的误差限()().fxfxx（）（）0ln1(ln),1(ln(x))()()1-4rxxxxxxxx设,的相对误差限为，求的误差限。解：有例即机动上页下页首页结束工科研究生公共课程数学系列121212121121212(,,),,,,((,,))()()((,,))((,,))(,,)()nnnnnkkknrrnnkfxxxxxxxxxffxxxxxfxxxfxxxfxxxfx当，是多元函数，，的近似值为，，则误差限，，相对误差限，，112()(,,)=110m-0.2-0.1-nkknxfxxxlldllmddmsld，已测得某场地长的值为，宽的值为d=80m,已知，。试求面积的绝对误（同例13差限与相对误差。）限例1-52,,,()()()800.21100.127()()()27()0.31%8800rssslddlldsssldldmssslds解：因由误差限函数公式故机动上页下页首页结束工科研究生公共课程数学系列210.12sgtgtsts设，假定是准确的，而对的测量有的误差，证明当增加时，的绝对误差增加，而相对误差却减少。例1-62210.1,(),2()()0.1()0.10.2()1()2rtsgtsgtsgttgtsgtsstgttt证明：由题意知（）再由函数绝对误差限公式得，故，相对误差限由于增加时，其近似值也增加，因此结论成立。机动上页下页首页结束工科研究生公共课程数学系列1.3误差定性分析与避免误差危害一、算法的稳定性用一个算法进行计算，由于初始数据误差在计算中传播使计算结果误差增长很快，