• 【矩阵论】1.准备知识(下)


    矩阵论的所有文章,主要内容参考北航赵迪老师的课件

    [注]由于矩阵论对计算机比较重要,所以选修了这门课,但不是专业搞数学的,所以存在很多口语化描述,而且对很多东西理解不是很正确与透彻,欢迎大家指正。我可能间歇性忙,但有空一定会回复修改的。

    矩阵论
    1. 准备知识——复数域上矩阵,Hermite变换)
    1.准备知识——复数域上的内积域正交阵
    1.准备知识——Hermite阵,二次型,矩阵合同,正定阵,幂0阵,幂等阵,矩阵的秩
    2. 矩阵分解——SVD准备知识——奇异值
    2. 矩阵分解——SVD
    2. 矩阵分解——QR分解
    2. 矩阵分解——正定阵分解
    2. 矩阵分解——单阵谱分解
    2. 矩阵分解——正规分解——正规阵
    2. 矩阵分解——正规谱分解
    2. 矩阵分解——高低分解
    3. 矩阵函数——常见解析函数
    3. 矩阵函数——谱公式,幂0与泰勒计算矩阵函数
    3. 矩阵函数——矩阵函数求导
    4. 矩阵运算——观察法求矩阵特征值特征向量
    4. 矩阵运算——张量积
    4. 矩阵运算——矩阵拉直
    4.矩阵运算——广义逆——加号逆定义性质与特殊矩阵的加号逆
    4. 矩阵运算——广义逆——加号逆的计算
    4. 矩阵运算——广义逆——加号逆应用
    4. 矩阵运算——广义逆——减号逆
    5. 线性空间与线性变换——线性空间
    5. 线性空间与线性变换——生成子空间
    5. 线性空间与线性变换——线性映射与自然基分解,线性变换
    6. 正规方程与矩阵方程求解
    7. 范数理论——基本概念——向量范数与矩阵范数
    7.范数理论——基本概念——矩阵范数生成向量范数&谱范不等式
    7. 矩阵理论——算子范数
    7.范数理论——范数估计——许尔估计&谱估计
    7. 范数理论——非负/正矩阵
    8. 常用矩阵总结——秩1矩阵,优阵(单位正交阵),Hermite阵
    8. 常用矩阵总结——镜面阵,正定阵
    8. 常用矩阵总结——单阵,正规阵,幂0阵,幂等阵,循环阵


    在这里插入图片描述

    1.6 内积

    X = ( x 1 x 2 ⋮ x n ) , Y = ( y 1 y 2 ⋮ y n ) ∈ C n X=\left(

    x1x2xn" role="presentation">x1x2xn
    \right),Y=\left(
    y1y2yn" role="presentation">y1y2yn
    \right)\in C^n X= x1x2xn ,Y= y1y2yn Cn

    1.6.1 复向量内积

    ( X , Y ) = Δ Y H X = x 1 y 1 ‾ + x 2 y 2 ‾ + ⋯ + x n y n ‾ = ∑ i = 1 n x i y i ‾ = t r ( Y H X ) (X,Y)\overset{\Delta}{=}Y^HX=x_1\overline{y_1}+x_2\overline{y_2}+\cdots+x_n\overline{y_n}=\sum\limits_{i=1}\limits^{n}x_i\overline{y_i}=tr(Y^HX) (X,Y)=ΔYHX=x1y1+x2y2++xnyn=i=1nxiyi=tr(YHX)

    ( Y , X ) = X H Y = x 1 ‾ y 1 + x 2 ‾ y 2 + ⋯ + x n ‾ y 1 = ∑ i = 1 n x i ‾ y i = t r ( X H Y ) (Y,X)=X^HY=\overline{x_1}y_1+\overline{x_2}y_2+\cdots+\overline{x_n}y_1=\sum\limits_{i=1}\limits^{n}\overline{x_i}y_i=tr(X^HY) (Y,X)=XHY=x1y1+x2y2++xny1=i=1nxiyi=tr(XHY)

    若取 Y = X Y=X Y=X ,则其内积
    ( X , X ) = X H X = x 1 x 1 ‾ + x 2 x 2 ‾ + ⋯ + x n x n ‾ = ∑ i = 1 n x i x i ‾ = ∣ x 1 ∣ 2 + ∣ x 2 ∣ 2 + ⋯ + ∣ x n ∣ 2 = ∣ X ∣ 2 = t r ( X H X ) = t r ( X X H )

    (X,X)=XHX=x1x1¯+x2x2¯++xnxn¯=i=1nxixi¯=|x1|2+|x2|2++|xn|2=|X|2=tr(XHX)=tr(XXH)" role="presentation">(X,X)=XHX=x1x1¯+x2x2¯++xnxn¯=i=1nxixi¯=|x1|2+|x2|2++|xn|2=|X|2=tr(XHX)=tr(XXH)
    (X,X)=XHX=x1x1+x2x2++xnxn=i=1nxixi=x12+x22++xn2=X2=tr(XHX)=tr(XXH)

    向量内积性质
    1. ( X , X ) ≥ 0 (X,X)\ge 0 (X,X)0 ;若 X ≠ 0 , ( X , X ) > 0 X\neq 0 ,(X,X)> 0 X=0,(X,X)>0
    2. ( X , Y ) = ( Y , X ) ‾ (X,Y)=\overline{(Y,X)} (X,Y)=(Y,X)
    3. ( k X , Y ) = k ( X , Y ) , ( X , k Y ) = k ‾ ( X , Y ) (kX,Y)=k(X,Y),(X,kY)=\overline{k}(X,Y) (kX,Y)=k(X,Y),(X,kY)=k(X,Y)
    4. ( X + Y , W ) = ( X , W ) + ( Y , W ) (X+Y,W)=(X,W)+(Y,W) (X+Y,W)=(X,W)+(Y,W)
      ( X , Y + W ) = ( X , Y ) + ( X , W ) (X,Y+W)=(X,Y)+(X,W) (X,Y+W)=(X,Y)+(X,W)
    5. ∣ ( X , Y ) ∣ 2 ≤ ∣ X ∣ ⋅ ∣ Y ∣ \vert (X,Y) \vert^2\le \vert X \vert\cdot\vert Y \vert (X,Y)2XY

    1.6.2 复矩阵内积

    a. 定义

    ( A , B ) = Δ t r ( B H A ) = t r ( A H B ) = ∑ a i j b i j ‾ , A , B ∈ C m , n ( A , A ) = Δ t r ( A H A ) = t r ( A A H ) = ∑ a i j a i j ‾ = ∑ ∣ a i j ∣ 2

    (A,B)=Δtr(BHA)=tr(AHB)=aijbij¯,A,BCm,n(A,A)\xlongequalΔtr(AHA)=tr(AAH)=aijaij¯=|aij|2" role="presentation">(A,B)=Δtr(BHA)=tr(AHB)=aijbij¯,A,BCm,n(A,A)\xlongequalΔtr(AHA)=tr(AAH)=aijaij¯=|aij|2
    (A,B)=Δtr(BHA)=tr(AHB)=aijbij,A,BCm,n(A,A)Δ tr(AHA)=tr(AAH)=aijaij=aij2

    矩阵A的模长: ∣ ∣ A ∣ ∣ = ( A , A ) = t r ( A A H ) = ∑ ∣ a i j ∣ 2 \vert \vert A \vert \vert=\sqrt{(A,A)}=\sqrt{tr(AA^H)}=\sqrt{\sum\vert a_{ij} \vert^2} ∣∣A∣∣=(A,A) =tr(AAH) =aij2

    在这里插入图片描述

    b. 性质
    1. ( A , A ) = t r ( A H A ) = ∑ ∣ a i j ∣ 2 ≥ 0 (A,A)=tr(A^HA)=\sum\vert a_{ij} \vert^2 \ge 0 (A,A)=tr(AHA)=aij20 ;若 A ≠ 0 A\neq 0 A=0 ,则 ( A , A ) > 0 (A,A)>0 (A,A)>0

    2. ( A , B ) = ( B , A ) ‾ (A,B)=\overline{(B,A)} (A,B)=(B,A)

    3. ( k A , B ) = k ( A , B ) , ( A , k B ) = k ‾ ( A , B ) (kA,B)=k(A,B),(A,kB)=\overline{k}(A,B) (kA,B)=k(A,B),(A,kB)=k(A,B)

    4. (A+B,D)=(A,D)+(B,D),(D,A+B)=(D,A)+(D,B)

    5. ∣ ( A , B ) ∣ 2 ≤ ∣ A ∣ ⋅ ∣ B ∣ \vert (A,B) \vert^2 \le \vert A\vert\cdot\vert B \vert (A,B)2AB

    c. 矩阵的内积形式
    列分块(常用)

    A = ( a 11 ⋯ a 1 p ⋮ ⋱ ⋮ a n 1 ⋯ a n p ) ∈ C n × p = ( α 1 , ⋯   , α p ) , 其中 α i 为 n 维列向量 ( n × 1 阶矩阵 ) A H = ( a 11 ‾ ⋯ a n 1 ‾ ⋮ ⋱ ⋮ a 1 p ‾ ⋯ a n p ‾ ) ∈ C p × n = ( α 1 ‾ T ⋮ α p ‾ T ) ,其中 α 1 ‾ T 是 n 维行向量 ( 1 × n 阶矩阵 )

    A=(a11a1pan1anp)Cn×p=(α1,,αp),αin(n×1)AH=(a11¯an1¯a1p¯anp¯)Cp×n=(α1¯Tαp¯T)α1¯Tn(1×n)" role="presentation">A=(a11a1pan1anp)Cn×p=(α1,,αp),αin(n×1)AH=(a11¯an1¯a1p¯anp¯)Cp×n=(α1¯Tαp¯T)α1¯Tn(1×n)
    AAH= a11an1a1panp Cn×p=(α1,,αp),其中αin维列向量(n×1阶矩阵)= a11a1pan1anp Cp×n= α1TαpT ,其中α1Tn维行向量(1×n阶矩阵)

    A H A = ( α 1 ‾ T ⋮ α p ‾ T ) ( α 1 , ⋯   , α p ) = ( α 1 ‾ T α 1 α 1 ‾ T α 2 ⋯ α 1 ‾ T α p α 2 ‾ T α 1 α 2 ‾ T α 2 ⋯ α 2 ‾ T α p ⋮ ⋮ ⋱ ⋮ α p ‾ T α 1 α p ‾ T α 2 ⋯ α p ‾ T α p ) = ( ( α 1 , α 1 ) ( α 2 , α 1 ) ⋯ ( α p , α 1 ) ( α 1 , α 2 ) ( α 2 , α 2 ) ⋯ ( α p , α 2 ) ⋮ ⋮ ⋱ ⋮ ( α 1 , α p ) ( α 2 , α p ) ⋯ ( α p , α p ) ) = ( ( α 1 , α 1 ) ‾ ( α 1 , α 2 ) ‾ ⋯ ( α 1 , α p ) ‾ ( α 2 , α 1 ) ‾ ( α 2 , α 2 ) ‾ ⋯ ( α 2 , α p ) ‾ ⋮ ⋮ ⋱ ⋮ ( α p , α 1 ) ‾ ( α p , α 2 ) ‾ ⋯ ( α p , α p ) ‾ )

    AHA=(α1¯Tαp¯T)(α1,,αp)=(α1¯Tα1α1¯Tα2α1¯Tαpα2¯Tα1α2¯Tα2α2¯Tαpαp¯Tα1αp¯Tα2αp¯Tαp)=((α1,α1)(α2,α1)(αp,α1)(α1,α2)(α2,α2)(αp,α2)(α1,αp)(α2,αp)(αp,αp))=((α1,α1)¯(α1,α2)¯(α1,αp)¯(α2,α1)¯(α2,α2)¯(α2,αp)¯(αp,α1)¯(αp,α2)¯(αp,αp)¯)" role="presentation">AHA=(α1¯Tαp¯T)(α1,,αp)=(α1¯Tα1α1¯Tα2α1¯Tαpα2¯Tα1α2¯Tα2α2¯Tαpαp¯Tα1αp¯Tα2αp¯Tαp)=((α1,α1)(α2,α1)(αp,α1)(α1,α2)(α2,α2)(αp,α2)(α1,αp)(α2,αp)(αp,αp))=((α1,α1)¯(α1,α2)¯(α1,αp)¯(α2,α1)¯(α2,α2)¯(α2,αp)¯(αp,α1)¯(αp,α2)¯(αp,αp)¯)
    AHA= α1TαpT (α1,,αp)= α1Tα1α2Tα1αpTα1α1Tα2α2Tα2αpTα2α1Tαpα2TαpαpTαp = (α1,α1)(α1,α2)(α1,αp)(α2,α1)(α2,α2)(α2,αp)(αp,α1)(αp,α2)(αp,αp) = (α1,α1)(α2,α1)(αp,α1)(α1,α2)(α2,α2)(αp,α2)(α1,αp)(α2,αp)(αp,αp)

    行分块

    A = ( a 11 ⋯ a 1 p ⋮ ⋱ ⋮ a n 1 ⋯ a n p ) ∈ C n × p = ( α 1 α 2 ⋮ α n ) , 其中 α i 为 p 维行向量 ( 1 × p 阶矩阵 ) A H = ( α 1 ‾ T α 2 ‾ T ⋯ α n ‾ T ) , 其中 α i ‾ T 为 p 维列向量 ( p × 1 阶矩阵 )

    A=(a11a1pan1anp)Cn×p=(α1α2αn),αip(1×p)AH=(α1¯Tα2¯Tαn¯T),αi¯Tp(p×1)" role="presentation">A=(a11a1pan1anp)Cn×p=(α1α2αn),αip(1×p)AH=(α1¯Tα2¯Tαn¯T),αi¯Tp(p×1)
    AAH= a11an1a1panp Cn×p= α1α2αn ,其中αip维行向量(1×p阶矩阵)=(α1Tα2TαnT),其中αiTp维列向量(p×1阶矩阵)

    A A H = ( α 1 α 2 ⋮ α n ) ( α 1 ‾ T α 2 ‾ T ⋯ α n ‾ T ) = ( α 1 α 1 ‾ T α 1 α 2 ‾ T ⋯ α 1 α n ‾ T α 2 α 1 ‾ T α 2 α 2 ‾ T ⋯ α 2 α n ‾ T ⋮ ⋮ ⋱ ⋮ α n α 1 ‾ T α n α 2 ‾ T ⋯ α n α n ‾ T ) = ( ( α 1 , α 1 ) ( α 1 , α 2 ) ⋯ ( α 1 , α n ) ( α 2 , α 1 ) ( α 2 , α 2 ) ⋯ ( α 2 , α n ) ⋮ ⋮ ⋱ ⋮ ( α n , α 1 ) ( α n , α 2 ) ⋯ ( α n , α n ) )

    AAH=(α1α2αn)(α1¯Tα2¯Tαn¯T)=(α1α1¯Tα1α2¯Tα1αn¯Tα2α1¯Tα2α2¯Tα2αn¯Tαnα1¯Tαnα2¯Tαnαn¯T)=((α1,α1)(α1,α2)(α1,αn)(α2,α1)(α2,α2)(α2,αn)(αn,α1)(αn,α2)(αn,αn))" role="presentation">AAH=(α1α2αn)(α1¯Tα2¯Tαn¯T)=(α1α1¯Tα1α2¯Tα1αn¯Tα2α1¯Tα2α2¯Tα2αn¯Tαnα1¯Tαnα2¯Tαnαn¯T)=((α1,α1)(α1,α2)(α1,αn)(α2,α1)(α2,α2)(α2,αn)(αn,α1)(αn,α2)(αn,αn))
    AAH= α1α2αn (α1Tα2TαnT)= α1α1Tα2α1Tαnα1Tα1α2Tα2α2Tαnα2Tα1αnTα2αnTαnαnT = (α1,α1)(α2,α1)(αn,α1)(α1,α2)(α2,α2)(αn,α2)(α1,αn)(α2,αn)(αn,αn)

    1.6.3 模长

    a. 复向量模长
    列向量模长

    X = ( x 1 x 2 ⋮ x n ) ∈ C , 则其模长 ∣ X ∣ = ∣ x 1 ∣ 2 + ∣ x 2 ∣ 2 + ⋯ ∣ x n ∣ 2

    X=(x1x2xn)C,X∣=x12+x22+xn2" role="presentation">X=(x1x2xn)C,X∣=x12+x22+xn2
    X= x1x2xn C,则其模长X∣=x12+x22+xn2

    模长性质

    ∣ k ⋅ X ∣ = ∣ k ∣ ⋅ ∣ X ∣ , k ∈ C \mid k\cdot X \mid = \mid k \mid \cdot \mid X \mid ,k\in C kX∣=∣kX,kC

    ∣ X k ∣ = ∣ X ∣ ∣ k ∣ , ( k ≠ 0 ⃗ ) \vert \frac{X}{k} \vert=\frac{\vert X\vert}{\vert k \vert},(k\neq \vec{0}) kX=kX,(k=0 )

    ∣ X ± Y ∣ ≤ ∣ X ∣ + ∣ Y ∣ \vert X\pm Y\vert \le \vert X\vert + \vert Y \vert X±YX+Y

    模平方公式

    令 X = ( x 1 x 2 ⋮ x n ) ∈ C , 则 { ① X H X = ∣ X ∣ 2 ② t r ( X H X ) = t r ( X X H ) = ∣ X ∣ 2 , 其中 ∣ X ∣ 2 = ∣ x 1 ∣ 2 + ∣ x 2 ∣ 2 + ⋯ + ∣ x n ∣ 2 = ∑ i = 1 n ∣ x i ∣ 2

    \begin{aligned} &令X=\left( \begin{matrix} x_1\\ x_2\\ \vdots\\ x_n \end{matrix} \right)\in C,则 \left\{ \begin{aligned} &① X^HX=\mid X \mid^2\\\\ &②tr(X^HX)=tr(XX^H)=\mid X \mid^2 \end{aligned}" role="presentation">\begin{aligned} &令X=\left( \begin{matrix} x_1\\ x_2\\ \vdots\\ x_n \end{matrix} \right)\in C,则 \left\{ \begin{aligned} &① X^HX=\mid X \mid^2\\\\ &②tr(X^HX)=tr(XX^H)=\mid X \mid^2 \end{aligned}
    \right.,其中\\ &\mid X \mid^2 = \mid x_1 \mid ^2 + \mid x_2 \mid^2+\cdots+\mid x_n \mid^2=\sum\limits_{i=1}\limits^{n}\mid x_i \mid^2 \end{aligned} X= x1x2xn C, XHX=∣X2tr(XHX)=tr(XXH)=∣X2,其中X2=∣x12+x22++xn2=i=1nxi2

    区分 :复数的模平方和复数平方的模

    • 复数的模平方 ∣ x 1 ∣ 2 = ( a + b i ) ( a − b i ) = a 2 + b 2 ≥ 0 \mid x_1 \mid^2=(a+bi)(a-bi)=a^2+b^2 \ge 0 x12=(a+bi)(abi)=a2+b20
    • 复数平方的模: ∣ x 1 2 ∣ = ( a + b i ) 2 = a 2 − b 2 + 2 a b i \mid x_1^2 \mid=(a+bi)^2= a^2-b^2+2abi x12∣=(a+bi)2=a2b2+2abi

    eg
    X = ( 1 i 1 ) ∈ C 3 , X X H = ( 1 i 1 ) ( 1 − i 1 ) = ( 1 ∗ 1 ∗ 1 ) ∴ t r ( X H X ) = t r ( X X H ) = ∣ 1 ∣ 2 + ∣ − i 2 ∣ 2 + ∣ 1 ∣ 2 = 3

    X=(1i1)C3,XXH=(1i1)(1i1)=(111)tr(XHX)=tr(XXH)=|1|2+|i2|2+|1|2=3" role="presentation">X=(1i1)C3,XXH=(1i1)(1i1)=(111)tr(XHX)=tr(XXH)=|1|2+|i2|2+|1|2=3
    X= 1i1 C3,XXH= 1i1 (1i1)= 111 tr(XHX)=tr(XXH)=∣12+i22+∣12=3

    b. 复矩阵模长

    令 A = ( a i j ) ∈ C , 则模长 ( 矩阵的 F 范数 ) ∣ ∣ A ∣ ∣ = ∑ i = 1 n ∣ a i j ∣ 2

    A=(aij)C,(F)∣∣A∣∣=i=1naij2" role="presentation">A=(aij)C,(F)∣∣A∣∣=i=1naij2
    A=(aij)C,则模长(矩阵的F范数)∣∣A∣∣=i=1naij2

    复矩阵的模平方公式

    设 A = ( a i j ) n × p , 则 t r ( A H A ) = t r ( A A H ) = ∣ ∣ A ∣ ∣ 2 = ∑ i = 1 n ∣ a i j ∣ 2

    A=(aij)n×p,tr(AHA)=tr(AAH)=∣∣A2=i=1naij2" role="presentation">A=(aij)n×p,tr(AHA)=tr(AAH)=∣∣A2=i=1naij2
    A=(aij)n×p,tr(AHA)=tr(AAH)=∣∣A2=i=1naij2

    eg
    X = ( 1 i 1 ) ∈ C 3 × 1 , 则 ∣ X ∣ 2 = X H X = ( 1 , − i , 1 ) ( 1 i 1 ) = 1 + − i 2 + 1 = 3

    X=(1i1)C3×1,X2=XHX=(1,i,1)(1i1)=1+i2+1=3" role="presentation">X=(1i1)C3×1,X2=XHX=(1,i,1)(1i1)=1+i2+1=3
    X= 1i1 C3×1,X2=XHX=(1,i,1) 1i1 =1+i2+1=3
    对于方阵 A = ( a i j ) n × n A=(a_{ij})_{n\times n} A=(aij)n×n ,有
    { t r ( A ) = Δ a 11 + a 22 + ⋯ + a n n = ∑ i = 1 n λ i d e t ( A ) = ∣ A ∣ = ∏ i = 1 n λ i
    \begin{aligned} \left\{ \begin{aligned} &tr(A)\overset{\Delta}{=} a_{11}+a_{22}+\cdots+a_{nn}=\sum\limits_{i=1}\limits^{n}\lambda_i\\ &det(A) = \mid A \mid = \prod\limits_{i=1}\limits^{n}\lambda_i \end{aligned}" role="presentation">\begin{aligned} \left\{ \begin{aligned} &tr(A)\overset{\Delta}{=} a_{11}+a_{22}+\cdots+a_{nn}=\sum\limits_{i=1}\limits^{n}\lambda_i\\ &det(A) = \mid A \mid = \prod\limits_{i=1}\limits^{n}\lambda_i \end{aligned}
    \right. \end{aligned}
    tr(A)=Δa11+a22++ann=i=1nλidet(A)=∣A∣=i=1nλi

    1.7 正交

    1.7.1 向量正交

    X = ( x 1 x 2 ⋮ x n ) , Y = ( y 1 y 2 ⋮ y n ) ∈ C n

    X=(x1x2xn),Y=(y1y2yn)Cn" role="presentation">X=(x1x2xn),Y=(y1y2yn)Cn
    X= x1x2xn ,Y= y1y2yn Cn

    X ⊥ Y    ⟺    ( X , Y ) = 0 = x 1 y 1 ‾ + x 2 y 2 ‾ + ⋯ + x n y n ‾ = x 1 ‾ y 1 + x 2 ‾ y 2 + ⋯ + x n ‾ y n ‾ = ( Y , X )

    XY(X,Y)=0=x1y1¯+x2y2¯++xnyn¯=x1¯y1+x2¯y2++xn¯yn¯=(Y,X)" role="presentation">XY(X,Y)=0=x1y1¯+x2y2¯++xnyn¯=x1¯y1+x2¯y2++xn¯yn¯=(Y,X)
    XY(X,Y)=0=(Y,X)=x1y1+x2y2++xnyn=x1y1+x2y2++xnyn

    正交性质
    1. X ⊥ Y ⇒ a X ⊥ b Y X\bot Y\Rightarrow aX\bot bY XYaXbY

      证: ( a X , b Y ) = b ‾ Y H a X = a b ‾ Y H X = a b ‾ ( X , Y ) = 0 (aX,bY)=\overline{b}Y^HaX=a\overline{b}Y^HX=a\overline{b}(X,Y)=0 (aX,bY)=bYHaX=abYHX=ab(X,Y)=0

    2. 勾股定理: X 1 ⊥ X 2 ⊥ ⋯ ⊥ X n ⇒ ∣ c 1 X 1 ± c 2 X 2 ± ⋯ ± c n X n ∣ 2 = ∣ c 1 X 1 ∣ 2 + ∣ c 2 X 2 ∣ 2 + ⋯ + ∣ c n X n ∣ 2 X_1\bot X_2\bot \cdots \bot X_n\Rightarrow \vert c_1X_1\pm c_2X_2\pm \cdots \pm c_nX_n\vert^2=\vert c_1X_1\vert^2+\vert c_2X_2\vert^2+\cdots+\vert c_nX_n\vert^2 X1X2Xnc1X1±c2X2±±cnXn2=c1X12+c2X22++cnXn2

      此时, X 1 , X 2 , ⋯   , x n X_1,X_2,\cdots,x_n X1,X2,,xn 称为一个正交组

    1.7.2 单位向量

    X ≠ 0 ⃗ X\neq \vec{0} X=0 X ∣ X ∣ \frac{X}{\vert X \vert} XX 是一个单位向量( ∣ X ∣ X ∣ ∣ = 1 \vert \frac{X}{\vert X \vert} \vert=1 XX=1

    在这里插入图片描述

    1.7.3 优阵(正交阵)

    预:非单位列向量

    半:p个n维列向量(p

    a. 预-半优阵(预-半正交阵)

    α 1 , α 2 , ⋯   , α p \alpha_1,\alpha_2,\cdots,\alpha_p α1,α2,,αp n n n 维列向量,且 p ≤ n p\le n pn ,且 α 1 ⊥ α 2 ⊥ ⋯ ⊥ α p \alpha_1\bot\alpha_2\bot \cdots\bot\alpha_p α1α2αp
    ,则称 A = ( α 1 , α 2 , ⋯   , α p ) A=(\alpha_1,\alpha_2,\cdots,\alpha_p) A=(α1,α2,,αp) 为预半优阵

    判定

    A = ( α 1 , α 2 , ⋯   , α p ) A=(\alpha_1,\alpha_2,\cdots,\alpha_p) A=(α1,α2,,αp) 是预半优阵    ⟺    A H A = ( ( α 1 , α 1 ) ⋯ 0 ⋮ ⋱ 0 0 ⋯ ( α p , α p ) ) \iff A^HA=\left(

    (α1,α1)000(αp,αp)" role="presentation">(α1,α1)000(αp,αp)
    \right) AHA= (α1,α1)000(αp,αp) 是对角阵,其中 α 1 , α 2 , ⋯   , α p \alpha_1,\alpha_2,\cdots,\alpha_p α1,α2,,αp n n n 维列向量

    在这里插入图片描述

    区分 A H A A^HA AHA p × p p \times p p×p 阶满秩方阵,而 A A H AA^H AAH n × n n\times n n×n 不满秩方阵

    b. 半优阵(半正交阵)

    A = ( α 1 , α 2 , ⋯   , α p ) A=(\alpha_1,\alpha_2,\cdots,\alpha_p) A=(α1,α2,,αp) 是预半优阵,其中 α i \alpha_i αi n n n 维列向量,若满足 ∣ α 1 ∣ = ∣ α 2 ∣ = ⋯ = ∣ α p ∣ = 1 \vert \alpha_1 \vert=\vert \alpha_2 \vert=\cdots=\vert \alpha_p \vert = 1 α1=α2==αp=1 ,则A为半优阵

    判定

    A = ( α 1 , ⋯   , α p ) A=(\alpha_1,\cdots,\alpha_p) A=(α1,,αp) 是半优阵    ⟺    α 1 ⊥ ⋯ ⊥ α p \iff \alpha_1\bot\cdots\bot\alpha_p α1αp ,且 ∣ α 1 ∣ = ⋯ = ∣ α p ∣ = 1 \vert \alpha_1 \vert=\cdots=\vert \alpha_p \vert=1 α1==αp=1    ⟺    A H A = I p \iff A^HA=I_{p} AHA=Ip

    性质
    1. 保模长 A为半U阵,则 ∣ A x ∣ 2 = ∣ x ∣ 2 \vert Ax \vert^2=\vert x \vert^2 Ax2=x2

      ∣ A x ∣ 2 = ( A x ) H ( A x ) = x H A H A x = ∣ X ∣ 2 \vert Ax \vert^2=(Ax)^H(Ax)=x^HA^HAx=\vert X\vert^2 Ax2=(Ax)H(Ax)=xHAHAx=X2

    2. 保正交 A为半优阵, x ⊥ y x\bot y xy ,则 A x ⊥ A y Ax\bot Ay AxAy

    c. 预-优阵(预-单位正交阵)

    α 1 , α 2 , ⋯   , α n \alpha_1,\alpha_2,\cdots,\alpha_n α1,α2,,αn n n n 维列向量,且 α 1 ⊥ α 2 ⊥ ⋯ ⊥ α n \alpha_1\bot\alpha_2\bot\cdots\bot\alpha_n α1α2αn ,则 A = ( α 1 , α 2 , ⋯   , α n ) A=(\alpha_1,\alpha_2,\cdots,\alpha_n) A=(α1,α2,,αn) 是预优阵

    eg
    X 1 = ( 1 i i ) , X 2 = ( 2 i 1 1 ) , X 3 = ( 0 1 − 1 ) ( X 1 , X 2 ) = 0 , ( X 2 , X 3 ) = 0 , ( X 1 , X 3 ) = 0 , ∴ X 1 ⊥ X 2 ⊥ X 3 , A = ( X 1 , X 2 , X 3 ) 是预 − 优阵

    X1=(1ii),X2=(2i11),X3=(011)(X1,X2)=0,(X2,X3)=0,(X1,X3)=0,X1X2X3,A=(X1,X2,X3)" role="presentation">X1=(1ii),X2=(2i11),X3=(011)(X1,X2)=0,(X2,X3)=0,(X1,X3)=0,X1X2X3,A=(X1,X2,X3)
    X1= 1ii ,X2= 2i11 ,X3= 011 (X1,X2)=0,(X2,X3)=0,(X1,X3)=0,X1X2X3,A=(X1,X2,X3)是预优阵

    判定

    A = ( α 1 , α 2 , ⋯   , α n )    ⟺    A H A = ( ( α 1 , α 1 ) ⋯ 0 ⋮ ⋱ 0 0 ⋯ ( α n , α n ) ) 是对角阵 其中, α 1 , α 2 , ⋯   , α p 是 n 维列向量

    A=(α1,α2,,αn)AHA=((α1,α1)000(αn,αn))α1,α2,,αpn" role="presentation">A=(α1,α2,,αn)AHA=((α1,α1)000(αn,αn))α1,α2,,αpn
    A=(α1,α2,,αn)AHA= (α1,α1)000(αn,αn) 是对角阵其中,α1,α2,,αpn维列向量

    在这里插入图片描述

    d. 优阵(正交阵)

    α 1 , α 2 , ⋯   , α n \alpha_1,\alpha_2,\cdots,\alpha_n α1,α2,,αn n n n 维列向量, α 1 ⊥ α 2 ⊥ ⋯ ⊥ α n \alpha_1\bot\alpha_2\bot\cdots\bot\alpha_n α1α2αn ∣ α 1 ∣ = ⋯ = ∣ α n ∣ = 1 \vert \alpha_1 \vert=\cdots=\vert \alpha_n \vert=1 α1==αn=1 ,则 A A A 是一个优阵(正交阵)

    性质
    • A = A n × n A=A_{n\times n} A=An×n 为优阵( A H A = I A^HA=I AHA=I),即 A A A 的列向量 α 1 , α 2 , ⋯   , α n \alpha_1,\alpha_2,\cdots,\alpha_n α1,α2,,αn 为单位正交向量组
    • A − 1 = A H A^{-1}=A^H A1=AH
    • A H A = I , 且 A A H = I A^HA=I,且AA^H=I AHA=I,AAH=I
    判定
    1. A A A 是优阵    ⟺    A H A = I    ⟺    A − 1 A = I    ⟺    A A H = I \iff A^HA=I\iff A^{-1}A=I\iff AA^H=I AHA=IA1A=IAAH=I    ⟺    A = ( α 1 , α 2 , ⋯   , α n ) \iff A=(\alpha_1,\alpha_2,\cdots,\alpha_n) A=(α1,α2,,αn) ,且 α 1 ⊥ α 2 , ⋯ ⊥ α n \alpha_1\bot\alpha_2,\cdots\bot\alpha_n α1α2,αn ∣ α 1 ∣ = ⋯ = ∣ α n ∣ = 1 \vert \alpha_1\vert=\cdots=\vert\alpha_n\vert=1 α1==αn=1

    2. ∣ A x ∣ 2 = ∣ x ∣ 2 \vert Ax\vert^2=\vert x \vert^2 Ax2=x2 A A A 是优阵
      ∵ ∣ A x ∣ 2 = ( A x ) H ( A x ) = x H A H A x = x H I x = ( x , x ) = ∣ x ∣ 2 \because \vert Ax\vert^2 = (Ax)^H(Ax)=x^HA^HAx=x^HIx=(x,x)=\vert x \vert^2 Ax2=(Ax)H(Ax)=xHAHAx=xHIx=(x,x)=x2

    3. x ⊥ y ⇒ A x ⊥ A y x\bot y \Rightarrow Ax\bot Ay xyAxAy A A A 是优阵
      ∵ ( A x , A y ) = ( A y ) H A x = y H A H A x = ( x , y ) = 0    ⟺    A x ⊥ A y \because (Ax,Ay)=(Ay)^HAx=y^HA^HAx=(x,y)=0\iff Ax\bot Ay (Ax,Ay)=(Ay)HAx=yHAHAx=(x,y)=0AxAy

    4. ( A x , A y ) = ( x , y ) (Ax,Ay)=(x,y) (Ax,Ay)=(x,y) A A A 是优阵

    优阵构造

    预优阵到优阵
    A = ( α 1 , ⋯   , α n ) 是预优阵 ⇒ A = ( α 1 ∣ α 1 ∣ , α 2 ∣ α 2 ∣ , ⋯   , α n ∣ α n ∣ ) 是优阵

    A=(α1,,αn)A=(α1|α1|,α2|α2|,,αn|αn|)" role="presentation">A=(α1,,αn)A=(α1|α1|,α2|α2|,,αn|αn|)
    A=(α1,,αn)是预优阵A=(α1α1,α2α2,,αnαn)是优阵

    在这里插入图片描述

    优阵到优阵

    A = ( α 1 , α 2 , ⋯   , α n ) A=(\alpha_1,\alpha_2,\cdots,\alpha_n) A=(α1,α2,,αn) 为优阵,则

    1. k = ± 1 , k A = ( k α 1 , k α 2 , ⋯   , k α n ) k=\pm1,kA=(k\alpha_1,k\alpha_2,\cdots,k\alpha_n) k=±1,kA=(kα1,kα2,,kαn) 为优阵
    2. B = ( β 1 , β 2 , ⋯   , β n ) B=(\beta_1,\beta_2,\cdots,\beta_n) B=(β1,β2,,βn) 为优阵,其中 β \beta β 组为 α \alpha α 组的重排
    3. (封闭性)若 A A A B B B 为同阶优阵,则 A B AB AB 也为优阵

    向量构造优阵

    将向量作为镜面阵的法向量,构造镜面阵(优阵+H阵)

    α = ( a 1 a 2 ⋮ a n ) ∈ C \alpha=\left(

    a1a2an" role="presentation">a1a2an
    \right)\in C α= a1a2an C A = I n − 2 α α H ∣ α ∣ 2 A=I_n-\frac{2\alpha\alpha^H}{\vert \alpha \vert^2} A=Inα22ααH 是一个优阵

    1. A H = A A^H=A AH=A A 2 = I ( A − 1 = A ) A^2=I(A^{-1}=A) A2=I(A1=A)

    2. A A A 为优阵 ( A H A = I ) (A^HA=I) (AHA=I)
      1. A 2 = ( I n − 2 α α H ∣ α ∣ 2 ) ( I n − 2 α α H ∣ α ∣ 2 ) = I n 2 − 4 α α H ∣ α ∣ 2 + 4 ( α α H ) ( α α H ) ∣ α ∣ 4 = I n − 4 α α H ∣ α ∣ 2 + 4 α ( α H α ) α H ∣ α ∣ 4 = I n − 4 α α H ∣ α ∣ 2 + 4 α ( ∣ α ∣ 2 ) α H ∣ α ∣ 4 = I n − 4 α α H ∣ α ∣ 2 + 4 α α H ∣ α ∣ 2 = I n 2. A H A = ( I n − 2 α α H ∣ α ∣ 2 ) H ( I n − 2 α α H ∣ α ∣ 2 ) = ( I n − 2 α α H ∣ α ∣ 2 ) ( I n − 2 α α H ∣ α ∣ 2 ) = A 2 = I ∴ A 是 U 阵

      1.A2=(In2ααH|α|2)(In2ααH|α|2)=In24ααH|α|2+4(ααH)(ααH)|α|4=In4ααH|α|2+4α(αHα)αH|α|4=In4ααH|α|2+4α(|α|2)αH|α|4=In4ααH|α|2+4ααH|α|2=In2.AHA=(In2ααH|α|2)H(In2ααH|α|2)=(In2ααH|α|2)(In2ααH|α|2)=A2=IAU" role="presentation">1.A2=(In2ααH|α|2)(In2ααH|α|2)=In24ααH|α|2+4(ααH)(ααH)|α|4=In4ααH|α|2+4α(αHα)αH|α|4=In4ααH|α|2+4α(|α|2)αH|α|4=In4ααH|α|2+4ααH|α|2=In2.AHA=(In2ααH|α|2)H(In2ααH|α|2)=(In2ααH|α|2)(In2ααH|α|2)=A2=IAU
      1.A22.AHA=(Inα22ααH)(Inα22ααH)=In2α24ααH+α44(ααH)(ααH)=Inα24ααH+α44α(αHα)αH=Inα24ααH+α44α(α2)αH=Inα24ααH+α24ααH=In=(Inα22ααH)H(Inα22ααH)=(Inα22ααH)(Inα22ααH)=A2=IAU

    eg
    α = ( 1 1 1 ) , 其 U 阵为 I 3 − 2 α α H ∣ α ∣ 2 = ( 1 0 0 0 1 0 0 0 1 ) − 2 3 ( 1 1 1 1 1 1 1 1 1 ) = ( 1 3 − 2 3 − 2 3 − 2 3 1 3 − 2 3 − 2 3 − 2 3 1 3 )

    α=(111),UI32ααH|α|2=(100010001)23(111111111)=(132323231323232313)" role="presentation">α=(111),UI32ααH|α|2=(100010001)23(111111111)=(132323231323232313)
    α= 111 ,U阵为I3α22ααH= 100010001 32 111111111 = 313232323132323231

  • 相关阅读:
    使用python连接Linux服务器发送指定命令
    ChatGPT 控制机器人的基本框架
    TensorFlow搭建CNN-LSTM混合模型实现多变量多步长时间序列预测(负荷预测)
    【uniapp/uview1.x】u-upload 在 v-for 中的使用时, before-upload 如何传参
    react经典面试题解析
    ssm基于微信小程序的新生自助报到系统+ssm+uinapp+Mysql+计算机毕业设计
    基于C++的通讯发报应用课程设计
    vscode代码拼写错误检测插件
    MPEG算法及代码例程
    onps栈移植说明(1)——onps栈的配置及裁剪
  • 原文地址:https://blog.csdn.net/qq_40479037/article/details/126942347