练习8 插值拟合.docVIP

练习28 插 值 拟 合 数学知识背景 插值在机械加工等工程技术合数据处理等学科研究中有很多直接的运用，而且插值又是数值微分和常微分方程数值解计算的基础。主要的插值方法有拉各朗日、分段线性、三次样条三种插值对于一些数据进行拟合可以有很多方法，通常用最小二乘法对直线进行拟合，对于离散的数据，可以进行曲线拟合。 主要内容 【本练习考查知识点】 本练习介绍MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析。通过实例学习如何用插值方法与拟合方法解决实际问题。 练习过程 (1)函数y=exp(-,在n个节点上（n不要太大，如5～10）用拉格朗日、分段线性、三次样条三种插值方法，计算m个插值点的函数值（m要适中,如50～100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n，再作比较，由此作初步分析。 编制计算拉格朗日插值的名为lagr1的M文件： function 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 在命令区中运行如下程序： n=5;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=lagr1(x0,y0,x); y2=interp1(x0,y0,x); y3=interp1(x0,y0,x,spline); [xyy1y2y3] plot(x,z,r,x,y,k:,x,y1,r,x,y2,m,x,y3,b) gtext(Lagr.),gtext(Piece.-linear),gtext(y=exp(-x^2)) 运行后，所得程序结果如图28－1所是，得到各点节点和插值点的值如下：（因篇副问题，此处没有全部列出） ans = -2.0000 0.0183 NaN 0.0183 0.0183 -1.9333 0.0238 NaN 0.0416 -0.0137 -1.8667 0.0307 NaN 0.0649 -0.0348 -1.8000 0.0392 NaN 0.0882 -0.0457 -1.7333 0.0496 NaN 0.1115 -0.0472 -1.6667 0.0622 NaN 0.1348 -0.0398 -1.6000 0.0773 NaN 0.1581 -0.0243 -1.5333 0.0953 NaN 0.1814 -0.0013 -1.4667 0.1164 NaN 0.2047 0.0284 -1.4000 0.1409 NaN 0.2281 0.0642 -1.3333 0.1690 NaN 0.2514 0.1054 -1.2667 0.2010 NaN 0.2747 0.1512 -1.2000 0.2369 NaN 0.2980 0.2011 -1.1333 0.2768 NaN 0.3213 0.2543 -1.0667 0.3205 NaN 0.3446 0.3101 -1.0000 0.3679 NaN 0.3679 0.3679 -0.9333 0.4185 NaN 0.4100 0.4269 -0.8667 0.4718 NaN 0.4522 0.4865 -0.8000 0.5273 NaN 0.4943 0.5459 -0