的相对误差.ppt

* 求 的小正根. 解 只有一位有效数字. 则具有3位有效数字. 若改用 例7 由求根公式 * 例8 计算 （用四位数学用表）. 由于 ， 只有一位有效数字. 具有三位有效数字 （这里 ）. 则 若利用 直接计算 * 此例说明，可通过改变计算公式避免或减少有效数字 的损失. 类似地，如果 和 很接近时，由 用右边算式有效数字就不损失. 也应该用右端算式代替左端. 当 很大时， * 一般情况，当 时，可用泰勒展开 取右端的有限项近似左端. 如果无法改变算式，则采用增加有效位数进行运算； 在计算机上则采用双倍字长运算，但这要增加机器计算时间和多占内存单元. * 3. 要防止大数“吃掉”小数 在数值运算中参加运算的数有时数量级相差很大， 例9 其中 . 而计算机位数有限，如不注意运算次序就可能出现大数“吃掉”小数的现象，影响计算结果的可靠性. 在五位十进制计算机上，计算 * 把运算的数写成规格化形式 由于在计算机内计算时要对阶， 若取 ， 对阶时 ，在五位的计算机中表示为 结果显然不可靠，这是由于运算中出现了大数 52492 “吃 掉”小数 造成的. 机器 0 ， 因此 * 于是 如果计算时先把数量级相同的一千个 相加，最后 再加52492，就不会出现大数“吃”小数现象， 这时 * 4. 注意简化计算步骤，减少运算次数 同样一个计算问题，如果能减少运算次数，不但可节省计算机的计算时间，还能减少舍入误差. * 例10 的值， 若直接计算 再逐项相加，一共需做 次乘法和 n 次加法. 若采用秦九韶算法 （3.4） 计算多项式 只要 n 次乘法和n 次加法就可算出 Pn (x)的值. * (3) 松弛技术. (1) 以直代曲技术. (线性化) 1.4 算法设计的常用技术 (2) 迭代技术. * 1.4.1 以直代曲技术. (线性化) 【例】已知可导函数 在 的值，求 . 解： * 1.4.2 迭代技术 【例】解非线性方程 . 解：从某个初值猜测值 出发，构造迭代序列： 逐次逼近方程的解. 例如：对a0, 求a的平方根等价于解方程 x2-a = 0, 解方程 x=(x+a/x)/2 . 从某个初值x0出发，构造迭代序列： * 【例12】用迭代法（4.4）求 . 解：取初值猜测值 ，若要求精确到10-6，由（4.4）得迭代序列： 迭代终止. 注意到 ，可知，只要迭代3次误差即小于 【实验1.2】编写程序实现算法（4.4），并计算 , 相邻两次迭代的差不超过 时终止迭代. * 1.4.3 松弛技术（relaxation） 【例】求积分 . 解：用梯形公式得： 引入松弛因子 进行加权平均得： 取 得著名的Simpson公式： * 常用数学软件： (3) Matlab (Matrix Laboratory)：最流行的科学与工程计算软件. (1) Fortran: 经典的科学计算工具. 1.5 数学软件 (4) Maple：擅长符号计算. (5) SAS，SPSS: 最流行的统计分析软件. (2) Mathematica: 擅长符号计算. (6) Lindo,Lingo: 常用的优化计算工具. * * 由(2.1)′可得 当 有 位有效数字时 反之，由 证明 * 知 至少有 位有效数字. 定理说明，有效位数越多，相对误差限越小. * 由于 知 ， 故只要取 ， 即只要对 的近似值取4位有效数字，其相对误差限就 小于0.1%. 此时由开方表得 . 要使 的近似值的相对误差限小于0.1%，需取 设取 位有效数字， 例3 几