实验二 插值法.doc

实验二 插值法 实验2.1 实验程序： function charpt2_1 %输入插值结点 result=inputdlg({请输入插值结点数N：},charpt_2,1,{2}) N1=str2num(char(result)); result=inputdlg({请输入区间下界：},charpt_2,1,{-1}) N2=str2num(char(result)); result=inputdlg({请输入区间上界：},charpt_2,1,{1}) N3=str2num(char(result)); if(N11) errordlg(结点输入错误！);return;end %插值结点小于1时错误 f=inline(1./(1+25*x.^2));a=N2;b=N3; %标准函数f x1=linspace(a,b,N1+1); y1=feval(f,x1); x=a:0.001:b;inter1=Lagrange(x1,y1,x); fplot(f,[a,b],r-); %画标准函数f的图形 hold on; plot(x,inter1,b-); %画插值逼近函数 legend(原函数,插值多项式函数) xlabel(x);ylabel(y=f(x) o and y=Ln(x)--); function y=Lagrange(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if(j~=k) p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=s+p*y0(k); end y(i)=s; end 实验结果： n=2 n=3 n=5 n=10 n=15 n=20 实验分析： 多项式插值逼近结果如图所示，当结点数由3个结点增加为5个结点时，函数图形越接近，当插值结点增加到10时，曲线光滑，函数逼近效果在曲线的中间部分（-0.6，0.6）比较好，但是由上图在两个端点处x=-1和x=1附近时会出现在与原函数f(x)偏差会很大（龙格现象）。可以看出，适当提高插值多项式次数，可以提高逼近的精度，但是太高反而会产生不良的现象。 实验程序： function charpt2_2 %输入插值结点 result=inputdlg({请输入插值结点数N：},charpt_2,1,{2}) N1=str2num(char(result)); result=inputdlg({请输入区间下界：},charpt_2,1,{-5}) N2=str2num(char(result)); result=inputdlg({请输入区间上界：},charpt_2,1,{5}) N3=str2num(char(result)); if(N11) errordlg(结点输入错误！);return;end %插值结点小于1时错误 f=inline(x./(1+x.^4));a=N2;b=N3; %标准函数f x1=linspace(a,b,N1+1); y1=feval(f,x1); x=a:0.001:b;inter1=Lagrange(x1,y1,x); fplot(f,[a,b],r-); %画标准函数f的图形 hold on; plot(x,inter1,b-); %画插值逼近函数 legend(原函数,插值多项式函数) xlabel(x);ylabel(y=f(x) o and y=Ln(x)--); function y=Lagrange(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if(j~=k) p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=s+p*y0(k); end y(i)=s