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

《数值分析》课程实验报告 用拉格朗日和牛顿插值法求解函数值 算法名称 用拉格朗日和牛顿插值法求函数值 学科专业 xxxxx 作者姓名 xxxx 作者学号 xxxxx 作者班级 xxxxxx 大学 二〇一五年十二月 《数值分析》课程实验报告 实验名称  用拉格朗日和牛顿插值法求解函数值  成绩 一、问题背景 在工程技术与科学研究中， 常遇到考察两个变量间的相互关系问题。 两个变量间的关系可以 通过函数表示，若 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 的一些数值如表： x lnx 用拉格朗日插值法计算的近似值。 2.函数 f(x)= x 的一些数值如表： x 用牛顿插值法计算√的近似值，画出插值函数与原函数的图形做比较。 三、算法描述 1.拉格朗日插值法：设已知 x0 ， x ， x 2 ， ..., xn 及 yi =f( xi Ln ( x) 为不超过 n 次 1 )(i=0,1,.....,n), 多项式且满足 Ln ( xi ) yi （ i=0,1,...n） .易知 Ln ( x) l 0 （ x） y 0 l n ( x) y n . = +....+ 其中， li (x) 均为 n 次多项式且满足式 （ 3）（ i,j=0,1,...,n ）,再由 x j l i ( x) （ j i）为 n 次多项式 n x xj 的 n 个根知 l i (x) j 0 =c i i .最后，由 n n 1 li (x j ) c ( xi x j ) 1 ( xi x j ) j 0 j 0 j i c= j i ,i=0,1,...,n. n n x xj . 0 xi xj li ( x)yi 总之， Ln (x) ， l i ( x) j l i ( x) = i 0 = j i 式为 n 阶 Lagrange 插值公式，其中， 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) syms t; n=length(x); ny=length(y); if n~=ny error end L=; for k=1:n lk=1; for j=1:n if j~=k lk=lk*(t-x(j))/(x(k)-x(j)); end end; L=L+y(k)*lk; end simplify(L); L=collect(L); yt=subs(L,'t' ,xt); 2.牛顿插值 建立函数文件： function [yt,N]=NewtInterp(x,y,xt) syms t; n=length(x); ny=length(y); if n~=ny error end a=zeros(1,n); N=y(1); w=1; for k=1:n-1 yy=zeros(1,n); for j=k+1:n yy(j)=(y(j)-y(k))/(x(j)-x(k)); end a(k)=yy(k+1); w=w*(t-x(k)); N=N+a(k)*w; y=yy; end yt=subs(N,'t' ,xt); simplify(N); N=collect(N); N=vpa(N,6); 五、实验结果及分析 1.拉格朗日插值法 在命令窗口输入： x=[,,,,]; y=[,,,,]; xt=;