昆明理工大学数值分析上机报告2.docxVIP

第 第 PAGE 10 页 函数插值方法(课题五) 一、问题提出 对于给定的一元函数 y ? f j ? ?0,1, , n?。试用 Lagrange ?x?  的 n+1 个节点值 y j ? f ?x ? j Lagrange 插值多项式。数据如下： （1） x j y j 求五次 Lagrange 多项式 L f ?0.596?, f ?0.99? 的值。 （2） 公式求其插值多项式或分段二次 0.40.550.65 0.4 0.55 0.65 0.80 0.95 0.4107 0.5781 0.6967 0.90 1.00 5 5 5 1.2538 2 ?x ?， 和分段三次插值多项式， 计算 5 123 1 2 3 4 5 6 7 j 0.368 0.135 0.050 0.018 0.007 0.002 0.001 j y 试构造 Lagrange 多项式 L 结果 f ?1.8?? 0.165299 f ?6 ?x ?，计算 f ?1.8?的值。 ?? 0二、要求 1、利用 Lagrange 插值公式 6.15 ? ? ??L ?x?? ?n ??n x ? xi ? y ? ? 编写出插值多项式程序； n?k ?0? i?0 n ? ? i?k x ? x k k i ? 2、给出插值多项式或分段三次插值多项式的表达式； 3、根据节点选取原则，对问题（2）用三点插值或二点插值，其结果如何； 4、对此插值问题用 Newton 插值多项式其结果如何。 2．作业环境(包括选用的程序语言、运行环境) 本题中的插值多项式程序采用的编程语言为 c++，因此运行环境可以在装有 Microsoft VC++的 windows XP 或 2000 的系统下运行程序。3．数学（理论背景）描述 在生产实践和科学研究所遇到的大量函数中，相当一部分是通过测量或实 验得到的。虽然其函数关系 y=f(x)在某个区间[a，b]上是客观存在的，但是却不知道具体的解析表达式，只能通过观察、测量或实验得到函数在区间［a，b］ 上一些离散点上的函数值、导数值等， 因此，希望对这样的函数用一个比较简单的函数表达式来近似地给出整体上的描述。还有些函数，虽然有明确的解析表达式，但却过于复杂而不便于进行理论分析和数值计算，同样希望构造一个既能反映函数的特性又便于计算的简单函数，近似代替原来的函数。插值法就是寻求近似函数的方法之一。 在用插值法寻求近似函数的过程中，根据所讨论问题的特点，对简单函数的类型可有不同的 选取，如多项式、有理式、三角函数等，其中多项式结构简单，并有良好的性质，便于数值计算和理论分析，因此被广泛采用。 1 n设函数 y=f(x)在区间［a,b］上有定义, y0 , y ,...y 且已知函数在区间［a,b 1 n 1上 n+1 个互异点 x0 , x 1 ,...x n 上的函数值，若存在一个简单函数 y=p(x )，使其 1n n1经过 y=f(x)上的这 n+1 个已知点( x0 , y0 ),( x , y ),…， (x , y )，即 1 n n 1 i ip( x )= y , i i 1那么，函数 p(x)称为插值函数，点 x0 , x 1  ,...x n  称为插节点， 点 1n n1( x0 , y0 ),( x , y ),…，( x , y ) 称为插值点，包含插值节点的区间［a,b 1 n n 1 插值区间,求 p (x)的方法称为插值法，f(x)称为被插函数。若 p(x)是次数不超过n 的多项式,用 Pn(x)表示，即 p (x) ? a  n 0 a x ? a x 2 1 2 ? ... ? a xn n 则称 pn (x) 为 n 次插值多项式，相应的插值法称为多项式插值；若P(x)为分段多项式，称为分段插值，多项式插值和分段插值称为代数插值。 数值计算公式 Lagrange 插值公式： ? ? ? n ??n ? x ? x ? L x ? ? ?  i ? y ?n ? 分段三次插值公式 k ?0? i?0 ? i?k x ? x k k i ? I ( x)? ?n h i?0 [ ?f i (x?) ?f i i i (x ) ] , 其中 I h (x ) ? f i i , I h (x ) ? f i i (i ? 0,1,^ n) ；? i (x), ? i (x) 如下： ?? x ? x ?2 ? x ? x ? ????i?1 ? ?1? 2 i ? , x ? x ? x (i ? 0略去） ?? x ? x ? ? x ? x