• 検索結果がありません。

行列 固有値問題

N/A
N/A
Protected

Academic year: 2021

シェア "行列 固有値問題"

Copied!
26
0
0

読み込み中.... (全文を見る)

全文

(1)

数値計算

大阪大学基礎工学部 永原正章

年 月 日 限

(2)

休講 補講 知

本日 限目 補講

来週 月 日 木 休講 月 日 木 限目 補講

教室 限目 限目

(3)

行列 固有値問題

固有値

行列

R

× 固有値

λ

対応 固有

= λ

固有値問題

固有値 存在 領域 求 固有値 求

(4)

固有値 存在 領域

正方行列

ρ( ) ≤ ∥ ∥

成 立

ρ( ) = | λ ( ) |

半径

∥ ∥ =

| |

| |

行列

行列 固有値 全 原点 中心 半径

∥ ∥

円 内部 境界 存在

(5)

固有値 存在 領域

ρ( ) ≤ ∥ ∥

∥ ∥

固有値

(6)

固有値 存在 領域

半径

∥ ∥

円 広

精密 簡単 領域 求

定理

(7)

定理 ⋆ ⋆ ⋆

行列

= [ ] C

×

σ( ) =

=1

, :=

{

C : | − | ≤

=1̸=

| | }

中心

半径

|

1

| + · · · + |

(1)

| + |

( +1)

| + · · · + | |

円 内部 境界

(8)

定理 ⋆ ⋆ ⋆

行列

= [ ] C

×

σ( ) =

=1

, :=

{

C : | − | ≤

=1̸=

| | }

半径中心

|

1

| + · · · + |

(1)

| + |

( +1)

| + · · · + | |

円 内部 境界

(9)

定理 ⋆ ⋆ ⋆

定理 行列

C

×

σ( )

成 立 行列 固有値 領域

(10)

例題 ⋆ ⋆ ⋆

行列

=

[ 0 1

2 3 ]

行列 対 円 求 固有値 存在領域 図示

二 円 1 2

1

= { ∈ C : | | ≤ 1 } ,

2

= { ∈ C : | − 3 | ≤ 2 }

二 円 1 2

1

= { ∈ C : | | ≤ 2 } ,

2

= { ∈ C : | − 3 | ≤ 1 } ,

(11)

例題 ⋆ ⋆ ⋆

= (

1

2

) (

1

2

) λ = 1, 2

0 1 2 3 4 5

2 1

1 2

2 1

2 1 0 1

(12)

定理 証明

複素数

λ

固有値

λ σ( )

λ =

=1

{

C : | − | ≤

=1̸=

| | }

行列 固有値

λ

対角成分 11

,

22

, . . . ,

等 場合 明

λ

λ

対角成分 等

λ ̸ = = 1, . . . ,

(13)

定理 証明

行列 対角成分 非対角成分 次 分解

= +

=

 

 

11

0 . . . 0

0

22

0

0 . . . 0

 

 

, =

 

 

0

12

. . .

1

21

0

1, 1

. . .

, 1

0

 

 

.

行列 λ

= λ = ( λ ) +

複素数

λ

固有値 λ

= λ

正則

(14)

定理 証明

λ

= λ = ( λ ) +

正則

0 C

存在

λ

= 0

{ ( λ ) + } = 0 λ ̸ = = 1, . . . , λ

正則

= ( λ )

1

(15)

定理 証明

= ( λ )

1

∥ ∥

= ( λ )

1

≤ ∥ ( λ )

1

∥ ∥

.

= [

1

, . . . , ]

R

∥ ∥

=

=1,...,

| |

行列

= [ ] R

×

∥ ∥

=

∈R , ̸=0

∥ ∥

=

=1,...,

=1

| |

(16)

定理 証明

不等式

∥ ∥

≤ ∥ ( λ )

1

∥ ∥

1 ≤ ∥ ( λ )

1

=

=1,...,

=1̸=

1

| − λ | | |

番号

∈ { 1, 2, . . . , }

存在

1 1

| − λ |

=1

̸

=

| |,| − λ| ≤

=1

̸

=

| |

λ

λ

(17)

定理 証明

固有値 固有値 等

定理 使 定理

成 立

(18)

練習問題

定理 用 次 行列 固有値 存在 領域 図示

= [ 1 2

3 4

]

(19)

解答例

= [ 1 2

3 4 ]

行列 対 円 求 固有値 存在領域 図示

二 円 1 2

1

= { ∈ C : | − 1 | ≤ 2 } ,

2

= { ∈ C : | − 4 | ≤ 3 }

二 円 1 2

1

= { ∈ C : | − 1 | ≤ 3 } ,

2

= { ∈ C : | − 4 | ≤ 2 } ,

(20)

解答例

= (

1

2

) (

1

2

)

λ

±

= 5 ± 33

2 5.372, 0.372

(21)

固有値 数値計算

行列

C

× 固有値 全 求

以上 場合 一般 厳密解 求 不可能 現実的 問題 解 場合 固有値 求 数値計算 必要

固有値 求 分解法

(22)

分解

行列

C

× 正則 次 分解 一意 的 可能

=

行列

=

=

成 立 行列 対角成分 正 右上三角行列

(23)

分解法

正則行列 固有値 求 分解法

1

:=

= 1, 2, . . .

:=

分解

+1

:=

掛 算 入 替 分解法 行列 固有値 近似値 得

(24)

固有値 分解法

行列

C

× 正則行列 固有値

λ

1

, λ

2

, . . . , λ

相異

1

| >

2

| > · · · > | > 0

満足 1

=

分解法

作 行列

→ ∞

次 行列 収 束

=

 

 

λ

1

. . . 0 λ

2

0 . . . 0 λ

 

 

.

(25)

実験

A=rand(5,5); A=A'*A; eig=gsort(spec(A));

An=A;

N=20;

err=[];

for n=1:N

[Q,R]=qr(An); //QR��

An=R*Q; //��������

err=[err,norm(eig-gsort(diag(An)))];//����

end

figure;plot(log10(err))//�������

(26)

実験

参照

関連したドキュメント

[r]

[r]

I Samuel Fiorini, Serge Massar, Sebastian Pokutta, Hans Raj Tiwary, Ronald de Wolf: Exponential Lower Bounds for Polytopes in Combinatorial Optimization. Gerards: Compact systems for

難病対策は、特定疾患の問題、小児慢性 特定疾患の問題、介護の問題、就労の問題

⽉⽇ 時間 事象・対応内容

処理対象水に海水由来の塩分が含まれており,腐食