6月18日宿題用テキスト

Gauss-Legendre 求積公式の導出
千葉豪
平成 26 年 4 月 9 日
以下の (N − 1) 次の多項式 f (x) を考える。
f (x) = aN −1 xN −1 + aN −2 xN −2 + ... + a1 x + a0
(1)
この多項式の区間 [−1, 1] の積分を以下の求積公式で求めるとする。
Z
M
X
1
−1
f (x)dx =
f (xm )wm
(2)
m=1
上式右辺における xm が求積公式の離散点、wm がその重みに対応する。
式 (2) の左辺は、以下のように解析的に解くことが出来る。
Z
Z
1
−1
f (x)dx =
=
=
1
−1
aN −1 xN −1 + aN −2 xN −2 + ... + a1 x + a0 dx
1
1
1
aN −1 xN +
aN −2 xN −1 + ... + a1 x2 + a0 x
N
N −1
2
N
−1
X
1
−1
an αn
(3)
n=0
ここで、αn は定数であり、次のように書ける。


2
, (n = even),
αn =
n+1

0,
(n = odd)
(4)
一方、式 (2) の右辺も以下のように書き下せる。
M
X
f (xm )wm =
m=1
M X
−1
−2
aN −1 xN
+ aN −2 xN
+ ... + a1 xm + a0 wm
m
m
m=1
= aN −1
M
X
!
−1
xN
m wm
+ ... + a1
m=1
=
M
X
m=1
N
−1
X
an βn
!
xm wm + a0
M
X
!
wm
m=1
(5)
n=0
ここで、
βn =
M
X
m=1
1
xnm wm
(6)
である。
以上より、式 (2) が任意の多項式で成り立つ、すなわち任意の an の値について成り立つ
ためには、以下の条件が必要となることが分かる。
αn = βn ,
n = 0, 1, ..., N − 1
(7)
いま我々は求積セットとして、M 個の (xm , wm ) のセットを考えているが、式 (7) から明
らかなように方程式の数は N 本であるため、M = N/2 のときに一意的に xm 、wm の値が
得られる。すなわち、N/2 点の求積セットにより (N − 1) 次の多項式の求積が厳密に行われ
る（N 点の求積セットにより (2N − 1) 次の多項式の求積が厳密に行われる）ということが
分かる。
(2N − 1) 次多項式の積分が N 点の Gauss-Legendre 求積により厳密に求まることが知ら
れている。従って、上記で述べたような決め方に基づく求積セットは Gauss-Legendre 求積
セットとなるはずである。ここでは、(2N − 1) 次多項式の積分が N 点の Gauss-Legendre
求積により厳密に求まることを、文献 [1] の記述に従って示す。
(2N − 1) 次の任意の多項式 g2N −1 (x) が [−1, 1] で定義されているとする。そして、(N − 1)
次の多項式 GN −1 (x) が以下の条件を満たすとする。
GN −1 (xm ) = g2N −1 (xm ), xm ∈ [−1, 1],
Z
Z
1 1
1 1
(ii)
GN −1 (x)dx =
g2N −1 (x)dx
2 −1
2 −1
(i)
m = 1, 2, ..., N
(8)
(9)
すなわち、GN −1 (x) は、区間 [−1, 1] の積分が g2N −1 (x) の積分と同一となり、また、区間
[−1, 1] の N 個の点で g2N −1 (x) と交差する、g2N −1 (x) よりも低次の関数である。
xm が与えられている場合に GN −1 (x) の x ∈ (−1, 1) における積分を求積によって求める
ことを考えよう。任意の GN −1 (x) に対して求積が成立するためには式 (7) が成り立つ必要
がある。この方程式の本数は N なので、未知数 wm として N 個が必要となる。すなわち、
N 点の求積により、GN −1 (x) の積分は以下の式で厳密に行うことが出来る。
Z
1
−1
GN −1 (x)dx =
N
X
wm GN −1 (xm )
(10)
m=1
従って、条件 (i)(ii) を満足する GN −1 (x) を用いることにより、多項式 g2N −1 (x) の積分を
求積法によって以下のように厳密に計算できることが分かる。
1
2
Z
1
−1
g2N −1 (x)dx =
1
2
Z
1
−1
GN −1 (x)dx =
N
N
1X
1X
wm GN −1 (xm ) =
wm g2N −1 (xm ) (11)
2 m=1
2 m=1
さて、それでは、条件 (i)(ii) を満足する xm はどのように決められるであろうか。
ここで、(2N − 1) 次の多項式 f2N −1 (x) を、g2N −1 (x) と GN −1 (x) を用いて次のように定
義できるであろう。
f2N −1 (x) = g2N −1 (x) − GN −1 (x)
(12)
2
関数 f2N −1 (x) は GN −1 (x) と g2N −1 (x) の差異に対応しているが、条件 (i) より f2N −1 (xm ) = 0
となることから、f2N −1 (x) は N − 1 次の多項式 FN −1 (x) を用いて次のように書くことが出
来る。
f2N −1 (x) = (x − x1 )(x − x2 )...(x − xN )FN −1 (x)
(13)
さらに、条件 (ii) より
1
2
Z
1
−1
f2N −1 (x)dx = 0
(14)
であるため、FN −1 (x) について以下の式が成り立たなければならない。
Z
1
−1
(x − x1 )(x − x2 )...(x − xN )FN −1 (x)dx = 0
(15)
FN −1 (x) は (N − 1) 次の多項式であるため、(N − 1) 次までの Legendre 多項式で記述する
ことが出来る。従って、上式が任意の g2N −1 (x) について成り立つためには、(x − x1 )(x −
x2 )...(x − xN ) が PN (x) の定数倍でなければならないことになる。以上より、関数 GN −1 (x)
と g2N −1 (x) の交差する点 xm は、Legendre 多項式 PN (x) の零点に対応することが分かる。
3

Download Report