モーメント法とは

5.1　モーメント法とは

線形空間と演算子

　線形作用素あるいは演算子（linear operator）を

L

，既知の強制関数（known forcing function）を

g

，未知関数を

f

として，

\begin{matrix} (1) & L (f) = g \end{matrix}

で表される場合を考える．まず，内積（inner product）

⟨ f, g ⟩

あるいはモーメント（moment）は次の関係を満足する．

\begin{aligned} (2) & ⟨ f, g ⟩ = ⟨ g, f ⟩ \\ (3) & ⟨ α f + β g, h ⟩ = α ⟨ f, h ⟩ + β ⟨ g, h ⟩ \\ (4) & ⟨ f^{*}, f ⟩ {\begin{cases} > 0 & (f \neq 0) \\ = 0 & (f = 0) \end{cases} \end{aligned}

ただし，

α

，

β

はスカラー，

*

は複素共役を示す．また，随伴作用素（adjoint operator）

L^{a}

は次式によって定義される．

\begin{matrix} (5) & ⟨ L f, g ⟩ = ⟨ f, L^{a} g ⟩ \end{matrix}

ここで，

L^{a} = L

なら作用素は自己随伴（self-adjoint）であるという．

f

が実数ならば，

L f

もまた実数である．

f

が実数あるいは複素数のとき，

⟨ f^{*}, L f ⟩ > 0 (f \neq 0)

なら演算子は正定値（positive definite），

⟨ f^{*}, L f ⟩ \geq 0 (f \neq 0)

なら半正定値（positive semidefinite），

⟨ f^{*}, L f ⟩ < 0 (f \neq 0)

なら負定値（negative definite）であるという．いま，任意の

g

に対して

L (f) = g

の解

f

が唯一決まるとき，逆作用素（inverse operator）

L^{- 1}

より，

\begin{matrix} (6) & f = L^{- 1} (g) \end{matrix}

【例】変域 $0 \leq x \leq 1$ において定義される $g (x)$ ， $f (x)$ について， $\begin{matrix} (7) & - \frac{d^{2} f}{d x^{2}} = g (x) \\ (8) & f (0) = f (1) = 0 \end{matrix}$ とする境界値問題を考えよう．

　演算子

L

は，

\begin{matrix} (9) & L = - \frac{d^{2}}{d x^{2}} \end{matrix}

これより，与式は，

\begin{matrix} (10) & L (f) = f (x) \end{matrix}

この問題における内積（inner product）

⟨ f, g ⟩

は，例えば，

\begin{matrix} (11) & ⟨ f, g ⟩ = \int_{0}^{1} f (x) g (x) d x \end{matrix}

あるいは，

\begin{matrix} (12) & ⟨ f, g ⟩ = \int_{0}^{1} w (x) f (x) g (x) d x \end{matrix}

ただし，

x (x) > 0

は任意の重み関数（arbitrary weighting function）を示す．式

(11)

より，式

(10)

と

g

の内積は，

\begin{array}{rcl} ⟨ L f, g ⟩ & = & \int_{0}^{1} L (f) g (x) d x \\ (13) & = & \int_{0}^{1} (- \frac{d^{2} f}{d x^{2}}) g d x \end{array}

部分積分して，

\begin{array}{rcl} ⟨ L f, g ⟩ & = & {[(- \frac{d f}{d x}) g]}_{0}^{1} - \int_{0}^{1} (- \frac{d f}{d x}) \frac{d g}{d x} d x \\ (14) & = & - {[\frac{d f}{d x} g]}_{0}^{1} + \int_{0}^{1} \frac{d f}{d x} \frac{d g}{d x} d x \end{array}

上式の第2項をさらに部分積分して，

\begin{array}{rcl} (15) & \int_{0}^{1} \frac{d f}{d x} \frac{d g}{d x} d x & = & {[f \frac{d g}{d x}]}_{0}^{1} - \int_{0}^{1} f \frac{d^{2} g}{d x^{2}} d x \end{array}

これより，式

(14)

の内積は，

\begin{array}{rcl} ⟨ L f, g ⟩ & = & - {[\frac{d f}{d x} g]}_{0}^{1} + ({[f \frac{d g}{d x}]}_{0}^{1} - \int_{0}^{1} f \frac{d^{2} g}{d x^{2}} d x) \\ (16) & = & \int_{0}^{1} f (- \frac{d^{2} g}{d x^{2}}) d x + {[f \frac{d g}{d x} - \frac{d f}{d x} g]}_{0}^{1} \end{array}

上式の第2項は，境界条件の項である．いま，

\begin{matrix} (17) & g (0) = g (1) = 0 \end{matrix}

のとき，この項はゼロになり，

\begin{array}{rcl} ⟨ L f, g ⟩ & = & \int_{0}^{1} f (- \frac{d^{2} g}{d x^{2}}) d x \\ (18) & = & ⟨ f, L^{a} g ⟩ \end{array}

このとき，

\begin{matrix} (19) & L = L^{a} = - \frac{d^{2}}{d x^{2}} \end{matrix}

ゆえ，この作用素は自己随伴である．　また，

\begin{array}{rcl} ⟨ f^{*}, L f ⟩ & = & \int_{0}^{1} f^{*} (- \frac{d^{2} f}{d x^{2}}) d x \\ = & {[f^{*} (- \frac{d f}{d x})]}_{0}^{1} - \int_{0}^{1} \frac{d f^{*}}{d x} (- \frac{d f}{d x}) d x \\ = & \int_{0}^{1} \frac{d f^{*}}{d x} \frac{d f}{d x} d x \\ (20) & = & \int_{0}^{1} {| \frac{d f}{d x} |}^{2} d x \end{array}

したがって，

f

が複素数であっても

L

は正定値である．

重み付き残差法

　電磁界解析では，

L

は微積分演算子（integro-differential operator），

g

は既知の駆動関数（known driving function），例えば，与えられた入射波，

f

は解くべき未知電流などである．ここでは，

f

を基底関数（basis functions）

f_{n}

で次のように展開する．

\begin{matrix} (21) & f ≃ \sum_{n = 1}^{N} a_{n} f_{n} \equiv \bar{f} \end{matrix}

ただし，

a_{n}

は未知係数（unknown coefficients）を示す．

\bar{f}

は

f

の近似解ゆえ，次のように残差

R

が生じる．

\begin{matrix} (22) & R = g - L (\bar{f}) = g - \sum_{n = 1}^{N} a_{n} L (f_{n}) \end{matrix}

この残差

R

に重みをかけて領域全体で積分したものがゼロになるように，

\begin{array}{rcl} \int w_{m} R d Ω & = & \int w_{m} {g - L (\bar{f})} d Ω \\ (23) & = & 0 (m = 1, 2, \dots, N) \end{array}

とし，未知定数

a_{n} (n = 1, 2, \dots, N)

を定める方法が重み付き残差法（method of weighted residuals）である．ただし，

w_{m}

は重み関数（weighting function）を示す．

マトリクス方程式

　演算子

L

は線形ゆえ，

\begin{matrix} (24) & L (\bar{f}) = \sum_{n = 1}^{N} a_{n} L (f_{n}) \end{matrix}

これより，式

(23)

は，

\begin{matrix} (25) & \int w_{m} {g - \sum_{n = 1}^{N} a_{n} L (f_{n})} d Ω = 0 (m = 1, 2, \dots, N) \end{matrix}

あるいは，

\begin{matrix} (26) & \sum_{n = 1}^{N} a_{n} \int w_{m} L (f_{n}) d Ω = \int w_{m} g d Ω (m = 1, 2, \dots, N) \end{matrix}

内積を式

(11)

のように表すと，

\begin{matrix} (27) & \sum_{n = 1}^{N} a_{n} ⟨ w_{m}, L f_{n} ⟩ = ⟨ w_{m}, g ⟩ (m = 1, 2, \dots, N) \end{matrix}

これは次のように行列表示できる．

\begin{matrix} (28) & [Z] a = b \end{matrix}

ここで，

[Z]

は

N \times N

正方行列で，

\begin{array}{rcl} [Z] & = & (\begin{array}{c} z_{11} & z_{12} & \dots & z_{1 N} \\ z_{21} & z_{22} & \dots & z_{2 N} \\ ⋮ & ⋮ \\ z_{N 1} & z_{N 2} & \dots & z_{N N} \end{array}) \\ (29) & = & (\begin{array}{c} ⟨ w_{1}, L f_{1} ⟩ & ⟨ w_{1}, L f_{2} ⟩ & \dots & ⟨ w_{1}, L f_{N} ⟩ \\ ⟨ w_{2}, L f_{1} ⟩ & ⟨ w_{2}, L f_{2} ⟩ & \dots & ⟨ w_{2}, L f_{N} ⟩ \\ ⋮ & ⋮ \\ ⟨ w_{N}, L f_{1} ⟩ & ⟨ w_{N}, L f_{2} ⟩ & \dots & ⟨ w_{N}, L f_{N} ⟩ \end{array}) \end{array}

行列要素

z_{m n}

は，

\begin{matrix} (30) & z_{m n} = \int w_{m} L (f_{n}) d Ω = ⟨ w_{m}, L f_{n} ⟩ \end{matrix}

また，列ベクトル

a

，

b

は，

\begin{matrix} (31) & a = (\begin{matrix} a_{1} \\ a_{2} \\ ⋮ \\ a_{N} \end{matrix}), b = (\begin{matrix} b_{1} \\ b_{2} \\ ⋮ \\ b_{N} \end{matrix}) = (\begin{matrix} ⟨ w_{1}, g ⟩ \\ ⟨ w_{2}, g ⟩ \\ ⋮ \\ ⟨ w_{N}, g ⟩ \end{matrix}) \end{matrix}

列ベクトルの要素

b_{m}

は，

\begin{matrix} (32) & b_{m} = \int w_{m} g d Ω = ⟨ w_{m}, g ⟩ \end{matrix}

行列

[Z]

の逆行列

[Z]^{- 1}

が存在する場合，

\begin{matrix} (33) & a = [Z]^{- 1} b \end{matrix}

これより，

\bar{f}

は，

\begin{matrix} (34) & \bar{f} = \sum_{n = 1}^{N} a_{n} f_{n} = f_{T} a = f_{T} [Z]^{- 1} b \end{matrix}

ここで，

f_{T}

は

f

の転置を示し，

\begin{matrix} (35) & f = (\begin{matrix} f_{1} \\ f_{2} \\ ⋮ \\ f_{N} \end{matrix}) \end{matrix}

基底関数については，区分的基底関数（local or subsectional basis functions）と全領域基底関数（global or entire-domain basis functions）があり，スカラーだけでなくベクトルで表される場合もある．また，積分の次元については，線状素子の散乱問題などは1次元の線積分，境界面上の等価波源に対する問題では2次元の面積分，分布波源に対する電磁波問題では3次元の体積積分となる．電磁界解析の分野では，このように線形代数方程式に変換して行列計算によって統一的に解く方法の総称をモーメント法（method of moments: MoM）と呼んでいる

^{†}

．このとき，重み関数は，試行関数（testing function）ともいう．現在では，電磁界が満足する積分方程式を求めてこのように計算する方法を単に，積分方程式法（integral equation methods）ともいう．

†

Roger F. Harrington, “Matrix Methods for field problems,” Proc. of IEEE, vol.55, no.2, pp.136-149 (1967).
Roger F. Harrington, “Field Computation by Moment Methods,” Chapter 4. Wire antennas and scatters, Wiley-IEEE Press (1993).
“Method of Moments in Antennas and Scattering,” R. C. Hansen, Editor, Artech House (1990).

5.1 モーメント法とは

線形空間と演算子

重み付き残差法

マトリクス方程式

5.1　モーメント法とは