宁波大学-数值方法.pptVIP

应用： 讲解编程，加深对算法的理解，掌握算法编程 §1 Lagrange Polynomial 注：? 通常不能确定 ?x , 而是估计 , ?x?(a,b) 将 作为误差估计上限。 ?当 f(x) 为任一个次数? n 的多项式时， , 可知 ，即插值多项式对于次数? n 的多项式是精确的。 Quiz: 给定 xi = i +1, i = 0, 1, 2, 3, 4, 5. 下面哪个是 l2(x)的图像？ y 0 - - - 1 0.5 - 0.5 1 2 3 4 5 6 x y 0 - - - 1 0.5 - 0.5 1 2 3 4 5 6 x y 0 - - - 1 0.5 - 0.5 1 2 3 4 5 6 x A B C §1 Lagrange Polynomial 例：已知 分别利用 sin x 的1次、2次 Lagrange 插值计算 sin 50? 并估计误差。 解： n = 1 分别利用x0, x1 以及 x1, x2 计算 ?利用 这里 而 ?sin 50? = 0.7660444… ) 18 5 ( 50 sin 1 0 ? ? p L 0.77614 外推 /* extrapolation */ 的实际误差 ? ?0.01001 ?利用 sin 50? ? 0.76008, 内插 /* interpolation */ 的实际误差 ? 0.00596 内插通常优于外推。选择要计算的 x 所在的区间的端点，插值效果较好。 §1 Lagrange Polynomial n = 2 ) 18 5 ( 50 sin 2 0 ? ? p L 0.76543 ?sin 50? = 0.7660444… 2次插值的实际误差 ? 0.00061 高次插值通常优于低次插值 但绝对不是次数越高就越好，嘿嘿…… HW: p.112-113 #2，#4, #5 §1 Lagrange Polynomial Lab 10. Lagrange Polynomial Given a set of sample points of a certain function f , use Lagrange polynomial to approximate function values at some given points . Input There are several sets of inputs. For each set: The 1st line contains an integer 20 ? n ? 0 which is the degree of Lagrange polynomial. n = ?1 signals the end of file. The 2nd line contains n+1 distinct real numbers . The 3rd line contains n+1 real numbers . The last line of a test case consists of an integer m 0 and m real numbers . The numbers are separated by spaces and new lines. §1 Lagrange Polynomial Output (? represents a space) For each ai , you are supposed to print the function value at this point in the following format: fprintf(outfile, f(%6.3f)?=?%12.8e\n, a, f ); The outputs of two test cases must be seperated