数学实验三报告总结计划.docxVIP

实验 3 插值与数值积分 实验报告 一、实验目的 1、掌握用 Matlab 计算拉格朗日 、分段线性 、三次样条三种插值的方法 ，改变节点的数目 ， 对三种插值结果进行初步分析 ； 2、掌握用 Matlab 及梯形公式 、辛普森公式计算数值积分 ； 3、通过实例学习用插值和数值积分解决实际问题 。 二、实验内容 10、表 3.7 给出的 x,y 数据位于机翼剖面的轮廓线上 ， y1 和 y2 分别对应轮廓的上下线 。假设 需要得到 x 坐标每改变 0.1 时的 y 坐标 。 试完成加工所需数据 ，画出曲线 ，求机翼剖面的面 积。 表 3.7 机翼剖面轮廓线数据 x 0 3 5 7 9 11 12 13 14 15 y1 0 1.8 2.2 2.7 3.0 3.1 2.9 2.5 2.0 1.6 y2 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 曲线绘制 利用 Matlab 编程画出机翼轮廓线 ，内容如下 ： 三次样条插值 x=[0 3 5 7 9 11 12 13 14 15]; y1=[0 1.8 2.2 2.7 3.0 3.1 2.9 2.5 2.0 1.6]; y2=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]; u=0:0.1:15; v1=spline(x,y1,u); v2=spline(x,y2,u); plot(u,v1,u,v2);grid; xlabel( x );ylabel( y1 或 y2 ); gtext( y1(x) );gtext( y2(x) ); 其中 ， u 为插值设置了步长和范围 ， grid 命令可以为曲线图添加网格线 ， xlabel 、 ylabel 分 别为横坐标 、纵坐标添加标签 ，gtext 命令可以实现曲线名称的添加 。得到的结果如下 ： 分段线性插值 将 v1 、v2 部分代码改为 v1=interp1(x,y1,u); v2=interp1(x,y2,u); 得到的图形如下 ： 可见 ，用分段线性插值会使得曲线不够光滑 （特别是区间 [10,15] 的部分 ）。 拉格朗日插值 根据教材内容 ，用 Matlab 编程 ，内容如下 ： function y=lagr(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0; for k=1:n p=1; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 并相应地改变 v1 、 v2 ，输出的结果如下 ： 可见 ，曲线出现了严重的振荡 ，对于此题 ，拉格朗日插值法已不适用 。 面积计算 利用梯形公式计算图形面积 ，相关代码如下 ： v=v1-v2; trapz(u,v) 对于三次样条插值法 ，输出的结果为 ： ans = 11.3444 对于分段线性插值法 ，输出的结果为 ： ans = 10.7500 由于三次样条插值法得到的图形更接近实际图形 ，用它算得的面积应更准确 。 11 、图 3.13 是欧洲一个国家的地图 （图参考教材 ）， 为了算出它的国土面积 ，首先对地图 作如下测量 ，以由西向东方向为 x轴，由南向北方向为 y 轴，选择方便的原点 ，并将从最西 边界点到最东边界点在 x 轴上的区间适当地划分为若干段 ，在每个分点的 y 方向测出南边界 点和北边界点的 y 坐标 y1 和 y2，这样就得到了表 3.8 的数据 （单位 ： mm ）。 根据地图的比例我们知道 18mm 相当于 40km ，试由测量数据计算该国国土的近似面 积，与它的精确值 41288km 2 做比较 。 表 3.8 地图边界点数据 x 7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 61.0 68.5 76.5 80.5 91.0 y1 44 45 47 50 50 38 30 30 34 36 34 41 45 46 y2 44 59 70 72 93 100 110 110 110 117 118 116 118 118 x 96 101 104 106.5 111.5 118 123.5 136.5 142 146 150 157 158 y1 43 37 33 28 32 65 55 54 52 50 66 66 68 y2 121 124 121 121 121 122 116 83 81 82 86 85 68 曲线绘制 相关代码如下 ： x=[7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 61.0