计算方法-第4讲--拉格朗日插值(2012新).ppt

* * Hbld-jpk@163.com 河北联合大学 求任一插值点 处的插值 ? ? ? ? ? 或无封闭形式, 节点可视为由 产生， 表达式复杂, 或未知。 ? 已知 n+1个节点 其中 互不相同，不妨设 3. 插值问题的提法 ? ? ? ? ? 构造一个(相对简单的)函数 通过全部节点, 即 再用 计算插值，即 ? 5、求解插值问题的基本思路 二、插值法的一般定义 ？ 定理1 证明 设有n+1个互不相同的节点 则存在唯一的多项式： 使得 构造方程组 三、一般插值多项式的原理 令： 方程组的矩阵形式如下： 所以方程组（4）有唯一解。 证毕 此定理说明只要n+1个节点互异，满足上述插值条件的多项式是唯一存在的。 注意：如果不限制多项式的次数，插值多项式并不唯一。 程序设计 X={x0,x1,x2,x3}={10,11,12,13}; y={y0,y1,y2,y3}={2.3026,2.3979,2.4849,2.5649}; A=Transpose[Table[{x0^j,x1^j,x2^j,x3^j},{j,0,3}]]; MatrixForm[%]; AA=LinearSolve[A,y]//N X1={1,x,x^2,x^3}; X1.AA N[%/.x-11.75,10] 可产生向量、矩阵等 一般插值多项式的原理的应用程序 A={{0,-1},{1.5,4.25},{5.1,35.21}} g1=ListPlot[Table[A],Prolog-AbsolutePointSize[10]]; Interpolation[A,InterpolationOrder-2] g2=Plot[%[x],{x,0,5.1}]; Show[g1,g2] N[%%%[3.66],5] 绘制点图 点的绝对直径 插值、插入 程序设计