拉格朗日插值matlab程序.pdf

拉格朗日插值matlab程序

拉格朗日插值的调用函数

functiony=lagrange(x0,y0,x)

n=length(x0);m=length(x);

fori=1:m

z=x(i);

L=0.0;

forj=1:n

T=1.0;

fork=1:n

ifk~=j

T=T*(z-x0(k))/(x0(j)-x0(k));

end

end

L=T*y0(j)+L;

end

y(i)=L;

end

四个图在一起：

x=[-1:0.05:1];

y=1./(1+25*x.^2);

x0=[-1:0.4:1];

y0=1./(1+25*x0.^2);

y1=lagrange(x0,y0,x);

x0=[-1:0.2:1];

y0=1./(1+25*x0.^2);

y2=lagrange(x0,y0,x);

x0=[-1:0.1:1];

y0=1./(1+25*x0.^2);

y3=lagrange(x0,y0,x);

plot(x,y,-r)

holdon

plot(x,y1,-b,x,y2,-r,x,y3,-r)

1/4

拉格朗日插值matlab程序

l5和fx在一起：

x=[-1:0.05:1];

y=1./(1+25*x.^2);

x0=[-1:0.4:1];

y0=1./(1+25*x0.^2);

y1=lagrange(x0,y0,x);

plot(x,y,-r)

holdon

plot(x,y1,-b)

l10和fx在一起：

x=[-1:0.05:1];

y=1./(1+25*x.^2);

x0=[-1:0.2:1];

y0=1./(1+25*x0.^2);

y2=lagrange(x0,y0,x);

plot(x,y,-r)

holdon

plot(x,y2,-b)

l20和fx在一起：

x=[-1:0.05:1];

y=1./(1+25*x.^2);

x0=[-1:0.1:1];

y0=1./(1+25*x0.^2);

y3=lagrange(x0,y0,x);

plot(x,y,-r)

holdon

plot(x,y3,-b)

2/4

拉格朗日插值matlab程序

3/4

拉格朗日插值matlab程序

4/4