数模培训_matlab插值方法解析.ppt

function y=lagrange(x0,y0,x) ii=1:length(x0); y=zeros(size(x)); for i=ii ij=find(ii~=i); y1=1; for j=1:length(ij), y1=y1.*(x-x0(ij(j))); end y=y+y1*y0(i)/prod(x0(i)-x0(ij)); end 算例：给出f(x)=ln(x)的数值表，用Lagrange计算ln(0.54)的近似值。 x=[0.4:0.1:0.8]； y=[-0.916291,-0.693147,-0.510826,-0.356675,-0.223144]; lagrange(x,y,0.54) ans = -0.6161 （精确解-0.616143) Runge现象 问题的提出：根据区间[a,b]上给出的节点做插值多项式p(x)的近似值，一般总认为p(x)的次数越高则逼近f(x)的精度就越好，但事实并非如此。 反例： 在区间[-5,5]上的各阶导数存在，但在此区间上取n个节点所构成的Lagrange插值多项式在全区间内并非都收敛。 取n=10,用Lagrange插值法进行插值计算。 x=[-5:1:5]; y=1./(1+x.^2); x0=[-5:0.1:5]; y0=lagrange(x,y,x0); y1=1./(1+x0.^2); ％绘制图形 plot(x0,y0,--r) ％插值曲线 hold on plot(x0,y1,‘-b) ％原曲线 * 再输入以下命令: xi=1:0.2:5; yi=1:0.2:3; zi=interp2(x,y,temps,xi,yi,cubic); mesh(xi,yi,zi) 画出插值后的温度分布曲面图. To MATLAB Exam_3 * To MATLAB exam_5 返回 * 插值函数griddata格式为: cz =griddata（x，y，z，cx，cy，‘method’） 用MATLAB作散点数据的插值计算 要求cx取行向量，cy取为列向量。 被插值点 插值方法 插值节点 被插值点的函数值 ‘nearest’ 最邻近插值 ‘linear’ 双线性插值 ‘cubic’ 双三次插值 v4- Matlab提供的插值方法 缺省时, 双线性插值 * To MATLAB exam_6 返回 * 作业1：在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）*（-50，150）里的哪些地方船要避免进入。 * 3.做海底曲面图 To MATLAB hd1 返回 4.作出水深小于5的海域范围,即z=5的等高线. * 作业2：山区地貌：在某山区测得一些地点的高程如下表：(平面区域1200=x=4000,1200=y=3600)，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较。 返回 * 作业3： 返回 * 作业4： 返回 * 作业5： 返回 * 作业6： 返回 * * 数学建模与数学实验 插 值 * 实验目的 实验内容 2、掌握用数学软件包求解插值问题。 1、了解插值的基本内容。 [1]一维插值 [2]二维插值 [3]实验作业 * 拉格朗日插值 分段线性插值 三次样条插值 一 维 插 值 一、插值的定义 二、插值的方法 三、用Matlab解插值问题 返回 * 返回 二维插值 一、二维插值定义 二、网格节点插值法 三、用Matlab解插值问题 最邻近插值 分片线性插值 双线性插值 网格节点数据的插值 散点数据的插值 * 一维插值的定义 已知 n+1个节点 其中 互不相同，不妨设 求任一插值点 处的插值 ? ? ? ? ? 节点可视为由 产生,， 表达式复杂,， 或无封闭形式,， 或未知.。 ? * 构造一个(相对简单的)函数 通过全部节点, 即 再用 计算插值，即 ? ? ? ? ? ? 返回 * 称为拉格朗日插值基函数。 已知函数f(x)在n+1个点x0,x1,…,xn处的函数值为 y0,y1,…,yn 。求一n次多项式函数Pn(x)，使其满足： Pn(xi)=yi,i=0,1,…,n. 解决此问题的拉格朗日插值多项式公式如下 其中Li(x) 为n次多项式： 拉格朗日(Lagrange)插值 * 拉格朗日(Lagrange)插值 特别地: 两点一次(线性)插值多项式: 三点