レッスン1part

Transcription

1 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム レッスン 4 ノルム このレッスンでは行列解析理解の基礎となる事項を学ぶ これらは代数と解析の境界領域にある話題であり 行列計算の誤差解析を理解するための基礎知識である 主題は実または複素有限次元ノルム空間上のノルムと線形変換である このような線形変換の具体例として m 行列を念頭において頂くのがよい このレッスンのキーワードは 線形変換の連続性 ( 有界性 ) 展開係数の有界性 有限次元ノルム空間の完備性 ノルムの同値性 成分ごとの収束とノルム収束の同値性 演算子ノルム ハーン バナハの定理 である 4. 線形変換の有界性と連続性この節では線形変換の連続性の同値な定義を学ぶ これらの事項は次節以降で学ぶ 有限次元ノルム空間とそれらの間の線形変換に関する基礎的な性質を理解する上で必要となる (I) ベクトル列の収束与えられた有限または無限次元ノルム空間 X 内のベクトルの無限列 {,, } { } がa X に収束する to coverge とは実数列の収束 a をいう こ のときa を列 { } の極限 ( または極限値 )lmt といい aまたは lm = aと書く 詳し くいえば aとは 任意のε > に対して自然数 N( ε ) を十分大きくとれば > N( ε ) を 満たすすべての に対して a < ε となる ことである 極限は存在すればひとつしかない 実際 aかつ b なら ノルムの性質よ り a b = ( b) ( a) b + a であるから a b = すなわ ち a b= でなければならない (II) 線形変換の連続性有限または無限次元ノルム空間 X から同じ体上の第 のノルム空間 Y への線形変換 について 次の各主張は互いに同値である : () ( 有界性 ) すべての X に対して M を満たすような正定数 M が存在する ( こ こに は X 上のノルム は Y 上のノルムを表す ) このような は有界である bouded という () ( = における連続性 ) 任意に与えられたε > に対してδ > を十分小さくとれば < δ を満たすすべての X に対して < ε が成り立つ このような は = において連続であるという Coyrght 再履修線形代数研究会

2 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム (3) ( = における連続性 ) ならかならず 証明 () () (3) () を示す が成り立つ () () () (3) < ε / M δ とすれば () により M < Mε / M = ε とする このとき { } が に収束しないとすれば 適当なε > と適当な 自然数の列 < < をとれば ε ( = < < ) が成立するはずである () により < δ ならかならず < ε となるようなδ > がとれるはずである ゆ δ え ある番号から先の はすべて < を満たし < が無限に多くの = の値 ε に対して満たされなければならない しかしこれは ε ( = < < ) と矛盾する (3) () () が真でなければ (3) も真でないことを示せばよい 実際 () が真でなければ すべ ての自然数 に対して > すなわち > を満たす る y = とおけば = / y は真ではあり得ない がとれるはずであ y すなわち y しかし y (III) (II) の各項はさらに以下の各項と同値である ( 証明は練習問題とする ): (4) ( 一様連続性 ) a X を任意かつ特定の一点とすれば 任意 ε > に対してδ > 小さくとれば > ゆえ を十分 a < δ を満たすすべての X に対して a < ε が成り立つ ここにδ はε のみに関係し 特定のa には依存しない (5) a X を任意かつ特定の一点とすれば aは a を意味する 以上を総合するとこういえる : ノルム空間 X より Y への線形変換 が有界または X 内の 任意かつ特定の 点で連続なら X 上で一様連続である 有界性と連続性は同値である この節を終える前に今後の議論によく出てくる次の 3 種の集合を定義しておく : 点 a を中心とする半径 r > の開球 oe shere { X : a < r} 点 a を中心とする半径 r > の閉球 closed shere { X : a r} 点 a を中心とする半径 r > の球の表面 surface { X : a = r} Coyrght 再履修線形代数研究会

3 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム 4. 展開係数の有界性 本節で学ぶ ( 線形変換としての ) 展開係数の有界性 は 有限次元ノルム空間の基本的性 質を導く上でもっとも直接的で使いやすい事実である 実際 次節以降で順次示すように 有限次元ノルム空間内のコーシー列は収束する ( 完備性 ) 有限次元ノルム空間内の有界列は収束する部分列を含む ( 完備性 ) 有限次元ノルム空間上の 種のノルムは同値である ( ノルムの同値性 ) 有限次元ノルム空間上の線形変換は連続である ( 連続性 有界性 ) などの諸事実はすべて本節の事実から出てくる 本節で示すのは次の事実である : Y を与えられた 次元ノルム空間 { b,, b } をY の基底 任意ベクトル y Y をこの基底で展開したものを y = g( y) b+ + g ( y) b とすれば 適当な正定数 α, β に対して () α y g ( ) ( ) y + + g y β y がすべての y Y に対して成り立つ ここに展開係数 g ( ),, ( ) y g y は それぞれ Y からス カラー体への線形変換 ( すなわち 線形汎関数 lear fuctoal) を表し g ( y) + + g ( y) 自体 Y 上のノルムを表す 証明 g ( y),, g ( y) のそれぞれはY からスカラー体への線形変換を表し y ' ( y) + + ( y) g g は Y 上の一つのノルムを表す ことの証明は練習問題とする さて y = g( y) b+ + g ( y) b のノルムをとれば y g ( y) b + + g ( y) b ( g ( y) + + g ( y))ma{ b,, b } () α = /ma{ b,, b } とおけば () の前半が出る ここにα は基底 { b,, b } となる と Y 上のノルム には依存するが y にはまったく依存しないことは明らかである 不等式後半の証明に入る 証明はやや入り組んでいるが 必要となる解析学上の予備知識 は実数および複素数の完備性だけである さて X = R ( またはC ) とすれば 変換 : (3) [ ] = = b + + b は X よりY 上への 対 線形変換を表す ゆえに もY より X 上への 対 線形変換を表 す X 上に - ノルムを与え S = { X : = } ( X の単位球の表面 ) とすれば (4) 十分小さな δ > をとれば すべての S に対して δ Coyrght 再履修線形代数研究会 3

4 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム が真であることを示す 実際 (4) を否定すれば任意の自然数 に対して < / ( ) = を満たす S 上のベクトル列 { } がとれるはずである これより 他方 { } は有 界列だから 収束する部分列 { }( =,, ) を含む ( { } の第 成分の列から 収束列を抽出し それに対応する第 成分の列から収束列を抽出し を繰り返せばよい ) その極限を () と書けば ( =,, ) そこで ( ) ( ) ( ) () () () () = と書けば () = ( ) = ( ) b + + ( ) b () () ( ) () ( ) () () b b ( =,, ) ma{,, } ところが仮定により である ゆえに () = これは () = を意味する する と = = () () () ゆえに () ゆえに主張 (4) は真でなければならない つぎに (4) を成立させるδ > に対して (5) y < δ を満たすすべての Y y に対して = y = g ( y) + + g ( y) < これは矛盾である が成立することを示す かりに y < δ にもかかわらず y となるような y Y が存在 したとする y = とおけば = y かつ ゆえ y δ となる これは (4) と矛盾する ( / S ) ゆえに(5) = = < =δ は真でなければならない (5) は (6) すべての y Y に対して g ( y) + + g ( y) (/ δ ) y が成り立つことを意味する これは証明すべき () の後半に他ならない 4.3 有限次元ノルム空間内のコーシー列は収束する この節では前節の結果を応用し Coyrght 再履修線形代数研究会 4

5 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム 有限次元ノルム空間内のコーシー列は収束する ことを示す ここにノルム空間内の列 {,, } { } がコーシー列 Cauchy sequece であるとは どんな小さなε > を与えても自然数 ( ε ) を十分大きくとれば q, > ε ( ) を満たす すべての自然数, q に対して q < ε が成り立つことをいう 証明与えられた有限次元ノルム空間を X dm X = > { } を X 内のコーシー列とす ( ) ( ) る X の基底 { b,, b } を一つとり = b+ + b( =,, ) と書けば 任 意の自然数 q, に対して q = ( ) b + + ( ) b ( ) ( q ) ( ) ( q ) 前節の結果を適用すると適当な正定数 β に対して + + β ( ) ( q) ( ) ( q) q ( ) ( ) が成り立つはずである これは各成分の列 { },,{ } が実数または複素数のコーシー列 を表すことを示す ゆえにこれらの各列は収束し,, ( ) () ( ) () とすれば () () の展開形から b + + b となることは明らかである 4.4 有限次元ノルム空間内の有界列は収束する部分列を含むこの節では 4. 節の応用例として 列コンパクト性 sequetal comactess: 有限次元ノルム空間内の有界列は収束する部分列を含む を示す 証明問題の有限次元ノルム空間を X とし dm X = ( > ) とする いま {,, } { } を任意の有界列 すなわち 適当な正定数 α に対して α( =,, ) を満たすようなベクトル列 とする X の基底 { b,, b } を一つとり 任意の をこの基底 ( ) ( ) で展開したものを = b + + b ( =,, ) とする すると 展開係数の有界 性 (4. 節 ) より ( ) ( + + ) β βα ( =,, ) を満たす正定数 β が存在す る 各列 () ( ) { },,{ } は有界数列ゆえ,, ような自然数の部分列 ( ) () ( ) () < < がとれる ゆえに対応する { } ( =,, ) となる の部分列は Coyrght 再履修線形代数研究会 5

6 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム b + + b () () に収束する 4.5 有限次元ノルム空間上のノルムはどの つも同値である 4. 節の結果の応用として以下 (I)(II)(III) を証明する : (I) と ' を同じ 次元ノルム空間 X 上の 種のノルムとすれば すべての X に対し て α ' β が成り立つような正定数 α, β が存在する これをノルムの同値性 equvalece という (II) 与えられたベクトル列 { } に対して a が成立すれば a ' も成り立 つ すなわち X 内のベクトル列が一つのノルムに関して収束すれば 他のすべてのノルムに 関して同一の極限に収束する (III) いま { } を X 内のベクトル列 { b,, b } を X の任意の基底とすれば X 上の任 意かつ特定のノルム に関する収束 a と成分ごとの収束 comoetwse covergece ( ),, ( ) ( ) ( ) a aとは同値である ここに = b+ + b a= ab + + a b としている 証明 (I) { b,, b } を X の基底とし 任意の X をこの基底で展開したものを (*) = a+ + aとする 4. 節の結果を適用すれば すべての X に対して + + c を成立させるような正定数 c が存在する 一方 (*) の ' - ノルムをとれば ' ma{ a ',, a '}( + + ) c'( + + ) c' c β が出る 両ノル ムの役割を交換すれば適当な正定数 c '' に対して c'' ' がすべて X に対して成り立つ ことがわかる ゆえに / c'' = α とおけば α ' β がすべての X に対して成 り立つ (II) (I) から簡単に従う (III) 練習問題とする 例 R ( または C ) 上の,, - ノルム + + ( / + + ) ma Coyrght 再履修線形代数研究会 6

7 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム に対して次の不等式が成り立つ ( ここに = [ ],, ): 最初二つの不等式は簡単に証明できる 最後の不等式中の は明らか また j < j = ( ) より が従う 例 本節の結果 (III) により R ( または C ) 上の任意のノルムに関する収束 a ( ) ( ) と成分ごとの収束,, a aとは同値であることがわかる ここに [ ] ( ) ( ) =, a= a a としている 4.6 有限次元ノルム空間上の線形変換 4. 節からの応用として次の事実を示す : (I) ( 有界性 ) 次元ノルム空間 X から有限または無限次元ノルム空間 Y への線形変換 は有界である (II) 次元ノルム空間 X から有限または無限次元ノルム空間 Y への任意の線形変換 に対し て = かつsu{ : = } = を満たす X が存在する ゆえに su{ : = } = ma{ : = } と書いてよい 証明 (I) {,, } を X の基底とし X b b をこの基底によって展開し = b + + b とすれば 4. 節により + + β を満たす正定数 β が存在し = b + + b γ( ++ ) γβ ( γ ma{ b,, b } ) { : = } S は実数の有界集合を表す ゆえに実数の完備性により (II) (I) により集合 su S α は確かに存在する すると 各自然数 に対して α (/ ) < α =を満たすベクトル列 {,, } がとれる {,, } は有限次元空間 X 内の有界列を 表すから 4.5 節の結果により収束する部分列 { } ( = < < ) をもつ 極限を かつ とす Coyrght 再履修線形代数研究会 7

8 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム = ( =,, ) ゆえ = である しかも れば 同じ部分列に対して ゆえに 他方 α は明らか 極限 の一意性より = α = su{ : = } 4.7 演算子ノルム いま を有限または無限次元ノルム空間 X から同じ体上の有限または無限次元ノルム空 間 Y への有界線形変換とする 次の量 を の演算子ノルム oerator orm という : () = su = su = f{ K : K } = X ここに最後の量はすべての に対して K が真であるような K の値全体の下限 を表す () 中の等号成立の証明は練習問題とする ゆえに とは X が単位球の表面上 ( または単位閉球内 ) をくまなく動いたときの の値の上限を表す とくに X が有限次元ならsu 記号は ma 記号で置換可能であることは前節において証明済みである () 式が実際にノルムを定義していることは ( すなわち ノルムの公理を満たす ) ことは次節において示す 例 (I) R ( または C ) 上のノルム / + + ( + + ) ma ( [ ] = R ) / は使いやすいノルムとして実務計算上重要であるが ( + + ) ( ) もノルムを表すことが知られている ( 腕試し問題 参照) 例を挙げると = [ ] なら = + = / = ( + ) = = ma{, } = 参考のため R における原点を中心とする単位閉球の表面の図を付す : Coyrght 再履修線形代数研究会 8

9 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム y ( 上図の説明 ) 平面は R を現す 一番外側の正方形は = ma{, y} = y のグラフの 中間の円は y / = ( + y ) = のグラフ 一番内側の菱形は y = + y = のグラ フを表す,, m -ベクトルノルム対応する A = a j C の演算子ノルムは次式によって与えら れる : (a) (b) (c) A su A = ma aj ( 最大列和ノルム ) = j = su = ma a j = m j = m A A ( 最大行和ノルム ) A su A = σ ( A ) = * = ma AA の最大固有値の平方根 ( 最大特異値ノルム ) 証明 (a) m m j j a とする 簡単な計算で = = α ma a = A m = a α j j = j= が出る = C ) これより A α とくに = e = 第 単位ベクトルをとれば ( [ ] e = かつ Ae = α が成り立つ ゆえに α A Coyrght 再履修線形代数研究会 9

10 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム (b) j j m j a とする すると簡単に = j = β ma a = A = ma a j j β が従う m j = = C ) これより A β とくに として = a / a ( a ) ( [ ] j j j j j = ( a = ) によって定義されるベクトルをとれば かつ A β が成り立つ j = = ゆえに A β (c) レッスン から A の特異値分解を UAV * = Σ ( U : m 次ユニタリ行列 V : 次ユ ニタリ行列 σ Σ= σ ( σ σ ) ) とすれば Av σ ( v = ) かつ Av = σ ( v = ) となっている ここ v は V の第 列を表す 一言注意すると m A R の場合 X = Y = m R, R としも (a)(b)(c) はそのままの形 で成り立つ すなわち 実行列の演算子ノルムの値はそれを実空間間の変換と考えても 複素 空間間の変換と見なしても同じ値となる これはベクトルノルム が絶対ノル ム absolute orm を表す すなわち = ( ここに は の各成分をその絶対値で置き換 えたものを表す ) を満たすことに起因する また = [ ],, 自体を X = C より Y = C への変換とみなせば 演算子ノルム,,, は それぞれ ベクトルノ ルム と全く同一となる 例 与えられた自然数 m, に対して R m またはC m は m 次元ベクトル空間を作る 明ら { B,, Bq, Bm} は基底の一例である ただし B q は ( B q ) j = ( (, j) = (, q) の かに とき ) ( B q ) j = ( それ以外のとき ) とする そして A = a j m, m, / () aj, ( ), ma F aj a ( ) j, j, j=, j= A A A に対して Coyrght 再履修線形代数研究会

11 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム はすべて R m またはC m 上のノルムを表す とくに F はフロベニウス ノルム Frobeus orm と呼ばれ 実務計算によく使われる これはベクトルノルムの単純な拡張に過ぎない こ れらはすべて同値であるから (4.5 節 ) 与えられた行列の列{ A } の任意かつ特定の行列ノ ルム に関する収束 A A と成分ごとの収束 ( ) a a ( =,, m, j =,, ) と は同値である j j 4.8 演算子ノルムの性質 演算子ノルムは応用上多用される この節では演算子ノルムの重要な性質を学ぶ (I) AB, を有限次元ノルム空間 X からノルム空間 Y への線形変換とし 与えられたベクトルノ ルム に対応する演算子ノルムを同じ記号 で表せば 次の関係が成り立つ : () A ; A = A= () ca = c A ( c は任意のスカラーを表す ) (3) A+ B A + B (4) A A ( X ) ()()(3) は演算子ノルムが実際にノルムの公理を満たすことを示す (II) X, Y, Z をノルム空間とし X, Y は有限次元とする B: X Y, A: Y Z を線形変換 とすれば (5) AB A B (III) A を有限次元ノルム空間 X からそれ自身への線形変換とすれば A の任意の固有値 λ に 対して λ A が成り立つ ここに A は X 上の任意のベクトルノルムに対応する演算子ノル ムを表す また A の固有値とは A の ( 任意の ) 行列表現の固有値を表す ( これは行列表現によ らない量を表す ) (IV) B を有限次元ノルム空間 X からそれ自身への線形変換とする このとき B < なら ( I B) が存在し ( ) ( ) I B B および ( + ) /( ) I B I B B が成り立つ Coyrght 再履修線形代数研究会

12 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム ここに は演算子ノルムを表す 証明 :(I) () 定義より A は明らか A= なら当然 A = である 逆に A = なら 演算子ノルムの定義より すべての に対して A ( / ) = が成り立つはずである こ れはベクトルノルムの性質より すべての を表す () これも演算子ノルムの定義から直ちに従う に対して A = であること すなわち A= (3) 任意の X ( = ) に対して ( A+ B) = A+ B A + B A + B ここで = を満たすすべての について su をとれば A+ B A + B が得られる (4) = の場合は問題の不等式は明らかに成り立つ なら / = ゆえ A ( / ) Aが成り立つはずである これより A A が出る (II) 任意の X ( = ) に対して (I)(4) により AB = A( B) A B A B ここで = を満たすすべての に対して su をとれば AB A B が得られる (III) ( 必ずしも行列とは限らない ) 線形変換の固有値の復習から始める dm X = > とし X の任意の基底 { b,, b } をとり 任意の X をこの基底で展開したものを = b + + b A = Ab + + Ab となる ここで Ab,, Ab X ゆ と書けば え これらも同じ基底で展開し Ab j = ajb ( j =,, ) と書けば j j j j j= = = j= = A = a b = ( a ) b となる 以上の計算は形式的な行列積の形に書ける : () = [ b b ] A = [ b b ] a a a a [ ] b b A{ b} { b} A の固有値とは 特定の基底 { b,, b } に関する行列表現 A{ } の固有値と定義する このい い方が許されるためには これらの固有値が基底 す必要がある 実際 他の任意の基底 { b ',, b '} { b,, b } b の選び方に無関係であることを示 をとれば これらの基底の間には適当な可 Coyrght 再履修線形代数研究会

13 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム b b b b V ( これも形式的な行列積形 ) なる関係がある こ 逆行列 V を介して [ ' '] = [ ] れより, A の 基底 { b ',, b '} () = [ b b ] V = [ b ' b '] ' ' に関する行列表現は ' ' a a ' A = [ b ' b ']( V ) V a a ' [ ' '] b b A{ b'} { b'} これより 二つの行列表現は互いに相似の関係にあることがわかる : () A = V A V { b'} { b} 従って二つの行列表現は全く同一の固有値を共有する さて 行列の固有値 - 固有ベクトルの関係から A の任意の固有値 λ に対して A = λ を満たす固有ベクトル ( X ) が存在することがわかる λ に対応する固有ベクトル を =を満たすようにとれば ( これは常に可能 ) A A = λ = λ = λ が出る (IV) B < なら B の任意の固有値 λ は前項により λ B < を満たす I B の固有値は その行列表現からかならず λ の形であるから とはならない ゆえに ( I B) は確かに存 在する ( I B)( I B) = Iより ( ) ( ) ( ) I B = I+ B I B だから ノルムをとれば ( ) + ( ) + ( ) I B I B I B B I B が得られる これを ( I B) について解 けば I B B が出る また ( ) 式を ( I B) I= B( I B) と変形し ノル ( ) /( ) ムをとれば ( ) /( ) I B I B B が得られる 固有値を使わない証明法についてはレ ッスン末 腕試し問題 参照 4.9 演算子ノルムの応用例 ノルムは データに一定の変動を与えたとき答えはどう変動するか を定量的に評価する ための尺度として多用される 例えば A を可逆行列とするとき B がどの程度小さければ A+ B の可逆性が保証されるか A+ B の固有値は A の固有値はどれくらい近くにあるか 行列方程式 A = b( ただし A は可逆行列 ) の解はデータ Ab, が変動するとどう変動するか Coyrght 再履修線形代数研究会 3

14 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム などの評価にノルムが使われる このような解析は摂動論 erturbato theory 安定性解析 stablty aalyss 感度解析 sestvty aalyss などと呼ばれている 本節では前節までの結果を利用し このようなノルムの応用例を学ぶ 例, AB R ( または 逆行列を表し C ) とし A は可逆行列とする A B < なら A+ B も可 () ( A+ B) A /( A B ) () ( A+ B) A B A /( A B ) が成り立つ ここに は任意の演算子ノルムを表す 証明 A = I の場合は前節 (IV) で証明済みである これを応用すれば ()() が得られる 実際 A+ B= A( I+ A B) ゆえ ( A+ B) = ( I+ A B) A ( 仮定 A B < と前節 (IV) によ I+ A B が存在する ) ノルムをとり前節(III)(IV) を利用すれば () 式が出る り ( ) つぎに A = ( I+ A B)( A+ B) I= ( A+ B)( A+ B) = A( I+ A B)( A+ B) より A ( A+ B) = A B( A+ B) = ( + ) + ( + ) A B A B A B これより = + (III)(IV) を使えば () 式が得られる 例 バワー ファイクの定理 Bauer-Fe theorem ルダン標準形が対角行列 ) X AX かつ特定の固有値を μ とすれば A B( I A B) A が出る ノルムをとり 再び前節 λ = D λ A C を対角化可能な行列とし ( ジョ とする B C とし A+ B の任意 (3) m μ λ X X B が成り立つ ここに は対角行列のノルムが対角成分の絶対値最大値を与えるような演算子ノ ルム ( 例 :,, ) を表すものとする (3) 式の右辺は B のみならず X にも依存す ることに注意 数 cod( ) X X X は X の条件数 codto umber と呼ばれ 誤差解析に よく出てくる量である ( レッスン 特異値分解 で出てきている ) また (3) はエルミート 行列の固有値単調定理 エルミート行列 ABC,, が A+ B= Cを満たせば Coyrght 再履修線形代数研究会 4

15 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム α + β γ α + β( =,, ) ここにα α, β β, γ γ は,, ABC の固有値を表す ( レッスン 8 シュール分解 Part II) が固有値ごとの比較を行ってい るのにくらべると弱い結果である 証明 : 前節 (IV) を使うと簡単である まず μ λ ( =,, ) と仮定してよい すると = A+ B μi = X A+ B μi X= D μi + X BX ところが det( ) det ( ) det{( ) } μ λ ( =,, ) だから D μ I+ D μi X BX = I は可逆行列を表す ゆえに det( ( ) ) これを前節 (IV) に照らすと これは ( D μi) X BX を意味する 前節 (III) とノルムに関 する仮定を使うと ( D μi) X BX ( D μi) X B X = X X B m μ λ 分母を払えば証明すべき不等式が得られる 例 3 行列方程式の解の摂動問題 ( レッスン.7 節の結果の一般化 ), たはC ) b, b R ( または C A A R ( ま ) を与えられた行列とし A は可逆行列 与えられた ベクトルノルムと演算子ノルムを同じ記号 で表すものとし A A < とする すると ( ) も可逆行列となる ( A A A A < ) そこで方程式 A+ A= A I+ A A A = b と ( A+ A)( + ) = b+ bの解の差 を評価してみよう ここに, + R ( またはC ) である 辺々相引けば単純な計算で これまでにもよく使ったノルムの性質を使うと = ( I+ A A) A ( b A) が出る A A b ( ) ( ) + A A A b + A A A A A A A b A ( ここに A = b より A b ) すなわち (4) cod( A) A b ( + ) A A b cod( A) A ここに cod( A) = A A は例 にも出てきた条件数である この式を見ると cod( A) A / A が に比べて小さければ データの相対的変動の和のほぼ A の条件数倍が 解の相対的変動として表れ得る ことがわかる ゆえに 条件数が大きければ データの小さな変動が大きな解の相対的変動を起こし得ることになる このようなわけで 係数行列の条件 Coyrght 再履修線形代数研究会 5

16 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム 数が大きい方程式は悪条件である ll-codtoed と呼ばれる ただし次の例のように例もある ことに注意 : A = + の場合は A の -ノルム( =,, ) 条件数はいずれも だから 大きければいくらでも大きくなるが A ( + ) = b+ b( A= としている ) は厳密に解 けるから 上の誤差評価式は悲観的過ぎることになる 条件数の推定は数値計算上大事な話題であるがここでは 参考のため 行列の条件数の a b 例を挙げるにとどめよう : 実際 A = ( ad bc ) c d の - 条件数は簡単な計算によっ て ma{ a + b, c + d}ma{ b + d, a + c} cod( A) = A A = ad bc 4. ハーン バナハの定理 (I) 次にのべる線形汎関数拡大定理をハーン バナハの定理 Hah-Baach theorem という : 与えられた実または複素有限次元ノルム空間 X の部分空間 M 上で定義された ( 有界 ) 線形汎関数 f すなわち M から対応するスカラー体への線形変換 f はその演算子ノルム の値を不変に保ちつつ X 上で定義された線形汎関数 f にまで拡大できる すなわち M 上 で f = f かつ f = f を満たす X 上の線形汎関数 f が存在する ここに f = ma{ f ( ) : =, X} f = ma{ f( ) : =, M} ( 注意 ) この定理は X が無限次元であってもこのままの形で成立するが 特殊な論法 ( ゾーン の補題 Zor s lemma) を必要とするため ここでは X を有限次元としている このため X 上 の線形汎関数は自動的に有界となる (4.6 節 ) X が無限次元ならこの定理が意味をもつため には f の有界性を陽に仮定する必要がある 証明 以下の証明は G. F. Smmos, Itroducto to oology ad Moder Aalyss, McGraw-Hll, 963, 48, 6 9 をもとにしている 解析学からの予備知識は これまで 通り 実数の完備性のみである M = {} なら M 上の線形汎関数は f = のみだから f = とすればよい そこで M を {} でない真部分空間 f を M 上で定義された線形汎関数とし f = と仮定しておく Coyrght 再履修線形代数研究会 6

17 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム こう仮定しても一般性が失われないことは明らかである 証明の核心部は次の補題である : ( 補題 ) M とすれば f は M M + sa{ } 上で定義された f = を満たす線形 汎関数 f に拡大できる X は有限次元としているから 補題の手続きを有限回繰り返せば f を X 上全体まで拡大 できることは明らかである ゆえに 補題を証明すれば定理の証明が済むことになる ( X が無 限次元の場合はこの補題と先ほど言及した ゾーンの補題 が拡大手続きの手段となる ) ( 補題の証明 ) 以下の証明は X が有限次元でも無限次元でも成り立つ X が実ノルム空間であ る場合をまず扱い その結果を複素ノルム空間の場合に応用する (A) X が実ノルム空間の場合 M 内の任意ベクトルは一意的に +α (α は実数 M ) の形に表せる ゆえに f が f を M 上に拡大した線形汎関数を表すための必要十分 f ( + α ) = f ( ) + α f ( ) = f( ) + αr の形をもつことである ここに 条件は r f( ) ( 実数!) と書いている 残る問題は r の値をどう選べば f が満たされるか だけである ( すでに f = であるから f は f = を意味する ) さて f の形から ( ) f 任意の M 任意の実数 α に対して f( ) + αr + α 任意の M 任意の実数 α に対して + α f( ) + αr + α 任意の M 任意の実数 α に対して f( ) + α αr f( ) + + α 任意の M 任意の実数 α に対して f( ) + r f( ) + + α α α α ( 最後の同値性は α > α < の場合に分けて検算されよ ), M この最後の条件を満たす r が存在することを次に示す 実際 任意の f( ) f( ) = f( ) f( ) ( f = ) に対して = ( + ) ( + ) f( ) + f( ) + + ゆえに を固定しすべての M に対応する右辺の値の集合の下限をとれば f ( ) + f{ f( ) + + : M} なる この関係はすべての M に対 Coyrght 再履修線形代数研究会 7

18 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム して成り立つから 左辺の上限をとれば a su{ f( ) + : M} f{ f( ) + + : M} b が得られる ゆえに r を ることがわかる a r bを満たすように選べば同値関係 ( ) の最後の条件が成立す (B) X が複素ノルム空間の場合 f は複素数値をとる関数であるから 実部と虚部を gh, と し f ( ) = g( ) + h( ) と書くことにする λ, μ を任意の実数とすれば 任意の uv, M 対して g( λu+ μv) + h( λu+ μv) = f( λu+ μv) = λf( u) + μf( v ) = λg( u) + μg( v) + { λh( u) + μh( v )} 実部を等置すれば (, λ μ は実数ゆえ ) g( λu+ μv) = λg( u) + μg( v) が得られる すなわち スカラーを実数に限定すれば g はノルム空間 X 上で定義された実線形汎関数と見なせる そ こで (A) の結果を利用して g を M 上に拡大し これを g ( 実線形汎関数 ) と呼ぼう ただし に g = g g( ) f( ) f = ゆえ 確かに g 従って g つぎに g ( ) + h ( ) = f( ) = f( ) = g { ( ) + h( )} = g( ) h( ) から h( ) = g( ) が出る 従って f ( ) = g( ) g( ) そこで f( ') = g( ') g( ') ( ' M) によって f を定義すると これが f = を満たしつつ f を M 上に拡大した線形汎関数となっていることを示そう 実際 直前の計算から f と f の 値は M 上で一致していることがわかる しかも g の性質から 任意の ', '' M に対して f ( ' + '') = f ( ') + f ( '') f ( α') = α f ( ') 任意の実数 α に対して が成り立つことがわかる そして後者の関係はたとえ α が複素数であっても成り立つことは次 の計算で検証できる : まず f( ') = g( ') g( ') = { g( ') g( ')} = f( ') ( g( ') = g( ')) ゆえに 任意の実数 λ, μ に対して f(( λ + μ) ') = f( λ') + f( μ') = λf( ') + f( μ ') = λ f( ') + μf( ') = ( λ+ μ) f( ') 以上により f は f を M 上に拡大した線形汎関数であることがわかる ' M 残るは f = を示すのみである これを示すには ' = を満たすすべての に 対して f ( ') が成立することを示せば十分である ( f は f を M 上に拡大した線形汎 関数であるから f f = はすでにわかっている ) 実際 ' M ' = f ( ') = re θ Coyrght 再履修線形代数研究会 8

19 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム ( f ( e θ ') ( 極表示 ) とすれば f ( ') = r = e θ f ( ') = f ( e θ ') = g ( e θ ') 数!) ゆえ f ( ') = g ( e ') g e ' e ' = ' = ( この計算から θ θ θ をとり f のノルム値はその実部のノルム値に等しい ことがわかる!) は実 f はすべての実数値 以上の証明をみると 各ステップは単純だが 構成は巧妙である 複素汎関数の扱いに複素数の性質が巧みに使われていて 彼らのいう slc roof の好例といえる 証明の複雑さはこの定理のパワーを暗示している ハーン バナハの定理の応用は以下 (II)(III)) の形をとることが多い (II) X を有限次元ノルム空間とすれば 任意の X に対して f ( ) = かつ f = を満たす線形汎関数 f が存在する 証明 (I) において M = sa{ } ととり f ( α) = α によって線形汎関数 f を定義すれ ば 明らかに f ( ) = かつ f = が成り立つ すると ハーン バナハの定理により f は要求される性質をもつ X 上の線形汎関数に拡大できる ( 注意 ) X を有限次元としているから f = を満たす任意の線形汎関数 f に対して f ( ) = を満たす が存在することはわかっている (4.6 節 ) (II) はこの双対問題 も真であることをいっている (II) は実は X が無限次元であっても成り立つ (III) X = R または C の場合与えられた に対して a = かつ a = を 満たすa X が存在する ここに a = ma{ a : = } 証明 (II) により f ( ) = かつ f = を満たす X 上の線形汎関数 f が存在する そして [ f f ] a= ( e ) ( e ) ( e, は単位ベクトル ) とすれば すべての X に対して ( ) f = a となる 例 (III) における X = R または C 上のノルムを -ノルム( ) とする b X を与えられたベクトルとし ( をここではb と書いている ) (III) で存在を保証されている Coyrght 再履修線形代数研究会 9

20 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム ab= b, a = を満たす a = [ a a ] X を示そう (a) = の場合 a = ma a = a は既知である ab= a b = a b を満たす a の例は a = b / b, a = b / bである ( ただし b = なら a = とおく ) 検算 : ab= ab + + ab = ( b/ b) b+ + ( b/ b) b= b + + b= b a = ma{,,} = (b) =の場合 a = a + + a = a は既知である ab = a b = a b を 満たすa の例 : b = ma b = b となる を一つとり a = ( ), a = b / b とする 検算 : ab= ( b / b ) b = b = b そして a = = (c) < <の場合 腕試し問題 中で示すように <, q / + / q = を満たす正数, q に対して q q / q ma{ : = } = = ( a + + a ) q a a a が成り立つ ab= a b q = a b ( q = + q ) を満たすa の例 : a = b / b ( =,, ただし b = なら a = とおく ) とすればよい 検算 : ab = ( b / b) b = b = b / ( ) / / / ( ( q q / ) ) ( q q ) ( q ) q = = b q b = b = b = = a a b b ゆえに ab= b = b b = a b 4. ハーン バナハの定理の応用例 本節の内容は W. Kaha, Numercal lear algebra, Caada Mathematcal Bullet 9 (966), 中の定理 (775 ページ ) をもとにしている つぎの事実 (I)(II) が成り立つ : Coyrght 再履修線形代数研究会

21 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム (I) m A = a j C を与えられた行列とし m C C 上に任意のノルムを与える す m ると = y = を満たす適当な C y C をとれば A = y A が 成り立つ ここに y A は与えられたベクトルノルムに対応する演算子ノルムを表す : y = y y y C y = ma{ : m, } A = ma{ A : C, = } 証明まず A = A を成立させる C = を適当にとる (4.6 節により可能 ) m ハーン バナハの定理により y A = A かつ y = を満たす y C がとれる A = A ゆえ 結局 y A = A となる (II) 与えられた可逆行列 A = a j C から非可逆行列までの最短距離は / A に等しい すなわち A < / A なら A+ A は必ず可逆行列であり A = / A を満たす A の 中に A+ Aを非可逆行列とするようなものが存在する ここに行列ノルム A A は与 えられたベクトルノルムに対応する演算子ノルムを表す 別の述べ方をすれば A/ A から非 可逆行列までの最短距離は/( ) = / cod( ) A A A によって与えられる cod( A) は A の条件数と呼ばれることはすでに述べた ベクトルノルムの選び方に無関係に A AA = I = が成り立つことに改めて注意 cod( ) ( 注意 )- ノルムに限定した場合はレッスン 特異値分解 で証明済みである 証明 A+ Aを非可逆行列とすれば ( A+ A) = を満たす C がとれる する と A = A A = A A A A A A これより A / A ( A ) が得られる 対偶をとれば問題の主張の前半が証明される つぎに (I) を A に適用すれば ( m= ) A = y A ただし y = = Coyrght 再履修線形代数研究会

22 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム を満たす, y C がとれるはずである A= y A / とおけば A = / A かつ A+ Aが非可逆行列となることを示そう 実際 A = y A ゆえ ( A+ A) A = ( A y / A ) A = ( y A )/ A = = が得られる A ゆえ これは A+ Aが非可逆行列であることを示す 前半の結果から A / A のはずである 他方 A = y / A y / A = / A 結局 A = / A が結論される 例 A = 3 4 より - 演算子ノルムに関して最短距離にある非可逆行列 4 B= A y / A を構築してみよう ( 証明参照 ) A = 3 ゆえに A = 3( 最 大行和ノルム ) としては A A および = を満たすベクトルをとればよい そのような = の一例は = [ ] である そして = [ ] A 3 となる y としては y A = A および y = を満たすベクトルをとればよい そのようなベクトルの一 例は y = [ ] である ゆえに B= A y / A = [ ]/3= = 3 8 また A より B までの距離は証明中で示したように/ = /3 A であるはずである 検算 : B A y / A = = = 3 3 A A の条件数 cod( A) を計算すると cod( ) = = 7 3= A A A 最後にひとこと : このレッスンの最も基礎的な結果は 線形変換の連続性 ( 有界性 ) 展開 係数の有界性 有限次元ノルム空間の完備性 ノルムの同値性 演算子ノルムと性質 ハ Coyrght 再履修線形代数研究会

23 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム ーン バナハの定理 である 証明のために解析学からの既知事項として仮定したのは実数 複素数の完備性だけである ハーン バナハの定理は一般のノルムに関して 与えられた可逆行列から非可逆行列までの最短距離 を算出するために必要であった この定理は双対問題の解を保証することを思えばその重要性を理解できよう 腕試し問題 問題 4. 与えられたベクトル空間 X 上のノルム ' の同値性 α ' β ( X α, β は正定数 ) を記号で ~ ' で表すことにすれば "~" は同値関係を表すことを示せ ( 略解 ( 回帰性 ) X 上の任意のノルム に対して ~ ( 対称性 ) ' に対して ~ ' なら '~ ( 推移性 ) ' '' に対して ~ ' かつ '~ '' なら ~ '' が成り立つことを示せばよい ) 問題 4.(4.8 節 (IV) の別証明 ) X を有限次元ノルム空間とし A を X からそれ自身への 線形変換とする 演算子ノルム A に対して 不等式 A < が満たされれば ( I A) が存 在し ( ) /( ) I A A が成り立つことを次の手順に従って示せ () 自然数 に対して B = I+ A+ + A とすれば ( I A) B = B ( I A) = I A () { B } はコーシー列を表すことを示せ ゆえに 4.3 節により { B } B とする : B B は収束する その極限を (3) () において の極限をとり A A を利用して ( I A) B = B( I A) = I すなわち B= ( I A ) を示せ (4) A = = B I A A A A A Coyrght 再履修線形代数研究会 3

24 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム および B B B B を利用して ( ) /( ) I A A を導け ( 略証 :() 単純な計算で出る () q q q = + + ( + + ) B B A A A A A ( ただし q ) これに 与えられた条件 A < を適用すればよい (3) 問題文中の指針に従えばよい ) 問題 4.3 (4.4 節の逆 ) 無限次元ノルム空間の単位球表面上から収束する部分列を全くもた ないベクトル列を選ぶことができる ことを次の手順によって示せ X を与えられた無限次元 ノルム空間とする (I) X を任意の 次元部分空間 ( =,, ) { b,, b } を X の基底とし X 外より任意のベクトルb を一つ選び X sa{ b, +, b, b } とすれば X + は + 次元部分空間を 表す このとき d dst( b, X) f{ b : X} ( b より X までの最短距離 ) = X () b を満たす () が存在することを示せ (II) ( )/ = ( )/d X + かつ () () () () y b b b を定義すれば 明らかに y () = X () y すると 任意の に対して () y が成り立つことを示せ (III) 単位球表面 S = { X : = } 上から 収束する部分列を全く含まない列を抽出できる ことを示せ 以上と 4.5 節の結果を総合すれば すべての有界列が収束する部分列を含むための必要十分条件は その空間が有限次元であることである ( 注意 ) (I)(II) は A. E. aylor ad D. C. Lay, Itroducto to Fuctoal Aalyss, Secod Edto, Kreger, 98, heorem 3.5( Resz s Lemma ),. 64 を多少変更したものである aylor ad Lay では X に相当する部分空間を 有限または無限閉部分空間 X としているた め 結論は < θ < を満たす各実数 θ に対して dst y X θ を満たす y () X () (, ) が存在 する こととなっている (III) は aylor ad Lay herem 3.6 (.65) と本質的に同じである ( 略証 (I) 各自然数 に対して (*) d b < d + (/ ) を満たす X を選べば { } は明らかに X 内の有界列を表すから 4.5 節により X 内のベクトル () に収束する部分列 Coyrght 再履修線形代数研究会 4

25 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム () を含む これを (*) 式に使えば d = b > を得る (II) 次を検算せよ : () () () () y = (/ d)( b ) = (/ d) d+ b (/ d) d = ( d+ X ) (III) b S を任意に取り sa{ b} X とする (II) により b S を適当にとれば dst( b, X ) = となる sa{ b, b} X とする 以下同様の手続きにより S 上のベクト { b, b, } { b } が得られる ここに ル列 + X = dst( b, ) ( =,, ) これより任意の自然数 m に対して b b が成 m り立つ このような S 上の列は収束する部分列を全く含まないことは明らか ) 問題 4.4 ( 行列ノルム計算問題 ) 次の行列の -ノルムと -ノルムを計算せよ: P= P P Q= Q Q * * [ ],,,,, ( 答 : P P P = = * = P P P 3 = = * = + Q Q Q = = * = 4 Q Q Q ) = = * = + 問題 4.5 ( 行列ノルムに関する問題 ) m A C B を A の任意小行列 ( ブロック ) とすれ ば B A が成り立つ ここに =,, のいずれかとする ( 略証 : -ノルムの定義より 任意の列ベクトル y の数個の成分を で置き換えたものを y ' で 表せば 明らかに y' y が成り立つ そこで B を A より第,, 行 第 j,, jq 削除して得られる小行列とすれば B ma{ A : = & =, j = j,, j } A ) j q B 問題 4.6 A = C または B A = なら A = ma{ B, C } が成り立つことを C 示せ ここに =,, のいずれかとする 列を ( 略証 : A = 最大列和 A = 最大行和 A = 最大特異値 = AA の最大固有値の平方 根 ここに BB CC AA= C C または B B を考慮すれば直ちに出る ) Coyrght 再履修線形代数研究会 5

26 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム 問題 4.7 c R, B R とすれば ( cc ) / BI = B Bc cc F cc F ( 注意 ) I cc ( 略証一般に cc は平面 = c 上への正射影を表す F はフロベニウス ノルムを表す X ( ) = ( tr( X X )) ここに tr( ) はトレースを表す 一般に F, j= / / j tr( P+ Q) = tr( P) + tr( Q), tr( cp) = c tr( P), tr( PQ) = tr( QP) が成り立つ ( ここに PQ, は正方行列 c はスカラー ) 以上から cc cc cc cc ( ) = tr(( ) ( )) = tr( ( ) ) cc cc cc cc F BI I BBI BBI ( ) ( ) = tr( BBI ( cc )) = tr( BB) tr( B cc B) = tr( BB) tr( Bc Bc ) cc cc cc F = B Bc cc ここに Bc は ゆえ Bc = Bc ) F 問題 4.8 任意の m a R, b R に対して ab ab a b および = = F ab a b が成立することを示せ m ( 注意 ) すでに学んだように =,, a R ( またはC m ) に対して a を演算子ノルムと見なしてもベクトルノルムと見なしても 値は一致する ( 略証まず ab を計算すると ( 前問略解参照 ) F = tr(( ) ) = tr( ) = tr( ) = = F ab ab ab aa bb aa bb aabb a b ab = ( ) ( ) ab ab の最大固有値 = ( aa )( 値は bb,,,であることより bb の最大固有値 ) ここでbb の固有 bb の最大固有値はbb である ゆえに ab = ( )( ) aa bb= a b = ab F ab a b = 後半の証明 : 演算子ノルムの性質から a b ) 問題 4.9 m a R, b R とするとき Xa = b を満たすすべての m X R のうち Coyrght 再履修線形代数研究会 6

27 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム X= X ba / a aは最小の -ノルム / b a をもつことを示せ ( を単に と書く ) ( 略証 Xa = b を満たす任意の m X R をとる ノルムをとれば X a b が得られる こ れより X b / a 一方 X = ba / a a R m は確かに Xa = b を満たし 前問の結果 を利用すれば X = ba / a a = b a / a = b / a が得られる ) 問題 4. 任意の m A R に対して A F ra( A) A A F が成立することを示せ ( 略証 A = trace( A A ) A = A の最大特異値 を使う A の特異値分解を F A = UΣV とする ここに, UVは直交行列 Σ = dag{ σ,, σ r,,,}( m 行列 ) ( σ σ r >, r ra( A)) σ,, σ r,,,は A の特異値を表す 簡単な計算によっ て = ( ) = + + A trace A A σ σ F r また A = σ はすでに知っている 問題の不 等式はこれらから得られる ) R ( またはC a= a a 問題 4. ( ベクトル -ノルム) この問題では [ ] ) の - ノルム / a = ( a + + a ) が実際にノルムの公理を満たすことを示す =,, の場合はすでにわかっているから > の場合のみを考える 三角不等式以外の要 請が満足されることは単純な計算で確認できるから この問題では三角不等式 ( 別名 ) a+ b a + b のみの証明を考える さらに ノルムの単調性および極限定理 a a ( ) の証明を行う (I) ab, なら a a b b q + q ( ただし q, >, + q = ) を示せ これは 相乗平均 / / 相加平均 すなわち ab ( a + b)/ の一般化を表す (II) ヘルダーの不等式 Hoelder s equalty: ab a b を示せ q = Coyrght 再履修線形代数研究会 7

28 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム ( a= [ a a ], b= [ b b ] ) = q = はコーシー シュワルツの不等式に還元する (III) 三角不等式 ( ミンコフスキーの不等式 Mows s equalty): a+ b a + b (IV) 単調性 : s t が成り立つことを示せ < なら 任意の [ ] (V) a a = ma a ( ) を示せ a= a a R ( またはC ) に対して a t a s ( 略証 (I) ab>, としてよい < < とし 関数 f( ) = + ( > ) を定義す ると のとき f ( ) ( f () = なら f '( ) = ( ) ) 直線 y = + は点 (,) における曲線 y = への接線となっていることに注意 つぎに a b なら = a/ b( ), = / とし a< bなら = b/ a, = / qとせよ q (II) ab, としてよい a ' = ( a / a ), b ' = ( b / b ) とおいて (I) の結果を利用すれば q t s ab a b + a b a b q q q (III) a+ b としてよい q q a+ b = a + b = a + b a + b が出る について和をとれば問題の不等式が出る = = a s + b s ( s a + b) = = ( ) q / q ( ) q / q s s = = a a + b + b a + b = = a ( ) + b ( ) ( ヘルダーの不等式による ) / q a b s ( ( ) = ( + )( ) ( a + b ) a+ b (/ q) q = ) ( a + b ) a+ b ( = / q ) 両辺を a+ b で割れば問題の不等式が出る (IV) a,, a を与えられた正数とすれば 関数 / f( ) ( a + + a ) は > で単調減少 関数を表すこと すなわち > なら f '( ) < を示せば十分である 実際 log f ( ) = (/ ) log( a + + a ) ( log は自然対数を表す ) を微分し 整理すれば Coyrght 再履修線形代数研究会 8

29 再履修線形代数 分解定理を主軸に整理整頓レッスン 4 ノルム ( f '( )/ f( )) = ulogu+ + ulogu られる ただしu /( a a + + a ), =,, が得 > であるから < u,, u < ゆえに右辺は負の実数を表す ゆえに f '( ) < + ここにb, ε ε ( ) は自然数値 / (V) ( 補題 ) ( b ε ) ( ) / のみを採るものとする 実際 ( b+ ε ) = とおけば ゆえに 項定理により + ( + ) = b + ε これより ( b + ε )/ ゆえに すなわち / ( ) b + ε ( 補題証了 ) つぎに a は a / = ( b+ ε ) ma a (b は高々 の自然数 ε は のとき ε となるような量 ) の形に書けることを確かめよ すると 補題 と (IV) から a ma a = a が得られる ) = a a R ( または 問題 4. ( 行ベクトルの - 演算子ノルム ) 与えられた [ ] の ノルムは 既知の行列演算子ノルムに関する結果を利用すれば a C ) a = ma a = a a = ma a = a a = a a = a によって与えられる 一般に q>, (/ ) + (/ q) = なら a a が成立することを示せ ここに a の - q 演算子ノルムは a = ma{ a : = } によって定義される ( 略証前問で証明したヘルダーの不等式を使う 実際 [ ] = = = を満たす任 意の に対して a = a a a ( ヘルダーの不等式による ) q これより q q = a a a が得られる ここで = a /( a a )( a ) = ( a = ) とすれ q q ば + q = q により = a = a が成り立つ 上の結果と合わせると a a q q が出る ) = q Coyrght 再履修線形代数研究会 9