大學線性代數初步
大學數學
數學 大 學 線性代數 , 大
數 .
大 學 線性代數學 ,
學 線性代數 . 大學線性代數
. 學 大 大 ( ) 線性
代數 , . 大學 線性代數 , 大
學 數學 . 大 大 數學
, 學 , 線性代數
. , 線性代數
數學 . , .
, .
學 , 數學 學 , 線性
代數 . 線性代數 .
, . 學
. , (Question).
, 大
. , 線性代數 . ,
學 線性代數 學 線性代數 , .
, ,
代. , .
, . , 性
, . , .
v
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
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
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.
性 .
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.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) ,
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 ,
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−w′∥2=⟨v−w′, v−w′⟩ = ⟨v−w+w−w′, v−w+w−w′⟩ = ⟨v−w,v−w⟩+⟨w−w′, w−w′⟩.
∥v − w′∥2=∥v − w∥2+∥w − w′∥2. w̸= w′ ∥w − w′∥ > 0. ∥v − w′∥2>
∥v − w∥2, ∥v − w′∥ > ∥v − w∥.
A∈ Mm×n, Ax = b , b∈ C(A). ,
, Ax = b inconsistent , b
, b0 Ax = b′ b′ b
. 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
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, y 線 a 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. ,
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, y′i = axi+ b. (xi, yi) 線 y = ax + b
. b′=
y′1 y′2 ... y′m
, b′=
x1 1 x2 1 ... ... xm 1
[a
b ]
∈ C(A). εi= y′i− yi,
線 y = ax + b 代 xi y′i yi . 代 xi
, εi , , ,
ε12+··· +εm2 . ε12+··· +εm2 b′ b ,
∥b − b′∥2. , 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 y′i= axi+ b. 性 :
(1) y1+··· + ym= y′1+··· + y′m.
7.3. Orthogonal Basis and Gram-Schmidt Process 153
(2) x1y1+··· + xmym= x1y′1+··· + xmy′m.
(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′=
y′1 y′2 ... y′m
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− y′1 y2− y′2
... ym− y′m
=
[x1(y1− y′1) +··· + xm(ym− y′m) (y1− y′1) +··· + (ym− y′m)
]
= [0
0 ]
.
(1), (2).
y′= (y′1+··· + y′m)/m, (1) y = y′. i = 1, . . . , m, y′i= 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.
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⟩ = ∥wi∥2̸= 0, ci=⟨w,wi⟩/⟨wi, wi⟩ = ⟨w,wi⟩/∥wi∥2. ,
∥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⟩
∥w1∥2 w1+··· +⟨w,wi⟩
∥wi∥2 wi+··· +⟨w,wn⟩
∥wn∥2 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
∥wi∥wi, 1
∥wj∥wj⟩ = 1
∥wi∥ 1
∥wj∥⟨wi, wj⟩.
i̸= j, ⟨ui, uj⟩ = ⟨wi, wj⟩ = 0; i = j, ⟨ui, ui⟩ = ⟨wi, wi⟩/∥wi∥2= 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⟩
∥w1∥2w1+··· +⟨v,wn⟩
∥wn∥2wn. w1, . . . , wn W orthonormal basis,
ProjW(v) =⟨v,w1⟩w1+··· + ⟨v,wn⟩wn.
Theorem 7.3.3, W orthogonal basis,
W⊥ basis, . , W
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 =
∥w1∥2
∥w2∥2 . ..
∥wn∥2
, (ATA)−1=
∥w11∥2
∥w12∥2
. ..
∥w1n∥2
.
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
∥w1∥2w1wT1+ 1
∥w2∥2w2wT2+··· + 1
∥wn∥2wnwTn.
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 ··· wn
∥w11∥2
∥w12∥2
. ..
∥w1n∥2
— wT1 —
— wT2 — ...
— wTn —
(7.3)
m× m matrix 1
∥w1∥2w1wT1+ 1
∥w2∥2w2wT2+··· + 1
∥wn∥2wnwTn (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 ··· wn
∥w11∥2
∥w12∥2
. ..
∥w1n∥2
⟨w1, ei⟩
⟨w2, ei⟩ ...
⟨wn, ei⟩
=
w 1 w2 ··· wn
⟨w1,ei⟩
∥w1∥2
⟨w2,ei⟩
∥w2∥2
...
⟨wn,ei⟩
∥wn∥2
(7.4) ei
⟨w1, ei⟩
∥w1∥2 w1+⟨w2, ei⟩
∥w2∥2 w2+··· +⟨wn, ei⟩
∥wn∥2 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∥wi∥2. ∥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⟩/∥w1∥2)w1,
w2= v2−⟨v2, w1⟩
∥w1∥2 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