この項目「レイリー・リッツ法」は翻訳されたばかりのものです。不自然あるいは曖昧な表現などが含まれる可能性があり、このままでは読みづらいかもしれません。(原文:en: Rayleigh?Ritz method
)レイリー・リッツ法(レイリー・リッツほう、英: Rayleigh?Ritz method)は、固有値問題に対する数値的近似解法の一つ。レイリー卿とヴァルター・リッツに名をちなむ。物理学上の境界値問題の解法として考案された。
固有値と固有ベクトルの近似をともなう全ての問題に応用でき、分野によってしばしば別名で呼ばれる。量子力学では、系を構成する粒子はハミルトニアンを用いて記述されるが、リッツ法(英語版)では試行波動関数を用いて最低エネルギー固有値に対応する固有関数を近似する。有限要素法の文脈では、数学的に等価なアルゴリズムが一般にリッツ・ガラーキン法(英語版)と呼ばれる。機械工学および構造工学では固有振動モードおよび共鳴周波数を近似する手法としてレイリー・リッツ法およびリッツ法という用語が用いられることが多い。 本手法は1908年から1909年にかけてヴァルター・リッツが発表したもので、リッツ法と呼ぶべきであるという主張もある[1][2]。A. W. Leissa[1]によれば、レイリー卿は1911年にリッツの業績を顕彰する論文を書いたが、彼自身が書籍他の刊行物において本手法をすでに何度も用いていたと述べている。後に異論も出たもののこの主張にくわえ、射影に単一ベクトルを用いる自明な場合、本手法はレイリー商の計算に帰着するという事実もあり、異論もあるもののレイリー・リッツ法という名称が現在まで用いられている。S.Ilanko[2]はリヒャルト・クーラントを引いて、レイリー卿とヴァルター・リッツがそれぞれ独立に、偏微分方程式の境界値問題と変分問題の等価性を活用し、有限のパラメータを決定すればよい極値問題で変分法を置き換えるというアイデアを独自に考案したとする。詳細については、リッツ法 数値線形代数において、レイリー・リッツ法は一般的[3]にサイズ N {\displaystyle N} の正方行列 A ∈ C N × N {\displaystyle A\in \mathbb {C} ^{N\times N}} についての固有値問題 A x = λ x {\displaystyle A\mathbf {x} =\lambda \mathbf {x} } の近似解を得るために用いられる。まず、行列 A {\displaystyle A} はより小さなサイズの行列へと射影される。射影は、その列が正規直交系をなす射影行列 V ∈ C N × m {\displaystyle V\in \mathbb {C} ^{N\times m}} により行われる。行列版のレイリー・リッツ法は最も単純で、以下のように書き下せる。 もし、行列 V ∈ C N × m {\displaystyle V\in \mathbb {C} ^{N\times m}} の列を正規直交基底として張られる線型部分空間が行列 A {\displaystyle A} の固有ベクトルに近い k ≤ m {\displaystyle k\leq m} 個のベクトルを含んでいれば、上記のレイリー・リッツ法はそれら固有ベクトルをよく近似する k {\displaystyle k} 個のリッツベクトルを与える。各リッツ対の精度は、容易に計算できる量 ‖ A x ~ i − λ ~ i x ~ i ‖ {\displaystyle \|A{\tilde {\mathbf {x} }}_{i}-{\tilde {\lambda }}_{i}{\tilde {\mathbf {x} }}_{i}\|} により評価できる。 最も簡単な m = 1 {\displaystyle m=1} の場合、 N × m {\displaystyle N\times m} 行列 V {\displaystyle V} は単位列ベクトル v {\displaystyle v} 、行列 V ∗ A V {\displaystyle V^{*}AV} はレイリー商 ρ ( v ) = v ∗ A v / v ∗ v {\displaystyle \rho (v)=v^{*}Av/v^{*}v} と一致するスカラーとなり、固有値問題の唯一の解は y 1 = 1 , μ 1 = ρ ( v ) {\displaystyle y_{1}=1,\mu _{1}=\rho (v)} 、唯一のリッツベクトルは v {\displaystyle v} それ自体となる。したがって、 m = 1 {\displaystyle m=1} の場合、レイリー・リッツ法はレイリー商の計算に帰着する。 別の有用なレイリー商とのつながりとして、各レイリー対 ( λ ~ i , x ~ i ) {\displaystyle ({\tilde {\lambda }}_{i},{\tilde {\mathbf {x} }}_{i})} に対して μ i = ρ ( v i ) {\displaystyle \mu _{i}=\rho (v_{i})} が成り立ち、したがってリッツ値は対応するレイリー商の理論から導かれるいくつかの性質をもつことが上げられる。
名称
行列の固有値問題への適用
m × m {\displaystyle m\times m} 行列 V ∗ A V {\displaystyle V^{*}AV} を計算する。ここで、 V ∗ {\displaystyle V^{*}} は V {\displaystyle V} の複素共役転置行列とする。
固有値問題 V ∗ A V y i = μ i y i {\displaystyle V^{*}AV\mathbf {y} _{i}=\mu _{i}\mathbf {y} _{i}} を解く。
リッツベクトル x ~ i = V y i {\displaystyle {\tilde {\mathbf {x} }}_{i}=V\mathbf {y} _{i}} およびリッツ値 λ ~ i = μ i {\displaystyle {\tilde {\lambda }}_{i}=\mu _{i}} を計算する
リッツ対と呼ばれる ( λ ~ i , x ~ i ) {\displaystyle ({\tilde {\lambda }}_{i},{\tilde {\mathbf {x} }}_{i})} をを元の行列 A {\displaystyle A} の固有値問題の近似解として出力する。