数值分析-插值法(共13页).doc

精选优质文档-倾情为你奉上第二章 插值法2.在区间-1,1上分别取n=10,20用两组等距节点对龙哥函数 f(x)=1/(1+25*x2)做多项式插值及三次样条插值，对每个n值，分别画出插值函数及f（x）的图形。（1）多项式插值先建立一个多项式插值的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); endendC=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;