201321080108-张勇杰-数值分析.pdfVIP

分段插值思想解决龙格振荡现象 201321080108 张勇杰 一、问题叙述 1       设f x  2 5x5 ，取不同插值结点x k0,1,,10k ，计算对应函数 x 1    10 值y f x k0,1,,10 ，构造 次插值函数多项式L (x)。通过与被插值函 k k 10 数比较发现插值函数出现振荡现象，L (x)对曲线局部特性的逼近很不理想。故 10 而寻求分段插值方法解决龙格振荡现象。 二、问题分析 x 5k k0,1,,10  10 首先分析问题产生原因，取插值点 k ，构造 次插值 1   多项式L (x)，如图 所示。图中实线为被插值函数f x ，虚线为插值函数L (x)。 10 10 空心点为插值结点处的函数值。图形显示插值函数出现振荡现象（在5和-5附近 误差很大）。 图 1 整体拉格朗日插值 整体插值会造成插值函数与被插值函数局部误差过大，从而影响插值效果。 可以通过分段插值思想来保证局部特性的逼近。通过多种分段插值手段，对比分 析插值效果。 分段线性插值是在相邻结点构成的小区间上分别做线性插值。分段三次埃尔 米特插值是在相邻结点构成的小区间上做三次埃尔米特插值，需要提供每个插值 结点的对应的函数一阶倒数值。三次样条插值不需要知道每个结点对应的一阶导 数值，通过给定的两个边界条件，逼近各结点上的一阶导数值从而得到相对光滑 插值函数。 三、实验程序及注释 1. 编写拉格朗日插值函数，保存为lagrange.m文件： functiony=lagrange(x0,y0,x); n=length(x0);m=length(x); fori=1:m z=x(i); s=0.0; fork=1:n p=1.0; forj=1:n ifj~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 调用拉格朗日函数，进行拉格朗日插值： x=-5:5;y=1./(x.^2+1); t=-5:0.01:5;f=1./(t.^2+1); l=lagrange(x,y,t); plot(t,l,k--,t,f,k-,x,y,ko); 2. 编写线性插值函数，保存为mylinear.m functiony=mylinear(x0,y0,x1,y1,x); m=length(x); fori=1:m y(i)=(x(i)-x1)/(x0-x1)*y0+(x(i)-x0)/(x1-x0)*y1; end 调用线性插值函数，进行分段线性插值： x=-5:5;y=1./(x.^2