龍愛芳，胡軍浩

（中南民族大學 數(shù)學與統(tǒng)計學學院，湖北 武漢430074）

隨著科學的進步，計算機技術(shù)的發(fā)展，很多領(lǐng)域涉及到定積分的計算，因此研究高精度的數(shù)值積分公式是有實際意義的.數(shù)值積分常見的有梯形公式和Simpson公式，它們的計算雖無需提供導數(shù)值，但代數(shù)精度不高.梯形公式有一次代數(shù)精度，Simpson公式有3次代數(shù)精度［1-6］.文獻［7-8］雖給出一個高精度的數(shù)值求積公式，但必須提供求積節(jié)點的一階導數(shù)值.文獻［9］給出了Newton-cotes求積公式的漸近性，雖可大大提高數(shù)值求積公式的代數(shù)精度，但同樣必須提供n+1階導數(shù)值，文獻［10］的求積公式?jīng)]有承襲性.本文從Hermite插值多項式出發(fā)，構(gòu)造了具有誤差量級為O（h5），且不需要計算導數(shù)值，只需要提供求積節(jié)點函數(shù)值的高精度數(shù)值求積公式.

1 高精度數(shù)值求積公式

1.1 Hermite插值多項式的構(gòu)造

構(gòu)造滿足插值條件f（xk）＝p（xk），f（xk+1）＝p（xk+1），f′（xk）＝p′（xk），f′（xk+1）＝p′（xk+1）次數(shù)不超過3的Hermite插值多項式的p（x）.即

上式中：αk（x），αk+1（x），βk（x），βk+1（x）為插值基函數(shù)，它們分別滿足

αk（x），αk+1（x），βk（x），βk+1（x）的表達式分別為

其中：h＝xk+1-xk.

構(gòu)造的插值多項式余項表達式為

式（1）中：ξ在xk與xk+1之間，并且與x有關(guān).

1.2 數(shù)值求積公式的構(gòu)造

對式（1）兩邊求積分，應用廣義積分中值定理，可得到積分中值定理.

定理1 設(shè)函數(shù)f（x）在區(qū)間［x，xk+1］有4階連續(xù)導函數(shù)，則

成立，η在xk與xk+1之間.

由式（2）得到帶有一階導數(shù)的數(shù)值求積分公式，即

由式（2）可知，數(shù)值求積公式（3）具有3次代數(shù)精度.在式（2）中，令xk＝a，xk+1＝x，則有

為了提高數(shù)值求積公式的代數(shù)精度，分析定理1中間點η的漸近性，則可得到定理2.

證明 令

應用3次洛必達法可得

應用式（4）可得

應用式（5），（6）可得

由式（7），可得到具有5次代數(shù)精度的第2個數(shù)值求積公式，記為

應用復化求積得

為了避免求導數(shù)，對上式的數(shù)值求積公式進行修正，應用公式

得到只需計算求積節(jié)點函數(shù)值，無須提供求積節(jié)點導數(shù)值的兩個數(shù)值求積公式，分別為

公式（9），（10）的誤差量級分別為O（h4），O（h5）.

2 數(shù)值試驗

應用梯形公式、求積公式（3），（8）計算，計算結(jié)果如表1所示.應用復化梯形公式、求積公式（9），（10）計算，計算結(jié)果如表2所示.

表1 梯形公式、求積公式（3），（8）的計算結(jié)果Tab.1 Numerical experiment of trapeziod formula and formula（3），（8）

表2 復化梯形公式，求積公式（9），（10）的計算結(jié)果Tab.2 Numerical experiment of compound trapeziod formula and formula（9），（10）

從計算結(jié)果看，公式（3）與公式（8）比梯形公式的精度高很多，但必須提供一階導數(shù)或四階導數(shù)；而公式（9）和公式（10），卻不用計算導數(shù)，計算的節(jié)點函數(shù)值的個數(shù)與復化梯形公式一樣.因此，計算量與復化梯形公式相當.復化梯形公式的誤差量級為O（h2），而公式（9）的誤差量級為O（h4），公式（10）的誤差量級為O（h5）.因此，公式（9）和公式（10）是非常有效的，無須計算導數(shù)的兩個數(shù)值積分公式.

［1］ 李慶楊.數(shù)值分析［M］.武漢：華中科技大學出版社，1986：83-106.

［2］ 李慶楊.數(shù)值計算原理［M］.北京：清華大學出版社，2000：251-282.

［3］ 李毅夫.Simpson公式余項中間點漸進性定理及Simpson公式的改進［J］.貴州師范大學學報：自然科學版，2007，25（4）：67-69.

［4］ 程海來.關(guān)于Simpson公式的兩點注記［J］.數(shù)學的實踐與認識，2007，27（21）：91-93.

［5］ 杜躍鵬，肖澤昌.改進Simpson公式及誤差分析［J］.高師理科學刊，2008，28（4）：27-29.

［6］ 劉彬清.一類高斯求積公式的極限性質(zhì)［J］.工程數(shù)學學報，2003，20（4）：137-139.

［7］ 吳新元.一個高精度數(shù)值積分公式［J］.計算物理，1988，5（4）：473-477.

［8］ 鄭華盛.高精度數(shù)值積分公式的構(gòu)造及其應用［J］.數(shù)學的實踐與認識，2007，37（15）：141-148.

［9］ 劉彬清.Newton-cotes數(shù)值求積公式的漸近性［J］.上海大學學報：自然科學版，2002，8（6）：503-506.

［10］ 劉彬清.一類高斯求積公式的極限性質(zhì)［J］.工程數(shù)學學報，2003，20（4）：137-139.