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

XIAN JIAOTONG UNIVERSITY 计算方法上 机报告 姓 名： 学 号： 班 级：能动 上课班级： 题目及求解： c_V 1 （ 4 2 1 1 ） 一、对以下和式计算：初TFkrrr亦万一需T?一需轲丿，要求: ① 若只需保留11个有效数字，该如何进行计算； ②若要保留30个有效数字，则又将如何进行计算； 1算法思想 （1）根据精度要求估计所加的项数，可以使用后骑误差估计，通项为: 16 3 + 12 18/? + 4 8n 16 3 + 1 2 1 8/? + 4 8n + 5 M 8/? + 6 1 4 16 8/? + 1 ve ； （2） 为了保证计算结果的准确性，写程序时，从后向前计算; （3） 使用Matlab吋，可以使用以下函数控制位数： digits（位数）或vpa（变量，精度为数） 2算法结构 5 = 0； 2 1 1 、 8^ + 4 8/2 + 5 8〃+ 6丿 for力=0丄2, if t 10-w end; for n = i」一hi — 2j?、0 3 Matlab源程序 clear;%清除工作空间变量 clear; %清除命令 口命令 m二in put（请输入有效数字的位数m = %输入有效数字的位数 s=0; for n=0:50 仁(1/165广(4/(85+1)?2/(8切+4)?1/(8勺+5)?1/(8勺+6)); 讦t = 10A(-m) %判断通项与精度的关系 break; end end; fprintf(儒要将n值加到n=%d\n：n?l); %需要将n值加到的数值 for i=n-l:-l:0 t=(l/16Ai)*(4/(8*i+l)-2/(8*i+4)-l/(8*i+5)-l/(8*i+6)); s=s+t; %求和运算 end S=vpa(s,m) %控制s的精度 4结果与分析 若保留11位有效数字，则27,此时求解得： s =3.1415926536； 若保留30位有效数字时，则n=22,此时求解得： s =3.14159265358979323846264338328。 通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结 果要求保留的准确数字位数的要求。 二、某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一 条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的 长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所 示： 分点 0 1 2 3 4 5 6 深度 9.01 96 7. 96 7. 97 8. 02 9. 05 10. 13 分点 7 8 9 10 11 12 13 深度 11. 18 12. 26 13. 28 13. 32 12.61 11.29 10. 22 分点 1 ■ 1 15 16 17 18 19 20 深度 9. 15 7. 90 7. 95 8.86 9.81 10. 80 10. 93 请用合适的曲线拟合所测数据点； 预测所需光缆长度的近似值，作出铺设河底光缆的曲线图； 1算法思想 如果使用多项式差值，则由于龙格现象，误差较大，因此，用相对较少的插值数据 点作插值，可以避免大的误差，但是如果又希望将所得数据点都用上，且所用数据点 越多越好，可以采用分段插值方式，即用分段多项式代替单个多项式作插值。分段多项 式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条 插值方法是一种具有较好“光滑性”的分段插值方法。 在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑，紧贴地 面，并且忽略水流对光缆的冲击。海底光缆线的长度预测模型如图2?1所示，光缆从A 点铺至B点，在某点处的深度为h 。 图2?1海底光缆线的长度预测模型 计算光缆长度时，用如下公式: r 20 厶=Jo = ￡°/(x)71 + /(x)2rfr 2算法结构 For i = 0,1,2,?…，n y戶 For k = 1,2 For i = n,n-\,???，k (Mz -x^j,) = Mj 卯一州二说 For i = 1,2,???，〃?1 xi+i —x{ 肉+i 1^ + hM) = ci ;1 一 ci = q. ;2=Z? 6Mj+] =4 % =Mo；〃“ 入 de。 2 =6?zz ;2=Z?n A =□，%=% 获取M的矩阵元素个数，存入m For k = 23,…、m 1 ak! =4 2 bk - lk - ck_x = dk - lk - 7k—\ =* 7k 9) For ￡ =加一1,加一2,???,1 （兀—％叽）/儿=＞叫 11) 12) 获取x的元素个数存入 Ink 13) For「