1 固有値問題の基礎
1.1 定義と標準形
1.1.1 行列の固有値問題
行列 \(A\) に対し、非零ベクトル \(x\) が \(Ax=\lambda x\) を満たすとき、\(\lambda\) を \(A\) の固有値、\(x\) を固有ベクトルという。ここで \(x\) は作用後も同じ方向を保ち、倍率 \(\lambda\) の分だけ長さ(厳密にはベクトルのスケール)が変わるという幾何学的意味をもつ。固有値問題は線形写像の「方向ごとの増減率」を取り出す作業として理解できる。
1.1.2 作用素の固有値問題
実数・複素数のベクトル空間から別のベクトル空間への線形作用素 \(T\) に対し、非零要素 \(f\) が \(Tf=\lambda f\) を満たすとき、\(\lambda\) を作用素の固有値と呼び、\(f\) を固有ベクトル(固有関数)に相当する対象として扱う。有限次元の場合は行列と同様に記述できるが、無限次元では位相・完備性・閉包などの構造が重要になり、同じ式でも解集合の振る舞いが大きく変わる。
1.2 固有値・固有ベクトルの性質
1.2.1 固有多項式と解の関係
固有値は特性方程式 \(\det(A-\lambda I)=0\) の解として与えられる。これは \(\lambda\) が \(A-\lambda I\) の非自明な核を作ることと等価であり、固有多項式の零点に対応する。多項式の係数は行列の要素から決まり、固有値の集合には代数的制約がかかるため、変数変換や行列の積の性質などが固有値の推定に影響する。
1.2.2 固有空間と幾何学的重複度
固有値 \(\lambda\) に対する固有空間は \(\ker(A-\lambda I)\) として定義され、次元により幾何学的重複度が測られる。幾何学的重複度は同じ固有値に属する独立な固有ベクトルの数に対応し、一般に代数的重複度(特性多項式での重さ)以下となる。不変部分空間がどれだけ豊富かを知る指標として有用であり、対角化可能性とも関係する。
1.3 線形代数的背景
1.3.1 既約分解と不変部分空間
線形作用素(行列)に対し、不変部分空間とは \(A\) の作用でその部分集合が保たれる部分空間のことを指す。既約分解は、作用がそれ以上分解できないブロック構造の単位に分かれる様子を体系化する考え方である。有限次元では、適切な基底変換により作用の構造がブロック化され、そのブロックごとの固有値や振る舞いを追跡できる。
1.3.2 対角化とジョルダン標準形
対角化とは、基底変換により \(A\) が対角行列 \(D\) に移ることを言う。これは十分条件として、各固有値について幾何学的重複度が代数的重複度に一致することが挙げられる。対角化できない場合にはジョルダン標準形が用いられ、固有値に加えて「鎖状の一般化固有ベクトル」が現れる。ジョルダンブロックは、冪に関する振る舞い(累乗や指数作用素の計算)に特有の多項式係数をもたらすため、理論と計算の両面で重要になる。
2 種類と拡張
2.1 特殊な行列・作用素
2.1.1 エルミート(自己共役)問題
エルミート(複素内積空間での自己共役)な行列 \(A=A^*\) では、固有値は実数となり、固有ベクトルは互いに直交して選べる。さらに十分条件が揃うと、基底全体を固有ベクトルで構成する直交対角化が成立する。これは数値計算における安定性や、物理での観測量に対応する理論と整合的である。
2.1.1.1 実固有値と直交固有ベクトル
実固有値の性質はスペクトルの解釈を明確にし、直交性はモード分解を容易にする。具体的には、異なる固有値に属する固有ベクトルの内積がゼロになり、係数の計算が独立に行えるため、エネルギーのような量の評価に直結する。
2.1.2 歪エルミート(反自己共役)問題
反自己共役な行列 \(A^*=-A\) では、固有値は純虚数(零を含む)になる。したがって指数型の時間発展では回転的な挙動と結びつきやすく、ノルムの変化と関連する議論が組み立てやすい。符号や虚部の扱いが安定性の議論に直結し、物理モデルの対称性とも整合することが多い。
2.1.3 一般化固有値問題(随伴行列を含む)
一般化固有値問題は、行列 \(A,B\) について \( Ax=\lambda Bx \) の形で定式化される。特に \(B\) が正則でないと分母的な解釈が難しくなるが、適切な仮定のもとでは有限次元でも体系的に解析できる。随伴(エルミート共役)を含むモデルは、対象の対称性を反映する形で現れ、スペクトルの位置や直交性が制約される。
2.2 連続体への拡張
2.2.1 ストゥルム=リウヴィル型
ストゥルム=リウヴィル型は、微分作用素に境界条件を組み合わせた形で固有値問題を扱う枠組みである。典型的には二階微分方程式を対象とし、固有値は実数であり、固有関数はある意味で直交性を満たすように設計される。さらに固有値列が離散的に並ぶ条件が整理され、有限次元の直交展開の連続版のように振る舞う。
2.2.2 偏微分方程式への固有値問題の導出
偏微分方程式の多くは、時間変数を分離して空間部分の固有値問題に帰着させることがある。たとえば拡散・波動・熱などの線形方程式では、空間のモードが固有値とともに時間方向へ指数的または振動的な因子として現れる。境界条件(固定、自由、周期など)は固有値を決める主要因であり、形状や媒質の性質がスペクトルに反映される。
2.3 非線形・一般化固有値問題
2.3.1 固有値が式の中で非線形に現れる場合
非線形固有値問題では \(\lambda\) が単に倍率として現れるのではなく、例えば \(F(\lambda)x=0\) のように方程式の中へ非線形に入る。すると、分岐や多重度、解の連続的変形が線形の場合より複雑になり、解析的に扱うには追加の仮定(正則性、単調性、解析接続など)が必要になる。計算では、線形化の近似や反復法の工夫が重要になる。
2.3.2 標数や特異性が関わる場合
数理的には「解の存在・一意性」を左右する要素として、作用素の標数、方程式の退化点、あるいは特異点近傍での挙動が挙げられる。行列的には特定の基底でランクが落ちる状況が現れ、固有空間の次元が突然変化することがある。連続体では係数の不連続や境界の角部などが特異性を生み、固有値の分布が標準形から外れる。
3 理論的性質
3.1 スペクトル理論の概観
3.1.1 固有値・連続スペクトル・残余スペクトル
無限次元では「固有値」だけでスペクトル全体を表せないことがある。作用素 \(T\) のスペクトルは、\(T-\lambda I\) が有界な逆作用素をもたない値全体として定義され、そこには固有値以外の領域が含まれる。連続スペクトルや残余スペクトルといった分類は、核の有無だけでなく像の閉性、可逆性の失敗の様式に基づく。
3.1.2 スペクトル半径と有界性
行列ではスペクトル半径 \(\rho(A)\) が固有値の絶対値の最大として定義できる。作用素の文脈でも、スペクトル半径は冪や指数の成長率と結びつき、どの程度の増幅が起こり得るかの目安になる。さらに、ノルムによる上界との関係により、有界性や収束性の見通しが得られる。
3.2 安定性とスペクトル
3.2.1 ダイナミクスにおける安定性判定
線形常微分方程式 \(\dot{x}=Ax\) では解が指数関数を通じて固有値により記述され、実部の符号が挙動(増大・減衰・中立)を左右する。離散時間の更新則では固有値の絶対値が重要となり、スペクトル半径が 1 未満かどうかが安定性の判断材料になる。固有値の位置は、初期条件によらない漸近挙動を捉えるための中心的指標である。
3.2.2 境界条件が与える影響
連続体モデルでは境界条件が固有値問題を規定し、結果として安定性にまで波及する。たとえば固定境界と自由境界は固有関数の形を変えるだけでなく、固有値の列の大小関係にも影響し、減衰率や振動周波数を変化させる。幾何形状の差(境界の曲率や連結性)もスペクトル分布を変える要因として現れる。
3.3 数学的性質の要点
3.3.1 存在定理の考え方
存在定理は、固有値や固有関数がどのような条件下で存在するか、またスペクトルがどのように分解されるかを保証する。有限次元では特性多項式の零点により固有値の存在が自動的に満たされるが、無限次元ではコンパクト性や可積分性、エネルギー汎関数の評価などが鍵になる。理論的には「有界列から収束部分列を取り出す」タイプの議論が多く見られ、固有関数の構成に繋がる。
3.3.2 正則性・冪性(線形化との関係)
固有関数がどれだけ滑らか(微分可能)かは、方程式の係数と境界条件に依存する。正則性が高いほど、固有値の数値近似やモードの評価が精度良く行える。さらに非線形問題を線形化する際、ヤコビアンに対応する線形作用素の固有値が局所的な振る舞いを決めることがある。ここで得られる冪性(冪に関する構造)や一般化固有ベクトルの有無は、単純な指数近似では捉えきれない補正項の存在につながる。
4 計算方法と応用
4.1 数値計算の基本方針
4.1.1 パワー法と逆反復法
パワー法は、適当な初期ベクトルを繰り返し \(x_{k+1}=Ax_k\) のように更新し、支配的な固有値に対応する方向へ収束させる方法である。最大絶対値の固有値が支配的であるほど速くなる。逆反復法は \((A-\mu I)^{-1}\) を用いる発想で、狙った位置 \(\mu\) に近い固有値を抽出しやすくする。実装では連立一次方程式の解を繰り返すため、行列の性質(疎性や前処理)が計算効率を左右する。
4.1.2 シフト付き反復法
シフト付き手法では \(A-\mu I\) を解く操作を組み合わせ、固有値の近傍へ反復を誘導する。シフトの選び方によって収束速度や安定性が変わり、誤差の伝播も影響を受ける。数値的には丸め誤差の蓄積を抑えるため、変換の設計と反復中の正規化が重要になる。
4.2 ブレイキングダウン(分解)に基づく手法
4.2.1 QR法
QR法は反復により行列を上三角(あるいはより簡単な形)へ近づけ、対角成分の収束から固有値を推定する代表的手法である。行列の分解 \(A=QR\) を繰り返し適用し、類似変換を通じてスペクトルを保存する点が本質である。エルミート行列に対しては特に安定性が良く、実務で広く利用される。
4.2.2 ハウスホルダー変換と三重対角化
三重対角化は、固有値計算の前処理として用いられる変換で、一般行列をより構造のある形へ写す。ハウスホルダー変換は直交変換に基づくため、数値的に誤差管理がしやすい。三重対角の形になるとQR反復などのコストが大幅に下がり、疎性がない場合でも計算が現実的になる。
4.3 大規模問題への工夫
4.3.1 反復解法(サブスペース法)
大規模行列では全固有値を直接求めることが難しく、サブスペース法がよく使われる。代表例として、複数ベクトルからなる部分空間上で小行列の固有問題を解き、その結果を元の空間へ戻して反復を行う。精度と計算量のバランスを制御でき、必要な固有値のみ(例えば最小・最大、あるいはスペクトル中の特定領域)を狙って計算できる。
4.3.2 疎行列・前処理の考え方
疎行列では非零要素の数に比例して計算できる操作が望ましいため、行列積や連立一次方程式の解き方が設計の中心になる。前処理は反復法の収束を改善するために、元の作用素を別の作用素で「条件数が良い形に見せる」工夫である。前処理の設計が適切だと、固有値探索や線形方程式の反復が安定し、反復回数が減る。
4.4 応用領域
4.4.1 振動・モード解析
機械・構造の振動問題では、系の自由度に対応する作用素の固有値が固有振動数、固有ベクトルがモード形状を与える。これにより複雑な応答が、独立なモード成分の重ね合わせとして理解できる。境界条件や材料特性の変化はスペクトルに反映され、設計段階での性能評価に利用される。
4.4.2 安定性解析
制御や力学系では、平衡点近傍の挙動を線形化して得られる固有値で安定・不安定の見通しを立てる。特に境界付近では一般化固有ベクトルや多重度が効き、単純な指数モデルでは現れない項が無視できなくなることがある。ゆえに固有値の幾何学的・代数的情報が実務判断に影響する。
4.4.3 画像処理・機械学習での利用(スペクトル手法)
画像処理や機械学習では、データやグラフの構造を行列(あるいは作用素)として表し、固有分解により低次元表現を得る手法がある。例えば共分散に基づく分解では主要成分が固有ベクトルとして現れ、特徴抽出や圧縮に役立つ。グラフ上の拡散やクラスタリングでもラプラシアンのスペクトルが中心となり、類似度の連結性を反映する形でクラスタの境界が形成される。