南大数值分析课件§2 复合求积.pptVIP

* §2 复合求积 /* Composite Quadrature */ Haven’t we had enough formulae? What’s up now? Oh come on, you don’t seriously consider h=(b?a)/2 acceptable, do you? Why can’t you simply refine the partition if you have to be so picky? Don’t you forget the oscillatory nature of high- degree polynomials! Uh-oh 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复合求积公式。 ? 复合梯形公式： 在每个 上用梯形公式： = Tn /*中值定理*/ §2 Composite Quadrature ? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2n 为偶数， 这时 ，有 §2 Composite Quadrature ? 收敛速度与误差估计： 定义 　　 若一个积分公式的误差满足 且C ? 0，则称该公式是 p 阶收敛的。 ~ ~ ~ 例：计算 解： 其中 = 3.138988494 其中 = 3.141592502 运算量基本相同 例： 用复化梯形公式、复化Simpson公式和复化 Cotes公式计算 要使精度要求达到 问N各取多少？ HW: p.273 1，3，4，5，6 §2 Composite Quadrature Q: 给定精度 ?，如何取 n ? 例如：要求 ，如何判断 n = ？ ? 上例中若要求 ，则 即：取 n = 409 通常采取将区间不断对分的方法，即取 n = 2k 上例中2k ? 409 ? k = 9 时，T512 = 3 S4 = 3.141592502 注意到区间再次对分时 可用来判断迭代 是否停止。 §2 Composite Quadrature Lab 13. Composite Trapezoidal Rule Use the Composite Trapezoidal rule with a given n 0 to approximate a given integral . You are supposed to write a function double CTR ( int n, double a, double b, double (*f)( ) ) to approximate the integral from a to b of the function f using the trapezoidal rule on n equal-length subintervals. Input There is no input file. Instead, you must hand in your function in a *.h file. The rule of naming the *.h file is the same as that of naming the *.c or *.cpp files. Output For each test case, you are supposed to return the approximation of the integral. §2 Composite Quadrature Sample Judge Program #include stdio.h #include math.h #include13.h ? double f1 ( double x ) { return (1.0/(1.0+sin(x)*sin(x))); } ? double f2 ( double x ) { return (x*log(x)); } ? void main( ) {