MATLAB插值实验报告数学实验.docVIP

新乡学院数学与信息科学系 实验报告 实验项目名称 插值实验 所属课程名称 数学实验 实 验 类 型 综合性实验 实 验 日 期 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 掌握用Matlab插值的方法，了解拉格朗日插值、线性插值、样条插值的基本思想，了解三种网格节点数据的插值方法的基本思想，了解掌握用Matlab计算一维差值和二维插值的方法。 【实验原理】 拉格朗日(Lagrange)插值。 已知函数f(x)在n+1个点x0,x1,…,xn处的函数值为 y0,y1,…,yn ．求一n次多项式函数Pn(x)，使其满足： Pn(xi)=yi,i=0,1,…,n. 解决此问题的拉格朗日插值多项式公式如下 其中Li(x) 为n次多项式： 称为拉格朗日插值基函数． 【实验环境 】 Matlab 7.0.1 Windows 8 专业版 二、 实验内容： 【实验方案】 1.利用interp2命令，以及画图命令。 2.会用各种插值命令。 3.熟悉并掌握拉格朗日插值。 【实验过程】（实验步骤、记录、数据、分析） 1.山区地貌：在某山区测得一些地点的高程如下表：(平面区域1200≤x ≤4000,1200≤y ≤3600)，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较． 3600 3200 2800 2400 2000 1600 1200 1480 1500 1550 1510 1430 1300 1200 980 1500 1550 1600 1550 1600 1600 1600 1550 1500 1200 1100 1550 1600 1550 1380 1070 1500 1200 1100 1350 1450 1200 1150 1010 1390 1500 1500 1400 900 1100 1060 950 1320 1450 1420 1400 1300 700 900 850 1130 1250 1280 1230 1040 900 500 700 y/x 1200 1600 2000 2400 2800 3200 3600 4000 解：原始图程序： x=1200:400:4000; y=1200:400:3600; z=[1480 1500 1550 1510 1430 1300 1200 980; 1500 1550 1600 1550 1600 1600 1600 1550; 1500 1200 1100 1550 1600 1550 1380 1070; 1500 1200 1100 1350 1450 1200 1150 1010; 1390 1500 1500 1400 900 1100 1060 950; 1320 1450 1420 1400 1300 700 900 850; 1130 1250 1280 1230 1040 900 500 700]; mesh(x,y,z); title(原始图),xlabel(x),ylabel(y),zlabel(z) 　　　　 最邻近插值图：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600); zi=interp2(x,y,z,xi,yi,nearest); mesh(xi,yi,zi) title(最邻近插值图),xlabel(x),ylabel(y),zlabel(z) 线性插值：[xi,yi]=meshgrid(1200:10:4000,1200:10:3600); zi=interp2(x,y,z,xi,yi,linear); mesh(xi,yi,zi) title(线性插值),xlabel(x),ylabel(y),zlabel(z)