拉格朗日插值matlab程序

《拉格朗日插值matlab程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、拉格朗日插值的调用函数functiony=lagrange(xO,yO,x)n=length(xO);m=length(x);fori=l:mz=x(i);L=0.0;forj=l:nT=1.0;fork=l:nifk~=jT=T*(z-xO(k))/(xO(j)-xO(k));endenclL=T*yO(j)+L;endy(i)=L;end四个图在一起：x=[-l:0.05:1];y=l./(l+25*x."2);x0=[-l:0.4:1];yO二l./(l+25*x0.‘2);yl=lagrange(xO,yO,x);x0=[-l:0.2:l];y0=

2、l./(l+25*x0."2);y2=lagrange(x0,yO,x);x0=[-1:0.1:1];y0=l./(l+25*x0.⑵;y3=lagrange(xO,yO,x);plot(x,y,'-r)holdonplot(x,yl,'一b‘,x,y2,'-r",x,y3,'-F)15和fx在一起:x二卜1:0.05:1];y=l./(l+25*x.A2);x0=[-l:0.4:l];y0=l./(l+25*x0.A2);yl=lagrange(xO/yO,x);plot(xzy；-r')holdonplot(x,yU・b‘)110和fx在一起：x二卜1

3、:0.05:1];y=l./(l+25*x.A2);x0=[-l:0.2:l];yO=l./(l+25*xO.A2);y2=lagrange(xO,yO,x);plot(x,y；-r')holdonplot(x,y2,'・bj120和fx在一起:x二卜1:0.05:1];y=l./(l+25*x.A2);x0=[-l:0.1:l];yO=l./(l+25*xO.A2);y3=lagrange(x0,y0,x);plot(x,y/-r')holdonplot(x,y3〉b‘)E3Figure：[o丨回「1=1回^