计算方法(B)上机作业.docx

《计算方法B》上机实习报告 学院：电信学院 姓名：许梦芸 班级：硕5026 学号：3115032012 代课老师：苏剑 1.对以下和式计算：，要求：(1)若只需保留11个有效数字，该如何进行计算；(2)若要保留30个有效数字，则又将如何进行计算； 解答： A:算法实现的思想及依据  = 1 \* GB3 ①首先对问题分析可知：当一般项116n(48n+1-28n+4-18n+5-18n+6)?时，截断误差不会对保留数字位数和精度产生影响；  = 2 \* GB3 ②当n=0时，估计出10S01,这表明S要保留m个有效数字，那么小数点后位数应该取到m-1位，因此误差上限?=0.1×10-（m-1）。 由题可得S(1/16^n)*(4/(8*n+1)) 保留11个有效数字，即误差ε0.5*10^（1-11） 保留30个有效数字，即误差ε0.5*10^（1-30） B:算法实现的结构 分为两个部分：顺序结构和循环结构。由于不知道循环次数，因此采用if循环，利用通项是否大于等于误差上限来判断循环条件的真伪。 Matalab中可以利用vpa函数控制最终结果的有效数位。 C:源程序及相关的注释说明 有效数字为11的情况 digita=11; %有效数字位数 for i=1:1:1000 S0=(1/(16^(i-1)))*(4/(8*i-7)); S1=(1/16^i)*(4/(8*i+1)); if((S0-S1)0.5*10^(1-digita)) n=i; break; end end S=0; digits(digita); for i=0:1:n S=S+vpa(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); End （1）运行结果： 算法同上所述 %有效数字为30的情况 digita=30;% 有效数字位数 for i=1:1:1000 S0=(1/(16^(i-1)))*(4/(8*i-7)); S1=(1/16^i)*(4/(8*i+1)); if((S0-S1)0.5*10^(1-digita)) n=i; break; end end S=0; digits(digita); for i=0:1:n S=S+vpa(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6)); end （2）运行结果： D:结果分析及总结 为保持计算精度，计算中对误差上取得较为严格：?=0.1×10-（m-1）；由于事先不知道循环次数，采用if循环。 编程中出现过错误：对vpa函数理解不当，vpa函数是计算精度控制函数，一般可与digits位控制函数配合使用，如果没有声明，系统默认digits控制精度20位；使用单独使用vpa函数，只控制显示结果的精度，函数内部变量的运算精度不受影响，matlab浮点运算按8字节双精度进行，本题的精度完全可以得到满足。 最后，由于本例matlab运算中的精度很高，不会出现大数逐次与小数相加减而导致误差，也没有出现相近数相减而导致精度降低，因此没有进行特别的处理，直接运算，结果精度满足要求。 2.某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示： 分点0123456深度9.018.967.967.978.029.0510.13分点78910111213深度11.1812.2613.2813.3212.6111.2910.22分点14151617181920深度9.157.907.958.869.8110.8010.93 (1)请用合适的曲线拟合所测数据点； (2)预测所需光缆长度的近似值，并作出铺设河底光缆的曲线图；解答： A:算法实现的思想及依据 多项式插值，随着插值数据点的数目增多，误差也会随之增大，因而要用较少的数据点才可避免大的误差，但我们总希望将所得的数据点都用上，分段插值可以解决这一矛盾。分段插值低阶次多项式不具有“光滑性”，因此选用分段三次多项式样条插值。拟合数据点使用多项式插值时由于龙格现象误差较大，故考虑使用三次样条插值。 B:算法实现的结构 参考教材给出的SPLINEM算法和TTS算法，在选定边界条件和选定插值点等距分布后，可以先将数据点的二阶差商求