2
応用数学解説
文献あり

ゲルシュゴリンの定理,あるいは線形代数祭の講評

273
0

はじめに

この記事は,友人が企画した「 線形代数祭 」というイベントに寄せた問題の解説である.

問題

A=(aij)n次複素正方行列とする.0t1とし,F(t)=(fij(t))
fij(t)={aij(i=j),taij(ij)
と定める.

  1. 複素数λF(t)の固有値であるとき,あるi{1,2,,n}が存在して
    |λaii|tji|aij|
    を満たすことを示せ.

  2. 複素数平面上の閉円板Di(t)
    Di(t)={zC||zaii|tji|aij|}
    で定める.D1(1),D2(1),,Dn(1)がどの2つも共通部分を持たないとき,Aは対角化可能であることを示せ.

なお,次の事実は証明せずに使ってよい.多項式関数w=f(z)が,複素数平面の円周C上に零点を持たないとき,Cの内側にあるf(z)の零点の個数Nは重複度を込めて
N=12πiC1wdwdzdz
である.ただし,積分路の向きは反時計回りにとる.

解答 (1)

x=(x1x2xn)TF(t)の固有値λに属する固有ベクトルとすると,λx=F(t)xより,すべてのi{1,2,,n}について
λxi=j=1nfij(t)xj=aiixi+jitaijxj
である.よって,iの値を|xi|=max{|x1|,|x2|,,|xn|}となるように選ぶと
|(λaii)xi|tji|aijxj|t|xi|ji|aij|
となる.x0よりxi0だから,不等式を|xi|で割れば示す式が得られる.

解答 (2)

Di(1)の半径をriとおく.仮定から,ijのときri+rj<|aiiajj|なので
min1i<jn(|aiiajj|rirj4)
は正数である.この値をδとおくと,i<jのとき
(ri+δ)+(rj+δ)ri+rj+|aiiajj|rirj2=|aiiajj|+ri+rj2<|aiiajj|
だから,円板Ei={zC|zaii|ri+δ}はどの2つも共通部分を持たない.

境界Eiはどの円板D1(1),D2(1),,Dn(1)とも交わらないので,小問1よりEiF(t)の固有値は属さない.すなわち,F(t)の固有多項式
p(t;z)=det(zIF(t))
Ei上に零点を持たない.よって,Eiに属するp(t;z)の零点の個数N
N=12πiEi(1p(t;z)ddzp(t;z))dz
である.右辺はtの連続関数であり,Nは整数なので,Nの値はtによらず一定である.よって
N=12πiEi(1p(0;z)ddzp(0;z))dz
だが,p(0;z)=det(zIF(0))の零点でEiに属するものはaiiだけだから,N=1である.

以上により,p(1;z)の零点,すなわちAの固有値は,すべてのEiにちょうど1つ属する.よって,Aは対角化可能である.

講評

この問題はゲルシュゴリンの定理 (Gershgorin circle theorem) と呼ばれる定理を証明させる問題である.ゲルシュゴリンの定理は,固有値の精度保証付き数値計算の基本となる定理である( 大石,2010 ).

ゲルシュゴリンの定理は非常に有名だが,にもかかわらず,厳密な論証が載っている文献はあまり多くない.本問の誘導とは異なる証明も含めて,網羅的なサーベイが ( Chi-Kwong Li & Fuzhen Zhang, 2019 ) にある.

参考文献

[1]
大石進一, 精度保証付き数値計算法の最近の到達点とMATLAB上のツールボックス, 計測と制御, 2010, pp. 273–278
[2]
Chi-Kwong Li & Fuzhen Zhang, Eigenvalue continuity and Gersgorin's theorem, Electron. J. Linear Algebra, 2019, pp. 619–625
投稿日:202424
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

数学、特に応用数学が好きです。Mathlogでは主に、数学とプログラミングを絡めたようなことを書けたらいいなと思っています。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 問題
  3. 解答 (1)
  4. 解答 (2)
  5. 講評
  6. 参考文献