数值分析实验题(华科)汇总.doc

数值分析实验题(华科)汇总 数值分析实验题(华科)汇总 PAGE / NUMPAGES 数值分析实验题(华科)汇总 数值分析实验作业 专业： 姓名： 学号： 实验 2.1 多项式插值的振荡现象 [问题提出 ]： 考虑在一个固定的区间上用插值逼近一个函数，显然 Lagrange 插值中使用的节点越多，插值多项式的 次数就越高，我们自然关心插值多项的次数增加时， Ln(x) 是否也更加靠近逼近的函数， Runge 给出的例子 是极著名并富有启发性的，设区间 [-1,1] 上函数 f (x) 1 1 25 x2 [实验内容 ]： 考虑区间 [-1 ， 1]的一个等距离划分，分点为 xi 1 2i , i 0,1,2,..., n n 则拉格朗日插值多项式为 n 1 Ln ( x) i 0 1 25x2 li (x) i 其中， li (x) ， i=0,1,2, ,n 是 n 次 Lagrange 插值函数。 [实验要求 ]： （ 1）选择不断增大的分点数目 n=2 ，3， 画出原函数 f(x) 及插值多项式函数 Ln(x) 在 [-1， 1]上的图像，比较并分析实验结果。 （ 2）选择其他的函数，例如定义在区间 [-5， 5]上的函数， h( x) x , g( x) arctan x 1 x4 重复上述的实验看其结果如何。 解：以下的 f(x) 、h(x) 、 g(x) 的为插值点用“ * ”表示，朗格朗日拟合曲线用连续曲线表示。通过三个函数的 拉格朗日拟合可以看到，随着插值点的增加，产生 Rung 现象。 （1） f(x) 多项式求值的振荡现象 n=2 多项式求值的振荡现象 n=3 1 1 0.9 f(x) 0.9 f(x) lagrange(x) lagrange(x) 0.8 0.8 0.7 0.7 0.6 0.6 y 0.5 y 0.5 0.4 0.4 0.3 0.3 0.2 0.2 0.1 0.1 0 0.40.60.81 0 0.40.60.81 -1-0.8-0.6-0.4 -0.200.2 -1-0.8 -0.6-0.4-0.200.2 x x 第1页 共13页 1 多项式求值的振荡现象 n=4 1.2 多项式求值的振荡现象 n=5 f(x) f(x) 0.8 lagrange(x) 1 lagrange(x) 0.6 0.8 0.4 0.6 y y 0.2 0.4 0 0.2 -0.2 0 -0.4 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -0.2 -0.8 -0.6 -0.4-0.2 0 0.2 0.4 0.6 0.8 1 -1 -1 x x 1.2 多项式求值的振荡现象 n=6 1 多项式求值的振荡现象 n=7 f(x) 0.9 f(x) 1 lagrange(x) lagrange(x) 0.8 0.8 0.7 0.6 0.6 y y 0.5 0.4 0.4 0.2 0.3 0.2 0 0.1 -0.2 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 0 -0.8 -0.6 -0.4-0.2 0 0.2 0.4 0.6 0.8 1 -1 -1 x x 1 多项式求值的振荡现象 n=8 1 多项式求值的振荡现象 n=9 0.8 f(x) f(x) lagrange(x) 0.8 lagrange(x) 0.6 0.4 0.6 0.2 0.4 y 0 y -0.2 0.2 -0.4 0 -0.6 -0.8 -0.2 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -0.4 -0.8 -0.6 -0.4-0.2 0 0.2 0.4 0.6 0.8 1 -1 -1 x x 1.6 多项式求值的振荡现象 n=10 1 多项式求值的振荡现象 n=11 1.4 f(x) 0.9 f(x) lagrange(x) lagrange(x) 1.2 0.8 1 0.7 0.8 0.6 y 0.6 y 0.5 0.4 0.4 0.2 0.3 0 0.2 -0.2 0.1 -0.4 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 0 -0.8 -0.6 -0.4-0.2 0 0.2 0.4 0.6 0.8 1 -1 -1 x x 第2页 共13页 (2) h(x) 0.6 多项式求值的振荡现象 n=2 0.6 多项式求值的振荡现象 n=3 h(x) h(x) 0.4 lagrange(x) 0.4 lagrange(x) 0.2 0.2 0 0 y y -0.2 -0.2 -0.4 -0.4 -0.6 -0.6 -0.8 -4 -3 -2