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' 立方插值。
所有的插值
用多项式函数(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' 立方插值。
所有的插值