计算方法之数值微分与积分.pptVIP

* （1）由I(f )与T2n(f )、Tn(f )之间的关系 既然如此，计算出T2n(f )后，不妨取： 显然比T2n(f )更接近于精确积分值。 而事实上，上式就是辛普森积分计算式Sn(f )。 下面以n=1为例，说明Sn(f )与T2n(f )、Tn(f )之间的关系： * 由于： o a=x0 x1=(a+b)/2 x2=b T2(f ) T1(f ) S1(f ) y=f(x) y x * 不难证明： 同理可得： （2）由I(f )与S2n(f )、Sn(f )之间的关系 * 既然如此，计算出S2n(f )后，不妨取： 显然也比S2n(f )更接近于精确积分值。 而事实上，上式就是柯特斯积分计算式Cn(f )。 （3）同理，在计算了C2n(f )后，不妨取 显然将比C2n(f )更接近于精确积分值。 * 实际使用Romberg积分法： （1）计算积分区间等分数分别为1,2,4,8,16,…时的序列T1 ,T2 ,T4 ,T8 ,T16 ,…， （2）计算S1 , S2 , S4 , S8 , …， （3）计算C1 , C2 , C4 , …， （4）计算R1 , R2 , …,直到│R2k-Rk│ε为止。 将各序列统一标记为T，可将它们列于下表： * k 梯形 辛普森 柯特斯 龙贝格 0 1 2 3 4 … … … … … 为便于编程计算，现将序列的计算式改写如下： * 龙贝格积分算法描述 * 一、数据说明 a , b , h — 区间端点和积分步长 MaxD — 区间最大二分次数 T[MaxD][4] — 存放积分序列 EPS — 表示误差控制量 e — Romberg序列中的|R2k-Rk| 二、操作步骤 Step1 输入区间端点a , b Step2 h=b-a ; T[0][0]=h[f(a)+f(b)]/2.0 ； k=1 , e=EPS+1.0 * Step3 While ( kMaxD And eEPS ) Do Step4 Step5 For m=1 To k Do If m3 Then 跳出Step5 EndFor m Step6 If k4 Then e=|T[k][3]-T[k-1][3]| Step7 k=k+1 ; h=h/2.0 Step8 If kMaxD Then Output(“Method Failed”) Else Output(T[k][3]) * 例6．7 利用Romberg积分算法重新计算 使误差不超过0.5?10-6。 k 0 0.9207355 1 0.9397933 0.9461459 2 0.9445135 0.9460869 0.9460830 3 0.9456909 0.9460833 0.9460831 0.9460831 结果表明：积分区间只需二分3次，涉及9个积分节点处函数值的计算，即能满足精度要求。 * 一、知识回顾 n次插值型求积公式至少具有n阶代数精度 二、问题提出 在节点数目固定为n+1（x0 , x1 , … xn）的条件下能否使数值积分公式 高 斯 型 求 积 公 式 具有最大的代数精度？ * 以n=1为例，节点数为2，做线性插值 b a y=L1(x) o y x y=f(x) x0 x1 并非必须选取端点a和b，选择适当位置的x0，x1,可使数值积分A0f(x0)+A1f(x1)比T(f)更精确. * 即未知数变为4个：A0，A1，x0，x1 所以最大代数精度可达到3阶，因为方程组： 有解的条件是m=3，即方程数量等于未知数个数。 同理：n+1个积分节点的插值型数值积分公式最大代数精度可达2n+1。 定义：按此构造的积分公式即Gauss积分公式，对应节点x0,…,xn称Gauss点。 * 三、两点Gauss求积公式 不失一般性，考虑将任何积分作如下变换 换元方法是： 解方程组 得A0=A1=1 * * 四、n+1点的Gauss型求积公式 可以证明n+1次勒让德（Legendre）多项式 例如： 的零点即Gauss点，所得公式又称Gauss-Legendre求积公式。 积分系数Ai仍按下式求解 * n xk Ak 2 ?0.5773503 1 3 ?0.7745967 0 0.5555556 0.8888889 4 ?0.8611363 ?0.3399810 0.3478548 0.6521452 5 ?0.9061798 ?0.538