安徽大学数值分析实验报告一.doc

数值分析实验报告一 姓 名 张骐 学 号 专业 金融数学 实验日期 2017年12月6日 实验性质 验证性实验 实验名称： 实验一 数据的插值与拟合分析 龙格现象的发生与防止，插值效果的比较 实验目的： 观察拉格朗日插值的龙格现象，探索避免此现象发生的方法，比较不同方法的插值结果。 实验内容： 将区间[-5,5]作10等分，计算函数在每个节点的值，做出插值函数的图形并与的图形比较。 实验要求三： 对函数做拉格朗日插值。在MATLAB中用内部函数plot利用插值点绘制函数的图形。 对函数做牛顿插值。在MATLAB中用内部函数plot利用插值点绘制函数的图形。 对函数做三次样条插值。在MATLAB中用内部函数plot利用插值点绘制函数的图形。 在MATLAB中用内部函数ezplot直接绘制函数的图形，并与以上方法做出的插值函数的图形进行比较（自编程序，用不同的颜色，不同节点符号将（1）--（3）的结果画在一张图上）。 程序： function [C,L]=lagran(x,y) （编写lagran插值函数文件） w=length(x); n=w-1; L=zeros(w,w); for k=1:n+1 V=1; for j=1:n+1 if k~=j V=conv(V,poly(x(j))/(x(k)-x(j))); end end L(k,:)=V; end C=y*L 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 function s=csfit(x,y,dx0,dxn) （编写三次样条插值函数文件） n=length(x)-1; h=diff(x); d=diff(y)./h; a=h(2:n-1); b=2*(h(1:n-1)+h(2:n)); c=h(2:n); u=6*diff(d); b(1)=b(1)-h(1)/2; u(1)=u(1)-3*(d(1)); b(n-1)=b(n-1)-h(n)/2; u(n-1)=u(n-1)-3*(-d(n)); for k=2:n-1 temp=a(k-1)/b(k-1); b(k)=b(k)-temp*c(k-1); u(k)=u(k)-temp*u(k-1); end m(n)=u(n-1)/b(n-1); for k=n-2:-1:1 m(k+1)=(u(k)-c(k)*m(k+2))/b(k); end m(1)=3*(d(1)-dx0)/h(1)-m(2)/2; m(n+1)=3*(dxn-d(n))/h(n)-m(n)/2; for k=0:n-1 s(k+1,1)=(m(k+2)-m(k+1))/(6*h(k+1)); s(k+1,2)=m(k+1)/2; s(k+1,3)=d(k+1)-h(k+1)*(2*m(k+1)+m(k+2))/6; s(k+1,4)=y(k+1); End x=-5:0.01:5; z=5./(1+x.^2); plot(x,z) x1=-5:1:5; y=5./(1+x1.^2); [C,L]=lagran(x1,y); xx=-5:0.1:5; yy=polyval(C,xx); hold on plot(xx,yy,x1,y) （描绘lagran函数插值图像以及插值点） [C,D]=newpoly(x1,y) x2=-5:0.01:5; y2=polyval(C,x2); plot(x2,y2,r:) （作牛顿插值图像） grid on x=-5:1:5; y=5./(1+x.^2); dx0=0.0739645; dxn=-0.0739645; S=csfit(x,y,dx0,dxn) x1=-5:0.01:-4;y1=polyval(S(1,:),x1-x(1)); x2=-4:0.01:-3;y2=polyval(S(2,:),x2-x(2)); x3=-3:0.01:-2;y3=polyval(S(3,:),x3-x(3)); x4=-2:0.01:-1;y4=po