插值和拟合资料讲解.docx

插值和拟合插值和拟合都是函数逼近或者数值逼近的重要组成部分他们的共同点都是通过已知一些离散点集M上的约束，求取一个定义在连续集合S（M包含于S）的未知连续函数，从而达到获取整体规律的目的，即通过"窥几斑"来达到"知全豹"简单的讲，所谓拟合是指已知某函数的若干离散函数值{f1,f2,…，f n},通 过调整该函数中若干待定系数f（入1,入2,…，入3）,使得该函数与已知点集的 差别（最小二乘意义）最小如果待定函数是线性，就叫线性拟合或者 线性回归（主要在统计中），否则叫作非线性拟合或者非线性回归表 达式也可以是分段函数，这种情况下叫作样条拟合而插值是指已知某函数的在若干离散点上的函数值或者导数信息，通过求解该函数中待定形式的插值函数以及待定系数，使得该函数在给定离散点上满足约束插值函数又叫作基函数，如果该基函数定义在整个定义域上，叫作全域基，否则叫作分域基如果约束条件中只有函数值的约束，叫作Lagrange插值，否则叫作Hermite插值从几何意义上将，拟合是给定了空间中的一些点，找到一个已知形式 未知参数的连续曲面来最大限度地逼近这些点；而插值是找到一个( 或几个分片光滑的)连续曲面来穿过这些点。

一、 概念的引入1. 插值与拟合在现实生活中的应用l 机械制造：汽车外观设计l采样数据的重新建构：电脑游戏中场景的显示，地质勘探，医学领域(CT)2. 概念的定义l插值：基于［a,b］区间上的n个互异点，给定函数f (x),寻找某个函数去逼近f (x)若要求e(x)在xi处与f (xi)相等，这类的函数逼近问题称为插值问题，xi即是插值点l 逼近： 当取值点过多时，构造通过所有点的难度非常大此时选择一个次数较低的函数最佳逼近这些点，一般采用最小二乘法l 光顾： 曲线的拐点不能太多，条件：①二阶几何连续②不存在多余拐点③曲率变化较小l 拟合：曲线设计过程中用插值或通过逼近方法是生成的曲线光滑(切变量连续) 光顾二、 插值理论设函数y=f(x)在区间［a,b］上连续，在［a,b］上有互异点x ,x ,…，x处取值0 1 ny ,y ,y如果函数e(x)在点x上满足e(x )二y (i=0,1,2,…，n),则称0 1 n i i i仅供学习与交流，如有侵权请联系网站删除谢谢3^(x)是函数y=f(x)的插值函数，x,x,…，x是插值节点若此时e(x)是代数0 1 n多项式P(x),则称P(x)为插值多项式。

显然f(x)~e(x), xW[a,b]1. 拉格朗日插值构造n次多项式P (x)= y l (x)=y l (x)+y l (x)+“・+y l (x),这是不超过 n k k 0 0 1 1 n nn次的多项式，其中基函数l (x) =k显然丨(x)满足丨(x)二k k i此时 P (x)^f(x),误差 R (x)=f(x)-P (x) =n n n其中 £ (a,b)且依赖于 x, = (x-x ) (x-x ”・(x-x )0 1 n很显然，当n=1、插值节点只有两个x ,x时k k+1P(x)=yl(x)+y l (x)1 k k k+1 k+1其中基函数l (x)二l (x) =k k+12. 牛顿插值构造 n 次多项式 N (x)=f(x )+f(x ,x )(x-x )+f(x ,x ,x )(x-x )(x-x )+•••n 0 0 1 0 0 1 2 0 1+f(x ,x ,x ,…，x )(x-x )(x-x ”・(x-x )0 1 2 n 0 1 n称为牛顿插值多项式,其中(二个节点,一阶差商)(三个节点,二阶差商)(n+1个节点，n阶差商)注意：由于插值多项式的唯一性，有时为了避免拉格朗日余项R (x)中n+1阶导 n数的运算，用牛顿插值公式R (x)=f(x)-N (x)=f(x,x，—, x )w (x),n n 0 n n+1其中 3 (X)=(X-X )(X-X )…(X-X )n+1 0 1 n3. 分段插值 子区间内，避免函数在某些区间失真1) 线性插值已知n+1个不同节点x,x,…，x，构造分段一次线性多项式P(x),使之满足0 1 nl P(x)在［a,b］上连续l P(X)=ykkl P(x) 在 ［x ,x ］上是线性函数,P(x) =i i+12) 两点带导数插值---避免尖点、一阶连续区间［a,b］上两个互异节点x,x，已知实数y ,y ,m ,m ，为了构造次数i i+1 i i+1 i i+1不大于3的多项式 满足条件引入 , 使之满足可以求出此时 = + ，其中4. 三次样条插值 二阶可导对于给定n+1个不同节点x ,x ,…，x及函数值y,y,…，y，其中0 1 n 0 1 na=x1〈…n二b。

构造三次样条插值函数S(x)S(x)称为三次样条函数时需满足:0l S(x) 在 ［a,b］上二阶导数连续l S(x)=y (k=0,1,…，n)kk丨每个子区间［x,x ］上 S(x)是三次多项式(k=0,1,…，n)k k+15.例题已知函数y=f(x)的观测值X1234Y0-5-63求三次插值多项式e(x)及e(2.5).解：(1)拉格朗日插值P(x)=yl(x)+yl(x)+yl(x)+yl(x)3 0 0 1 1 2 2 3 3=(-5) + (-6) +3 =x3-4x2+3 e(x)~P (x)= X3-4X2+3 0(2.5)=2.53-4*2.52+3=-6.3753(2)牛顿插值xf(x)一阶差商二阶差商三阶差商o'2-5-53-6-1243951N (x) =f (x ) +f (x ,x ) (x-x ) +f (x ,x ,x ) (x-x ) (x-x ) +f (x ,x ,x ,x ) (x-x ) (x-3 0 0 1 0 0 1 2 0 1 0 1 2 3 0x)(x-x)12=0+(-5)(x-1)+2*(x-1)(x-2)+1*(x-1)(x-2)(x-3) =x3-4x2+3三、Mat lab在插值中的应用1. Lagrange 插值1 )方法回顾对给定的n个节点x ,x ,…，x及对应的函数值y ,y ,…，y，利用n次Lagrange1 2 n 1 2 n插值多项式公式，插值区间内任意x的函数值y可以通过下式求出：2) Matlab 实现函数 Lagrange.mfunction y=lagrange(x0,y0,x)n=length(x0);m=length(x);for i=1:mz=x(i);s=0.0;for k=1:np=1.0;for j=1:nif j~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;end3) 例题的Matlab实现x=[1 2 3 4];y=[0 -5 -6 3];lagrange(x,y,2.5)2. Runge现象及分段线性插值1) Runge 现象Runge在本世纪初发现：在[-1,1]上用n +1个等距结点作插值多项式P (x),使 n其在个结点的值与函数y(x)=1/(1+25x2)在结点的值相等。

但在nT®时，插值多项式P (x)在区间中部趋于y(x)但对于0.726W|x|W1的x, P (x)严重 nn发散通过下面的例子，以图形的方式体会Runge现象(f(x)=1/(1+x2))x=[-5:1:5];y=1./(1+x•⑵;x0=[-5:0.1:5];y0=lagrange(x,y,x0);y1=1./(1+xO•⑵；%绘制图形plot(x0,y0,'--r')hold onplot(x0,y1,'-b')2) Matlab实现分段插值 维插值in terpll yi = i nterp1(x,y,xi)对(x,y)进行插值，计算插值点xi的函数值l yi = interp1(y,xi)默认x=1:n, n是向量y的元素个数l yi = interp1(x,y,xi,'method')指定特定算法插值，method可以是如下字符串0 linear线性插值0 spli ne三次样条插值0 cubic 三次插值要求：x是单调，但不要求连续等距如果x连续等距，可以选用快速插值法调用函数时只需在method前加，如"*spline"3) 例题① 用一维线性插值解决Runge现象y2=interp1(x,y,x0);plot(x0,y2,'*m')② 正弦曲线的插值示例x=0:0.1:10;y=sin(x);xi=0:0.25:10;yi=interp1(x,y,xi);plot(x,y,'o',xi,yi)3. Hermite 插值1)方法介绍已知n个插值节点x ,x ,…，x及起对应的函数值y ,y ,…，y和一阶导数值1 2 n 1 2 ny',y',…，y '，则计算插值区域内任意X的函数值y的Hermite插值公式1 2 n其中： ,2) Matlab 实现 Hermite.mfunction y=hermite(x0,y0,y1,x)n=length(x0);m=length(x);for k=1:myy=0.0;for i=1:nh=1.0;a=0.0;for j=1:n•广 ・~ ・if j~=ih二h*((x(k)-xO(j))/(xO(i)-xO(j))厂2;a=1/(x0(i)-x0(j))+a;endendyy=yy+h*((x0(i)-x(k))*(2*a*y0(i)-y1(i))+y0(i));endy(k)=yy;end3) 例题如下数据表，构造Hermite多项式，并求出Sin0.34的近似值x0.30 |0.320.35Sin(x)0.29552|0.314570.34290Cos(x)0.95534|0.949240.93937x0=[0.3 0.32 0.35];y0=[0.29552 0.31457 0.34290];y1=[0.95534 0.94924 0.93937];x=[0.3:0.005:0.35];y=hermite(x0,y0,y1,x);plot(x,y)y=hermite(x0,y0,y1,0.34);ysin(0.34)y2=sin(x);hold onplot(x,y2,'--r')4. 三次样条插值y3=interp1(x,y,x0,'*spline');y3=spline(x,y,x0);plot(x0,y3,'-g')四、 数据拟合1. 方法介绍在实际生活中，往往需要从一组实验数据(x ,y )中寻找出变量x,y之间的函ii数关系。

由于观测数据不可避免出现误差，因此并不需要y二f(x) —定要经过所 有的点，而只要求在给定点x 上误差△ i=f(x)-y按某种标准达到最小。