《拉格朗日插值matlab程序》由会员分享,可在线阅读,更多相关《拉格朗日插值matlab程序(4页珍藏版)》请在金锄头文库上搜索。
1、拉格朗日插值的调用函数function y=lagrange(x0,y0,x)n=length(x0);m=length(x);for i=1:mz=x(i);L=0.0;for j=1:nT=1.0;for k=1:nif k=jT=T*(z-x0(k)/(x0(j)-x0(k);endendL=T*y0(j)+L;endy(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= lagrang
2、e(x0,y0,x);x0=-1:0.1:1;y0=1./(1+25*x0.2);y3= lagrange(x0,y0,x);plot(x,y,-r)hold onplot(x,y1,-b,x,y2,-r,x,y3,-r)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) hold on 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) hold on 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) hold on plot(x,y3,-b)