資源描述:
《拉格朗日插值matlab程序》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、拉格朗日插值的調(diào)用函數(shù)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四個(gè)圖在一起: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回^