拉格朗日插值法在工程-应用中的算法实现.docVIP

2010年第 1期 总第 154期 林区教学 Teach ing of Fo re stry Region  No. 1 2010 General No. 154 拉格朗日插值法在工程应用中的算法实现 徐 小 丽 (苏州工业职业技术学院 ,江苏 苏州 ) 　　摘 　要 :“学以致用 ”是每一门学科都致力追求的境界 ,数学自然也不例外 。可人们往往淹没于一堆高深的数 学理论及烦琐的分析 、论证之中 ,早已失去了数学的乐趣 ! 特别是用数学来解决实际的问题更是感到无所适从 。介 绍拉格朗日插值法在现实分析中的应用 ,并通过 C语言程序来实现这一数学分析法的自动化 ,为复杂的工程分析 研究提供了一条数学算法的捷径 。 关键词 :插值法 ;拉格朗日多项式 ;算法实现 ; C语言 中图分类号 : O29　　　文献标志码 : A 　　　文章编号 : 1008 - 6714 ( 2010) 01 - 0017 - 03 y1- y0 　　一 、插值算法原理 p ( x) = y0+ x1- x0 ( x - x0) ;由前式整理得 : 在科学研究和实际工程设计中 ,几乎所有的问题都可 以用函数 y = f ( x)来表示其某种内在规律的数量关系 。但  p ( x) = x - x1 x0- x1  y 0 + x - x0 x1- x0 y1; 记 : l0( x ) = x-x1 x- x  , l 1 ( x ) 理想化的函数关系在实际的工程应用中是很难奢求的 。 = x - x0 。由一次函数的性质有 : 对于那些没有明显解析式的函数关系表达式则只能通过 实验或观测得到 y = f ( x)在某区间 [ a, b ]上的一系列点的 函数值 yi=f ( xi ) , ( i = 0, 1, 2, …, n) , 即一张表 ( x , yi ) , ( i = 0, 1, 2, …, n) 。根据给定的函数表构造一个既能反映函 数 f ( x)的特性 , 又便于计算的简单函数 p ( x)来近似 f ( x) , 并使得 p ( xi) = f ( xi) , ( i = 0, 1, 2, …, n) ,这就是插值法 。 设函数 y = f ( x) 在区间 [ a, b ]上有定义 , x0, x1, …, xn 是 [ a, b ]上取定的 n + 1个互异节点 ,且在这些点处的函数 值 f ( x0) , f ( x1) , …, f ( xn)为已知 ,即 yi= f ( x ) , ( i = 0, 1, 2, …, n) 。若存在一个和 f ( x)近似的函数 φ( x) ,满足 x1- x0 l ( x0) = l ( x1) = 1, l0( x1) = l ( x0) = 0; l0( x) + l ( x) = 1 其中 l ( x) , l1( x)称为线性插值基函数 。 当 n = 2时即为抛物插值 , 这也是一种常用的代数插 值 。根据给定函数 f ( x ) 在三个互异节点 x0, x1, x2的值 f ( x0) , f ( x1) , f ( x2) ,构造次数不超过二次的多项式 P ( x) = ax+ ax + a, 使满足二次插值条件 P ( x) = f ( x ) , ( i = 0, 1, 2) 。 P ( x)的参数直接由插值条件决定 , 并满足下面 的代数方程组 : 2 φ( xi) = f ( xi ) , ( i = 0, 1, 2, …, n) ( 1 - 1) a0+ a1x0+ a2x0= y0 则称 φ( x)为 f ( x ) 的一个插值函数 , 点 xi 为插值节点 , 式 ( 1 - 1)称为插值条件 , 区间 [ a, b ]称为插值区间 , 而误差 函数 R ( x) = f ( x) - φ( x)称为插值余项 。 由于代数多项式具有数值计算和理论分析方便的特 点 ,故在工程实际应用中我们多选代数插值 。即求一个次 数不超过 n次的多项式 P ( x) = anxn+ an- 1xn - 1+… + a1x + a0( 1 - 2) 满足 p ( x) = f ( x) , ( i = 0, 1, 2, …, n) , 则称 P ( x )为 f ( x) a0+ a1x1+ a2x21=y1 a0+ a1x2+ a2x22=y2 仿线性插 值 , 用基 函 数 的 方法 求 解 方 程 组 。求 二 次 式 l ( x) = 1, l( x) = l ( x) = 0。因 x, x是 l( x)的两个零 点 ,于是 l ( x) = c ( x - x1) ( x - x2) ,且 l ( x0) = 1, 确定系 1 数