试验二插值法与数据拟合课案.doc

试验二 插值法 实验目的 学会Lagrange 插值和牛顿插值等基本插值方法； 讨论插值的Runge现象，掌握分段线性插值方法 学会Matlab提供的插值函数的使用方法，会用这些函数解决实际问题。 实验要求 按照题目要求完成实验内容； 写出相应的Matlab 程序； 给出实验结果(可以用表格展示实验结果)； 分析和讨论实验结果并提出可能的优化实验。 写出实验报告。 实验步骤 1、用编好的Lagrange 插值法程序计算书本P66 的例1、用牛顿插值法计算P77的例1。 2、已知函数在下列各点的值为： 0.2 0.4 0.6 0.8 1.0 0.98 0.92 0.81 0.64 0.38 试用4次牛顿插值多项式对数据进行插值，根据｛｝，画出图形。 3、在区间[-1,1]上分别取用两组等距节点对龙格函数作多项式插值，对不同值，分别画出插值函数及的图形。 4、下列数据点的插值 0 1 4 9 16 25 36 49 64 0 1 2 3 4 5 6 7 8 可以得到平方根函数的近似，在区间[0,64]上作图。 用这9个点作8次多项式插值。 附：试验报告格式样本（正式报告这行可删除） 佛山科学技术学院 实 验 报 告 课程名称 数值分析 实验项目 插值法 专业班级 12.数学与应用数学（师范） 姓名 叶楚欣 学号 2012214103 指导教师 黄国顺 成 绩 日 期 月 日 一、实验目的 1、学会Lagrange 插值、牛顿插值和 分段线性插值等基本插值方法； 2、讨论插值的Runge现象，掌握分段线性插值方法 3、学会Matlab提供的插值函数的使用方法，会用这些函数解决实际问题。 二、实验原理 1、拉格朗日插值多项式 2、牛顿插值多项式 3、分段线性插值 三、实验步骤 1、用MATLAB编写独立的拉格朗日插值多项式函数 2、用MATLAB编写独立的牛顿插值多项式函数 3、利用编写好的函数计算本章P66例1、P77例1的结果并比较。 4、已知函数在下列各点的值为： 0.2 0.4 0.6 0.8 1.0 0.98 0.92 0.81 0.64 0.38 试用4次牛顿插值多项式对数据进行插值，根据｛｝，画出图形。 5、在区间[-1,1]上分别取用两组等距节点对龙格函数作多项式插值，对不同值，分别画出插值函数及的图形。 6、下列数据点的插值 0 1 4 9 16 25 36 49 64 0 1 2 3 4 5 6 7 8 (1)可以得到平方根函数的近似，在区间[0,64]上作图。 用这9个点作8次多项式插值。 源程序与解答： 1、Lagrange 插值多项式源代码 function yi=Lagrange(x, y, xi) % Lagrange 插值多项式，其中 % x --- 向量，全部的插值节点 % y --- 向量，插值节点处的函数值 % xi --- 标量，自变量x % yi --- xi 处的函数估计值 n=length(x); m=length(y); if n~=m error(The lengths of X and Y must be equal); return; end p=zeros(1,n); % 对向量p赋初值0 for k=1:n t=ones(1,n); for j=1:n if j~=k if abs(x(k)-x(j))eps error(the DATA is error!); return; end t(j)=(xi-x(j))/(x(k)-x(j)); end end p(k)=prod(t); end yi=sum(y.*p); （2）function yi=New_Int(x, y, xi) % Newton 基本插值公式，其中 % x --- 向量，全部的插值节点，按行输入 % y --- 向量，插值节点处的函数值，按行输入 % xi --- 标量，自变量x % yi --- xi 处的函数估计值 n=length(x); m=leng