计算方法word教案.doc

计算方法 第一章 绪 论 §1.0 引 言 §1.1 数值算法概论 (1) 计算方法的研究内容、对象与特点 (2) 基本求解步骤 §1.2 预备知识、误差 (1) 误差的来源 (2) 误差分析、数值稳定性的分析和说明 (3) 误差的基本概念——绝对误差 相对误差 有效数字 (4) 数值算法 §1.0 引 言 现代科学的三个重要组成部分: 科学理论, 科学实验, 科学计算。 它们相辅相成，互相独立，可以互相补充又都不可缺少，作为三种科学研究手段之一的科学计算是一门工具性、方法性、边缘性的新学科，发展迅速，它的物质基础是计算机(包括其软硬件系统)，其理论基础主要是计算数学。 科学计算的核心内容是以现代化计算机以及数学软件为工具，以数学模型为基础进行模拟研究。 出现了形如:计算数学,计算物理学,计算力学,计算化学, 计算生物学,计算地质学, 计算经济学等许多新学科及其发展。并已形成一系列专门研究数学问题的数值解法的算法软件，如目前流行的数学软件主要有以下几种： 符号运算软件： Mathematica, Maple 矩阵处理软件： Matlab Matlab简介 统计处理软件： SAS, Spss, Origin 数学CAD软件： MathCAD 等功能强大的著名数学软件。 §1.1 数值算法概论 §1.1.1 计算方法的研究内容、对象与特点 §1.1.2 基本求解步骤 说明： (1) 数学模型是通过科学实验或者观察分析一系列数据后，用数学作为工具近似地描述客观事物的一种数学表达式。 在数学模型中，往往包含了若干参量如物体比重、阻力系数、热交换系数等，这些物理参数通常由实验仪器测得，根据仪器的精密程度，物理参数的确定也会产生一定的误差。 (2) 在建立了数学模型之后，并不能立刻用计算机直接求解，还必须寻找用计算机计算这些数学模型的数值方法，即将数学模型中的连续变量离散化，转化成一系列相应的算法步骤，编制出正确的计算程序，再上机计算得出满意的数值结果。 (3) 算法：从给定的已知量出发，经过有限次四则运算及规定的运算顺序，最后求出未知量的数值解，这样构成的完整计算步骤称为算法。 评价算法的两个主要标准：计算速度和计算精度，此外，还有计算存贮量等。 一个面向计算机,计算复杂性好,又有可靠理论分析的算法就是一个好算法. 计算复杂性是算法好坏的标志，它包括时间复杂性(指计算时间多少)和空间复杂性(指占用存储单元多少)。对很多数值问题使用不同算法，其计算复杂性将会大不一样，例如对20阶的线性方程组若用代数中的Cramer法则作为算法求解，其乘除法运算次数需要，若用每秒运算1亿次的计算机计算也要30万年，这是无法实现的，而用计算方法中介绍的Gauss消元法求解，其乘除法运算次数只需3 060次，这说明选择算法的重要性。当然有很多数值方法事先不可能知道其计算量，故对数值方法除理论分析外，还必须通过数值试验检验其计算复杂性。作为基本要求希望读者能适当做一些计算机上的数值试验，对加深算法的理解是极有好处的。 例1.1：计算多项式的值。 算法1 算法例2：计算n次多项式的值。 算法1247) (又称为Horner算法(1819)) 计算。 说明：需乘法n次，加法n次，存储单元n+3个。 上述秦九韶算法的结构是递归的，它通过一次式的反复计算，逐步降低多项式的次数，直到归结为零次式为止。若以多项式的次数（或项数）定义为求值问题的规模，则秦九韶算法的特点是，在递归计算的过程中问题的规模逐次减1。 例：§1.2 预备知识和误差 §1.2.1 误差的来源 实际问题(建立数学模型(研究计算方法(编程上机计算(解结果 模型误差: 在建立数学模型过程中,不可能将所有因素均考虑,必然要进行必要的简化,这就带来了与实际问题的误差。 测量误差: 测量已知参数时,数据带来的误差。 截断误差: 在设计算法时,必然要近似处理,寻求一些简化。 舍入误差: 计算机的字长是有限的,每一步运算均需四舍五入，由此产出的误差称舍入误差。如：π、1／3，……取小数点8位、16位。 [截断误差的实例] 例1.4： 当很小时，可用作为的近似值，其截断误差小于。 例1.5： 已知： 求的近似值，并估计误差。 分析：对函数用Taylor展开，用多项式 近似代替，则数值方法的截断误差为 。 解：利用展开式的前三项，取n=2， 截断误差为： 数值计算方法主要讨论截断误差和舍入误差的影响,不讨论模型误差和测量误差。 §1.2.2 误差分析的重要性以及数值稳定性 一个数值方法进行计算时，由于原始数据有误差，在计算中这些误差会传