数值分析插值法.docxVIP

第二章 插值法 2.在区间[-1,1]上分别取 n=10,20 用两组等距节点对龙哥函数 f(x)=1/(1+25*x^2)做多项式插值及三次样条插值，对每个n 值，分别画出插值函数及f（x）的图形。 多项式插值 ①先建立一个多项式插值的 M-file； 输入如下的命令（如牛顿插值公式）： function [C,D]=newpoly(X,Y) n=length(X); D=zeros(n,n) D(:,1)=Y for j=2:n for k=j:n D(k,j)=(D(k,j-1)- D(k-1,j-1))/(X(k)-X(k-j+1)); end end C=D(n,n); for k=(n-1):-1:1 C=conv(C,poly(X(k))) m=length(C); C(m)= C(m)+D(k,k); end ②当 n=10 时，我们在命令窗口中输入以下的命令： clear,clf,hold on; X=-1:0.2:1; Y=1./(1+25*X.^2); [C,D]=newpoly(X,Y); x=-1:0.01:1; y=polyval(C,x); plot(x,y,X,Y,.); grid on; xp=-1:0.2:1; z=1./(1+25*xp.^2); plot(xp,z,r) 得到插值函数和f（x）图形： ③当 n=20 时，我们在命令窗口中输入