LA@标准正交化@施密特正交化
文章目录
基本定理👺
n n n维欧式空间中任意一个正交向量组能够扩充为正交基
- n n n维欧式空间 V V V中任意一个正交向量组能够扩充为正交基
- 证明:设
A
0
:
α
1
,
⋯
,
α
m
A_0:\alpha_1,\cdots,\alpha_m
A0:α1,⋯,αm是一组正交向量组,对
p
=
n
−
m
p=n-m
p=n−m作数学归纳法
- 首先, n n n维空间 V V V的基包含的向量个数是 n n n
- p = 0 p=0 p=0时, A A A本身就是 V V V的正交基
- 设 p = k p=k p=k时定理成立,即存在 B : β 1 , ⋯ , β k B:\beta_1,\cdots,\beta_k B:β1,⋯,βk使得: α 1 , ⋯ , α m , β 1 , ⋯ , β k \alpha_1,\cdots,\alpha_m,\beta_1,\cdots,\beta_k α1,⋯,αm,β1,⋯,βk是 V V V的正交基
- 因为
p
=
k
+
1
p=k+1
p=k+1时,有
m
<
n
m<n
m<n,一定存在向量
β
\beta
β不能被
A
0
A_0
A0线性表示
- 显然 ∑ j = 1 m k j α j = β \sum_{j=1}^{m}k_j\alpha_j=\beta ∑j=1mkjαj=β无解,即 β ≠ ∑ j = 1 m k j α j \beta\neq{\sum_{j=1}^{m}k_j\alpha_j} β=∑j=1mkjαj
- 可见 β − ∑ j = 1 m k j α j ≠ 0 \beta-\sum_{j=1}^{m}k_j\alpha_j\neq{0} β−∑j=1mkjαj=0
- 构造向量 α m + 1 \bold{\alpha}_{m+1} αm+1= β − ∑ j = 1 m k j α j \beta-\sum_{j=1}^{m}k_j\alpha_j β−∑j=1mkjαj,则 α m + 1 ≠ 0 \alpha_{m+1}\neq{\bold{0}} αm+1=0
- 则
(
α
i
,
α
m
+
1
)
(\alpha_i,\alpha_{m+1})
(αi,αm+1)=
(
β
,
α
i
)
−
k
i
(
α
i
,
α
i
)
(\beta,\alpha_i)-k_i(\alpha_i,\alpha_i)
(β,αi)−ki(αi,αi),
(
i
=
1
,
2
,
⋯
,
m
)
(i=1,2,\cdots,m)
(i=1,2,⋯,m)
- 令 ( α i , α m + 1 ) (\alpha_i,\alpha_{m+1}) (αi,αm+1)= 0 0 0,则 k i = ( β , α i ) ( α i , α i ) k_i=\frac{(\beta,\alpha_i)}{(\alpha_i,\alpha_i)} ki=(αi,αi)(β,αi)(1), ( i = 1 , 2 , ⋯ , m ) (i=1,2,\cdots,m) (i=1,2,⋯,m)
- 这就是说,待定系数 k i k_i ki按照式(1)的方式计算,可得到满足 ( α i , α m + 1 ) = 0 (\alpha_i,\alpha_{m+1})=0 (αi,αm+1)=0, ( i = 1 , 2 , ⋯ , m ) (i=1,2,\cdots,m) (i=1,2,⋯,m),的 α m + 1 \alpha_{m+1} αm+1
- 即 α m + 1 \alpha_{m+1} αm+1和向量组 A 0 A_0 A0的所有向量正交,所以 A 1 : α 1 , ⋯ , α m , α m + 1 A_1:\alpha_1,\cdots,\alpha_m,\alpha_{m+1} A1:α1,⋯,αm,αm+1构成正交向量组
- 由归纳法假设, A 1 A_1 A1可以扩充成正交基
- 所以定理得证
扩充正交向量组到正交基
- 上述定理的证明过程给出了一个具体的扩充正交向量组的方法:
- 即从任意一个非零向量出发,扩充出向量空间的一个正交基
- 只需要按照证明中的步骤逐个扩充,最后得到的正交向量组就是一个正交基
- 对正交基做单位化,就得到一组标准正交基
标准正交基存在定理
- 对于 n n n维欧式空间中任意一组基 A n : a 1 , ⋯ , a n A_n:\bold{a}_1,\cdots,\bold{a}_n An:a1,⋯,an,总可以找到一组标准正交基 B n : b 1 , ⋯ , b n B_n:\bold{b}_1,\cdots,\bold{b}_n Bn:b1,⋯,bn,使得 L ( A i ) = L ( B i ) L(A_i)=L(B_i) L(Ai)=L(Bi), i = 1 , ⋯ , n i=1,\cdots,n i=1,⋯,n
- 证明:
- 首先,取 b 1 = 1 ∣ a 1 ∣ a \bold{b}_1=\frac{1}{|\bold{a}_1|}\bold{a} b1=∣a1∣1a,假设已经求出 B m B_m Bm,它是单位正交的,且: L ( A i ) L(A_i) L(Ai)= L ( B i ) L(B_i) L(Bi), i = 1 , ⋯ , m i=1,\cdots,m i=1,⋯,m
- 再求 b m + 1 \bold{b}_{m+1} bm+1;因为 L ( A m ) = L ( B m ) L(A_m)=L(B_m) L(Am)=L(Bm),且对于 a m + 1 ∈ A n \bold{a}_{m+1}\in{A_n} am+1∈An但 a m + 1 ∉ A m \bold{a}_{m+1}\not\in{A_m} am+1∈Am,所以 a m + 1 \bold{a}_{m+1} am+1不能由 A m , B m A_m,B_m Am,Bm线性表示
- 构造向量
ξ
m
+
1
=
a
m
+
1
−
∑
j
=
1
m
k
j
b
j
\xi_{m+1}=\bold{a}_{m+1}-\sum_{j=1}^{m}k_j\bold{b}_{j}
ξm+1=am+1−∑j=1mkjbj;
- 这里 k j , j = 1 , ⋯ , m k_j,j=1,\cdots,m kj,j=1,⋯,m是待定系数
- 对等式两边同时取内积运算: ( ξ m + 1 , b i ) (\xi_{m+1},\bold{b}_i) (ξm+1,bi)= ( a m + 1 , b i ) − ∑ j = 1 m k j ( b i , b j ) (\bold{a}_{m+1},\bold{b}_{i})-\sum_{j=1}^{m}k_j(\bold{b}_i,\bold{b}_j) (am+1,bi)−∑j=1mkj(bi,bj)= ( a m + 1 , b i ) − k i ( b i , b i ) (\bold{a}_{m+1},\bold{b}_{i})-k_i(\bold{b}_i,\bold{b}_i) (am+1,bi)−ki(bi,bi)
- 令上式为0,可得: k i = ( a m + 1 , b i ) ( b i , b i ) k_i=\frac{(\bold{a}_{m+1},\bold{b}_i)}{(\bold{b}_i,\bold{b}_i)} ki=(bi,bi)(am+1,bi), i = 1 , ⋯ , m i=1,\cdots,m i=1,⋯,m
- 所以
ξ
m
+
1
=
a
m
+
1
−
∑
j
=
1
m
(
a
m
+
1
,
b
j
)
(
b
j
,
b
j
)
b
j
\xi_{m+1}=\bold{a}_{m+1}-\sum_{j=1}^{m}\frac{(\bold{a}_{m+1},\bold{b}_j)}{(\bold{b}_j,\bold{b}_j)}\bold{b}_{j}
ξm+1=am+1−∑j=1m(bj,bj)(am+1,bj)bj时:
- 显然 ξ m + 1 ≠ 0 \xi_{m+1}\neq{0} ξm+1=0且 ( ξ m + 1 , b i ) (\xi_{m+1},\bold{b}_i) (ξm+1,bi)=0, i = 1 , ⋯ , m i=1,\cdots,m i=1,⋯,m
- 这说明 ξ m + 1 \bold{\xi}_{m+1} ξm+1与 b 1 , ⋯ , b m \bold{b}_1,\cdots,\bold{b}_m b1,⋯,bm都正交
- 再作标准化,得到单位向量 b m + 1 \bold{b}_{m+1} bm+1:令 b m + 1 = ξ m + 1 ∣ ξ m + 1 ∣ \bold{b}_{m+1}=\frac{\xi_{m+1}}{|\xi_{m+1}|} bm+1=∣ξm+1∣ξm+1
- b 1 , ⋯ , b m , b m + 1 \bold{b}_1,\cdots,\bold{b}_{m},\bold{b}_{m+1} b1,⋯,bm,bm+1是一单位正交向量组,同时 L ( A m + 1 ) = L ( B m + 1 ) L(A_{m+1})=L(B_{m+1}) L(Am+1)=L(Bm+1)
- 由归纳法原理,定理成立
标准正交化
- 基的标准正交化:设 V 1 V_1 V1是 V V V的一个基,寻找 V V V的另一个与 V 1 V_1 V1等价的标准正交基 V 2 V_2 V2的过程称为基的标准正交化
正交化(schmidt正交化)
-
一个向量组线性无关是该向量组成为正交向量组的必要条件(不充分条件)
- 将一个基中的各个向量单位化是容易的,利用公式 a = a ∣ ∣ a ∣ ∣ \bold{a}=\frac{\bold{a}}{||\bold a||} a=∣∣a∣∣a即可
- 但是将它们正交化比较复杂,这里介绍一种方法能够对一组正交向量组进行正交化
- 本方法的原理参考标准正交基存在定理👺
-
对于一个线性无关组 A A A,可以通过施密特正交化方法,求出一个等价的正交向量组 B B B,
- schmidt方法是一种递推计算的方法,可以将一组基正交化
- 标准化只需要利用正交公式计算即可
- 从几何的角度,正交化的过程是被正交化的向量被正交分解(投影到各个相互垂直的轴上)的过程
- 被正交化的向量维数越高,或向量越多(尤其是向量数量),正交分解的过程就越繁琐,计算量越大
-
设 A : α 1 , α 2 , ⋯ , α n A:\alpha_1,\alpha_2,\cdots,\alpha_n A:α1,α2,⋯,αn
-
令 β 1 = α 1 \beta_1=\alpha_1 β1=α1
- β i = α i − ∑ k = 1 i − 1 ( α i , β k ) ( β k , β k ) β k = α i − ∑ k = 1 i − 1 ( α i , β k ) ∣ ∣ β k ∣ ∣ 2 β k ( i = 2 , ⋯ , n ) \begin{aligned} \beta_{i} =&\alpha_i-\sum_{k=1}^{i-1} \frac{(\alpha_i,\beta_{k})}{(\beta_{k},\beta_{k})}\beta_{k} \\ =&\alpha_i-\sum_{k=1}^{i-1}\frac{(\alpha_i,\beta_{k})}{||\beta_{k}||^2}\beta_{k} \quad (i=2,\cdots,n) \\ \end{aligned} βi==αi−k=1∑i−1(βk,βk)(αi,βk)βkαi−k=1∑i−1∣∣βk∣∣2(αi,βk)βk(i=2,⋯,n)
-
令 B = β 1 , ⋯ , β n B=\beta_1,\cdots,\beta_n B=β1,⋯,βn,则 B B B正交向量组
-
-
正交性证明:(todo)
- β j \beta_j βj= β j = α j − ∑ k = 1 j − 1 ( α j , β k ) ( β k , β k ) β k \beta_{j} =\alpha_j-\sum_{k=1}^{j-1} \frac{(\alpha_j,\beta_{k})}{(\beta_{k},\beta_{k})}\beta_{k} βj=αj−∑k=1j−1(βk,βk)(αj,βk)βk
- 为了便于讨论和书写,令:
- F i = ∑ k = 1 i − 1 ( α i , β k ) ∣ ∣ β k ∣ ∣ 2 β k ( i = 2 , ⋯ , n ) F_i=\sum_{k=1}^{i-1}\frac{(\alpha_i,\beta_{k})}{||\beta_{k}||^2}\beta_{k} \quad (i=2,\cdots,n) Fi=∑k=1i−1∣∣βk∣∣2(αi,βk)βk(i=2,⋯,n)
- β i = α i − F i \beta_i=\alpha_i-F_i βi=αi−Fi
- ( β i , β j ) (\beta_i,\beta_j) (βi,βj)= ( α i , α j ) − ( α i , F j ) − ( F i , α j ) + ( F i , F j ) (\alpha_i,\alpha_j)-(\alpha_i,F_j)-(F_i,\alpha_j)+(F_i,F_j) (αi,αj)−(αi,Fj)−(Fi,αj)+(Fi,Fj)=0
- 或者,令 B m : β 1 , ⋯ , β m B_m:\beta_1,\cdots,\beta_m Bm:β1,⋯,βm,依次证明 β m + 1 \beta_{m+1} βm+1与正交向量组 B m B_m Bm正交
例
- 设
a
1
,
a
2
,
a
3
\bold{a}_1,\bold{a}_2,\bold{a}_3
a1,a2,a3分别是
(
1
,
2
,
−
1
)
T
(1,2,-1)^T
(1,2,−1)T,
(
−
1
,
3
,
1
)
T
(-1,3,1)^T
(−1,3,1)T,
(
4
,
−
1
,
0
)
T
(4,-1,0)^T
(4,−1,0)T,试用施密特正交化它们
- 令 F i = ∑ k = 1 i − 1 ( a i , b k ) ∣ b k ∣ 2 b k F_i=\sum_{k=1}^{i-1}\frac{(\bold{a}_i,\bold{b}_k)}{|\bold{b}_k|^2}\bold{b}_k Fi=∑k=1i−1∣bk∣2(ai,bk)bk, i = 2 , ⋯ , r i=2,\cdots,r i=2,⋯,r
- 取 b 1 = a 1 \bold{b}_1=\bold{a}_1 b1=a1= ( 1 , 2 , − 1 ) T (1,2,-1)^T (1,2,−1)T
-
b
2
=
a
2
−
F
2
\bold{b}_2=\bold{a_2}-F_2
b2=a2−F2
- ∣ b 1 ∣ 2 = 6 |\bold{b}_1|^2=6 ∣b1∣2=6
- ( a 2 , b 1 ) = − 1 + 6 − 1 = 4 (\bold{a_2,b_1})=-1+6-1=4 (a2,b1)=−1+6−1=4
- b 2 = ( − 1 , 3 , 1 ) T − − 1 + 6 − 1 1 + 4 + 1 ( 1 , 2 , − 1 ) T \bold{b}_2=(-1,3,1)^T-\frac{-1+6-1}{1+4+1}(1,2,-1)^T b2=(−1,3,1)T−1+4+1−1+6−1(1,2,−1)T= 5 3 ( − 1 , 1 , 1 ) T \frac{5}{3}(-1,1,1)^{T} 35(−1,1,1)T
-
b
3
=
a
3
−
F
3
\bold{b}_3=\bold{a}_3-F_3
b3=a3−F3
- ∣ b 2 ∣ 2 = 25 9 ( 1 + 1 + 1 ) |\bold{b}_2|^2=\frac{25}{9}(1+1+1) ∣b2∣2=925(1+1+1)= 25 3 \frac{25}{3} 325
- ( a 3 , b 1 ) = 4 − 2 = 2 (\bold{a_3,b_1})=4-2=2 (a3,b1)=4−2=2, ( a 3 , b 2 ) = 5 3 ( − 4 − 1 + 0 ) = − 25 3 (\bold{a_3,b_2})=\frac{5}{3}(-4-1+0)=-\frac{25}{3} (a3,b2)=35(−4−1+0)=−325
- F 3 = 2 6 ( 1 , 2 , − 1 ) T F_3=\frac{2}{6}(1,2,-1)^T F3=62(1,2,−1)T+ − 25 3 ⋅ 3 25 ⋅ 5 3 ( − 1 , 1 , 1 ) T -\frac{25}{3}\cdot\frac{3}{25}\cdot\frac{5}{3}(-1,1,1)^{T} −325⋅253⋅35(−1,1,1)T= 1 3 ( ( 1 , 2 , − 1 ) T − 5 ( − 1 , 1 , 1 ) T ) \frac{1}{3}((1,2,-1)^T-5(-1,1,1)^T) 31((1,2,−1)T−5(−1,1,1)T)= 1 3 ( 6 , − 3 , − 6 ) \frac{1}{3}(6,-3,-6) 31(6,−3,−6)= ( 2 , − 1 , − 2 ) (2,-1,-2) (2,−1,−2)
- b 3 = ( 4 , − 1 , 0 ) T − ( 2 , − 1 , − 2 ) T = ( 2 , 0 , 2 ) T \bold{b}_3=(4,-1,0)^T-(2,-1,-2)^T=(2,0,2)^T b3=(4,−1,0)T−(2,−1,−2)T=(2,0,2)T= 2 ( 1 , 0 , 1 ) T 2(1,0,1)^T 2(1,0,1)T
- 标准化(单位化): e 1 = b 1 ∣ b ∣ \bold{e}_1=\frac{\bold{b}_1}{|\bold{b}|} e1=∣b∣b1= 1 6 ( 1 , 2 , − 1 ) T \frac{1}{\sqrt{6}}(1,2,-1)^T 61(1,2,−1)T; e 2 = 3 5 5 3 ( − 1 , 1 , 1 ) T \bold{e}_2=\frac{\sqrt{3}}{5}\frac{5}{3}(-1,1,1)^{T} e2=5335(−1,1,1)T= 3 3 ( − 1 , 1 , 1 ) T \frac{\sqrt{3}}{3}(-1,1,1)^{T} 33(−1,1,1)T; e 3 \bold{e}_3 e3= 1 2 2 2 ( 1 , 0 , 1 ) T \frac{1}{2\sqrt{2}}2(1,0,1)^T 2212(1,0,1)T= 2 2 ( 1 , 0 , 1 ) T \frac{\sqrt{2}}{2}(1,0,1)^T 22(1,0,1)T
解析几何角度解释正交化过程
-
以三维向量为例
-
-
a 1 = b 1 \bold{a_1=b_1} a1=b1
-
b 2 = a 2 − c 2 \bold{b_2=a_2-c_2} b2=a2−c2, c 2 \bold{c}_2 c2是 a 2 \bold{a}_2 a2在 b 1 \bold{b}_1 b1上的投影向量(分量),由分量公式: c 2 \bold{c}_2 c2= Prj b 1 a 2 ⋅ b 1 ∣ b 1 ∣ \text{Prj}_{\bold{b}_1}\bold{a}_2\cdot{\bold{\frac{b_1}{|b_1|}}} Prjb1a2⋅∣b1∣b1= ( a 2 , b 1 ) ∣ b 1 ∣ 2 b 1 \bold{\frac{(a_2,b_1)}{|b_1|^2}b_1} ∣b1∣2(a2,b1)b1
-
b 3 = a 3 − c 3 \bold{b_3=a_3-c_3} b3=a3−c3,而 c 3 \bold{c}_3 c3是 a 3 \bold{a}_3 a3在平行于 b 1 , b 2 \bold{b_1,b_2} b1,b2的平面上的投影向量,由于 b 1 ⊥ b 2 \bold{b_1\perp{b_2}} b1⊥b2,所以 c 3 \bold{c}_3 c3等于分别在 b 1 , b 2 \bold{b_1,b_2} b1,b2上的投影向量 c 31 , c 32 \bold{c_{31},c_{32}} c31,c32之和: c 3 = c 31 + c 32 = ( a 3 , b 1 ) ∣ b 1 ∣ 2 b 1 \bold{c_3=c_{31}+c_{32}=\frac{(a_3,b_1)}{|b_1|^2}b_1} c3=c31+c32=∣b1∣2(a3,b1)b1+ ( a 3 , b 2 ) ∣ b 2 ∣ 2 b 2 \bold{\frac{(a_3,b_2)}{|b_2|^2}b_2} ∣b2∣2(a3,b2)b2
例:扩充正交基
- 已知 a 1 = ( 1 , 1 , 1 ) T \bold{a}_1=(1,1,1)^T a1=(1,1,1)T;求一组非零向量 a 2 , a 3 \bold{a_2,a_3} a2,a3,使 a 1 , a 2 , a 3 \bold{a_1,a_2,a_3} a1,a2,a3两两正交
- 解:
- 两两正交等价于: ( a 1 , a 2 ) = 0 (\bold{a_1,a_2})=0 (a1,a2)=0, ( a 1 , a 3 ) = 0 (\bold{a_1,a_3})=0 (a1,a3)=0, ( a 2 , a 3 ) = 0 (\bold{a_2,a_3})=0 (a2,a3)=0
- 先处理前两个方程,它们可以用同一个线性方程建模: a 1 T x = 0 \bold{a_1^{T}x=0} a1Tx=0,即 ( x 1 + x 2 + x 3 ) = 0 (x_1+x_2+x_3)=0 (x1+x2+x3)=0
- 显然该方程组仅包含一个线性方程,并且系数矩阵的秩为1,基础解系包含的向量个数为 n − r = 3 − 1 = 2 n-r=3-1=2 n−r=3−1=2,自由未知数的个数也是2
- 为自由未知数数组取2组线性无关的值:
q
1
=
(
1
,
0
)
T
,
q
2
=
(
0
,
1
)
T
q_1=(1,0)^T,q_2=(0,1)^T
q1=(1,0)T,q2=(0,1)T,对应的,基础解系取值为:
- ξ 1 = ( − 1 , 1 , 0 ) T \xi_1=(-1,1,0)^T ξ1=(−1,1,0)T; ξ 2 = ( − 1 , 0 , 1 ) T \xi_2=(-1,0,1)^T ξ2=(−1,0,1)T
- 然而
ξ
1
,
ξ
2
\xi_1,\xi_2
ξ1,ξ2此时不是正交的,需要对其正交化:
- 取 c 1 = ξ 1 \bold{c}_1=\xi_1 c1=ξ1
- c 2 \bold{c}_2 c2= ξ 2 − ( ξ 2 , c 1 ) ∣ c 1 ∣ 2 c 1 \xi_2-\frac{(\xi_2,c_1)}{|c_1|^2}c_1 ξ2−∣c1∣2(ξ2,c1)c1= ( − 1 , 0 , 1 ) T − 1 2 ( − 1 , 1 , 0 ) T (-1,0,1)^T-\frac{1}{2}(-1,1,0)^T (−1,0,1)T−21(−1,1,0)T= ( − 1 2 , − 1 2 , 1 ) T (-\frac{1}{2},-\frac{1}{2},1)^T (−21,−21,1)T= 1 2 ( − 1 , − 1 , 2 ) T \frac{1}{2}(-1,-1,2)^T 21(−1,−1,2)T
- 显然, c 1 , c 2 \bold{c_1,c_2} c1,c2是基础解系的线性组合,它们仍然满足 a 1 T x = 0 \bold{a_1^{T}x=0} a1Tx=0
- 所以可以取 a 2 = ( − 1 , 1 , 0 ) T \bold{a}_2=(-1,1,0)^T a2=(−1,1,0)T, a 3 = 1 2 ( − 1 , − 1 , 2 ) T \bold{a}_3=\frac{1}{2}(-1,-1,2)^T a3=21(−1,−1,2)T
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了