计算方法290712.pptVIP

第一章总结 一、误差来源 模型误差/描述误差 观测误差/测量误差 截断误差：把无穷数用有限数表示。 舍入误差：受字长限制，取有限位数字运算。 三、 应选用数值稳定性好的计算公式 例： 计算积分 解： 采用分部积分法 得到In的递推公式 假设用4位小数计算9个积分值（I0I1…I8) 0.1008 I8 0.1124 I7 0.1268 I6 0.1455 I5 0.1709 I4 0.2073 I3 0.2642 I2 0.3699 I1 0.6321 I0 真值 In 用 递推公式推出的值I1, I2,…,I8 0.1008 -0.7280 I8 0.1124 0.2160 I7 0.1268 0.1120 I6 0.1455 0.1408 I5 0.1709 0.1704 I4 0.2073 0.2074 I3 0.2642 0.2642 I2 0.3699 0.3679 I1 0.6321 0.6321 I0 真值 In 改进递推公式： 当n=8时， 0.0409I80.1111 取初值 递推公式n=8,7,…,1倒推计算，误差缩小n!倍。 用 递推公式推出的值I1, I2,…,I8 0.1008 0.0760 -0.7280 I8 0.1124 0.1155 0.2160 I7 0.1268 0.1264 0.1120 I6 0.1455 0.1456 0.1408 I5 0.1709 0.1709 0.1704 I4 0.2073 0.2073 0.2074 I3 0.2642 0.2642 0.2642 I2 0.3699 0.3699 0.3679 I1 0.6321 0.6321 0.6321 I0 真值 In 1.1 科学与工程计算含义 1.2 误差的产生和分类 1.3 误差和有效数字 1.3.1 绝对误差和绝对误差限 1.3.2 有效数字 1.3.3 相对误差和相对误差限 1.3.4 有效数字与相对误差 1.4 运算误差分析 1.4.1 函数运算误差 1.4.2 算数运算误差 1.5 数值稳定性和减小运算误差 1.5.1 数值稳定性 1.5.2 减小运算误差 第一章 绪论 1.5.2 减小运算误差 一、 两个相近的数相减，会严重损失有效数字 例： 1、 设y=x-A，其中A在运算时不发生误差 2、 解决相近的数相减的问题 ⑴ 多保留有效数字 ⑵ 变换计算公式 常见的变换公式： 当x1和x2接近时， 当x→0时， 当x充分大时， 当f(x)和f(x*)很接近，又需作f (x)-f (x*)运 算时，为避免有效数字丢失，可采用泰勒 展开式： 取右端有限项近似左端。 3、 举例 ⑴ 解： ⑵ 解： ⑶ 解： ⑷ 解： 二、 要防止大数“吃掉”小数 例1： 1、计算机字长有限，在进行加减法运算时，要 对阶和规格化，对阶是以大数为基准的。 例2： 2、当参加运算的两个数数量级相差很大时，要 注意运算次序，可能避免大数吃掉小数。 例3： 例4： 在5位10进制计算机上计算 解： 应先算 例5： 在5位10进制计算机上求解 可用大数吃掉小数 不能用大数吃掉小数 例6： 用消元法解方程组 若只用3位数计算，结果是否可靠？ 若只用3位数计算，1和2会被1015吃掉，结果x2=1不准确。 解： 三、 绝对值太小的数不宜做除数 四、 简化计算步骤，减少运算次数 例1： 解1： 计算254次 解2： 计算14次 例2： 计算多项式 法1： 直接计算akxk(k=1,2,…,n) 再逐项相加，需计算 次乘法和n次加法。 * 一、 绝对误差 二、 绝对误差限 m—整数 若四舍五入保留n位得到绝对误差限 x1,x2,…xn是x*的有效数字。 1、如果近似数x*的误差限是某一位的半个单位， 由该位到x*的第一位非零数有n位，就有n位有 效数字。 例： 30.120 ε*= n=5 2、若x* =± 0.x1x2…xn×10m是四舍五入得到的， 则n为有效数字。 例： x1=3.14159 x1*=3.141 n≠4 |e1* | = |x1-x1*| = 0.00059 0.0005 x2=3.14159 x2*=3.142 n=4 3、若x* =± 0.x1x2…xn×10m不是四舍五入得到的， 则n不是有效数字。 4、 n相同的两个近似数x1*、x2* ，其e1 *和e2*不 一定相等。 例1： x=0.9460831… x1*=0.9460832 x2*=0.946