插值实验报告.docxVIP

新乡学院数学与信息科学系 实验报告 实验项目名称 插值实验 所属课程名称 数学实验 实 验 类 型 综合性实验 实 验 日 期 2013-6-6 班 级 数学与应用数学一班 学 号 11111011013 姓 名 高亚丹 成 绩 PAGE PAGE 10 一、实验概述： 【实验目的】 了解插值的基本内容以及插值方法，如：一维插值（拉格朗日插值、分段线性插值、三次样条插值），二维插值（最临近插值、分段线性插值、双线性插值）。 【实验原理】 拉格朗日(Lagrange)插值: 已知函数 f(x)在 n+1 个点 x0,x1,…,xn 处的函数值为 y0,y1,…,yn ．求一 n 次多项式函数 Pn(x)，使其满足： Pn(xi)=yi,i=0,1,…,n. 解决此问题的拉格朗日插值多项式公式如下 P (x) ? ?n L (x) ? y n i i i?0 其中 Li(x) 为 n 次多项式： L (x) ? (x ? x 0 )(x ? x 1 ) (x ? x ??i?1 ? ? )(x ? x i?1 ) (x ? x ) ??n ? ? i (x i x )(x i x ) (x i x i?1 )(x i x i?1 ) (x i x ) n 特别地: 两点一次(线性)插值多项式: L ?x?? x ? x 1 y x ? x 0 y 1 x ? x 0 0 1 x ? x 1 1 0 三点二次(抛物)插值多项式: ? ? ?x ? x ?? ?x ? x ? ?x ? x ?? ?x ? x ? ?x ? x ?? ?x ? x ? 2L x ? ?x 2 1 xx? x ? x x x ? x ? y0 ? ?x 0 ? ?? ?2? ? ? ? ? 2 ? 2 ? y ? ?x 0 1? ?2? x ? 1 ? ? 2 ? 1 ? y 0 1 0 2 1 0 1 2 2 0 2 1 分段线性插值: L (x) ? ?n n j ?0 y l (x) j j ? x ? x ? j ?1 , x ? x ? x ? x j x j ?1 j ?1 j ? x ? x l (x) ? ? j ?1 , x ? x ? x j ? x ? x j j ?1 ? j j ?1 ?0, 其他 ?? 计算量与 n 无关; n 越大，误差越小. lim L (x) ? g (x), x ? x ? x n?? n 0 n 一 维 插 值 函 数 ： yi=interp1(x,y,xi,’method’) yi 为 xi 处得插值结果 x、y 为插值节点xi 为被插值点 method 为插值方法，其中包括 nearest（最邻近插值）、linear（线条插值）、spline （三次样条插值）、cubic（立方插值） 4. 二 维 插 值 函 数 ： z=interp2(x0,y0,z0,x,y,’method’) z 为被插值点函数值 x0、y0 为插值节点 x、y 为被插值点 method 为插值方法，其中包括 nearest（最邻近插值）、linear（双线性插值）、cubic （双三次插值） 【实验环境】 Matlab 7.0 Microsoft Windows7 Professional 版本 2002 Service Pack 3 二、实验内容： 【实验方案】 用不同的方法对已知函数进行插值，用 matlab 画出相应的图形。 【实验过程】（实验步骤、记录、数据、分析） 1 山区地貌：在某山区测得一些地点的高程如下表： (平面区域 1200 ≤x ≤ 4000,1200≤y ≤3600)，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较． 3600 1480 1500 1550 1510 1430 1300 1200 980 3200 1500 1550 1600 1550 1600 1600 1600 1550 2800 1500 1200 1100 1550 1600 1550 1380 1070 2400 1500 1200 1100 1350 1450 1200 1150 1010 2000 1390 1500 1500 1400 900 1100 1060 950 1600 1320 1450 1420 1400 1300 700 900 850 1200 1130 1250 1280 1230 1040 900 500 700 y/x 1200 1600 2000 2400 2800 3200 3600 4000 通过此例对最近邻点插值、双线性插值方法和双三次插值方法的插值效果进行比较． 解：最邻近插值： z=[1