计算方法B上机报告材料.docxVIP

实用标准文案 实用标准文案 精彩文档 精彩文档 计算方法 B 上机报告 第 1 题 03某通信公司在一次施工中，需要在水面宽度为 20 米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长 0 3 分点 0 1 2 3 4 5 6 深度 9.01 8.96 7.96 7.97 8.02 10.1 9.05 3 分点 7 8 9 10 11 12 13 11.1 12.2 13.2 13.3 12.6 11.2 10.2 深度 8 6 8 2 1 9 2 分点 14 15 16 17 18 19 20 10.8 10.9 深度 9.15 7.90 7.95 8.86 9.81 请用合适的曲线拟合所测数据点； 估算所需光缆长度的近似值，并作出铺设河底光缆的曲线图； 问题分析和算法思想： 本题的主要目的是对 21 个测量数据进行拟合，同时对拟合曲线进行线积分即可得到河底光缆长度的近似值，可以用的插值方法很多：多项式插值、Lagrange 插值、Newton 插值、三次样条插值等。由于数值点较多时，采用高次多项式插值将产生很大的误差， 用拉格朗日插值多项式会出现龙格现象。故为了将所有的数据点都用上，且题中光缆为柔性，可光滑铺设于水底，鉴于此特性，采用三次样条插值的方法较为合适。 计算光缆长度近似值，只需将每两点之间的距离算出，然后依次相加，所得的折线长度，即为光缆长度的近似值。 光缆长度计算公式： l ? ?20 1? ?f ??x??2 dx ? ?19 ?k ?1 ?20 1? ?f ??x ??2 dx 0 算法结构： k ?0 k 0 三次样条算法结构见《计算方法教程》P110。源程序： clear;clc; x=0:20; 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]; d=y; plot(x,y,k.,markersize,15) hold on %%%计算二阶差商for k=1:2 for i=21:-1:(k+1) d(i)=(d(i)-d(i-1))/(x(i)-x(i-k)); end end %%%假定d的边界条件，采用自然三次样条for i=2:20 d(i)=6*d(i+1); end d(1)=0; d(21)=0; %%%追赶法求解带状矩阵的m值a=0.5*ones(1,21); b=2*ones(1,21); c=0.5*ones(1,21); a(1)=0;c(21)=0; u=ones(1,21); u(1)=b(1); r=c; yy(1)=d(1); %%%追的过程for k=2:21 l(k)=a(k)/u(k-1); u(k)=b(k)-l(k)*r(k-1); yy(k)=d(k)-l(k)*yy(k-1); end %%%赶的过程m(21)=yy(21)/u(21); for k=20:-1:1 m(k)=(yy(k)-r(k)*m(k+1))/u(k); end %%%利用插值点画出拟合曲线k=1; nn=100; xx=linspace(0,20,nn); l=0; for j=1:nn for i=2:20 if xx(j)=x(i) k=i; break; else k=i+1; end end h=1; xbar=x(k)-xx(j); xmao=xx(j)-x(k-1); s(j)=(m(k-1)*xbar^3/6+m(k)*xmao^3/6+(y(k-1)-m(k-1)*h^2/6)*xbar+(y(k)-m(k) *h^2/6)*xmao)/h; sp(j)=-m(k-1)*(x(k)-xx(j))^2/(2*h)+m(k)*(xx(j)-x(k-1))^2/(2*h)+(y(k)-y(k- 1))/h-(m(k)-m(k-1))*h/6; l(j+1)=(1+sp(j)^2)^0.5*(20/nn)+l(j);%利用第一类线积分求河底光缆的长度end %%%绘图 title(光缆插值曲线) xlabel(分点) ylabel(深度) plot(xx,s,r-,linewidth,1.5) grid disp([所需光缆长度为,num2str(l(nn+1)),米]) 运行结果： 光缆插值曲线 14 13 12 11 度深 10 9 8 7 0 2 4 6 8 10 分点  12 14 16 18 20 实用标准文案 实用标准文案 精彩文档 精彩文档 第 2 题 假定某天的气温变化记录如下表所示，试用数据拟合的方法找出这一