数值分析插值法.docx

第二章 插值法

2.在区间[-1,1上]分别取n=10,20用两组等距节点对龙哥函数f(x)=/1(1+25*x^2做)及三次样条插值，对每个n值，分别画出插值函数及f（x）的图形。

多项式插值

①先建立一个多项式插值的M-fil；e输入如下的命令（如牛顿插值公式）：function[C,D]=newpoly(X,Y)n=length(X);

D=zeros(n,n)D(:,1)=Y

forj=2:n

fork=j:n

D(k,j)=(D(k,j-1)-D(k-1,j-1))/(X(k)-X(k-j+1));

end

多项式插值

endC=D(n,n);

fork=(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,holdon;

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);gridon;

xp=-1:0.2:1;

z=1./(1+25*xp.^2);plot(xp,z,r)

得到插值函数和f（x）图形：

＇飞IIIII