计算方法_论.ppt

计 算 方 法 本课程的基本要求 掌握数值方法的基本原理 掌握常用的科学与工程计算的基本方法 能用所学方法在计算机上算出正确结果 本课程的任务： ● 建立各种数学问题的数值计算算法的方法和理论。 通俗地讲，就是为各种实际问题提供有效的数值近似解方法。 ● 提供在计算机上实际可行的、理论可靠的、计算复杂 性好的各种常用算法。 参考书： 1.《计算方法引论》（第二版），徐萃薇，孙绳武编著。 高等教育出版社，2002  2. 《数值分析》, 李庆扬，王能超，易大义编。 清华大 学出版社 ,2001 3. 《数值计算方法》, 关治，陈景良。清华大学出版社， 1990 绪 论 绪 论 本章内容 §1 引言 §2 误差的来源及分类 §3 误差的度量 §4 误差的传播 §5 减少运算误差的原则 本章要求掌握的内容 1 引言 计算方法又称数值分析, 它是研究各种数 学问题的数值解法及其理论的一门学科。 基本的数学问题: 1.大型线性代数方程组Ax=b求解; 2.矩阵A的特征值和特征向量计算; 3.非线性方程 求解(求根); 4.积分 计算; 5.常微分方程初值问题求解; 6.其它。 2 误差的来源及分类 在用数值方法解题过程中可能产生的误差归纳起来有如下几类： 1. 模型误差 2. 观测误差 3. 截断误差 4. 舍入误差 用数学方法解决一个具体的实际问题，首先要建立数学模型，这就要对实际问题进行抽象、简化，因而数学模型本身总含有误差，这种误差叫做模型误差 数学模型是指那些利用数学语言模拟现实而建立起来的有关量的描述 数学模型的准确解与实际问题的真解不同 在数学模型中通常包含各种各样的参变量，如温度、长度、电压等，这些参数往往是通过观测得到的，因此也带来了误差，这种误差叫观测误差 数学模型中的参数和原始数据，是由观测和试验得到的 由于测量工具的精度、观测方法或客观条件的限制,使数据含有测量误差,这类误差叫做观测误差或数据误差 根据实际情况可以得到误差上下界 数值方法中需要了解观测误差,以便选择合理的数值方法与之适应 精确公式用近似公式代替时,所产生的误差叫截断误差 例如, 函数f(x)用泰勒(Taylor)多项式 在数值计算中只能对有限位字长的数值进行运算 需要对参数、中间结果、最终结果作有限位字长的处理工作，这种处理工作称作舍入处理 用有限位数字代替精确数，这种误差叫做舍入误差，是数值计算中必须考虑的一类误差 3 误差的度量 4 算术运算误差 5 减少运算误差若干原则 （2）防止大数“吃掉”小数 例 求二次方程x2-105x+1=0的根 解：按二次方程求根公式 x1=(105+(1010-4)1/2)/2 x2=(105-(1010-4)1/2)/2 在8位浮点数计算得 x1=(105+105 )/2=105 (正确）, x2=(105-105 )/2=0 (错误） 产生错误的原因 ① 出现大数1010吃掉小数4的情况 ② 分子部分出现两个相近数相减而丧失有 效数位常称为灾难性的抵消 （3）绝对值太小的数不宜做除数 当分母为两个相近数相减时,会丧失有效数字 例1.8 计算 算法2。分成三组因子。每组只取六位小数计算 a=0.0005/0.0003=1.666667(有舍入) b=0.0143/0.0125=1.144000 c=0.0012/0.0135=0.088889 (有舍入) D=a*b*c=1. 666667* 1.144000* 0.088889 =0.169482,准确到小数后5位。 b c a （4）简化计算步骤，减少运算次数 x255=xx2x4x8x16x32x64x128 原先要做254次乘法现只需14次即可 又如计算多项式 p(x)=anxn? an-1xn-1 ? … ? a1x ? a0 的值 若直接计算akxk,再逐项相加，一共要做 n+(n-1)+…+2+1=n(n+1)/2次乘法和n次加法 如果将前n项提出x，则有 p(x)=（anxn-1? an-1xn-2 ? … ? a1 ）x ? a0 =((anxn-2?an-1xn-3?…? a2)x?a1）x ? a0 =(…(anx ? an-1)x?…?a2)x ? a1)x ?