数值分析的一章.docVIP

绪论 数学是研究数与形的科学。其中研究怎样利用手指、算盘、算尺、计算器、计算机等工具，来求出数学问题数值解答的学问，就是计算方法，或称计算数学。它是数学中最古老的部分，但只是在电子计算机出现以后，人们获得了高速度、自动化的计算工具，才为众多浩繁的数值问题的解决展现了光明的前景。从此，科学研究与工程设计的手段，发生了由模型试验向数值计算的巨大转变。 在农业科学研究中，数值计算方法已经成为不可缺少的有力工具。我校开设数值计算方法课已有二十几年的历史，许多研究生通过学习，在毕业论文中引用了数值计算方法解决应用问题，提高了论文水平，也有许多在职教师和科研人员学习这门课程后，将数值方法引用到科研课题中，取得了较好的结果。本讲义旨在为农业院校高年级学生、研究生提供一本学习数值计算方法的入门工具。 §1 数值计算方法的意义、内容与方法 由计算机解决科学计算问题需经历几个过程：由实际问题建立数学模型，根据数学模型提出求解的数值计算方法，编出程序，上机求出结果。通过以上过程，可以看出：数值计算方法是计算机、数学和应用科学之间的桥梁，是程序设计对数值结果进行分析的依据，是用计算机进行科学计算全过程的一个重要环节。 数值计算方法的研究对象就是构造计算机能用的算法。 §2 算法 一、算法的概念 计算机实质上只会做加减乘除等基本运算。研究怎样通过计算机所执行的基本运算，求得各类问题的数值解或近似数值解，就是计算数学的根本课题。 由基本运算及运算顺序的规定所构成的完整的解析步骤，称为算法。 描述算法可以有不同的方式。例如，可以用日常语言和数学语言加以叙述，也可以借助形式语言（算法语言）给出精确的说明，也可以用框图直观地显示算法的全貌。下面举个简单的例子。 例：求解二元一次联立方程组 EMBED Equation.3 这个方程组的行列式解法可表述如下： 首先判别 是否为零，存在两种可能： （1）如果，则令计算机计算 然后输出计算的结果x1，x2。 （2）如果D= 0，则或是无解，或有无穷多组解，这是奇异的情形。 下面用框图（也称流程图）来形象地描述上面的算法。 开 始 开 始 输 入 输 入 D=a11 D=a11a22-a12 Y D=0 D=0 N 输出无解信息 输出无解信息 输出 x1 输出 x1, x2 结 束 结 束 框图 这里我们使用了两种形式的框，一种是矩形框 ，称为叙述框，计算公式就填在这种框内；另一种是棱形框 ，称为判断框，表示算法的判断检查部分。 二、算法的优劣 我们知道，计算机的特点是运算速度快，存贮的信息量大，并能自动完成极其复杂的计算过程。计算机功能虽然很强，但是否可降低对算法的要求呢？许多事实说明，如果算法选择不当，计算机的利用率就得不到充分发挥，有时甚至不能得到满意的解答。一个好的算法，要求有以下几个优点： 1．计算量小 例如求解线性方程组，原则上可以用行列式解法。如解n阶方程组，要计算n + 1个n阶行列式的值，总共需要做n! (n - 1) (n + 1) 次乘法运算。设n = 20，采用每秒一亿次的计算机，要连续计算三十万年才能完成。这当然是没有实际意义的。因此，计算量大小是衡量算法优劣的一项重要标准。 在估计计算量时，我们将区分主次抓住计算过程中费时较多的环节。比如，由于加减操作的机器时间比乘除少得多，对和式可以忽略加法而只统计乘法的次数。又如算式中需要多次计算函数值，每求一次 时，通常需要进行多次加减乘除运算，因此对算式只要统计调用f 的次数就可以了。 2．存贮量小 计算程序所占用的工作单元的数量称为算法的内存量。尽管计算机能贮存大量信息，但计算大型算题时有些微机仍不能使用。因此，尽量节约存贮量有经济价值，是衡量算法质量的又一标准。 3．逻辑结构简单 设计算法时应考虑的另一个因素是逻辑结构问题。虽然计算机能够执行极其复杂的计算方案，但是计算方案的每个细节都需要编程人员制定。因此算法的逻辑结构应尽量简单，才能使编制程度、维修程序和使用程序比较方便。 由此可见，虽然计算机是一种强有力的计算工具，但不能因此忽略对算法的研究。应该以计算量大小、存贮量多少、逻辑结构是否简单作为评定算法优劣的标准。 §3 数值计算中的误差 用数值方法解决科学研究或工程技术中的实际问题，一般来说，产生误差是不可避免的，根本不存在绝对的严格和精确。但是我们可以认识误差，从而控制误差，使之局限于最小（或尽量小）的范围内。 一、误差的来源 运用数学工具解决实际问题，可能产生的误差主要