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

行列の演算

N/A
N/A
Protected

Academic year: 2025

シェア "行列の演算"

Copied!
22
0
0

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

全文

(1)

行列の演算

戸瀬 信之

ITOSE PROJECT

2011

04

19,26

at

駒場

2020

06

22

日 経済数学入門
(2)

ベクトルの演算の性質(復習)

定理

1

(1) ~ a , ~ b , ~ c ∈ K n

に対して

~ a + ~ b = ~ b + ~ a (1)

(~ a + ~ b ) + ~ c = ~ a + ( ~ b + ~ c ) (2) (2) ~ a , ~ b ∈ K n

λ, µ ∈ K

に対して

λ (µ~ a ) = (λµ)~ a (3)

(λ + µ) ~ a = λ~ a + µ~ a (4)

λ(~ a + ~ b ) = λ~ a + λ~ b (5)

(3)

和と定数倍

m × n

行列,すなわち

m

n

列の行列

A = (~ a 1 · · ·~ a j · · · ~ a n ) =

 a 1

a ... i

a ... m

= (a ij )

B = ( ~ b 1 · · · ~ b j · · · ~ b n ) =

 b 1

b ... i

b ... m

= ( b ij )

に対して

(4)

和と定数倍 (2)

その和と差を

A + B =

~ a 1 + ~ b 1 · · · ~ a j + ~ b j · · · ~ a n + ~ b n

=

a 1 + b 1

a i + ... b i

a m + ... b m

= (a ij + b ij )

A − B =

~ a 1 − ~ b 1 · · · ~ a j − ~ b j · · · ~ a n − ~ b n

=

a 1 − b 1

a i − ... b i

a m − ... b m

= ( a ij − b ij )

(5)

和と定数倍 (3)

定数

α ∈ K

に対して

A

α

倍を

α A = (α~ a 1 · · · α~ a j · · · α~ a n ) =

 α a 1

αa ... i

...

α a m

= (α · a ij )

(6)

定理 2

m

n

列の行列

A,B,C

に対して,以下が成立します.

(1) A + B = B + A

(2) (A + B) + C = A + (B + C ) (3) α(β A ) = (αβ) A

(4) (α + β) A = α A + β A (5) α(A + B) = αA + αB

A = (~ a j ) , B = ( ~ b j ) , C = (~ c j )

と行列

A , B , C

j

列を用いて定理

2

を示します.
(7)

証明

(1) A + B = (~ a j + ~ b j ) = ( ~ b j + ~ a j ) = B + A

から証明できます.ここで定理

1

(1)

を用いました.

(2) ( A + B ) + C = (~ a j + ~ b j ) + (~ c j ) =

(~ a j + ~ b j ) + ~ c j

= ~ a j + ( ~ b j + ~ c j )

= A + ( B + C )

から証明できます.ここで定理

1

(2)

を用いました.
(8)

行列 × 列ベクトル (1)

A = ~ a 1 · · · ~ a j · · · ~ a n

=

 a 1

a ... i

a ... m

~ x =

 x 1

x ... j

x ... n

 (6)

との積を

A ~ x = x 1 ~ a 1 + · · · + x j ~ a j + · · · + x n ~ a n ∈ K m

(9)

行列 × 列ベクトル (2)

A~ x = x 1 ~ a 1 + . . . + x j ~ a j + . . . + x n ~ a n

=

x 1 a i1 + . . . + x j ... a ij + . . . + x n a in

...

=

 a 1 ~ x

a ... i ~ x a m ... ~ x

(10)

行列 × 行列

X = (~ x 1 . . . ~ x k . . . ~ x ` )

n

`

列とします.

A ~ x 1 , . . . , A ~ x k , . . . , A ~ x ` ∈ K m

が定義される.

AX = A( ~ x 1 · · · ~ x k · · · ~ x ` )

= ( A ~ x 1 · · · A ~ x k · · · A ~ x ` )

=

a 1 ~ x 1 · · · a 1 ~ x k · · · a 1 ~ x `

... ... ...

a i ~ x 1 · · · a i ~ x k · · · a i ~ x `

... ... ...

a m ~ x 1 · · · a m ~ x k · · · a m ~ x `

(11)

行列 × 行列 (2)

特に

m = 1,

すなわち

A = a = ( a 1 · · · a n )

のとき

aX = ( a ~ x 1 · · · a ~ x ` )

これを用いると

AX =

 a 1 X

a i ... X a m ... X

(12)

行列 × 行列 (3)

aX = ( a ~ x 1 · · · a ~ x k · · · a ~ x ` )

= (· · · a 1 x 1k + a 2 x 2k + · · · + a n x nk · · · )

=

a 1 (· · · x 1k · · · ) + a 2 (· · · x 2k · · · )

... ...

+ a n (· · · x nk · · · )

= a 1 x 1 + · · · + a n x n

(13)

行列 × 行列 (4)— 行ベクトルを中心にして

まとめ

( a 1 · · · a j · · · a n )

 x 1

x ... j

x ... n

= a 1 x 1 + · · · + a j x j + · · · + a n x n

 a 1

a ... i

a ... m

 X =

 a 1 X a i ... X a m ... X

(14)

線型写像 (1)

F A ( ~ x ) = A ~ x

によって定る写像

F A : K n −→ K m ~ x 7→ A ~ x (7)

について定理

3

が成立します(F

A

の線型性)

定理

3

F A (~ x + ~ y ) = F A (~ x ) + F A (~ y ) (~ x , ~ y ∈ K n ) (8) F A (λ~ x ) = λF A (~ x ) (~ x ∈ K n , λ ∈ K) (9)

(15)

線型写像 (2)— 定理 3 の証明

~ x =

 x 1

x ... n

 , ~ y =

 y 1

y ... n

~ x

~ y

の成分表示をします.すると

F A (~ x + ~ y) = (x 1 + y 1 )~ a 1 + · · · + (x j + y j )~ a j + · · · + (x n + y n )~ a n

= x 1 ~ a 1 + · · · + x j ~ a j + · · · + x n ~ a n

+ y 1 ~ a 1 + · · · + y j ~ a j + · · · + y n ~ a n

= F A (~ x ) + F A (~ y )

F A (λ~ x ) = (λ x 1 ) ~ a 1 + · · · + (λ x j ) ~ a j + · · · + (λ x n ) ~ a n

= λ(x 1 ~ a 1 + · · · + x j ~ a j + · · · + x n ~ a n ) = λF A ( ~ x )

(16)

線型写像 (3)— 定理 3 の拡張

行列の積の形では

A(~ x + ~ y ) = A ~ x + A ~ y , A(λ~ x ) = λ (A ~ x ) (10)

と表されます.これをまとめて得られる

A (λ~ x + µ~ y ) = λ A ~ x + µ A ~ y (11)

および

(10)

を繰り返して得られる

A ( c 1 ~ x 1 + · · · + c ` ~ x ` ) = c 1 A ~ x 1 + · · · + c ` A ~ x ` (12)

も有用です.この公式

(12)

X = (~ x 1 · · · ~ x ` )

~ c = t ( c 1 · · · c ` )

用いて

A(X ~ c) = (AX )~ c (13)

(17)

写像の合成

2

つの写像

f : X → Y

g : Y → Z

があるとします.このとき

g ◦ f : X → Z x 7→ g (f (x ))

が定義できます.これを

f

g

の合成写像と呼びます.

さらに写像

h : Z → W

がある場合は

h ◦ g ◦ f : X → W

が定義できますが,これは

h ◦ (g ◦ f ) = (h ◦ g ) ◦ f

と,合成をどの順序で行っても変わりません.
(18)

行列の積の性質 (1)

定理

4 A

B

m × n

行列,

X

Y

n × `

行列,

Q

` × g

行 列とします.このとき次が成立します.

(1) (A + B)X = AX + BX (2) A ( X + Y ) = AX + AY (3) A (α X ) = (α A ) X = α( AX ) (4)(結合法則) (AX ) Q = A (XQ)

(19)

行列の積の性質 (2)— 証明の準備(その 1 )

定理

4

の証明の準備として

A = (~ a 1 · · · ~ a n ) , B =

~ b 1 · · · ~ b n

~ x ∈ K n

に対して

( A + B ) ~ x = A ~ x + B ~ x (14)

を示します.実際

(A + B)~ x =

~ a 1 + ~ b 1 · · · ~ a n + ~ b n

~ x

= x 1 ( ~ a 1 + ~ b 1 ) + · · · + x n (~ a n + ~ b n )

= x 1 ~ a 1 + · · · + x n ~ a n + x 1 ~ b 1 + · · · + x n ~ b n = A~ x + B~ x

と示すことができます.
(20)

行列の積の性質 (3)— 証明の準備(その 2 )

(10)

で以下の最初の等号を示しましたが

A (α~ x ) = α ( A ~ x ) = (α A )~ x (15)

が成立することも定理

4

(3)

の証明で使います.実際,2番目の 等号は

~ x = t (x 1 · · · x n )

として

(α A ) ~ x = α~ a 1 · · · α~ a n

~ x

= x 1 α~ a 1 + · · · + x n α~ a n = α( x 1 ~ a 1 + · · · + x n ~ a n ) = α( A ~ x )

(21)

行列の積の性質 (4)— 証明 (1)

(1)

両辺の

k

列を比較します.

( A + B ) ~ x k = A ~ x k + B ~ x k

から分かります((14)参照).

(2)

両辺の

k

列を比較します.

A (~ x k + ~ y k ) = A ~ x k + A ~ y k

から分かります((10)参照).

(3)

各辺の

k

列を比較しますが

(15)

A (α~ x k ) = (α A )~ x k = α( A ~ x k )

を導きます.
(22)

行列の積の性質 (5)— 証明 (2)

(4) (13)

を用いると

~ q ∈ K `

に対して

A (X ~ q) = (AX ) ~ q

が成立します.Q

t

列を

~ q t

とすると

A (X ~ q t ) = (AX ) ~ q t

ですが,これは示すべき式の

t

列が等しいことを意味します.これ を用いると

A(XQ) = A(X ~ q 1 · · · X ~ q g ) = (A(X ~ q 1 ) · · · A(X ~ q g ))

= ((AX )~ q 1 · · · (AX )~ q g ) = (AX )(~ q 1 · · · ~ q g ) = (AX )Q

参照

関連したドキュメント

モデルの形態と作成の考え方

参考文献 43 参考文献 [1] 数値計算の手順と実際 高田勝 春海佳三朗 コロナ社 [2] FOR TRAN77 数値計算プログラミング 森正武 岩波書店 [3] C

概要:大規模疫学データを解析する状況で数万人のデータに対し統計解析モデルを用いる場合,一般に条

[r]

Lecture Guidance and Review or the Matrix

[r]

オーバラップを行った SUMMA のアルゴリ ズムを実装せよ。また、マルチキャスト版 SUMMA と、性能を比較せよ。.. [L20]

オーバラップを行った SUMMA のアルゴリ ズムを実装せよ。また、マルチキャスト版 SUMMA と、性能を比較せよ。.. [L20]