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

1、数值分析课程实验报告用拉格朗日和牛顿插值法求解函数值算法名称用拉格朗日和牛顿插值法求函数值学科专业xxxxx作者姓名xxxx作者学号xxxxx作者班级xxxxxxxxx 大学二一五年十二月数值分析课程实验报告实验名称用拉格朗日和牛顿插值法求解函数值成绩一、问题背景在工程技术与科学研究中，常遇到考察两个变量间的相互关系问题。两个变量间的关系可以通过函数表示，若 x为自变量，y为因变量，则函数关系可描述为y=f(x)。大多数问题中，函数表达式y=f(x)未知，人们通常采用逼近的方法处理：取得一组数据点(xi,yi ) (i=0,1,2,n ),数据点可由不同方式取得(例如，可根据工程设计要求得到，

2、也可通过采样或实验取得)，然后构造一个简单函数 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的一些数值如表：x1.41.51.61.71.8Inx0.33650.40550.47000.53060.5878用拉格朗日插值法计算In 1.54的近似值。2函数f(x)=, x的一些数值如表：x0.40.50.60.70.80.63250.70710.77460.83670.8944用牛顿插值法计算0.54的近似值，画出插值函数与原函数的

3、图形做比较。三、算法描述1.拉格朗日插值法：设已知X。,X1,X2 ,., Xn 及yi=f( x )(i=0,i,.,n),Ln(x)为不超过n次多项式且满足Ln(X)* (i=0,1,.n).易知 Ln(x) =(x)y+.+ln(x) yn.其中,li(x)均为n次多项式且满足式(3)(i,j=0,1,.,n),再由Xj(ji )为n次多nX xj项式li(x)的n个根知li(x) =c ij i0.最后，1由nli(Xj) c(Xij 0j iXj) 1n(XXjj 0c= j i),i=0,1,., n.n X Xjn总之，Ln(X)=ih(x)yij 0 X Xj 0, i ( )

4、 = j i式为n阶Lagrange插值公式，其中，i()(i=0,1,.n)称为n阶Lagrange插值的基函数。2.牛顿插值法：插值法是利用函数f(X)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f(X)的近似值。如果这特定函数是多项式，就称它为插值多项式。当插值节点增减时全部插值基函数均要随之变化，这在实际计算中很不方便。为了克服这一缺点，提出了牛顿插值。牛顿插值通过求各阶差商，递推得 到的一个公式：f(x)=fx0+fx0,x1(x-x0)+fx0,x1,x2(x-x0)(x-x1)+.fx0,.x n(x-xO).(x-x

5、 n-1)+Rn(x)。四、主要代码1.拉格朗日插值建立函数文件：fun cti onyt,L=Lagl nterpl(x,y,xt)syms t ;n=len gth(x);ny=le ngth(y);if n=nyerrorendL=0.0;for k=1:nlk=1;for j=1: nif j=k lk=lk*(t-x(j)/(x(k)-x(j);enden d;L=L+y(k)*lk;endsimplify(L);L=collect(L);yt=subs(L, t ,xt);2.牛顿插值建立函数文件：fun cti onyt,N=Newt In terp(x,y,xt)syms t

6、;n=len gth(x);ny=le ngth(y);if n=nyerrorenda=zeros(1, n);N=y(1);w=1;for k=1:n-1yy=zeros(1, n);for j=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=1.4,1.5,1.6,1.7,1.8;y=0.3365,0.4055,

7、0.4700,0.5306,0.5878; xt=1.54;yt,L=Lagl nterpl(x,y,xt);z=1:0.05:4;yz=subs(L,t,z);figure;plot(z,log(z),-r,z,yz,-b)hold onplot(x,y,marker,+)hold onplot(xt,yt,marker,o),(x_k,y_k),x=1.54)lege nd(l n(x),拉格朗日插值多项式xlabel(x)ylabel(y)yt得到结果及图像如下：yt =0.4318得到In 1.54的近似值为0.4318。拉格朗日插值模型简单，结构紧凑，是经典的插值法。但是由于拉格朗日

8、的插值多项式和每个节点都有关，当改变节点个数时，需要重新计算。且当增大插值阶数时容易出现龙格现象。2.牛顿插值法在命令窗口输入：x=0.4 0.5 0.6 0.7 0.8;y=0.6325 0.7071 0.7746 0.8367 0.8944;xt=0.54;yt,N=Newtl nterp(x,y,xt)z=0.1:0.05:2;yz=subs(N,t,z);figure;plot(z,sqrt(z),-r,z,yz,-b)hold onplot(x,y,marker,+)hold onplot(xt,yt,marker,o)h=legend($sqrtx$,牛顿,$(x_k,y_k)$,

9、$x=0.54$);set(h,I nterpreter,latex)xlabel(x)ylabel(y)得到结果及图像如下：yt =0.7348N =-0.291667*tA4 + 0.925*tA3 - 1.30208*tA2 + 1.46125*t + 0.2046;边Figure 1 |Fil? dir 里 i户 w Insert IoclsHelp |O 3 U | o x- a|BQ_得到0.54的近似值为0.7348，插值函数为N =- 0.291667*tA4 + 0.925*tA3 - 1.30208*tA2 + 1.46125*t + 0.2046,其计算精度是相当高的。Lagrange插值法和Newton插值法解决实际问题中关于只提供复杂的离散数据的函数求值问