• 沒有找到結果。

大學線性代數初步

N/A
N/A
Protected

Academic year: 2022

Share "大學線性代數初步"

Copied!
19
0
0

加載中.... (立即查看全文)

全文

(1)

大學線性代數初步

大學數學

(2)

數學 大 學 線性代數 , 大

數 .

大 學 線性代數學 ,

學 線性代數 . 大學線性代數

. 學 大 大 ( ) 線性

代數 , . 大學 線性代數 , 大

學 數學 . 大 大 數學

, 學 , 線性代數

. , 線性代數

數學 . , .

, .

學 , 數學 學 , 線性

代數 . 線性代數 .

, . 學

. , (Question).

, 大

. , 線性代數 . ,

學 線性代數 學 線性代數 , .

, ,

代. , .

, . , 性

, . , .

v

(3)

Chapter 7

Inner Product Space

Rm dot product 性 , vector space. vector

space, dot product inner product . vector space V

v, w∈ V ⟨v,w⟩ ∈ R Proposition 1.4.2 性 ,

(1) ⟨v,w⟩ = ⟨w,v⟩, ∀v,w ∈ V.

(2) ⟨v,v⟩ ≥ 0 ⟨v,v⟩ = 0 v = O.

(3) r∈ R ⟨rv,v⟩ = ⟨v,rw⟩ = r⟨v,w⟩.

(4) ⟨u,v + w⟩ = ⟨u,v⟩ + ⟨u,w⟩, ∀u,v,w ∈ V.

V inner product space. inner product space 性 .

性 inner product 性 , finite dimensional

inner product space 性 . 大 inner product

, 大 Rm . ,

, v, w∈ Rm, v· w v, w ,

, ⟨v,w⟩ v, w . , v, w column vectors

, ⟨v,w⟩ = vTw.

7.1. Projection

R2 linear transformation projection. ,

w∈ R2, v∈ R2 w Projw(v) w (

Span(w)), v− Projw(v) w , Span(w) . ,

Rm subspace W . , W Rm

subspace, v∈ Rm, v W projection W w, v− w

W . .

Definition 7.1.1. W Rm subspace.

W={v ∈ Rm| ⟨v,w⟩ = 0, ∀w ∈ W}.

143

(4)

W Rm W . W orthogonal complement of W .

orthogonal complement projection .

Definition 7.1.2. W Rm subspace. v∈ Rm, w∈ W v−w ∈ W, w the projection of v on W .

v∈ Rm, the projection of v on W , .

ProjW(v) . 性, 性,

Rm w ProjW(v), w∈ W v− w ∈ W .

orthogonal complement 性 . W Rm subspace, v, v∈ W r, s∈ R. 性 (3)(4), w∈ W, ⟨rv+sv, w⟩ = r⟨v,w⟩ + s⟨v, w⟩. v, v∈ W, ⟨v,w⟩ = ⟨v, w⟩ = 0, ⟨rv + sv, w⟩ = 0.

v, v∈ W r, s∈ R, rv + sv∈ W, W Rm subspace.

orthogonal complement complement . W orthogonal

complement W W . W∩W . W∩W

subspace, O . .

Lemma 7.1.3. W Rm subspace. W∩W={O}.

Proof. W∩W subspace, O∈ W ∩W. w∈ W ∩W. w∈ W,

w∈ W ⟨w,w⟩ = 0. w∈ W, ⟨w,w⟩ = 0. 性 (2) w = O,

W∩W={O}. 

projection 性.

Proposition 7.1.4. W Rm subspace. v∈ Rm, v W projection .

Proof. w, w∈ W v W projection. w, w v−w ∈ W v− w∈ W. W subspace, (v− w) − (v − w)∈ W, w− w∈ W.

W subspace, w− w∈ W. Lemma 7.1.3 w− w= O, w = w,

性. 

projection 性. W subspace, 步

dimension. W W , dim(W) dim(W ) .

dim(W ) = n w1, . . . , wn W basis. A =



— w1 — ...

— wn

,

A i-th row wi ( row vector). wi∈ Rm, A n× m matrix.

v∈ W, wi∈ W, ⟨wi, v⟩ = 0, ∀i = 1,...,n. Av = O,

W A nullspace. , v∈ Rm A nullspace , Av = O

(5)

7.1. Projection 145

⟨wi, v⟩ = 0, ∀i = 1,...,n. w∈ W, w1, . . . , wn W basis, c1, . . . , cn∈ R w = c1w1+··· + cnwn.

⟨w,v⟩ = ⟨c1w1+··· + cnwn, v⟩ = c1⟨w1, v⟩ + ··· + cn⟨wn, v⟩ = 0,

A nullspace W. 性 .

Lemma 7.1.5. W Rm subspace. w1, . . . , wn W basis. A w1, . . . , wn ( row vector) row vector n× m matrix. W A nullspace.

Example 7.1.6. W = Span(



 1 2 2 1



,



 3 4 2 3



) W orthogonal complement W.

A =

[ 1 2 2 1 3 4 2 3

]

. Lemma 7.1.5, A nullspace{x ∈ R4| Ax = O} W. elementary row operation A 1-st row −3 2-nd row

[ 1 2 2 1

0 −2 −4 0

]

. 2-nd row −1/2

[ 1 2 2 1 0 1 2 0

] .



−1 0 0 1



,



 2

−2 1 0



A

nullspace, W basis.

Question 7.1. Lemma 7.1.5 w1, . . . , wn Span(w1, . . . , wn) = W linearly

independent, W A nullspace?

Rm column vector, Lemma 7.1.5 w1, . . . , wn

column vector A column vector ( A m×n matrix), Lemma

7.1.5 .

Corollary 7.1.7. W Rm subspace. w1, . . . , wn W basis. A w1, . . . , wn column vector m× n matrix. W AT nullspace.

Corollary 7.1.7 column space W = C(A), Corollary 7.1.7 .

Corollary 7.1.8. A m× n matrix. C(A)= N(AT).

dim(W) dim(W ) . W Rm

subspace dim(W ) = n. W basis row vectors n× m matrix A,

Lemma 7.1.5 W A nullspace. dim(W) A nullspace

( A nullity). A rank A row space ( A row vectors space W ) , rank(A) = dim(W ) = n. Theorem 4.4.13 A nullity A column 數 m rank(A), null(A) = m− rank(A) = m − n.

性 .

(6)

Theorem 7.1.9. W Rm subspace. dim(W) = m− dim(W).

Theorem 7.1.9 .

Corollary 7.1.10. W Rm subspace. (W)= W .

Proof. Theorem 7.1.9 dim((W)) = m− dim(W) = m− (m − dim(W)) = dim(W ). w∈ W, w∈ W, ⟨w,w⟩ = 0, w∈ (W).

W⊆ (W). dim(W ) = dim((W)) Corollary 4.3.6 (W)= W . 

, Corollary 7.1.10 dimension . W

infinite dimensional vector space subspace , (W)= W .

W⊆ (W), (W)⊆ W . finite dimensional , dimension , Corollary 4.3.6 (W)= W .

Question 7.2. A m× n matrix. C(A) = N(AT).

Theorem 7.1.9 W Rm subspace w1, . . . , wn W basis, dim(W) = m− n, W basis u1, . . . , um−n.

w1, . . . , wn, u1, . . . , um−n linearly independent, .

Corollary 7.1.11. W Rm subspace w1, . . . , wn W basis u1, . . . , um−n W basis. w1, . . . , wm, u1, . . . , um−n Rm basis.

Proof. w1, . . . , wn, u1, . . . , um−n linearly independent. , c1, . . . , cn, cn+1, . . . , cm∈ R 0 c1w1+··· + cnwn+ cn+1u1+··· + cmum−n= O.

w = c1w1+··· + cnwn, W vector space, w∈ W. w =−(cn+1u1+

··· + cmum−n) W vector space, w∈ W. 言 , w∈ W ∩W.

Lemma 7.1.3 w = O. w1, . . . , wn linearly independent c1w1+··· + cnwn= w = O c1 =··· = cn = 0. cn+1 =··· = cm= 0. c1, . . . , cn, cn+1, . . . , cm

0 , w1, . . . , wn, u1, . . . , um−n linearly independent.

dim(Rm) = m, w1, . . . , wn, u1, . . . , um−n Rm m linearly independent vectors,

w1, . . . , wn, u1, . . . , um−n Rm basis. 

W basis{w1, . . . , wn} W basis{u1, . . . , um−n} Rm basis, v∈ Rm, c1, . . . , cn, d1, . . . , dm−n∈ R

v = c1w1+··· + cnwn+ d1u1+··· + dm−num−n.

w = c1w1+··· + cnwn, w∈ W v− w = d1u1+··· + dm−num−n∈ W.

w the projection of v on W . projection 性. ,

.

(7)

7.1. Projection 147

Theorem 7.1.12. W Rm subspace w1, . . . , wn W basis u1, . . . , um−n W basis. v∈ Rm, v = c1w1+··· + cnwn+ d1u1+··· + dm−num−n, v W projection

ProjW(v) = c1w1+··· + cnwn.

Example 7.1.13. Example 7.1.6 v =



 2 0 1 4



W = Span(



 1 2 2 1



,



 3 4 2 3



)

projection.



−1 0 0 1



,



 2

−2 1 0



W basis,



 2 0 1 4





 1 2 2 1



,



 3 4 2 3



,



−1 0 0 1



,



 2

−2 1 0



 linear combination



 2 0 1 4



 = −



 1 2 2 1



 +



 3 4 2 3



 + 2



−1 0 0 1



 +



 2

−2 1 0



.

v W projection



 1 2 2 1



 +



 3 4 2 3



 =



 2 2 0 2



.

Rm subspace W , v∈ Rm, ProjW(v) ,

projection on W 數. ProjW Rm→ Rm 數 v∈ Rm

ProjW(v). ProjW:Rm→ Rm linear transformation.

v, v∈ Rm, w =ProjW(v) w= ProjW(v).

w, w∈ W v− w ∈ W v− w∈ W. ProjW(v + v) = ProjW(v) + ProjW(v) = w + w. 言 , w + w∈ W (v + v)− (w + w)∈ W. W,W vector space, w, w∈ W v− w,v− w∈ W w + w∈ W

(v + v)− (w + w) = (v− w) + (v− w)∈ W. r∈ R,

ProjW(rv) = rProjW(v) = rw. rw∈ W rv− rw ∈ W. W,W

vector space, w∈ W,v − w ∈ W, .

Proposition 7.1.14. W Rm subspace, ProjW :Rm→ Rm linear trans- formation.

Question 7.3. ProjW:Rm→ Rm image kernel.

Theorem 7.1.12 v W . W

basis w1, . . . , wn W basis u1, . . . , um−n, v w1, . . . , wn, u1, . . . , um−n

線性 , ProjW(v). 步 . A w1, . . . , wn column

vectors m× n matrix, W A column space C(A). ProjW(v) ,

(8)

w∈ W v− w ∈ W, w ProjW(v) . w∈ W ? column space , w∈ W x∈ Rn Ax = w. v−w ∈ W= (C(A))

, Corollary 7.1.8 v− w = v − Ax ∈ N(AT). x∈ Rn

AT(v−Ax) = O. 性 , x∈ Rn (ATA)x = ATv.

ProjW(v) , x∈ Rn Ax = ProjW(v), x

v− Ax ∈ W, (ATA)x = ATv . (ATA)x = ATv, x

Ax = ProjW(v) . (ATA)x = ATv ? .

Lemma 7.1.15. A∈ Mm×n rank(A) = n, ATA n× n invertible matrix.

Proof. A∈ Mm×n AT∈ Mn×m, ATA∈ Mn×n. Theorem 3.5.9, ATA invertible. (ATA)x = O nontrivial solution

ATA invertible matrix. x = u (ATA)x = O , (ATA)u = O.

(ATA)u = AT(Au), Au∈ N(AT). Au∈ C(A), Au∈ C(A) ∩ N(AT) =

C(A)∩C(A). Lemma 7.1.3 Au = O. rank(A) = n, Corollary

2.3.5 Ax = O nontrivial solution, u = O. (ATA)x = O

nontrivial solution, ATA n× n invertible matrix. 

A W basis column vector m× n matrix,

rank(A) = dim(W ) = n. Lemma 7.1.15 ATA invertible.

(ATA)x = ATv ATA inverse, x = (ATA)−1ATv.

(ATA)x = ATv , A ProjW(v). .

Theorem 7.1.16. W Rm subspace w1, . . . , wn W basis. A w1, . . . , wn column vector m× n matrix. v∈ Rm, v W projection

ProjW(v) = A(ATA)−1ATv.

Example 7.1.17. Theorem 7.1.16 Example 7.1.13 ,

v =



 2 0 1 4



W = Span(



 1 2 2 1



,



 3 4 2 3



) . A =



 1 3 2 4 2 2 1 3



,

AT=

[ 1 2 2 1 3 4 2 3

]

, ATA =

[ 10 18 18 38

]

inverse (ATA)−1= (1/28)

[ 19 −9

−9 5

] . Theorem 7.1.16

ProjW(v) = 1 28



 1 3 2 4 2 2 1 3



[ 19 −9

−9 5

][ 1 2 2 1 3 4 2 3

]



 2 0 1 4



 =



 2 2 0 2



.

Example 7.1.13 .

Theorem 7.1.16 W Rm subspace, W =Rm,

dim(W ) = n m. W =Rm , W projection ,

(9)

7.2. Inconsistent Systems 149

W={O}, Rm W =Rm .

dim(W ) = n < m . , W basis column vector

A, m× n matrix . A AT invertible.

(ATA)−1 A−1(AT)−1. Theorem 7.1.16 A(ATA)−1AT A(A−1(AT)−1)AT, identity matrix.

Theorem 7.1.16 Theorem 7.1.12 projection . W

basis , W basis. A(ATA)−1AT Rm v,

ProjW(v). A(ATA)−1AT PW , W projection

matrix. , PW ProjW :Rm→ Rm linear transformation

standard matrix representation. linear transformation standard matrix

representation , W basis column vectors m× n

matrix A , A projection matrix .

Question 7.4. Example 7.1.17 W projective matrix ?

Theorem 7.1.16 projection , projection matrix

(ATA)−1 .

projection .

7.2. Inconsistent Systems

projection .

linear system inconsistent , .

大 學 線.

性 (2) O , v∈ V ⟨v,v⟩ > 0,

∥v∥ =

⟨v,v⟩. R3 .

. 性 inner product

space projection .

Proposition 7.2.1. W Rm subspace v∈ Rm. w v W , w∈ W w̸= w, ∥v − w∥ > ∥v − w∥.

Proof. v− w= v− w + w − w. w =ProjW(v), v− w ∈ W. W vector space, w− w∈ W.

∥v−w2=⟨v−w, v−w⟩ = ⟨v−w+w−w, v−w+w−w⟩ = ⟨v−w,v−w⟩+⟨w−w, w−w⟩.

∥v − w2=∥v − w∥2+∥w − w2. w̸= w ∥w − w∥ > 0. ∥v − w2>

∥v − w∥2, ∥v − w∥ > ∥v − w∥. 

A∈ Mm×n, Ax = b , b∈ C(A). ,

, Ax = b inconsistent , b

, b0 Ax = b b b

(10)

. b0 Ax = b0 , .

Ax = b b C(A) subspace, Proposition 7.2.1

b b b0 b C(A) projection. projection

, b− b0∈ C(A)= N(AT). AT(b− b0) = O,

ATb0= ATb. (7.1)

Ax = b0

ATAx = ATb. (7.2)

(7.2) (7.1) , b0 Ax = b0,

ATAx = ATb. (7.2) .

x = x0 ATAx = ATb , b0= Ax0 b0∈ C(A) b0−b ∈ N(AT) = C(A). b0 b C(A) projection, Proposition 7.2.1 b0

Ax = b b b . ATAx = ATb ,

. (7.2) ?

Proposition 7.2.2. A∈ Mm×n b∈ Rm. ATAx = ATb .

rank(A) = n, .

Proof. b0 b C(A) projection, b0∈ C(A) b− b0∈ C(A) . b0∈ C(A) x0∈ Rn Ax0= b0. b− b0∈ C(A) C(A)= N(AT) (Corollary 7.1.8)

O = AT(b− b0) = ATb− ATb0= ATb− AT(Ax0),

x0 ATAx0= ATb, ATAx = ATb .

rank(A) = n, Lemma 7.1.15 ATA invertible, ATAx = ATb

. 

rank(A) < n , ATAx = ATb

Ax = ProjC(A)(b) , A m×n matrix Theorem 3.4.6 (

). .

Definition 7.2.3. A∈ Mm×n b∈ Rm. Ax = b.

ATAx = ATb Ax = b normal equations. normal equations , Ax = b least squares solution.

, rank(A) = n , ATAx = ATb x = (ATA)−1ATb, Ax = b least squares solution.

, Ax = b least squares solution, normal equation

, b C(A) projection. least squares solution ,

A . Ax = b consistent, b∈ C(A), b

C(A) projection b . Ax = b normal equations ATAx = ATb

(11)

7.2. Inconsistent Systems 151

. Ax = b least squares solution ,

Ax = b , normal equations Ax = b .

Example 7.2.4. Ax = b A =



 1 3 2 4 2 2 1 3



 b =



 2 0 1 4



,

least squares solution. ATA =

[ 10 18 18 38

]

ATb = [8

20 ]

, Ax = b

normal equations

10x1+ 18x2 = 8 18x1+ 38x2 = 20 . [x1

x2 ]

= [−1

1 ]

Ax = b least squares solution. rank(A) = 2, ATA invertible inverse (ATA)−1= (1/28)

[ 19 −9

−9 5

] . [x1

x2

]

= (ATA)−1ATb = 1 28

[ 19 −9

−9 5

][8 20

]

= [−1

1 ]

.



 1 3 2 4 2 2 1 3



 [−1

1 ]

=



 2 2 0 2



 b C(A) projection ( Example 7.1.17).

Question 7.5. W Rm subspace W = Span(w1, . . . , wn). A w1, . . . , wn column vector m× n matrix. v∈ Rm x = x0∈ Rn Ax = v normal equations ATAx = ATv , v W projection Ax0. w1, . . . , wn

linearly independent , normal equations , ProjW(v)

? W projection matrix PW A(ATA)−1AT ?

least squares solution , 線.

(x1, y1), . . . , (xm, ym), 線 y = ax + b, (xi, yi) 線. x1, . . . , xm y1, . . . , ym 數,

x, ya y b. ,

a, b

ax1+ b = y1

ax2+ b = y2

... axm+ b = ym.

Ax = b, A =



 x1 1 x2 1 ... ... xm 1



, x = [a

b ]

, b =



 y1

y2 ... ym



. ,

(x1, y1), . . . , (xm, ym) Ax = b , Ax = b least squares solution, ATAx = ATb. ,

(12)

xi, yi , x1, . . . , xm ( (x1, y1), . . . , (xm, ym)

線 , ), rank(A) = 2, Proposition 7.2.2

normal equations ATAx = ATb . , 線

?

a, b∈ R, i = 1, . . . , m, yi = axi+ b. (xi, yi) 線 y = ax + b

. b=



 y1 y2 ... ym



, b=



 x1 1 x2 1 ... ... xm 1



 [a

b ]

∈ C(A). εi= yi− yi,

線 y = ax + b 代 xi yi yi . 代 xi

, εi , , ,

ε12+··· +εm2 . ε12+··· +εm2 b b ,

∥b − b2. , Ax = b least squares solution x = [a

b ]

, a, b

線 y = ax + b b b C(A) , ∥b − b

. ε12+··· +εm2 . least squares solution

線, least squares line ( 線, 線), 學

line of regression ( 線).

Example 7.2.5. (−1,0),(1,1),(2,3) least square line y = ax + b.

−1a + b = 0 1a + b = 1 2a + b = 3,

−1 1 1 1 2 1

[ a b ]

=

0 1 3

. normal equations

[ −1 1 2 1 1 1

] −1 1 1 1 2 1

[ a b ]

=

[ −1 1 2 1 1 1

]0 1 3

,

6a + 2b = 7 2a + 3b = 4,

least squares solution a = 13/14, b = 5/7 least squares line y =13

14x +5 7. least squares line 性 .

Proposition 7.2.6. (x1, y1), . . . , (xm, ym). y = ax + b least squares line i = 1, . . . , m yi= axi+ b.:

(1) y1+··· + ym= y1+··· + ym.

(13)

7.3. Orthogonal Basis and Gram-Schmidt Process 153

(2) x1y1+··· + xmym= x1y1+··· + xmym.

(3) x y x1, . . . , xm y1, . . . , ym 數, x = (x1+··· + xm)/m, y = (y1+··· + ym)/m. (x, y) 線 y = ax + b , y = ax + b.

Proof. b =



 y1 y2 ... ym



, b=



 y1 y2 ... ym



 A =



 x1 1 x2 1 ... ... xm 1



, x = [a

b ]

Ax = b

least squares solution b= A [a

b ]

b C(A) projection. b− b∈ C(A)= N(AT),

AT(b− b) =

[ x1 ··· xm

1 ··· 1 ]



 y1− y1 y2− y2

... ym− ym



=

[x1(y1− y1) +··· + xm(ym− ym) (y1− y1) +··· + (ym− ym)

]

= [0

0 ]

.

(1), (2).

y= (y1+··· + ym)/m, (1) y = y. i = 1, . . . , m, yi= axi+ b,

y= ax + b. (3), y = ax + b. 

, 線,

線. 線 ,

線 y = ax2+ bx + c,a, b, c 數 ,

least squares solution. 線

線 . , .

7.3. Orthogonal Basis and Gram-Schmidt Process

Rm v subspace W projection .

. , Theorem 7.1.12, W basis

W basis, v basis Rm basis linear combination,

W , v W projection . w

basis linear combination, .

linear combination , . , Theorem

7.1.16, W basis column vectors A, W

projection matrix. , ATA inverse.

W basis. basis, .

basis, ?

, basis.

W Rn subspace w1, . . . , wn basis ⟨wi, wj⟩ = 0,∀i ̸= j.

w∈ W, w1, . . . , wn W basis, c1, . . . , cn∈ R w = c1w1+···+cnwn.

(14)

c1, . . . , cn, wi

, ci. i = 1, . . . , n, ⟨w,wi⟩.

⟨w,wi⟩ = ⟨c1w1+··· + cnwn, wi⟩ = c1⟨w1, wi⟩ + ··· + cn⟨wn, wi⟩ = ci⟨w1, wi⟩.

wi̸= O, ⟨wi, wi⟩ = ∥wi2̸= 0, ci=⟨w,wi⟩/⟨wi, wi⟩ = ⟨w,wi⟩/∥wi2. ,

∥wi∥ = 1, ci=⟨w,wi⟩. .

Proposition 7.3.1. W Rn subspace w1, . . . , wn basis ⟨wi, wj⟩ = 0,∀i ̸= j. w∈ W,

w = ⟨w,w1

∥w12 w1+··· +⟨w,wi

∥wi2 wi+··· +⟨w,wn

∥wn2 wn.

basis, linear combination ,

.

Definition 7.3.2. W Rm subspace w1, . . . , wn W basis i̸= j ⟨wi, wj⟩ = 0. w1, . . . , wn W orthogonal basis.

⟨wi, wi⟩ = 1, ∀i = 1,...,n, w1, . . . , wn W orthonormal basis.

, w1, . . . , wn W orthogonal basis, i = 1, . . . , n, ui = wi/∥wi∥, u1, . . . , un W orthonormal basis.

⟨ui, uj⟩ = ⟨ 1

∥wiwi, 1

∥wjwj⟩ = 1

∥wi 1

∥wj∥⟨wi, wj⟩.

i̸= j, ⟨ui, uj⟩ = ⟨wi, wj⟩ = 0; i = j, ⟨ui, ui⟩ = ⟨wi, wi⟩/∥wi2= 1.

前 , W orthogonal basis w1, . . . , wn,

W w1, . . . , wn 線性 . 前 Theorem 7.1.12

. , W orthogonal basis w1, . . . , wn,

v∈ Rm W projection . w =ProjW(v) w = c1w1+···+cnwn. v− w ∈ W, i = 1, . . . , n, ⟨v,wi⟩ − ⟨w,wi⟩ = ⟨v − w,wi⟩ = 0.

⟨v,wi⟩ = ⟨w,wi⟩ = ci⟨wi, wi⟩, ci=⟨v,wi⟩/⟨wi, wi⟩. Theorem

7.1.12 projection .

Theorem 7.3.3. W Rm subspace w1, . . . , wn W orthogonal basis.

v∈ Rm, v W projection

ProjW(v) =⟨v,w1

∥w12w1+··· +⟨v,wn

∥wn2wn. w1, . . . , wn W orthonormal basis,

ProjW(v) =⟨v,w1⟩w1+··· + ⟨v,wn⟩wn.

Theorem 7.3.3, W orthogonal basis,

W basis, . , W

(15)

7.3. Orthogonal Basis and Gram-Schmidt Process 155

orthogonal basis w1, . . . , wn W projection matrix. A w1, . . . , wn column vectors m× n matrix,

ATA =





— wT1

— wT2 — ...

— wTn





 w1 w2 ··· wn

(i, j)-th entry wTiwj=⟨wi, wj⟩, w1, . . . , wn orthogonal, ATA diagonal

matrix, (ATA)−1,

ATA =





∥w12

∥w22 . ..

∥wn2



, (ATA)−1=





∥w112

∥w122

. ..

∥w1n2





.

Theorem 7.1.16 projection , .

Theorem 7.3.4. W Rm subspace w1, . . . , wn W orthogonal basis.

v∈ Rm, ProjW(v) = PWv, PW m× m matrix PW = 1

∥w12w1wT1+ 1

∥w22w2wT2+··· + 1

∥wn2wnwTn.

Proof. wi m× 1 matrix, wiwTi m× m matrix.

A w1, . . . , wn column vectors m× n matrix. Theorem 7.1.16 PW = A(ATA)−1AT. w1, . . . , wn orthogonal, 前

A(ATA)−1AT=

 w 1 w2 ··· w n





∥w112

∥w122

. ..

∥w1n2









— wT1

— wT2 — ...

— wTn



 (7.3)

m× m matrix 1

∥w12w1wT1+ 1

∥w22w2wT2+··· + 1

∥wn2wnwTn (7.4)

(7.3) m× m matrix .

i = 1, . . . , m, ei∈ Rm, i-th entry 1 entry 0 vector.

m× m matrix ei matrix i-th column. (7.3) (7.4)

ei , (7.3) (7.4) , .

(7.3) ei

 w 1 w2 ··· w n





∥w112

∥w122

. ..

∥w1n2









⟨w1, ei

⟨w2, ei ...

⟨wn, ei



=

 w 1 w2 ··· w n







⟨w1,ei

∥w12

⟨w2,ei

∥w22

...

⟨wn,ei

∥wn2







(16)

(7.4) ei

⟨w1, ei

∥w12 w1+⟨w2, ei

∥w22 w2+··· +⟨wn, ei

∥wn2 wn.

 Question 7.6. ProjW:Rm→ Rm linear transformation linear transformation

standard matrix representation PW, Theorem 7.3.3 Theorem 7.3.4.

W orthogonal basis, v∈ Rm W

projection. Rm nonzero subspace, orthogonal

basis, Rm nonzero subspace orthogonal basis (

orthonormal basis). Gram-Schmidt process.

Rm nonzero subspace W , dim(W ) = n. w1, . . . , wk W ⟨wi, wj⟩ = 0, ∀i ̸= j, w1, . . . , wk linearly independent.

linearly independent c1, . . . , ck ∈ R 0 c1w1+··· + ckwk = O.

i = 1, . . . , k 0 =⟨c1w1+··· + ckwk, wi⟩ = ci∥wi2. ∥wi∥ ̸= 0, ci= 0. c1, . . . , ck 0 , w1, . . . , wk linearly independent.

W orthogonal basis, W w1, . . . , wn ⟨wi, wj⟩ = 0,

∀i ̸= j , linearly independent dim(W ) = n, Corollary 4.3.5

W basis. W nonzero vectors.

W̸= {O}, W nonzero vector v1. w1= v1

W1= Span(v1) = Span(w1). dim(W ) = 1, W = W1 w1 W orthogonal basis. dim(W ) > 1, W1( W, nonzero vector v2∈ W v2̸∈ W1.

v2, w2 ∈ W w2 ̸= O ⟨w1, w2⟩ = 0. ,

w2= v2− ProjW1(v2), w2 ∈ W1, W1 = Span(w1), ⟨w1, w2⟩ = 0.

w2 ̸= O. w2= O, v2 = ProjW1(v2)∈ W1, 初

v2 ̸∈ W1 . W = Span(w1), Proposition 1.4.9

ProjW1(v2) = (⟨v2, w1⟩/∥w12)w1,

w2= v2−⟨v2, w1

∥w12 w1.

Span(w1, w2) = Span(v1, v2), w1, w2 , w1, w2 Span(v1, v2), Span(w1, w2)⊆ Span(v1, v2). v1, v2 linearly independent

w1, w2 linearly independent, dim(Span(w1, w2)) = dim(Span(v1, v2)) = 2 Span(w1, w2) = Span(v1, v2). , W2= Span(w1, w2) = Span(v1, v2).

dim(W ) = 2, W = W2, w1, w2 W orthogonal basis. dim(W ) > 2, W2( W, nonzero vector v3∈ W v3̸∈ W2. v3, w3∈ W

w3̸= O ⟨w1, w3⟩ = ⟨w2, w3⟩ = 0. 前, w3= v3− ProjW1(v2),

w3∈ W2, W2= Span(w1, w2), ⟨w1, w3⟩ = ⟨w2, w3⟩ = O. v3̸∈ W2, 前 w3̸= O. w1, w2 W2 orthogonal basis, Theorem 7.3.3

參考文獻

相關文件

崑山科技大學資訊工程系99 學年度 學年度 學年度 學年度 第一學期 第一學期 第一學期線性代數 第一學期 線性代數 線性代數 線性代數平時 平時 平時考試題 平時 考試題 考試題 考試題. 姓名 姓名

近世代數的重要內容是集合及這些集合上的代數運算。 集合本身和作為代數運算的載體的 集合是不加區分的, 故實質上研究的是代數運算本身。 說更仔細些, 考慮非空集合S上一個或幾

elementary row operations reduced echelon form,. echelon form Gauss

augmented matrix [A |I 4 ], elementary row operation A

column vector

Proposition 9.4.2, A orthogonal diagonalizable, Spectral Theorem.. Theorem 9.4.6

Theorem 8.2.6 (3) elementary column operation.. determinant elementary row

, A echelon form ( reduced echelon form) pivot column vectors.. elementary row operations column