MATLAB数值实验一(数据的插值运算及其应用完整版).docVIP

佛山科学技术学院 实 验 报 告 课程名称 数值分析 实验项目 插值法与数据拟合 专业班级 机械工程 姓 名 余红杰 学 号 2111505010 指导教师 陈剑 成 绩 日 期 月 日 一、实验目的 1、学会Lagrange 插值、牛顿插值和三次样条插值等基本插值方法； 2、讨论插值的Runge现象 3、学会Matlab提供的插值函数的使用方法，会用这些函数解决实际问题。 二、实验原理 1、拉格朗日插值多项式 2、牛顿插值多项式 3、三次样条插值 三、实验步骤 1、用MATLAB编写独立的拉格朗日插值多项式函数 2、用MATLAB编写独立的牛顿插值多项式函数 3、用MATLAB编写独立的三次样条函数（边界条件为第一、二种情形） 4、已知函数在下列各点的值为： 0.2 0.4 0.6 0.8 1.0 0.98 0.92 0.81 0.64 0.38 根据步骤1,2,3编好的程序，试分别用4次拉格朗日多项式、牛顿插值多项式以及三次样条函数(自然边界条件)对数据进行插值，并用图给出 ｛｝，、和。 5、在区间[-1,1]上分别取用两组等距节点对龙格函数作多项式插值，对不同值，分别画出插值函数及的图形。 6、下列数据点的插值 0 1 4 9 16 25 36 49 64 0 1 2 3 4 5 6 7 8 可以得到平方根函数的近似，在区间[0,64]上作图。 用这9个点作8次多项式插值。 用三次样条（第一边界条件）程序求。 7、对于给函数在区间[-1,1]上取，试求3次曲线拟合，试画出拟合曲线并打印出方程，与第5题的结果比较。 四、实验过程与结果： 1、Lagrange 插值多项式源代码： function ya=lag(x,y,xa) %x 所有已知插值点 %y 插值点对应函数值 %xa 所求点，自变量 %ya 所求点插值估计量 ya=0; mu=1; %初始化 %循环方式求L系数，并求和： for i = 1:length(y) for j = 1:length(x) if i ~= j mu = mu * (xa - x(j) ) / ( x(i) - x(j) ); else continue end end ya = ya + y(i) * mu ; mu = 1; end Newton源代码： function ya = newton(x,y,xa) %x 所有已知插值点 %y 插值点对应函数值 %xa 所求点，自变量 %ya 所求点插值估计量 %建立系数零矩阵D及初始化： D = zeros(length(x)-1); ya = y(1); xi = 1; %求出矩阵D，该矩阵第一行为牛顿插值多项式系数： for i=1:(length(x)-1) D(i,1) = (y(i+1) -y(i))/(x(i+1) -x(i)); end for j=2:(length(x)-1) for i=1:(length(x)-j) D(i,j) = (D(i+1,j-1) - D(i,j-1)) / (x(i+j) - x(i)); end end %xi为单个多项式（x-x(1))(x-x(2))...的值 for i=1:(length(x)-1) for j=1:i xi = xi*(xa - x(j)); end ya = ya + D(1,i)*xi; xi = 1; end 三次样条插值多项式 （1）（第一边界条件）源代码： function y=yt1(x0,y0,f_0,f_n,x) _____________(1) %第一类边界条件下三次样条插值； %xi 所求点； %yi 所求点函数值； %x 已知插值点； %y 已知插值点函数值； %f_0左端点一次导数值； %f_n右端点一次导数值； n = length(x0); z = length(y0); h = zeros(n-1,1); k=zeros(n-2,1); l=zeros(n-2,1); S=2*eye(n); for i=1:n-1 h(i)= x0