计算方法30655.ppt

§1.1计算方法的任务与特点 基本的数学问题: 1.大型线性代数方程组Ax=b求解; 2.矩阵A的特征值和特征向量计算; 3.非线性方程 求解(求根); 4.积分 计算; 5.常微分方程初值问题求解; 6.其它。 求精确解(值)一般非常困难。例如： 1. 方程组阶数n很大，例如n=20,计算机运算速度 1亿次/秒,用不好的方法,大约需算30多万年; 好方法不到一分钟。另外，有计算结果可靠性 问题。 2. 特征值定义 3. 形式复杂时求根和求积分很困难。 4.线性微分方程易解， 如 非线性方程难解，如 §1.2 误差基础知识 4. 舍入误差。计算机字长有限，一般实数不能精确存储，于是产生舍入误差。例如：在10位十进制数限制下： 舍入误差很小，本课程将研究它在运算过程中 是否能有效控制。 二．误差基本概念 1．绝对误差。设 ——准确值， ——近似值。 称 　为 的绝对误差。 为 的绝对误差限。 2．相对误差。称 　　为 的相对误差。 实用中，常用 表示 的相对误差。 称 为 的相对误差限。 例1.3：测得某桌面的长a的近似值a*=120cm,宽b的 近似值b*=60cm。若已知|e(a*)|≤0.2cm, |e(b*)|≤0.1cm。 试求近似面积s*=a*b* 的绝对误差限与相对误差限。 §1.3 选用算法应遵循的原则 1.尽量简化计算步骤，减少乘除运算的次数. 例如，计算多项式 通常运算的乘法次数为 若采用递推算法, 则乘法次数仅为n. 又如 2.防止大数“吃掉”小数 当|a||b|时，尽量避免a+b 。例如，假设计算机 只能存放10位尾数的十进制数，则 3.尽量避免相近数相减 例如，当x很大时，应 4.避免绝对值很小的数做分母 当|b||a|时，应尽量避免 。 5. 选用数值稳定性好的算法，以控制舍入误差高速 增长 例如 若 （误差 ）则计算 时误差扩大了 倍,而 基本要求： 1.熟悉计算方法在解决实际问题中所处的地位，熟悉计算方法是以计算机为工具求近似解的数值方法； 2.熟悉绝对误差（限），相对误差（限）及有效数字概念； 3.熟悉公式（1.2）--（1.9）； 4.熟悉选用算法应遵循的原则； 作业： 作业集（A）第一章 1，2，3，4，5，6，7，8 计算方法 例2.1 用二分法求 在(1,2)内的根，要求绝对误差不超过 解： f(1)=-50 有根区间 中点 f(2)=140 -(1,2)+ f(1.25)0 (1.25,1.5) f(1.375)0 (1.25,1.375) f(1.313)0 (1.313,1.375) f(1.344)0 (1.344,1.375) f(1.360)0 (1.360,1.375) f(1.368)0 (1.360,1.368) §2.2 迭代法 2.收敛定理（定理2.2） 二．反幂法——求A的按摸最小的特征值。 设A可逆，由 求实对称矩阵A的全部特征值及特征向量 ——Jacobi旋转法 基本思想： 求一般矩阵全部特征值和特征向量的QR方法