《MATLAB实现拉格朗日插值69840.doc》由会员分享,可在线阅读,更多相关《MATLAB实现拉格朗日插值69840.doc(2页珍藏版)》请在三一文库上搜索。
MATLAB实现拉格朗日插值建立如下拉格朗日插值函数:function y=lagrange(x0,y0,x);n=length(x0);m=length(x);for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j); end end s=p*y0(k)+s; end y(i)=s;end画图程序如下:x=-5:1:5;y=1./(1+x.2);x0=-5:0.001:5;y0=lagrange(x,y,x0);y1=1./(1+x0.2);plot(x0,y0,r)hold onplot(x0,y1,g)注:画出的图形为n =10的图形得到图形如下:n=10的图像