数值分析课程实验报告-拉格朗日和牛顿插值法

《数值分析》课程实验报告用拉格朗日和牛顿插值法求解函数值算法名称用拉格朗日和牛顿插值法求函数值学科专业xxxxx作者姓名xxxx作者学号xxxxx作者班级xxxxxxxxx大学二0一五年十二月《数值分析》课程实验报告实验名称 用拉格朗日和牛顿插值法求解函数值 成绩一、 问题背景在工程技术与科学研究中，常遇到考察两个变量间的相互关系问题。两个变量间的关系可以通过函数表示，若x为自变量，y为因变量，则函数关系可描述为y=f(x)。大多数问题中，函数表达式y=f(x)未知，人们通常采用逼近的方法处理：取得一组数据点(xi,yi)(i=0,1,2,...,n),数据点可由不同方式取得(例如，可根据工程设计要求得到，也可通过采样或实验取得)，然后构造一个简单函数P(x)作为y=f(x)的近似表达式，即y=f(x)^P(x),对于y=f(x)^P(x)，若满足P(xi)=f(xi)=yi,i=0,1,2...,n,这类问题成为插值问题。二、 数学模型1.函数f(x)=lnx的一些数值如表：xInx用拉格朗日插值法计算的近似值。2.函数f(x)=IX的一些数值如表：x用牛顿插值法计算”的近似值，画出插值函数与原函数的图形做比较。三、算法描述Y Vx L(x)1.拉格朗日插值法：设已知X0，x1，x2,...,七及七=f(七)(i=0,1,.....,n),/)为不超过n次L(x)=y L(x)I/、yl(x)y多项式且满足nii(i=0,1,.・.n),易知n=0(x) 0+...・+nn.l(x) x l(x)其中，卜'均为n次多项式且满足式(3)(i,j=0,1,...,n),再由j(产i)为n次多项式卜)的n的n个根知((x)=c么°叶x-xj.最后l(x)=c^n(x-x)=1nFf(x-x)ij ij ijj=° j=°j丰i c=j丰i ,i=0,1,...,n.X-X寸 n .z1(X)y x-x.弓士L(X) 1 1I(X)j-0ij I(X)总之，n =i=0 ,i=j& 式为n阶Lagrange插值公式，其中，i(i=0,1,...n)称为n阶Lagrange插值的基函数。2.牛顿插值法：插值法是利用函数f(x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f(x)的近似值。如果这特定函数是多项式，就称它为插值多项式。当插值节点增减时全部插值基函数均要随之变化，这在实际计算中很不方便。为了克服这一缺点，提出了牛顿插值。牛顿插值通过求各阶差商，递推得到的一个公式：f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x0)...(x-xn-1)+Rn(x)四、主要代码1.拉格朗日插值建立函数文件：function[yt,L]=LagInterpl(x,y,xt)symst;n=length(x);ny=length(y);ifn~=nyerrorendL=;fork=1:nlk=1;forj=1:nifj~=klk=lk*(t-x(j))/(x(k)-x(j));endend;L=L+y(k)*lk;endsimplify(L);L=collect(L);yt=subs(L,'t',xt);2.牛顿插值建立函数文件：function[yt,N]=NewtInterp(x,y,xt)symst;n=length(x);ny=length(y);ifn~=nyerrorenda=zeros(1,n);N=y(1);w=1;fork=1:n-1yy=zeros(1,n);forj=k+1:nyy(j)=(y(j)-y(k))/(x(j)-x(k));enda(k)=yy(k+1);w=w*(t-x(k));N=N+a(k)*w;y=yy;endyt=subs(N,'t',xt);simplify(N);N=collect(N);N=vpa(N,6);五、实验结果及分析1.拉格朗日插值法在命令窗口输入：x=[〃,,];y=[〃〃];xt=;[yt,L]=LagInterpl(x,y,xt);z=1::4;yz=subs(L,'t',z);figure;plot(z,log(z),'--r',z,yz,'-b')holdonplot(x,y,'marker','+')holdonplot(xt,yt,'marker','o')legend('ln(x)','拉格朗日插值多项式','(x_k,y_k)','x=')xlabel('x')ylabel('y')yt得到结果及图像如下：yt=得到的近似值为。拉格朗日插值模型简单，结构紧凑，是经典的插值法。但是由于拉格朗日的插值多项式和每个节点都有关，当改变节点个数时，需要重新计算。且当增大插值阶数时容易出现龙格现象。2.牛顿插值法在命令窗口输入：x=[ ]；y=[];xt=;[yt,N]=NewtInterp(x,y,xt)z=::2;yz=subs(N,'t',z);figure;plot(z,sqrt(z),'--r',z,yz,'-b')holdonplot(x,y,'marker','+')holdonplot(xt,yt,'marker','o')h=legend('$\sqrt{x}$','牛顿','$(x_k,y_k)$','$x=$');set(h,'Interpreter','latex')xlabel('x')ylabel('y')得到结果及图像如下：yt=n=-*tA4+*tA3-*tA2+*t+

得到V的近似值为，插值函数为N=_冲4+*tA3-*tA2+*t+，其计算精度是相当高的。Lagrange插值法和Newton插值法解决实际问题中关于只提供复杂的离散数据的函数求值问题，通过将所考察的