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=nm作数学归纳法
    • 首先, 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=a11a,假设已经求出 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+1An a m + 1 ∉ A m \bold{a}_{m+1}\not\in{A_m} am+1Am,所以 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+1j=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+1j=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==αik=1i1(βk,βk)(αi,βk)βkαik=1i1∣∣βk2(α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=αjk=1j1(β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=1i1∣∣βk2(αi,βk)βk(i=2,,n)
      • β i = α i − F i \beta_i=\alpha_i-F_i βi=αiFi
      • ( β 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=1i1bk2(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=a2F2
    • ∣ b 1 ∣ 2 = 6 |\bold{b}_1|^2=6 b12=6
    • ( a 2 , b 1 ) = − 1 + 6 − 1 = 4 (\bold{a_2,b_1})=-1+6-1=4 (a2,b1)=1+61=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)T1+4+11+61(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=a3F3
    • ∣ b 2 ∣ 2 = 25 9 ( 1 + 1 + 1 ) |\bold{b}_2|^2=\frac{25}{9}(1+1+1) b22=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)=42=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(41+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} 32525335(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)T5(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=bb1= 1 6 ( 1 , 2 , − 1 ) T \frac{1}{\sqrt{6}}(1,2,-1)^T 6 1(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=53 35(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 22 12(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=a2c2, 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|}}} Prjb1a2b1b1= ( a 2 , b 1 ) ∣ b 1 ∣ 2 b 1 \bold{\frac{(a_2,b_1)}{|b_1|^2}b_1} b12(a2,b1)b1

  • b 3 = a 3 − c 3 \bold{b_3=a_3-c_3} b3=a3c3,而 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}} b1b2,所以 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=b12(a3,b1)b1+ ( a 3 , b 2 ) ∣ b 2 ∣ 2 b 2 \bold{\frac{(a_3,b_2)}{|b_2|^2}b_2} b22(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 nr=31=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 ξ2c12(ξ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)T21(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
posted @   xuchaoxin1375  阅读(24)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
点击右上角即可分享
微信分享提示