拉格朗日插值1.docVIP

拉格朗日插值 　　2013-2014(1)专业课程实践论文 题目 :拉格朗日插值法 　　队员 1:李 赫, 0918180206, R 数学 09-2 　　队员 2:杨泓絮, 0918180225, R 数学 09-2 　　- 　　专业课程实践论文成绩考核表 　　指导教师签字: 　　一、算法理论 　　(一)问题的背景 　　在实际问题中常用到这样的函数 () y f x =,其在某个区间 [, ]a b 上是存在 的。但是,通过观察或测试或试验只能得到在区间 [, ]a b 上有限个离散点 　　01, , , n x x x … 上的函数值 ()i i y f x =()0,1, , i n =…… 。或者 ()f x 的函数 ()f x 表 　　达式是已知的,但却很复杂而不便于运算;希望用一个既能反映函数 ()f x 的 特性,又便于计算的简单函数来描述它。 (二)插值问题的数学提法 　　已 知 函 数 在 1n +个 点 01, , , n x x x … 上 的 函 数 值 ()i i y f x =, 　　()0,1, , i n =…… 。 求 一 个 简 单 函 数 ()y P x =, 使 其 满 足 :()i i y x P =, 　　()0,1, , i n =…… 。即要求该简单函数的曲线要经过 ()y f x =上已知的这个 1 　　n +个 点 :()()()n n y x y x y x , 1, 10, 0, , , ?, 同 时 在 其 它 [, ]x a b ∈上 要 估 计 误 差 : 　　()()()R x f x P x =- 　　其中 　　()P x 为 ()f x 的插值函数, 01, , , n x x x … 称为插值节 点,包含插值节点的区间 [, ]a b 称为插值区间,求插值函数 ()P x 的方法称为插 值法。若 ()P x 是次数不超过 n 的代数多项式,就称 ()P x 为插值多项式,相应 的插值法称为多项式插值。若 ()P x 是分段的多项式,就是分段插值。若 ()P x 是三角多项式,就称三角插值。 　　下面是对拉格朗日插值法的介绍: 　　已知函数 ()y f x =在 1n +个不同的点 01, , , n x x x … 上的函数值分别为: 　　01, , , n y y y …… ,求一个次数不超过 n 的多项式 ()n i i P x y =, ()0,1, , i n =…… ,即 　　1n +个不同的点可以唯一决定一个 n 次多项式。 　　(1)插值基函数 　　过 1n +个不同的点分别决定 1n +个 n 次插值基函数 ()()()01, , , n l x l x l x …… 　　每个插值基本多项式 ()i l x 满足: 1. ()i l x 是 n 次多项式; 　　2. ()1i i l x =,而在其它 n 个 ()0i k l x = , ()k i ≠。由于 ()0i k l x =, 　　()k i ≠,故有因子: 　　　　()()()()n i i x x x x x x x x -?--?-+-110 　　因其已经是 n 次多项式,故而仅相差一个常数因子。令: 　　()()()()() 　　n i i i x x x x x x x x a x l -?--?-=+-110 　　由于 ()1i i l x =,可以定出 a ,进而得到: 　　　　()()()()()n i i i i i i n i i i x x x x x x x x x x x x x x x x x I -?--?--?--?-= 　　+-+-110110 　　(2) n 次拉格朗日型插值多项式 ()n P x 　　()n P x 是 1n +个 n 次插值基本多项式 ()()()01, , , n l x l x l x …… 的线性 组合,相应的组合系数是 01, , , n y y y …… 。即: 　　()()()() 　　x l y x l y x l y x P n n n +?++=1100, 　　从而 ()n P x 是一个次数不超过 n 的多项式,且满足 　　()n i i P x y = ()0, 1 　　, , i n =…… 　　　　源程序: 　　#include 　　#include 　　#include 　　float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值算法 */ { 　　int i,j; 　　float *a,yy=0.0; /*a作为临时变量,记录拉格朗日插值多项式 */ a=(float *)mall