科学工程计算与matlab编程6-2讲述.ppt

6.2 插值6.2.1 Lagrange插值 方法介绍 算例：给出f(x)=ln(x)的数值表， 用Lagrange计算ln(0.54)的近似值。 x=[0.4:0.1:0.8]； % 取样本点 y=log(x) ans = 0.916291,-0.693147,-0.510826,-0.356675,-0.223144 y0=lagrange(x,y,[0.54,0.55,0.78]) ans = -0.6161 -0.5978 -0.2484 y0(1)-log(0.54); %计算误差 ans = 4.352891839787265e-005 5.2.2 Hermite插值 方法介绍 不少实际问题不但要求在节点上函数值相等，而且要求一阶、二阶甚至更高阶导数值也相等，满足这一要求的插值多项式就是Hermite插值多项式。 只讨论函数值与一阶导数值均相等的Hermite插值。 MATLAB实现 % hermite.m function y=hermite(x0,y0,y1,x) n=length(x0); m=length(x); for k=1:m yy=0.0; for i=1:n h=1.0; a=0.0;