[2018年最新整理]3数值积分数值微分.ppt

第三章 数值积分数值微分 3.1 引言 7.2 Newton-Cotes 求积公式 五. 变步长的求积法(区间逐次分半法) 例题 3.4 Romberg求积公式 Romberg算法 计算步骤 例题(P129) 3.5 Gauss积分公式 7.6 数值微分 结果比较 首先已经将区间 [a, b] 进行 n 等分，用复合的梯形公式计算积分： 再将区间 [a , b] 进行 2 n 等分，实际上是将每个小区间再 2 等分，如此有： 变步长的求积法是先对较大的步长进行计算，然后将步长逐步折半，直至达到精度要求。 于是有变步长的梯形公式： 用变步长的梯形公式计算： T1=0.5[f (0)+ f (1)] =0.9207355 T2=0.5 T1+0.5f (1/2)=0.9397933 ………… The interval of integrate [a,b]: a=0 b=1 error :e=0.00001 |Ti+1-Ti| i=1 T1=0.920735 i=2, T2=0.939793 0.019058 i=3, T4=0.944514 0.004720 i=4, T8=0.945691 0.001177 i=5, T16=0.945985 0.000294 i=6, T32=0.946059 0.000073 i=7, T64=0.946077 0.000018 i=8, T128=0.946082 0.000005 integrate=0.946082 准确值：0.946083 由前面的结论： 则， ——梯形加速公式 这说明，可以用低阶的公式组合后成为一个高阶的公式。 类似， ——抛物线加速公式 还可得， ——柯特斯加速公式(Romberg公式) … … … … … … ? T1 = ) 0 ( 0 T ? T8 = ) 3 ( 0 T ? T4 = ) 2 ( 0 T ? T2 = ) 1 ( 0 T ? S1 = ) 0 ( 1 T ? R1 = ) 0 ( 3 T ? S2 = ) 1 ( 1 T ? C1 = ) 0 ( 2 T ? C2 = ) 1 ( 2 T ? S4 = ) 2 ( 1 T T16 = ) 4 ( 0 T S8 = ) 3 ( 1 T C4 = ) 2 ( 2 T R2 = ) 1 ( 3 T 用Romberg公式计算： Newton-Cotes积分公式，可以知道n为偶数时，n＋1个点数值积分公式有n＋1阶精度。是否有更高的代数精度呢？ n +1个点的数值积分公式，最高可以到多少代数精度？ 一点数值积分 0次代数精度 1次代数精度 例： 1次代数精度 例：在两点数值积分公式中,如果节点也作为未知量,则有4个未知量可以列出4个方程:(以f(x)在[-1,1]为例) 可解出： 得数值积分公式 具有3次代数精 度，比梯形公式 1次代数精度高 n +1个求积节点的数值积分公式，最高2n+1次代数精度。 证明： 取 易知： 也就是说，数值积分公式，对一个2n+2次的多项式是有误差的， 所以，n＋1个节点的数值积分公式不超过2n＋1次代数精度。 定理 如何构造最高精度2n+1次的公式？ [a,b]上的n+1次正交多项式的n+1个零点为Gauss点，数值积分公式有2n+1次的代数精度 (2)求出pn+1(x)的n+1个零点x0 , x1 , x2 , … xn 即为Gauss点. (1)求出区间[a,b]上的n+1次正交多项式pn+1(x) . (3)计算求积系数 Gauss型求积公式的构造方法 区间[-1,1]上的Gauss型求积公式,称为Gauss-Legendre求积公式,其Gauss点为Legendre多项式的零点. Gauss-Legendre求积公式 ——中矩形公式 一次代数精度 一个节点： 两个节点： 0.5555555556 0.8888888889 0.1012285363 0.2223810345 0.3137066459 0.3626837834 ±0.9602898565 ±0.7966664774 ±0.5255324099 ±0.1834346425 8 0.1294849662 0.2797053915 0.3818300505 0.4179591837 ±