资源描述：
《插值运算的matlab函数》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、插值运算的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)在区间［010］的曲线，并在曲线上插值节点xk=k,k=0,1...1
2、0及函数值，画出分段线性插值折线图x=0:10;y=sin(x);xi=0:0.25:10;yil=interp1(x,y,xi,'nearest');yi2=interpl(x,y,xi,'linear');yi3=interpl(x,y,xi,'spline');yi4=interp1(x,y,xi,'pchip‘)；yi5=interpl(x,y,xi,'cubic*);subplot(l,5,l)plot(x,y,'o；xi,yi1,'k—；xi,sin(xi),k‘)；titleCbfNearest');su
3、bplot(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
4、),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为边界类型，缺省为非扭结边界条件；val
5、conds表示边界值。边界类型:Complete1为给定边界条件的一阶导数；5oZ・kiior为非扭结；，period©为周期边界条件；teconcT为给定边界条件的二阶导数;、ariational，为自由边界例子：(1)已知y二f(x)的函数表及端点条件S,,(xl)=S,,(x4)=0X1245F(x)1342求三次样条插值函数S(x),并计算f(3),f(4.5)的近似值clear;clc;x=[l245];y=[l342];s=csape(x,y/variation,)value=fnval(s,[34.5])再
6、输入s.coefs可以得到三次样条插值分段表示的系数注释：V=FNVAL(F,X)orFNVAL(X,F)providesthevalueatthepointsinXofthefunctiondescribedbyF.(1)已知函数y=l/(25xA2+l)在[01]上的值如下表X00.250.50.751Y10.39030.13790.06640.0385求三次样条插值函数S(x),使满足S，(0)=()S(1)=0.074x=[0:0.25:lJ;y=l./(l+25*x.A2);s=csape(x,y,'compl
7、ete',[0-0.074])fnpk(s,T)s.coefs2高维插值函数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，为双三次插值。例子：利用函数peak
8、s产生一个山顶曲而数据，通过插值做岀更为紧密的山顶曲面clear;clc;[x,y,z]=peaks(10);subplot(1,2,1)mesh(x,y,z),holdon;plot3(x,y,z,'r**);subplot(l,2,2)[xi,yi]=meshgrid(-3:0.1:3,-3:0.1:3);zi=
