计算机数值方法教案.doc

第O章绪论一、教学设计1教学内容：数值计算方法这门课程的形成背景及主要研究内容、研究方法、主要特点；算法的有关概念及要求；误差的来源、意义、及其有关概念。数值计算中应注意的一些问题。2重点难点：算法设计及其表达法；误差的基本概念。数值计算中应注意的一些问题。3教学目标：了解数值计算方法的基本概念；掌握误差的基本概念：误差、相对误差、误差限、相对误差限、有效数字；理解有效数字与误差的关系。学会选用相对较好的数值计算方法。4教学方法：介绍与讨论二、教学过程1。1引论1课程简介：数学科学的一个分支，它研究数值计算方法的设计、分析和有关的理论基础与软件实现问题。另外，有一个较常用的名词“数值分析”，其包含的内容属于计算数学的一个部分。2历史沿革：数学最初导源于计算，计算曾经是古代数学的最重要的组成部分。各个时期的大数学家，在发展基础数学的同时也都对计算方法作出了重要贡献。例如：牛顿、拉格朗日、高斯、秦九韶等。直到20世纪40年代，由于技术手段和计算工具条件的不足，发展比较缓慢，作用也比较有限。3计算方法的形成：20世纪下半叶，计算机极大地扩展了数学的应用范围与能力。如：天气预报计算能力的提高与所用计算方法的效能密切相关。以原来分散在数学各分支的计算方法为基础的一门新的数学科学“计算数学”开始形成并迅速发展。4作用与意义：科学实验、科学理论、科学计算已成为人类进行科学活动的三大方法。这是伽利略、牛顿以来在科学方法论方面取得的重大进展。5计算方法的任务：将计算机不能直接计算的运算，化成在计算机上可执行的运算。例：，针对数值问题研究可在计算机上执行且行之有效的新系列计算公式。例：解线性方程组，已有Cram法则，但不可行。(几十万年)误差分析，即研究数值问题的性态和数值方法的稳定性。6计算机数值方法的研究对象：(与科学计算有关的数学问题是多种多样的，最基本类型有：)利用计算机解决科学计算问题的全过程大致如下：实际问题构造数学模型设计数值计算方法程序设计上机求出结果回到实际问题。数学模型举例：例1：鸡兔同笼：（共10只，34只脚）导致方程组；例2：曲边梯形的面积。相应地，本课程主要研究的数值问题有：函数的插值与逼近方法；微分与积分计算方法；线性方程组与非线性方程组计算方法；微分方程数值解等。7主要特点既有纯数学的高度抽象性与严密科学性的特点，同时又具有应用广泛性与数值试验的高度技术性。（要求先掌握基本数学知识，以及计算机的基本操作）8学习目的：学习一些常用的数值方法，掌握数值方法的基本理论，为进一步研究新算法奠定基础。初步掌握一种软件包：Mathematic,Matlab等的使用方法。9参考书目：1袁蔚平等编.计算方法与实习.南京 ：东南大学出版社,2000年7月2李庆杨等编.数值分析 .武汉：华中工学院出版社,1982年1月3葛福生编.数值计算方法.南京 ：河海大学出版社,1996年4月1。2数值问题与数值算法1数值问题：指输入数据与输出数据之间函数关系的一个确定而无歧义的描述。例：求二次方程的根，可算作一个数值问题；求常微分方程的解，却不能称作数值问题，需离散化。2数值方法：求解数值问题的计算机上可执行的系列计算公式。例1：Cram法则，Gauss消去法例2：求根公式3数值算法：指有步骤地完成解数值问题的过程，数值方法是它的前提和基础，它是数值方法的具体化。具备以下四个特性：目的性；确定性；可执行性；有穷性。(有别于常规的思维)算法设计的目的：可靠性好、计算精度高；计算复杂性好；为程序设计作准备。4算法设计及其表达法表达方法：自然语言法和图示法。例：通过二次方程求根的例子，说明数值方法与数值算法的区别，并演示算法常用的表达方法之一：自然语言法（图示法不加介绍）。(首先要选择数值方法：公式法或迭代法)主要步骤：（阅读课本后，要求自己解释）1输入数据2若怎样？（若，否则）3若，计算，若怎样？（）若怎样？（）若怎样？4输出2误差2-1 误差的基本概念1误差来源及种类：模型误差（忽略次要因素）观测误差（测量工具的限制）截断误差（有限代替无限，如Taylor展开）舍入误差（计算机字长位数有限），主要讨论。2举例说明误差分析的重要性：计算。递推公式（A）：，；，；递推公式（B）：，计算。递推公式（A）：，；，；递推公式（B）：，或，故，取中值即得计算结果见表nIn(A)In(A)In(B)In00.6321205590.63210.6321205590.632120558810.3678794410.36790.3678794410.367879441220.2642411180.26420.2642411180.264241117730.2072766460.20740.2072766470.20727664740.1708934160.17040.1708934120.170893411950.145532920.14800.145532940.145532940660.126802480.11200.12680235660.126802356670.112382640.21600.1123835040.112383504180.10093888-0.72800.10093196740.100931967490.091550087.5520.0916122930.091612293100.0844991999-74.520.08387707010.08387707014110.0705088012820.70.07735222890.07735222896120.153894386-98480.07177325360.0717732538913-1.00062702128020.06694770260.06694770317140.06273216410.0627321639150.05901753840.05901754088160.05571938500.0557193459170.05277045530.0527711191180.05013180450.05011985496190.04749571470.0477227558200.03256855810.0455448841计算结果见表，差之毫厘，失之千里，试分析其中原因（误差被放大或缩小）。3定义3.1设为准确值，为的一个近似值，称为近似值的绝对误差，简称误差。（绝对）误差限：，或例：4定义3.2近似值的误差与准确值之比，称为近似值的相对误差。（相对）误差限：。但实际中常以下式计算：，相应地。例：；例：已知，其近似值为，求的绝对误差限和相对误差限。考察常用的四舍五入(为什么要四舍五入？) 所引起的误差，不超过某一位数字的半个单位（个位为1，十位为10，）5定义3.3如果的误差绝对值不超过某一位数字的半个单位，若该位数字到的第一位非零数字共有位，则称用近似时具有位有效数字，简称有位有效数字。例：分别具有6,5位有效数字。如何描述有效数字？（一般情况下在计算机中数往往规格化，故有必要考察规格化数。）6定义3.4如果（其中）是的近似值，且满足不等式，则称有位有效数字。例：设，它的两个近似值和分别有3,4位有效数字。一般地，有效数字位数多，相对误差小，但上例例外。下面讨论相对误差与有效数字的关系。7定理3.1设（其中）是的近似值，（1）如果有位有效数字，则的相对误差限为；（2）若的相对误差限为，则至少有位有效数字。例：设（1），（2），分别求的有效数字位数与相对误差限。（用此例说明定理3.1的不唯一）例：要使的近似值的相对误差限小于0.1%，要取几位有效数字？（4位，4.472）（用定理3.1的（1）来解，但并不能保证最好的结果：例：，2,3,4位有效数字都可以使相对误差小于0.1%）例：已知作为的近似值有位有效数字，问作为的近似值有几位有效数字？3-2 数值运算的误差估计（不加证明）1（可先讲一元的）设给定多元函数，且设为的近似值，以作为的近似值，其误差分析可利用Taylor展开，其绝对误差绝对误差限为：相对误差为：或：相对误差限为：或2基本运算的误差（1），（2），（3），（4），（5），例：见书本P28，T10，应用以上理论作分析。3浮点数运算的误差设，而机器字长为，则在机器中不同阶的数在运算时，需对阶，这也会导致误差。例：在只有四位字长的机器上作下列运算，结果如何？ 1000+0.1+0.1+0.1+0.1+0.1+0.1+(11项)(Matlab: digits(4);x=vpa(0.1),1000+x+x.) 0.1+0.1+0.1+0.1+1000（以下略）记为，其绝对误差和相对误差分别为(相对误差只与机器字长有关)四则运算的误差在一定条件下有以下误差估计：，连加和连加的误差例：例：0.368 467 6 + 1070.632 754 4 + 0.493 200 1 + 0.480 010 0 = 1070.632 754 4对于加法，先参加运算的数，在计算和中引起的误差也较大，故应先让绝对值小的数参加运算；对于乘法，误差与数参加运算的次序无关。3-3数值方法的稳定性与算法设计原则1数值稳定：对于一个数值方法，若其计算结果受计算过程舍入误差影响较小，则称此数值方法是数值稳定的；否则称它是不稳定的。如：前例中方法（A）和（B）。（可作分析）引例：Hilbert矩阵作为系数的方程组。求线性方程组的计算解（1）精确解为：；（2）取3位有效数字，用Gauss消去法求解，其解：；Matlab:A=,b=,x=Ab;若中间结果也取3位有效数字，则解的过程如下，解为：1 0.5 0.333 1.830.5 0.333 0.25 1.080.333 0.25 0.2 0.7831 0.5 0.333 1.830 0.083 0.0835 0.1650 0.0835 0.0891 0.1741 0 -0.17 0.8360 0.083 0.0835 0.1650 0 0.0051 0.008011 0 0 1.10 0.083 0 0.03390 0 0.0051 0.008011 0 0 1.10 1 0 0.4080 0 1 1.57（3）同样取3位有效数字，用Gauss消去法求解方程组得到计算解为，即为精确解。2 -1 3 14 2 5 42 0 2 62 -1 3 10 4 -1 20 1 -1 52 0 2.75 1.50 4 -1 20 0 -0.75 4.52 0 0 180 4 0 -40 0 -0.75 4.51 0 0 90 1 0 -10 0 1 -6数值问题的计算解的精度，不但与数值方法的稳定性有关，还与数值问题的性态好坏有关。引出概念2病态问题：如果问题的输入数据的微小误差，引起输出结果的很大误差，（即输出对输入的扰动很敏感）则称这种问题为病态问题或坏条件问题（ill-posed）；否则称为良态问题。3四则运算中的稳定性问题（1）防止大数吃小数由于计算机中数的位数有限（假设7位），两数相加需对阶，于是会发生大数吃小数的现象。例：0.368 467 6 + 1070.632 754 4 + 0.493 200 1 + 0.480 010 0 = 1070.632 754 4方法：将小数放在前面相加。0.368 467 6 + 0.493 200 1 + 0.480 010 0 + 1070.632 754 4= 1.34167770000000 + 1070.632 754 4 = 1070.632 754 5（2）要避免两个相近数相减（可用此原则优化一元二次方程的求根公式，见同济P6：例5）例：都有五位有效数字，但却只有两位有效数字。方法：改变计算公式。例：，计算。取，直接计算得结果为，但若取，例：当很接近时，例：当很大时，例：，其中，试在8位字长的机器上求解此一元二次方程。（3）要避免除数