MATLAB对一组数据进行插值的方法

1、拉格朗日插值
用多项式函数（10.2）作为插值函数时，希望通过解方程组（10.3）而得到待定系数


function y=lagrange(x0,y0,x);
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
en


2、分段线性插值
用Matlab实现分段线性插值不需要编制函数程序，Matlab中有现成的一维插值函数interp1。
y=interp1(x0,y0,x,'method')
method指定插值的方法，默认为线性插值。其值可为：
'nearest' 最近项插值
'linear' 线性插值
'spline' 立方样条插值'cubic' 立方插值。
所有的插值方法要求x0是单调的。
当x0为等距时可以用快速插值法，使用快速插值法的格式
为 '*nearest'、'*linear' 、'*spline' 、'*cubic'


3、三次样条曲线插值


Matlab中三次样条插值也有现成的函数：
y=interp1(x0,y0,x,'spline')；
y=spline(x0,y0,x)；
pp=csape(x0,y0,c