数值分析实验2(wangwei).docVIP

实验2.1（多项式插值的振荡现象） 问题提出：考虑在一个固定的区间上用插值逼近一个函数。显然拉格朗日插值中使用的节点越多，插值多项式的次数就越高。我们自然关心插值多项式的次数增加时，是否也更加靠近被逼近的函数，龙格给出的一个例子是极著名并富有启发性的，设区间[-1，1]上函数 实验内容:考虑区间[-1，1]的一个等距划分，分点为 则拉格朗日插值多项式为 其中的是次拉格朗日插值基函数。 实验要求： （1）选择不断增大的分点数目画出原函数及插值多项式函数在 [-1，1]上的图像，比较并分析实验结果。 （2）选择其他的函数，例如定义在区间[-5，5]上的函数 重复上述的实验看其结果如何。 （3）区间[a,b]上切比雪夫点的定义为 以为插值节点构造上述各函数的拉格朗日插值多项式，比较其结果。 程序清单： 被逼近函数的函数文件func1.m function y=func1(x,c) %直接用被逼近函数计算函数值，c用来选择函数 if c==1; y=1./(1+25*x.^2); end; if c==2; y=x./(1+x.^4); end; if c==3; y=atan(x); end; 2．拉格朗日插值函数文件lagr.m function yy=lagr(x,y,xx) %用拉格朗日插值多项式进行插值,x为插值点的自变量列阵，y为对应的函数值列阵，xx为 %待插值的自变量列阵，yy为求得的对应于xx的函数值列阵 yy=xx-xx;%初始化为0向量 n=length(x)-1; for i=0:n; z=yy; k=0; for j=0:n; if j~=i k=k+1; if k==1; z=(xx-x(j+1))/(x(i+1)-x(j+1)); else z=z.*(xx-x(j+1))/(x(i+1)-x(j+1)); end; end; end; yy=yy+z*y(i+1); end; 主程序main1.m m=input(请输入节点数目:); chof=input(请选择函数:(1:1/(1+25*x^2)) 2:x/(1+x^4) 3:arctanx)); div=input(请选择节点方式:(1.均布节点 2.切比雪夫点)); %生成用于画图的自变量序列 if chof==1; t=-1:0.01:1; else; t=-5:0.01:5; end; y1=func1(t,chof);%直接求t对应的函数值 %按要求构造节点的自变量序列 n=m-1; for i=0:n; if div==1; if chof==1; nod(i+1)=-1+2*i/n; else; nod(i+1)=-5+10*i/n; end; end; if div==2; if chof==1; nod(i+1)=cos((2*i+1)*pi/(2*(n+1))); else nod(i+1)=5*cos((2*i+1)*pi/(2*(n+1))); end; end; end; nodv=func1(nod,chof);%求节点处函数值 y=lagr(nod,nodv,t); plot(t,y1,t,y); t=0; y=0; y1=0; nod=0; nodv=0; 实验结果及其分析： 分点数分别为2，3，5，8，11，15时的函数和插值多项式函数的图像如下： n=2: n=3 n=5: n=8: n=11: n=15: 分析：从图上可以看出随着分点数目的增加区间端部的插值函数的值与原函数的误差越来达，而区间中部的误差越来越小，这就是所谓的龙格现象，这是等距节点的高次插值多项式的典型病态现象。 对于函数情况如下： n=4: n=8: n=11: n=15: 对于函数情况如下： n=3: