资源描述：
《西安交通大学计算方法b上机报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、XI'ANJIAOTONGUNIVERSITY计算方法上机报告姓名：学号：班级：能动上课班级：题目及求解：c_V1（4211）一、对以下和式计算："初TFkrrr亦万一需T?一需轲丿，要求:①若只需保留11个有效数字，该如何进行计算；②若要保留30个有效数字，则又将如何进行计算；1算法思想（1）根据精度要求估计所加的项数，可以使用后骑误差估计，通项为:16"3+1218/?+48n+5M<8/?+61416"8/?+1ve；（2）为了保证计算结果的准确性，写程序时，从后向前计算;（3）使用Matlab吋，可以使用以下函数控制位数：digits（位数）或vpa（变量，精度为
2、数）2算法结构5=0；211、8^+48/2+58〃+6丿for力=0丄2,ift<10-wend;forn=i」一hi—2j・、03Matlab源程序clear;%清除工作空间变量%清除命令口命令m二input（请输入有效数字的位数m=%输入有效数字的位数s=0;forn=0:50仁(1/165广(4/(85+1)・2/(8切+4)・1/(8勺+5)・1/(8勺+6));讦t<=10A(-m)%判断通项与精度的关系break;endend;fprintf(儒要将n值加到n=%d：n・l);%需要将n值加到的数值fori=n-l:-l:0t=(l/16Ai)*(4/(
3、8*i+l)-2/(8*i+4)-l/(8*i+5)-l/(8*i+6));s=s+t;%求和运算endS=vpa(s,m)%控制s的精度4结果与分析若保留11位有效数字，则27,此时求解得：s=3.1415926536；若保留30位有效数字时，则n=22,此时求解得：s=3.14159265358979323846264338328。通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。二、某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而
4、估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示：分点0123456深度9.01&967.967.978.029.0510.13分点78910111213深度11.1812.2613.2813.3212.6111.2910.22分点1■1151617181920深度9.157.907.958.869.8110.8010.93请用合适的曲线拟合所测数据点；预测所需光缆长度的近似值，作出铺设河底光缆的曲线图；1算法思想如果使用多项式差值，则由于龙格现象，误差较大，因此，用相对较少的插值数据点作插值，可以避免大的误差，但是如果又希望
5、将所得数据点都用上，且所用数据点越多越好，可以采用分段插值方式，即用分段多项式代替单个多项式作插值。分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑，紧贴地面，并且忽略水流对光缆的冲击。海底光缆线的长度预测模型如图2・1所示，光缆从A点铺至B点，在某点处的深度为h。图2・1海底光缆线的长度预测模型计算光缆长度时，用如下公式:r20厶=Jo=£°/(x)71+/'(x)2rfr2算法结构1)Fori=0,1,2,•…，n1.1y
6、戶2)Fork=1,22.1Fori=n,n-,・・・，k2.1.1(Mz-x^j,)=>Mj3)卯一州二>说4)Fori=1,2,・・・，〃・14.1xi+i—x{4.2肉+i1^+hM)=>ci;1一ci=>q.;2=>Z?4.36Mj+]=>45)%=>Mo；〃“入de。2=>6?zz;2=>Z?n6)A=>□，%=>%7)获取M的矩阵元素个数，存入m8)Fork=23,…、m&1ak!=>4&2bk-lk-ck_x=>dk-lk-7k—=*7k9)9)For£=加一1,加一2,・・・,19.1（兀—％叽）/儿=＞叫11)12)获取x的元素个数存入Ink13)F
7、or「=12・・沙一113.1ifxA14)xk-XkA二>/l；Xk-X^>X；丘一母_1x3x3h2h2++(几_i一哝_]三-)元+(几一哝三■戊]//?=>孑ooOO3Matlab源程序clear;%清除工作空间变量%清除命令窗口命令x=0:l:20;%产生从0到20含21个等分点的数组X=0:0.2:20;y=[9.01/8.96/7.96/7.97/8.02/9.05/10.13/ll.l8,1226,13.28,1332,12.61,11.29,10.22,9.1
