计算方法B上机报告.pdfVIP

计算方法 B 上机实习报告 计算方法 B 上机实习报告 1 4 2 1 1 1. 对以下和式计算： S ，要求： n n 0 16 8n 1 8n 4 8n 5 8n 6 (1) 若只需保留 11 个有效数字，该如何进行计算； (2) 若要保留 30 个有效数字，则又将如何进行计算； 问题分析 ： 1 4 2 1 1 在该题中 S 的每一项 sn n ( ) 存在两个相近的数相 16 8n 1 8n 4 8n 5 8n 6 减的问题， 因此为了避免有效数字损失， 最好是改变运算顺序，分别将正数和负数相加，然 后再将其和相加。另外， sn 中有多个负数相加，可以按照绝对值递增的顺序求和，以减少 舍入误差的影响。 同时， 为了避免大数吃小数的问题， 本题先计算出保留目标有效数字所需 要的迭代次数，然后采用倒序相加的方法实现。 程序实现： clear;clc; m=input( 请输入要保留的有效数字位数： ); s1=0; s2=0; k=0; s=1; %%%% 判断多需要的迭代次数 while s=0.5*10^-(m-1) s=4/(16^k*(8*k+1))-(2/(16^k*(8*k+4))+1/(16^k*(8*k+5))+1/(16^k*(8*k+6))); k=k+1; end %%%% 正负数分别按照绝对值递增的顺序倒序相加 for n=(k-1):-1:0 a1=4/(16^n*(8*n+1)); a2=2/(16^n*(8*n+4)); a3=1/(16^n*(8*n+5)); a4=1/(16^n*(8*n+6)); s1=a1+s1; s2=a4+a3+a2+s2; end S=s1-s2; S=vpa(S,m) 运算结果： 总结心得： 在计算求和问题中， 应特别注意相近数相减的问题， 这样会造成有效数字灾难性的损失。 另外在两个数量级相差较大的数字相加减时， 较小数的有效数字会被丧失， 因此要按照从小 到大的顺序相加。 在上题计算中分别对正负相采用倒序相加， 这样就有效的避免了 “大数吃 小数”的问题。 2. 某通信公司在一次施工中，需要在水面宽度为 20 米的河沟底部沿直线走向铺 设一条沟底光缆。 在铺设光缆之前需要对沟底的地形进行初步探测， 从而估计所 需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据 ( 单 位 : 米 ) 如下表所示： 分点 0 1 2 3 4 5 6 深度 9.01 8.96 7.96 7.97 8.02 9.05 10.13 分点 7 8 9 10 11 12 13