西安交大计算方法B2017大作业.docxVIP

实用标准 计算方法 B 上机报告 姓 名 ： 学 号 ： 班 级 ： 学 院 ： 任课教师： 2017 年 12 月 29 日 文案大全 实用标准 实用标准 文案大全 文案大全 题目一： 题目内容 某通信公司在一次施工中，需要在水面宽度为 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 深度 11.18 12.26 13.28 13.32 12.61 11.29 10.22 分点 14 15 16 17 18 19 20 深度 9.15 7.90 7.95 8.86 9.81 10.80 10.93 请用合适的曲线拟合所测数据点； 估算所需光缆长度的近似值，并作出铺设河底光缆的曲线图； 实现题目的思想及算法依据 首先在题目（1）中要实现的是数据的拟合，显然用到的是我们在第三章中数据近似的知识内容。多项式插值时，这里有 21 个数据点，则是一个 20 次的多项式，但是多项式插值随着数据点的增多，会导致误差也会随之增大，插值结果会出现龙格现象，所以不适用于该题目中点数较多的情况。为了避免结果出现大的误差，同时又希望尽可能多地使用所提供的数据点，提高数据点的有效使用率， 这里选择分段插值方法进行数据拟合。分段插值又可分为分段线性插值、分段二次插值和三次样条插值。由于题目中所求光缆的现实意义，而前两者在节点处的光滑性较差，因此在这里选择使用三次样条插值。 根据课本 SPLINEM 算法和 TSS 算法，采用第三种真正的自然边界条件，在选定边界条件和选定插值点等距分布后，可以先将数据点的二阶差商求出并赋值给右端向量d，再根据 TSS 解法求解三对角线线性方程组从而解得M 值。求出 M 后，对区间进行加密，计算 200 个点以便于绘图以及光缆长度计算。 对于问题（2），使用以下的公式: L = ? 20 0 f ( x )ds ? ? 20 0 f (x) 1? f (x)2 dx ? ?19 k ?0 ?k ?1 k f (x) 1? f (x)2 dx 算法结构 1. For i ? 0 ,1, 2 ,? ? ?, n 1.1 y ? M i i For k ? 1,2 For i ? n, n ? 1,?, k 2.1.1 (M ? M ) /( x ? x ) ? M i i ?1 i i ? k i x ? x ? h 1 0 1 4. For i ? 1,2,? , n -1 x ? x i?1 i ?h i?1 h /(h ? h ) ? c ;1? c ? a ;2 ? b i?1 i 6M ? d i?1 i i?1 i i i d 0 ? M ;d 0 n ? M ;? ?c n 0 0 2 ?b ;? 0 n ? a ;2 ? b n n b ?? ,d ?? 1 1 1 1 For k ? 2,3,?, m ! 获取 M 的矩阵元素个数，存入 m a k b / ? ?l k?1 k - l ? c ? ? k k k?1 k 7.3 d - l ?? ?? k k k?1 k 8. ? / ? ? M m m m 9. For k ? m ?1, m ? 2,?,1 9.1 (? ? c ? M ) / ? ? M k k k ?1 k k 10. 1? k ! 获取 x 的元素个数存入 s 11. For i ?1,2,?, s ?1 11.1 if ~x ? x i then i ? k ；break else i ?1? k 12. x k ? x k?1 ? h; x ? ~x ? x; ~x ? x k k?1 ? x? [M k ?1 x 3 ? M 6 k x?3 6 ( y  k ?1 M k ?1 h2 )x ? ( y 6 k ? M h2 )x?] / h ? ~y k 6 matlab 源程序 n=20; x=0:n; y=[9.01 8.96 7.96 7.97 8.02 9.05 10.13 11.18 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]; M=y; %用于存放差商，此时为零阶差商h=zeros(1,n+1); c=zeros(1,n+1); d=zeros(1,n+1); a=zeros(1,n+1); b=2*ones