7.对于非光滑曲线,在上选择不同的和,用等距节点做次拉.docVIP

7．对于非光滑曲线，在上选择不同的和，用等距节点做次拉格朗日插值，观察误差大小和收敛情况。 程序： x0为n维插值节点向量，y0为n维被插函数值向量 x为插值点（也可以是m维向量） y为插值点x处的函数值（近似值，也可以是与x同维的向量） function y=lagrangen(x0,y0,x) n=length(x0); %m=length(x); %for i=1:m z=x;%(i); s=0; for k=1:n L=1; for j=1:n if j~=k L=L.*(z-x0(j))/(x0(k)-x0(j)); end end s=s+L*y0(k); end y=s;%y(i)=s; %end 程序运行 k=7; x1=0:0.1:1; y1=abs(sin(k.*pi.*x1)); xx1=0.354; yy1=lagrangen(x1,y1,xx1); x2=0:0.01:1; y2=abs(sin(k.*pi.*x2)); xx2=0.635; yy2=lagrangen(x2,y2,xx2); plot(x1,y1,g,xx1,yy1,o,x2,y2,r,xx2,yy2,.) grid 图形输出为：