《西安交通大学_计算方法上机报告》由会员分享，可在线阅读，更多相关《西安交通大学_计算方法上机报告（34页珍藏版）》请在金锄头文库上搜索。
1、计算方法上机报告计算方法上机报告姓 名： 学 号：班 级：上课班级： 计算方法上机报告说明： 本次上机实验使用的编程语言是Matlab语言，编译环境为MATLAB 7.11.0，运行平台为Windows 7。1. 对以下和式计算：，要求： 若只需保留11个有效数字，该如何进行计算； 若要保留30个有效数字，则又将如何进行计算；（1） 算法思想1、根据精度要求估计所加的项数，可以使用后验误差估计，通项为： ；2、为了保证计算结果的准确性，写程序时，从后向前计算；3、使用Matlab时，可以使用以下函数控制位数： digits(位数)或vpa(变量，精度为数)（2）算法结构1.；2.for if
2、end;3.for （3）Matlab源程序clear; %清除工作空间变量clc; %清除命令窗口命令m=input(请输入有效数字的位数m=); %输入有效数字的位数s=0; for n=0:50 t=(1/16n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6); if t=10(-m) %判断通项与精度的关系 break; endend;fprintf(需要将n值加到n=%dn,n-1); %需要将n值加到的数值for i=n-1:-1:0 t=(1/16i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6); s=s+
3、t; %求和运算ends=vpa(s,m) %控制s的精度 （4）结果与分析 当保留11位有效数字时，需要将n值加到n=7， s =3.1415926536； 当保留30位有效数字时，需要将n值加到n=22, s =3.14159265358979323846264338328。 通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。2. 某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单
4、位:米)如下表所示：分点0123456深度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）算法思想 如果使用多项式差值，则由于龙格现象，误差较大，因此，用相对较少的插值数据点作插值，可以避免大的误差，但是如果又希望将所得数据点都用上，且所用数据点越多越好，可以采用分段插值方式，即用分段多项式代替单个多
5、项式作插值。分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑，紧贴地面，并且忽略水流对光缆的冲击。海底光缆线的长度预测模型如下所示，光缆从A点铺至B点，在某点处的深度为。海底光缆线的长度预测模型计算光缆长度时，用如下公式：（2）算法结构1.For 1.1 2.For 2.1 For 2.1.1 3.4.For 4.1 4.2 4.3 5.6.7.获取M的矩阵元素个数，存入m8.For 8.1 8.2 8.3 9.10.For 10.1 11.获取x
6、的元素个数存入s12.13.For 13.1 if then ；break else 14.（3）Matlab源程序clear;clc; x=0:1:20; %产生从0到20含21个等分点的数组X=0: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=length(X);% 求三次样条插值函数s(x) M=y; for k=2:3;
7、 %计算二阶差商并存放在M中 for i=n:-1:k; M(i)=(M(i)-M(i-1)/(x(i)-x(i-k+1); endendh(1)=x(2)-x(1); %计算三对角阵系数a,b,c及右端向量dfor i=2:n-1; h(i)=x(i+1)-x(i); c(i)=h(i)/(h(i)+h(i-1); a(i)=1-c(i); b(i)=2; d(i)=6*M(i+1);end M(1)=0; %选择自然边界条件M(n)=0;b(1)=2;b(n)=2;c(1)=0;a(n)=0; d(1)=0; d(n)=0; u(1)=b(1); %对三对角阵进行LU分解y1(1)=d(1
8、);for k=2:n; l(k)=a(k)/u(k-1); u(k)=b(k)-l(k)*c(k-1); y1(k)=d(k)-l(k)*y1(k-1);endM(n)=y1(n)/u(n); %追赶法求解样条参数M(i)for k=n-1:-1:1; M(k)=(y1(k)-c(k)*M(k+1)/u(k);ends=zeros(1,N);for m=1:N; k=1; for i=2:n-1 if X(m)=x(i); k=i-1; break; else k=i; endend H=x(k+1)-x(k); %在各区间用三次样条插值函数计算X点处的值 x1=x(k+1)-X(m); x
9、2=X(m)-x(k); s(m)=(M(k)*(x13)/6+M(k+1)*(x23)/6+(y(k)-(M(k)*(H2)/6)*x1+(y(k+1)-(M(k+1)*(H2)/6)*x2)/H;end% 计算所需光缆长度L=0; %计算所需光缆长度for i=2:N L=L+sqrt(X(i)-X(i-1)2+(s(i)-s(i-1)2);enddisp(所需光缆长度为 L=);disp(L);figureplot(x,y,*,X,s,-) %绘制铺设河底光缆的曲线图xlabel(位置,fontsize,16); %标注坐标轴含义ylabel(深度/m,fontsize,16);titl
10、e(铺设河底光缆的曲线图,fontsize,16);grid;（4）结果与分析 铺设海底光缆的曲线图如下图所示： 仿真结果表明，运用分段三次样条插值所得的拟合曲线能较准确地反映铺设光缆的走势图，计算出所需光缆的长度为 L=26.4844m。3. 假定某天的气温变化记录如下表所示，试用数据拟合的方法找出这一天的气温变化的规律；试计算这一天的平均气温，并试估计误差。时刻0123456789101112平均气温15141414141516182020232528时刻131415161718192021222324平均气温313431292725242220181716（1）算法思想在本题中，数据点的
11、数目较多。当数据点的数目很多时，用“多项式插值”方法做数据近似要用较高次的多项式，这不仅给计算带来困难，更主要的缺点是误差很大。用“插值样条函数”做数据近似，虽然有很好的数值性质，且计算量也不大，但存放参数的量很大，且没有一个统一的数学公式来表示，也带来了一些不便。另一方面，在有的实际问题中，用插值方法并不合适。当数据点的数目很大时，要求通过所有数据点，可能会失去原数据所表示的规律。如果数据点是由测量而来的，必然带有误差，插值法要求准确通过这些不准确的数据点是不合适的。在这种情况下，不用插值标准而用其他近似标准更加合理。通常情况下，是选取使最小，这就是最小二乘近似问题。在本题中，采用“最小二乘法”找出这一天的气温变化的规律，使用二次函数、三次函数、四次函数以及指数型函数，计算相应的系数，估算误差，并作图
