插值运算的matlab函数.docxVIP

插值运算的matlab函数 1一维插值函数 interpl () 命令格式： yi=interpl (x,y,xi,method) x为插值节点构成的向量，y为插值节点函数值构成的向量，yi是被插值点xi的 插值结果，Fethod，是采用的插值方法，缺省时表示分线段性插值，Nearest4为 最邻近插值；Tinea/为分线段性插值spline，为三次样条插值；，pchip，为分段 Hermite 插值；’cubic为分段 Hermite 插值 例子：I田i出y=sin(x)在区间［0 10］的曲线，并在曲线上插值节点xk=k, k=0,1...10 及函数值，画出分段线性插值折线图 x=0:10; y=sin(x); xi=0:0.25:10; yil =interp 1 (x,y,xi,nearest); yi2=interpl(x,y,xi,linear); yi3=interp l(x,y,xi,spline); yi4=interp 1 (x,y,xi,pchip‘)； yi5=interpl(x,y,xi,cubic*); subplot(l,5,l) plot(x,y,o；xi,yi 1 ,k—；xi,sin(xi),k‘)； titleC\bfNearest); subplot(l,5,2) plot(x,y,o,xi,yi2,k?-；xi,sin(xi),k:‘)； title(\bfLinear); subplot(l,5,3) plot(x,y,o；xi,yi3,k—；xi,sin(xi),k‘)； titleCAbfSpline*); subplot(l,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=interp 1 (x,y,xi,spline) 命令格式2： pp=spline(x,y),输出三次样条函数分段表示的结构 pchip()命令格式与spline()完全相同 csape()为可输入边界条件的三次样条函数 命令格式：pp=csape(x,y,conds,valconds),x为插值节点构成的向量，y为插值节点 函数值构成的向量；conds为边界类型，缺省为非扭结边界条件；valconds表示 边界值。 边界类型:Complete1为给定边界条件的一阶导数；5oZ?kiior为非扭 结；，period?为周期边界条件；teconcT为给定边界条件的二阶导数;、ariational， 为自由边界 例子： (1)已知y二f(x)的函数表及端点条件S,,(xl)=S,,(x4)=0 X 1 245 F(x) 1342 求三次样条插值函数S (x),并计算f (3) , f(4.5)的近似值 clear; clc; x=[l 2 4 5]; y=[l 3 4 2]; s=csape(x,y/variation,) value=fnval(s,[3 4.5]) 再输入s.coefs可以得到三次样条插值分段表示的系数 注释：V = FNVAL(F,X) or FNVAL(X,F) provides the value atthe points in X of the function described by F . (1)已知函数y=l/(25xA2+l)在[0 1]上的值如下表 X 0 0.25 0.5 0.75 1 Y 1 0.3903 0.1379 0.0664 0.0385 求三次样条插值函数S(x),使满足S，(0)=()S(1)=0.074 x=[0:0.25:lJ; y=l./(l+25*x.A2); s=csape(x,y,complete,[0 -0.074]) fnpk(s,T) s.coefs 2高维插值函数 interp2()二维插值函数 命令格式:zi=interp2(X,YZ,xi,yi,method) X为插值节点构成的矩阵(X二严x,e=(l,l,???,l))； Y为插值节点构成的矩阵 (Y=y%)； Z为插值节点函数值构成的矩阵，Zij=f(Xi,Yj) rmethod4^采用的插值 方法，缺省时表示分双线性插值，Nearest4为最邻近插值「lineal为分线段性插 值spline，为三次样条插值；’cubic，为双三次插值。 例子：利用函数peaks产生一个山顶曲而数据，通过插值做岀更为紧密的山顶曲 面 clear; clc; [x,y,z]=pe