数值计算的可靠性.pdf

athematics Education 数学教育 数值计算的可靠性 ř Michal K ížek 张智民 谨以此文祝贺Ivo Babu ška 教授90 岁寿辰 引言 首先，我们介绍一下使用计算机进行数值计算的几个基本术语。众所周知，计算 机使用2 进制表示数字，通常分为32 位 （bit ）和64 位两种。如果是32 位，有一位表 示符号，7 位用来记录方幂，剩下24 位用于分数，这样一来，每个计算机表达的数和 原先给的数之间的精度为 2 的负25 次方，大约是8 位十进制小数，这就是常说的 “单 精度”。如果是64 位，有一位表示符号，11 位用来记录方幂，剩下52 位用于分数，这 样以来，精度为2 的负53 次方，大约是 16 位十进制小数，这就是 “双精度”。所谓浮 点运算，都是将每一个参与计算的数字表示成32 或64 位的浮点数来进行的。不难看出， 使用单精度和双精度，计算过程产生的舍入误差是不一样的。 其次，读者还需了解目前广泛使用的两个商业数学和工程计算软件 ：MATLAB 和 Maple 。MATLAB 是矩阵实验室 （Matrix Laboratory ）的简称，由美国MathWorks 公 司出品，它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿 真等诸多强大功能集成在一个易于使用的视窗环境中，在工业设计、科学计算以及课 堂教学等方面都有非常广泛的应用。相比于MATLAB ，Maple 的特点是符号运算，其 功能覆盖诸多的数学领域，如微积分、线性代数、微分方程、积分变换、概率论和数 理统计、图论、张量分析、微分几何、线性规划、组合数学、抽象代数、泛函分析、 数论、复分析和实分析、特殊函数、编码和密码理论、优化等。与历史更久的Fortran 和C 等计算机语言相比，MATLAB 和Maple 是更高级、更加人性化的语言，其语法相 对简单，经常是直接输入数学公式即可达到计算目的。 有了这些准备，现在言归正传。对下面这个定积分 1 1 n x I = x e dx