计算方法上机作业.docx

计算方法上机报告姓名：学号：班级：上课班级：说明：本次上机实验使用的编程语言是Matlab语言，编译环境为MATLAB1.11.0,运行平台为WindOWS7。1.对以下和式计算：s=y(-i?1Zj16n8n÷18n+48n+58n+6)TTTr八，要求：若只需保留11个有效数字，该如何进行计算；若要保留30个有效数字，则又将如何进行计算；(1)算法思想1、根据精度要求估计所加的项数，可以使用后验误差估计，通项为:(211148n+48n+5-8n+6Jcc16n8n+l<f2、为了保证计算结果的准确性，写程序时，从后向前计算；3、使用MatIab时，可以使用以下函数控制位数：digits（位数）或vpa（变量，精度为数）（2）算法结构t=_L（_l21i_1s-0;16n8n+18n+48n+58n+6/12. for11-0,1,2,''Jif10end;3. for11=j-lj-2,0s-5÷t;O)MatIab源程序clear;%清除工作空间变量clc;%清除命令窗口命令m=input(谛输入有效数字的位数m三,);%输入有E数字的位数s=0;forn=0:50t=(l16n)*(4(8*n+l)-2(8*n+4)-l(8*n+5)-l(8*n+6);ift<=10-m)%判断通项与精度的关系break;endend;fprintf。需要将n值加到n=%dn,zn-l);%需要将n值加到的数值t=l16i)*(4(8*i+l)-2(8*i+4)-l(8*i+5)-l(8*i+6);s=s+t;%求和运算S=VPa（S,m）%控制S的精度（4）结果与分析当保留11位有效数字时，需要将n值加到n=7,s=3.1415926536；当保留30位有效数字时，需要将n值加到n=22,s=3.14159265358979323846264338328o通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。2.某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据（单位:米）如下表所示：分点O123456深度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)算法思想如果使用多项式差值，则由于龙格现象，误差较大，因此，用相对较少的插值数据点作插值，可以避免大的误差，但是如果又希望将所得数据点都用上，且所用数据点越多越好，可以采用分段插值方式，即用分段多项式代替单个多项式作插值。分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑,紧贴地面，并且忽略水流对光缆的冲击。海底光缆线的长度预测模型如下所示，光缆从A点铺至B点，在某点处的深度为海底光缆线的长度预测模型计算光缆长度时，用如下公式:尸20f(×)ds0=f /(x)Vl(x)2dxJO=Er%w+f() k=0 k2 dx =£ V(x)2+(y)2(2)算法结构1.Fori=0,12,"Lly尸也2.Fork=l,22.1For=n,n-l2.1.1(MLM-I)/(阳一&_。=弧3.XLXO=%4.For=l,2,n-l4.1勺+1一*产力h142力hi(A4L+i)=c1；lc产%；2=%4.36MH4=df5.do="。也=MfI认OnCo2=如甲产%;2=36.如=4/1=7.获取M的矩阵元素个数，存入m8.Fork2,3,，川8.1°1-="82%-4g="A83九山少女Tnyk9,%/"皿=”610.ForAm-l,印一2,/10.1"kfM*i)"A=%L获取X的元素个数存入S12.10413.Fori=1,2,，5113.1if*"，theni=k；break.，elsei+l=k14,XkTkr=h：XATnX；XTl=X33-AJ22YV>4/)4*M1不+如石+(JI-MI石)x+(A-Mlrr)x/力=yMatlab源程序clear;clc;x=0:l:20;%产生从0到20含21个等分点的数组X=O:0.2: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;%等分用位置的深度数据n=length(x);%等分点的数目N=Iength(X);%求三次样条插值函数S(X)M=y;fork=2:3;%计算二阶差商并存放在M中fori=n:-l:k;M(i)=(M(i)-M(i-l)(×(i)-x(i-k+l);endendh(l)=×(2)-×(l);%计算三对角阵系与a,b,c及右端向量dfori=2:n-l;h(i)=x(i+l)-x(i);c(i)=h(i)(h(i)+h(i-l);a()=l-c(i);b(i)=2;d(i)=6*M(i+l);endM(I)=O;%选择自然边界条件M(n)=0;b(l)=2;b(n)=2;c(D=0;a(n)=0;d(D=0;d(n)=0;u(l)=b(l);%对三对角阵进行LU分解yl(l)=d(l);fork=2:n;l(k)=a(k)u(k-l);u(k)=b(k)-l(k)*c(k-l);yl(k)=d(k)-l(k)*yl(k-l);endM(n)=yl(n)u(n);%追赶法求解样条参数M(i)fork=n-l:-l:l;M(k)=(yl(k)-c(k)*M(k+l)u(k);ends=zeros(l,N);form=l:N;k=l;fori=2:n-lifX(m)<=×(i);k=i-l;break;elsek=i;endendH=x(k+l)x(k);%在各区间用三次样条插值函数计算X点处的值xl=x(k+l)-×(m);x2=X(m)-x(k);s(m)=(M(k)*(xlA3)6+M伙+l)*(x27)6+(y(k)-(M(k)*(HA2)6)*xl+(y伙+1HM化+1)*(HA2)6)*x2)/H；end%计算所需光缆长度1.=0;%计算所需光缆长度fori=2:N1.=L+sqrt(X(i)-X(i-l)2+(s(i)-s(i-l)2);enddisp。所需光缆长度为L=');disp(L);figurePlOt(XMtX,s,U)%绘制铺设河底光缆的曲线图Xlabel('位置IfontSiZHl6);%标注坐标轴含义ylabel('m'fontsize,z16);tit底，铺设河底光缆的曲线图+OmSiZHI6)；grid;(4)结果与分析铺设海底光缆的曲线图如下图所示:68101214161820位置铺设河底光缆的曲线图Il10仿真结果表明，运用分段三次样条插值所得的拟合曲线能较准确地反映铺设光缆的走势图，计算出所需光缆的长度为L=26.4844mo3.假定某天的气温变化记录如下表所示，试用数据拟合的方法找出这一天的气温变化的规律；试计算这一天的平均气温，并试估计误差。（1）算法思想在本题中，数据点的数目较多。当数据点的数目很多时，用“多项式插值”方法做数据近似要用较高次的多项式，这不仅给计算带来困难，更主要的缺点是误差很大。用“插值样条函数”做数据近似，虽然有很好的数值性质，且计算量也不大，但存放参数MI的量很大，且没有一个统一的数学公式来表示，也带来了一些不便。另一方面，在有的实际问题中，用插值方法并不合适。当数据点的数目很大时，要求，（X）通过所有数据点，可能会失去原数据所表示的规律。如果数据点是由测量而来的，必然带有误差，插值法要求准确通过这些不准确的数据点是不合适的。在这种情况下，不用插值标准而用其他近似标准更加合理。通常情况下，是选取明使"二最小，这就是最小二乘近似问题。在本题中，采用“最小二乘法”找出这一天的气温变化的规律，使用二次函数、三次函数、四次函数以及指数型函数C二。e-Mb,计算相应的系数，估算误差，并作图比较各种函数之间的区别。（2）算法结构本算法用正交化方法求数据的最小二乘近似。假定数据以用来生成了G,并将y作为其最后一列（第列）存放。结果在。中，P是误差1.使用二次函数、三次函数、四次函数拟合时1 .将“时刻值”存入勺，数据点的个数存入输入拟合多项式函数P（X）的最高项次数小二一1,则拟合多项式函数为P（X）="（x）+孙（XHTa叫（X）,根据给定数据点确定GFor/=0,1,2,Ji-lFori=12，小2.1”尸。叱12.2，产。“+13FOrk=1，2,，3.1形成矩阵Q*3.1.1-SgMgn)(Eg 严=O3.1. 2N kk 3.1.3 For=fc+l,+2,113.1.3.10k3,2变换G'1到Gfc3. 2.1°ngkkFor/-k+l+2,nr11+13.2.24. 2.3For/-+1,m3.2.3.1处1“尸氏/4.解三角方程Rfl=h1j4.19升1/9。0=%4.2Fori二一L一2.，14.2.1nxg*/g“=s2EI5.计算误差%m/="+1"、使用指数函数拟合时现将指数函数进行变形：将C=y'=X代入C=Qei("靖得：y=Qe-W一靖对上式左右取对数得:Iny=Ina-bc22bcx-bx2令z=lny,0=11a-bc2,al=2bcfa2=-b则可得多项式：Z=04111X+2X2(3)MatIab源程序clear;%清除工作空间变量clc;%清除命令窗Ll命令x=0:24;%将时刻值存入数组y=15,14,14,14,14,15,16,18,20,20,23,25,28,31,34,31,29,27,25,24,22,20,18,0,16;,m=size(x);%将数据点的个数存入mT=sum(y(l:m)/m;fprint