数值分析第二次上机(共3页).docx

《数值分析第二次上机(共3页).docx》由会员分享，可在线阅读，更多相关《数值分析第二次上机(共3页).docx（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上数值分析第二次上机题目：在区间-1,1上分别取n=10、20用两组等距节点对龙格函数做多项式插值及三次样条插值，对每个n值，分别画出插值函数及的图形程序：function y1=lagrange(x0,y0,x1)n=length(x0);syms x;for k=1:n l(k)=x/x; for p=1:n if p=k l(k)=l(k)*(x-x0(p)/(x0(k)-x0(p); end endendz=0;for k=1:n z=z+l(k)*y0(k);endy1=subs(z,x,x1);程序：n=10时n=10;x0=-1:2/n:1;y0=1./(

2、1+25*x0.2);x=-1:.001:1;y1=lagrange(x0,y0,x);y=1./(1+25*x.2);plot(x,y,x,y1,-.,x0,y0,p);legend(Runge Function,插值函数,插值节点);title(n=10时的Lagrange插值的龙格现象);xlabel(x);ylabel(y);结果：程序：n=20时n=20;x0=-1:2/n:1;y0=1./(1+25*x0.2);x=-1:.001:1;y1=lagrange(x0,y0,x);y=1./(1+25*x.2);plot(x,y,x,y1,-.,x0,y0,p);legend(Runge Function,插值函数,插值节点);title(n=20时的Lagrange插值的龙格现象);xlabel(x);ylabel(y);结果：专心-专注-专业