spline插值【DOC精选】.docx

例5.6.1给定以下数据, 求出三次样条函数，并计算函数分别在-0.15，-0.05，0.05，0.18，0.25处的近似值，并作图。x0.10.20.150-0.20.3y0.950.840.861.060.500.72解：编程如下:clearx=[0.1,0.2,0.15,0,-0.2,0.3];y=[0.95,0.84,0.86,1.06,1.50,0.72];pp=spline(x,y);pp.coefsxx=[-0.15,-0.05,0.05,0.18,0.25];yy=ppval(pp,xx) %or:yy=spline(x,y,xx)fnplt(pp,k)hold on plot(x,y,o)hold on plot(xx,yy,r*)运行结果:ans = -36.3850 21.8592 -5.1164 1.5000 -36.3850 0.0282 -0.7390 1.0600227.6995 -10.8873 -1.8249 0.9500 -143.0047 23.2676 -1.2059 0.8600 -143.0047 1.8169 0.0484 0.8400yy = 1.2943 1.1016 1.0186 0.8409 0.8291写出三次样条函数为:练习: 给定函数，，节点，用三次样条插值求及在某些点处的近似值.MATLAB实现:yi=interp1(x,y,xi) %根据数据(x,y)给出分段线性插值函数的xi的值yi.yi=interp1(x,y,xi,’spline’) %使用三次样条插值.yi=interp1(x,y,xi,’cubic’) %使用分段三次插值.Yi=spline(x,y,xi) %同yi=interp1(x,y,xi,’spline’) pp=spline(x,y) %返回三次样条插值的分段多项式(pp)形式结构(使用“非扭结”端点条件).MATLAB中的spline函数可以对数据点进行三次插值，默认的边界条件强制使得插值第一段的三次项系数与第二段的三次项系数相同，最后一段与倒数第二段的多项式系数相同，这就所谓的非扭结（not-a-knot）条件.1一维插值函数?Interp1（）命令格式：yi=interp1（x,y,xi,’method’）x为插值节点构成的向量，y为插值节点函数值构成的向量，yi是被插值点xi的插值结果，‘method‘是采用的插值方法，缺省时表示分线段性插值，’nearest‘为最邻近插值；’linear‘为分线段性插值;’spline’为三次样条插值；’pchip’为分段Hermite插值；’cubic’为分段Hermite插值例子：画出y=sin(x)在区间[0?10]的曲线，并在曲线上插值节点xk=k，k=0,1…10及函数值，画出分段线性插值折线图x=0:10;y=sin(x);xi=0:0.25:10;yi1=interp1(x,y,xi,nearest);yi2=interp1(x,y,xi,linear);yi3=interp1(x,y,xi,spline);yi4=interp1(x,y,xi,pchip);yi5=interp1(x,y,xi,cubic);subplot(1,5,1)plot(x,y,o,xi,yi1,k--,xi,sin(xi),k:);title(\bfNearest);subplot(1,5,2)plot(x,y,o,xi,yi2,k--,xi,sin(xi),k:);title(\bfLinear);subplot(1,5,3)plot(x,y,o,xi,yi3,k--,xi,sin(xi),k:);title(\bfSpline);subplot(1,5,4)plot(x,y,o,xi,yi4,k--,xi,sin(xi),k:);title(\bfPchip);subplot(1,5,1)plot(x,y,o,xi,yi5,k--,xi,sin(xi),k:);title(\bfCubic);?spline()为三次样条函数命令格式1：yi=spline(x,y,xi),意义等同于yi=interp1(x,y,xi,spline)命令格式2：pp=spline(x,y)?，输出三次样条函数分段表示的结构pchip()命令格式与spline()完全相同?csape()为可输入边界条件的三次样条函数命令格式：pp=csape(x,y,conds,valconds),x为插值节点构成的向量，y为插值节点函数值构成的向量；conds为边界类型，缺省为非扭结边界条件；valconds表示边界