西安交通大学计算方法B上机报告.doc

《西安交通大学计算方法B上机报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、计算方法上机报告计算方法上机报告姓名：学号：班级：能动上课班级：计算方法（B）上机报告题目及求解：一、对以下和式计算：，要求：①若只需保留11个有效数字，该如何进行计算；②若要保留30个有效数字，则又将如何进行计算；1算法思想（1）根据精度要求估计所加的项数，可以使用后验误差估计，通项为：；（2）为了保证计算结果的准确性，写程序时，从后向前计算；（3）使用Matlab时，可以使用以下函数控制位数：digits(位数)或vpa(变量，精度为数)2算法结构；forifend;for3Matlab源程序clear;%清除工作空间变量cl

2、c;%清除命令窗口命令m=input('请输入有效数字的位数m=');%输入有效数字的位数s=0;forn=0:5032计算方法（B）上机报告t=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6));ift<=10^(-m)%判断通项与精度的关系break;endend;fprintf('需要将n值加到n=%d',n-1);%需要将n值加到的数值fori=n-1:-1:0t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6));s=

3、s+t;%求和运算ends=vpa(s,m)%控制s的精度4结果与分析若保留11位有效数字，则n=7，此时求解得：s=3.1415926536；若保留30位有效数字时，则n=22,此时求解得：s=3.14159265358979323846264338328。通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。二、某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依

4、据。已探测到一组等分点位置的深度数据(单位:米)如下表所示：分点012345632计算方法（B）上机报告深度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算法思想如果使用多项式差值，则由于龙格现象，误差较大，因此，用相对较少的插值数据点作插

5、值，可以避免大的误差，但是如果又希望将所得数据点都用上，且所用数据点越多越好，可以采用分段插值方式，即用分段多项式代替单个多项式作插值。分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑，紧贴地面，并且忽略水流对光缆的冲击。海底光缆线的长度预测模型如图2-1所示，光缆从A点铺至B点，在某点处的深度为。图2-1海底光缆线的长度预测模型计算光缆长度时，用如下公式：32计算方法（B）上机报告2

6、算法结构1）For1.12）For2.1For2.1.13）4）For4.14.24.35）6）7）获取M的矩阵元素个数，存入m8）For8.18.28.39）10）For10.111）获取x的元素个数存入s32计算方法（B）上机报告12）13）For13.1ifthen；breakelse14）3Matlab源程序clear;%清除工作空间变量clc;%清除命令窗口命令x=0:1:20;%产生从0到20含21个等分点的数组X=0:0.2:20;y=[9.01,8.96,7.96,7.97,8.02,9.05,10.13,11.1

7、8,12.26,13.28,13.32,12.61,11.29,10.22,9.15,7.90,7.95,8.86,9.81,10.80,10.93];%等分点位置的深度数据n=length(x);%等分点的数目N=length(X);%%求三次样条插值函数s(x)M=y;fork=2:3;%计算二阶差商并存放在M中fori=n:-1:k;M(i)=(M(i)-M(i-1))/(x(i)-x(i-k+1));endendh(1)=x(2)-x(1);%计算三对角阵系数a,b,c及右端向量d32计算方法（B）上机报告fori=2:n

8、-1;h(i)=x(i+1)-x(i);c(i)=h(i)/(h(i)+h(i-1));a(i)=1-c(i);b(i)=2;d(i)=6*M(i+1);endM(1)=0;%选择自然边界条件M(n)=0;b(1)=2;b(n)=2;c(1)=0;