################################################################### 線形代数の予備知識:固有値・固有ベクトルと逆問題 ################################################################### ========================================================= 固有値・固有ベクトル ========================================================= 行列 :math:`A` を考えた時,行列 :math:`A` が作る写像には,入力されたベクトル :math:`v` を定数倍するような方向が存在し得る. このある行列に対する特別な方向ベクトルのことを「固有ベクトル」,定数倍を「固有値」という.つまり, .. math:: A v = \lambda v を満たすベクトル :math:`v` を求めればよく,これを満たすn本の固有ベクトル :math:`v_k \ (k=1,2,\cdots,n)` は対応する固有値 :math:`\lambda_k` に対して, .. math:: A v_k = \lambda v_k \ \ \ \ (k=1,2,\cdots,n) を満たす. :math:`v_k \ \ (k=1,2,\cdots,n)` は互いに **線型独立** である. :math:`k=1,2,\cdots,n` 全ての :math:`k` でこの条件を満たすのでn本の固有ベクトルを集めた行列 :math:`V=[v_1,v_2,\cdots,v_n]` を用いて, .. math:: A V = V \Lambda と表現することができる.ここで対角行列 :math:`\Lambda` は, .. math:: \Lambda = \begin{bmatrix} \lambda_1 & 0 & \cdots & 0 \\ 0 & \lambda_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \lambda_n \end{bmatrix} である. 行列 :math:`A` の対角化,対角行列 :math:`\Lambda` は, .. math:: \Lambda = V^{-1} A V 逆に,行列 :math:`A` は対角行列 :math:`\Lambda` を用いて, .. math:: A = V \Lambda V^{-1} と表現することができる.行列 :math:`A` の対角化表示 :math:`A=V \Lambda V^{-1}` を求めることができた. ========================================================= 固有値・固有ベクトルと逆行列 ========================================================= 固有値・固有ベクトルが既知であれば,これらを利用して **行列反転** (つまりは **逆行列** )を計算することができる.まず,行列反転する対象となる連立方程式を .. math:: M v = w と書いておく. ここで, :math:`v, w` は :math:`M` の線型独立な固有ベクトル :math:`v_k` を用いた線型結合で表すことができて, .. math:: v = \sum_k^n a_k v_k .. math:: w = \sum_k^n b_k v_k と書ける. :math:`V` を用いて,行列-ベクトル積の形式で書けば, .. math:: v = V a .. math:: w = V b である. :math:`Mv=w` へ代入することで, .. math:: MVa &= Vb \\ b &= V^{-1} M V a = \Lambda a k成分を考えれば, .. math:: a_k = \dfrac{ b_k }{ \lambda_k } が得られる.もし, :math:`V^{-1}` が既知であれば, .. math:: b = V^{-1} w から, :math:`w` の線型結合の係数 :math:`b_k` をまとめた :math:`b` を求めることができ,これを用いて, .. math:: v = \sum _k^n a_k v_k = \sum _k^n \dfrac{b_k}{\lambda_k} v_k から行列反転( :math:`Mv=w` から :math:`v` を求める操作 )を行うことができる.線型結合を用いて,くどく書かれているが,これはつまり, .. math:: M ^{-1} = V \Lambda^{-1} V^{-1} を計算しているのと等価である. .. note:: **固有値,固有ベクトルが分かっていれば,逆行列を直接求めなくとも,** :math:`M ^{-1} = V \Lambda^{-1} V^{-1}` **により,逆行列を求めることができる** .行列反転に帰着されていた逆問題にとって, **固有値・固有ベクトルが計算上の強力なツールとなり得る.** .. note:: 単に逆行列を演算するのに比較し,対角表示した行列から逆問題を解くことには以下の利点が存在する. * 固有値・固有ベクトル毎の成分表示が行える.これは固有値 :math:`\lambda_k` に対応する成分がどの程度強いのかを表示する操作(主成分分析) に相当し,逆問題解析上の利点が生じる. * 非正方な行列の演算 *