插值运算matlab函数.docVIP

插值运算的matlab函数1一维插值函数nterp1（）命令格式：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表示边界值。边界类型：‘complete‘为给定边界条件的一阶导数；’not-a-knot‘为非扭结；’periodic‘为周期边界条件；’second’为给定边界条件的二阶导数;’variational’为自由边界例子：（1）已知y=f(x)的函数表及端点条件S’’(x1)=S’’(x4)=0X 1 2 4 5F(x) 1 3 4 2求三次样条插值函数S（x）,并计算f（3），f(4.5)的近似值clear;clc;x=[1 2 4 5];y=[1 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=1/(25x^2+1)在[0 1]上的值如下表X 0 0.25 0.5 0.75 1Y 1 0.3903 0.1379 0.0664 0.0385求三次样条插值函数S(x)，使满足S’(0)=0,S’(1)=-0.074x=[0:0.25:1];y=1./(1+25*x.^2);s=csape(x,y,complete,[0 -0.074])fnplt(s,r)s.coefs2高维插值函数interp2()二维插值函数命令格式：zi=interp2(X,Y,Z,xi,yi,’method’)X为插值节点构成的矩阵(X=e’*x,e=(1,1,…,1))；Y为插值节点构成的矩阵(Y=y’*e)；Z为插值节点函数值构成的矩阵，Zij=f(Xi,Yj) ;‘method‘是采用的插值方法，缺省时表示分双线性插值，’nearest‘为最邻近插值；’linear‘为分线段性插值;’spline’为三次样条插值；’cubic’为双三次插值。例子：利用函数peaks产生一个山顶曲面数据，通过插值做出更为紧密的山顶曲面clear;clc;[x,y,z]=peaks(10);s