实 验 2 插 值 与 拟 合 系 班 姓名 学号 【实验目的】 1、 掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法,改变节点的数目,对三种插值结果进行初步分析
2、 掌握用MATLAB作线性最小二乘的方法
3、 通过实例学习如何用插值方法与拟合方法解决实际问题,注意二者的联系和区别
【实验内容】 预备:编制计算拉格朗日插值的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'); -2-1
91y =ex p(-x2)Lagr
-linear
Spline012345678