next up previous
Next: About this document Up: システム数理IV 第3回 Previous: 5 刻み幅調節と埋め込み型公式

6 陰的公式

陰的公式については、以下の強力な結果が知られている:

s段陰的公式の到達可能次数は 2s である。」

これのため、いくらでも次数の高い公式をつくることができる。陽的な場合に は8次よりも高次のものはほとんど知られていないのとは対照的である。

これは単なる存在証明ではなく、実際に 2s 次の公式が構築されている。こ れは実はユニークに決まり、「陰的ガウス公式」と呼ばれるものがそれである。 以下、陰的ガウス公式と、それに関係したいくつかの公式についてその導出の 原理と計算法を簡単に説明しよう。

6.1 ガウス型数値積分

ガウス型公式について述べる前に、まずそれと非常に関係の深い数値積分公式 について説明しておく。

数値積分とは、一変数であれば

と、微分方程式ではなく単に関数として与えられたものに対して、ある固定区 間 での定積分を数値的に求めることである。一般に、解析的に 積分ができないようなものを計算しないといけないことはよくあるので、そう いうものが必要になることは多い。

数値積分法として、すでに皆さんが知っているであろうものは、たとえば

といったものであろう(知らない人は適当な数値計算の教科書で勉強したほう がいいかも)。なお、ロンバーグ加速については後でもう少し詳しく説明する。

例えば台形公式は、等間隔に区間を分けたとして、ある部分区間での積分を

と、両端を結んだ台形の面積で近似するということに相当する。これに対して、 もっとましな方法というのを理論的に考えられないであろうか?

一つの方法は、区間の中の点をいくつかとって、それらを通る補間多項式を作 り、それを積分することである。中の点を等間隔にとったものをニュートン・ コーツ型の公式という。

さて、一般に補間するための点をp個とった時、 p-1次の多項式を持って くればそのすべての点を通るように出来る。ニュートン・コーツ型の公式はそ のような最小次数多項式を使うものである。例えばもとの関数が実際にp-1 次以下の多項式ならば、厳密に正しいものが求まる。

この補間の原理を説明しておく。例えば補間に使う点が とあ るとすれば、以下のような多項式を考える

この多項式は、 で 1、それ以外のすべての で0になる。従って、 補間多項式は

とすればいい。

数値積分に使うには、あらかじめ の積分を求めておいて、それと を掛けて和をとるだけである。 等間隔でなく点をとることで精度をあげることはできないかというのが、ここ での問題である。

ここで、一般的に、不等間隔に点をとってなんとかするという観点から考える こともできるが、なかなか難しいので答のほうから考えていく。

答は、「ガウス・ルジャンドルの多項式の零点を点にとると、次数を最大にで きる」というものである。ガウス・ルジャンドルの多項式は、有限区間で帰納 的に定義できる直交関数系である。つまり、最初の項を定数として、k番目 の項は k次多項式であって、それまでの項すべてと直交するようにとるので ある。つまり、関数基底にグラム・シュミットの直 交化を適用したものがガウス・ルジャンドルの多項式である。

だいぶ疲れたので、今回の資料はここまで。次は、直交多項式をつかった積分 公式と、数値積分と陰的ルンゲ・クッタ公式の関係についてと、具体的な計算 法について。


Jun Makino
Mon May 18 20:15:02 JST 1998