插值matlab.ppt

插值问题;问题的提出;插值在数学发展史上是个老问题。它和Lagrange、Newton、Gauss等著名数学家的名字连在一起，它最初来源于天体计算——由若干观测点（即节点）计算任意时刻星球的位置（即插值点和插值）——的需要。 插值仍然在诸如图像处理、机械加工和数据处理等科学研究中有着许多直接的应用； 插值又是数值积分、数值微分、常微分方程数值解等数值计算的基础。; 例如98年美国赛A题，生物组织切片的三维插值处理，94年A题逢山开路，山体海拔高度的插值计算。2005年的雨量预报的评价的插值计算。;插值问题的提法：; 构造一个(相对简单的)函数;拉格朗日插值;2、抛物插值;3、n次拉格朗日插值;上述多项式称为 n次拉格朗日（Lagrange）插值多项式， ;（2）拉格朗日插值的matlab实现;（3） 函数的调用;（4）练习;随着节点的增加，多数情况时误差越来越小，但有时会出现很??的震荡，这种震荡称为龙格现象。高次插值多项式的这些缺陷，促使人们转而寻求简单的低次多项式插值。;龙格现象;;;分段线性插值;例;x=linspace(-6,6,100); y=1./(x.^2+1); x1=linspace(-6,6,5); y1=1./(x1.^2+1); plot(x,y,x1,y1,x1,y1,‘o’,‘LineWidth’,1.5)；;它具有良好的收敛性，用它计算插值时，只用到x左右两个节点，计算量与节点个数n无关，但n越大，分段越多，插值误差越小。实际上用函数表作插值计算时，分段线性插值就足够了，如数学、物理中用的特殊函数表，数理统计中用的概率分布表等。;实用插值方法;三次样条函数记作S(x)(a≤x≤b)，满足以下条件： 1、在每个小区间上是三次多项式； 2、在a≤x≤b上二阶导数连续； 3、 它也具有良好的收敛性。;输入:节点x0,y0, 插值点x (均为数组，长度自定义)； 输出:插值y (与x同长度数组)。 ;例; 例：在1-12的11小时内，每隔1小时测量一次温度，测得的温度依次为：5，8，9，15，25，29，31，30，22，25，27，24。试估计每隔1/10小时的温度值。; 二维插值; 要求x0,y0单调；x，y可取为矩阵，或x取行向量，y取为列向量，x,y的值分别不能超出x0,y0的范围．;例：测得平板表面3×5网格点处的温度分别为： 82 81 80 82 84 79 63 61 65 81 84 84 82 85 86 试作出平板表面的温度分布曲面z=f(x,y)的图形．;再输入以下命令: xi=1:0.2:5; yi=1:0.2:3; zi=interp2(x,y,temps,xi,yi,cubic); mesh(xi,yi,zi) 画出插值后的温度分布曲面图. ;山体地貌;程序设计： clear x0=[1200:400:4000]; y0=[1200:400:3600]; z0=[1130 1250 1280 1230 1040 900 500 700; 1320 1450 1420 1400 1300 700 900 850; 1390 1500 1500 1400 900 1100 1060 950; 1500 1200 1100 1350 1450 1200 1150 1010; 1500 1200 1100 1550 1600 1550 1180 1070; 1500 1550 1600 1550 1600 1600 1600 1550; 1480 1500 1550 1510 1430 1300 1200 980]; xi=1200:1:4000; %加密数据点 yi=1200:1