数值分析 插值法.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