资源预览内容
第1页 / 共4页
第2页 / 共4页
第3页 / 共4页
第4页 / 共4页
亲,该文档总共4页全部预览完了,如果喜欢就下载吧!
资源描述
拉格朗日插值的调用函数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= lagrange(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)
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号