清華大学_计算方法(数学实验)实验2插值与拟合.docVIP

实 验 2 插 值 与 拟 合 系 班 姓名 学号 【实验目的】 掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析。 掌握用MATLAB作线性最小二乘的方法。 通过实例学习如何用插值方法与拟合方法解决实际问题，注意二者的联系和区别。 【实验内容】 预备：编制计算拉格朗日插值的M文件： 以下是拉格朗日插值的名为y_lagrl的M文件： function y=y_lagr1(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=p*y0(k)+s; end y(i)=s; end 第1题（d） 选择函数y=exp(-x2) (-2≤x≤2),在n个节点上（n不要太大，如5~11）用拉格朗日、分段线性、三次样条三种插值方法，计算m个插值点的函数值（m要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n，在作比较，由此作初步分析。 运行如下程序： n=7;m=61;x=-2:4/(m-1):2; y=exp(-x.^2); z=0*x; x0=-2:4/(n-1):2; y0=exp(-x0.^2); y1=y_lagr1(x0,y0,x); y2=interp1(x0,y0,x); y3=interp1(x0,y0,x,spline); [xyy1y2y3] plot(x,z,w,x,y,r--,x,y1,b:,x,y2,m,x,y3,b) gtext(y=exp(-x^2)),gtext(Lagr.),gtext(Piece.-linear.),gtext(Spline), 运行后，得到各节点和插值点的值如下： X y y1 y2 y3 0 0.0667 0.1333 0.2000 … … 0.8667 0.9333 1.0000 … … 1.6000 1.6667 … … 1.9333 2.0000 1.0000 1.0000 1.0000 1.0000 0.9956 0.9958 0.9641 0.9947 0.9824 0.9831 0.9282 0.9797 0.9608 0.9624 0.8924 0.9561 … … 0.4718 0.4626 0.4995 0.4836 0.4185 0.4062 0.4523 0.4329 0.3679 0.3531 0.4051 0.3837 … … 0.0773 0.1292 0.1087 0.0534 0.0622 0.1271 0.0937 0.0347 … … 0.0238 0.0685 0.0334 0.0104 0.0183 0.0183 0.0183 0.0183 将三种插值结果y1,y2,y3与精确值y项比较，显然y1在节点处不光滑,拉格朗日插值出现较大的振荡,样条插值得结果是最好的.增加n值（使n=11），再运行以上程序，得到的图形如右图所示，比较这两个图可发现，节点增加后，三种插值方法结果的准确度均有所提高，因此可近似地认为：增加节点个数可以提高插值结果的准确程度。 第3题 用给定的多项式，如y=x3-6x2+5x-3，产生一组数据（xi,yi,i=1,2,…,n），再在yi上添加随机干扰（可用rand产生（0，1）均匀分布随机数，或用randn产生N（0，1）分布随机数），然后用xi和添加了随机干扰的yi作3次多项式拟合，与原系数比较。如果作2或4次多项式拟合，结果如何？ 解：2 编制y_2_3.m文件 n=15; x=0:8/(n-1):8; y=x.^3-6*x.^2+5*x-3; z=0*x; y0=y+rand(1,15); f=polyfit(x,y0,m); r=polyval(f,x) pl2ot(x,z,k,x,y,r:,x, r,b) 程序及运行结果如下： m=2 ,y_2_3 f =