数学[实验]之数值积分.pptVIP

一、 三种插值方法 插值多项式 根据n+1个节点得到唯一的n次多项式; 拉格朗日插值多项式 * 实验3 插值与数值积分 插值 可以理解为，在若干已知的函数值之间插入计算一些未知的函数值. 三种插值方法： 三次样条插值。 拉格朗日多项式插值、 分段线性插值、 教材P44----P64 参考书P219---P227 根据一组已知数据(xk,yk)(节点),认为来自于g(x),寻找一个靠近g(x)的简单函数f(x),使f(xk)=yk,并用f(x)求插值点的函数值. 数值积分根据一组已知数据(xk,yk),认为来自于f(x),近似地求一个定积分 . 四种数值积分方法: 矩形公式、 梯形公式、 辛普森 (Simpson)公式、 高斯(Gauss)公式 系数矩阵X为向量[x0,x1,…,xn]生成的凡德蒙阵，方程组简记作 XA=Y (4) 其中A=(an , an-1,...,a1, a0) 为多项式系数构成的列向量。方程组 有唯一的解，即根据n+1个节点可以确定唯一的 n 次插值多项式。 设 Ln (x)是n 次多项式： 对于节点(xj, yj )，应有 即多项式的系数满足方程组 （3） （1） （2） 向量A可以这样求出： A=X \ Y 这样我们就可以得到插值多项式函数了。 1 拉格朗日插值 x0=[0 2 4 5 7 8 10]; y0=[-3 -1 2 5 7 3 4]; X=vander(x0); Y=y0; A=X\Y x=[0:0.1:10]; y=polyval(A,x); plot(x,y,r-,x0,y0,b*) P=[ 0.0014 -0.0302 0.2151 -0.5135 0.1303 1.5119 -3.0000] 所求插值多项式为：p=A’ （5） li(x)是 n 次多项式，满足 （6） 令 （7） 显然(7)是满足(2)的多项式。由方程组(4)解的唯一性，（7）式的Ln(x)与(1)式相同。(5),(7)称为拉格朗日插值多项式，用Ln(x)计算插值称为拉格朗日插值。 为求满足 (2)的多项式Ln(x) ，作 计算机编程容易! 用另一种方法来构造Ln(x) 为作拉格朗日插值，先作 M文件定义一个函数lagr1，保 存在搜索路径内。 function y=lagr1(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 输入节点x0,y0及插值x 后，运行 y=lagr1(x0,y0,x)， plot(x,y,’r-’,x0,y0,’b+’) 即可作出插值函数的图形 最里层循环对一个插值点x构造li(x)，外一层循环对一个插值点构造Ln(x)，最外层循环求出每个插值点 x 处的多项式值 y。 x0=[0 2 4 5 7 8 10]; y0=[-3 -1 2 5 7 3 4]; x=0:0.1:10; y=lagr1(x0,y0,x); plot(x,y,r-,x0,y0,b+) 还是用前面的例子： Lagrange插值多项式震荡的很厉害,看19世纪range给出的一个 不收敛的反例. 2 分段线性插值 将每两个相邻的节点用直线连起来，如此形成的一条折线 就是分段线性插值函数，记作In(x)，它满足In(xj)=yj，且In(x)在 每个小区间[xj, xj+1]上是线性函数( j=0,1,…,n)。 (12) 分段线性插值函数In(x)有良好的收敛性，对于x∈[a, b]有 。 作分段函数 (13) 满足In(xj)=yj 的分段线性插值函数就是： 分段线性插值由Matlab中的库函数interp1(x0,y0,x)实现 x0 x1 x2 x3 x4 x5 l0(x) 1 l1(x) l3(x) x y o l5(x) x0 x1 x2 x3 x4