16广义逆的计算及应用

广义逆的计算及应用

  • 1.由 h e r m i t e hermite hermite标准形求 { 1 } \{1\} {1}
  • 2.由满秩分解求广义逆
  • 3.矩阵方程 A X B = D AXB=D AXB=D的相容性条件及通解
    • 3.1 定理1
    • 3.2 推论
  • 极小范数解

1.由 h e r m i t e hermite hermite标准形求 { 1 } \{1\} {1}

{1}逆即 A X A = A , AXA=A, AXA=A,先想办法把A矩阵表示出来
对任意矩阵,我们都有 E A = B , B EA=B,B EA=B,B H e r m i t e Hermite Hermite矩阵,现在再在B矩阵后乘上置换矩阵P,使得 E A P = [ I r K 0 0 ] EAP=\left[\begin{matrix} I_r& K\\ 0&0\end{matrix}\right] EAP=[Ir0K0] A = E − 1 [ I r K 0 0 ] P − 1 A=E^{-1}\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]P^{-1} A=E1[Ir0K0]P1试着假设 X X X矩阵是 X = P [ I r 0 0 L ] E , X=P\left[\begin{matrix} I_r& 0\\ 0&L \end{matrix}\right]E, X=P[Ir00L]E,接着来证明一下是否是这样。 A X A = A = E − 1 [ I r K 0 0 ] P − 1 P [ I r 0 0 L ] E E − 1 [ I r K 0 0 ] P − 1 AXA=A=E^{-1}\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]P^{-1}P\left[\begin{matrix} I_r& 0\\ 0&L \end{matrix}\right]EE^{-1}\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]P^{-1} AXA=A=E1[Ir0K0]P1P[Ir00L]EE1[Ir0K0]P1 = E − 1 [ I r K 0 0 ] [ I r 0 0 L ] [ I r K 0 0 ] P − 1 = E − 1 [ I r K L 0 0 ] [ I r K 0 0 ] P − 1 = E − 1 [ I r K 0 0 ] P − 1 = A =E^{-1}\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]\left[\begin{matrix} I_r& 0\\ 0&L \end{matrix}\right]\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]P^{-1}=E^{-1}\left[\begin{matrix} I_r& KL\\ 0&0 \end{matrix}\right]\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]P^{-1}=E^{-1}\left[\begin{matrix} I_r& K\\ 0&0 \end{matrix}\right]P^{-1}=A =E1[Ir0K0][Ir00L][Ir0K0]P1=E1[Ir0KL0][Ir0K0]P1=E1[Ir0K0]P1=A故,可知 A ( 1 ) = P [ I r 0 0 L ] E A^{(1)}=P\left[\begin{matrix} I_r& 0\\ 0&L \end{matrix}\right]E A(1)=P[Ir00L]E
16广义逆的计算及应用_第1张图片
简而言之,划重点啦!!!!
{ 1 } 逆 = P [ I r 0 0 L ] E \{1\}逆=P\left[\begin{matrix} I_r& 0\\ 0&L \end{matrix}\right]E {1}=P[Ir00L]E { 2 } 逆 = P [ I r 0 0 0 ] E \{2\}逆=P\left[\begin{matrix} I_r& 0\\ 0&0 \end{matrix}\right]E {2}=P[Ir000]E

2.由满秩分解求广义逆

像这种单纯的及复杂的概念的我觉得不会考(实在是记不住哇)
A A A进行满秩分解; A = F G , A ∈ C r m × n , F ∈ C r m × r , G ∈ C r r × n A=FG,A\in C^{m\times n}_{r},F\in C^{m\times r}_r,G\in C^{r\times n}_r A=FG,ACrm×n,FCrm×r,GCrr×n

  1. G ( i ) F ( 1 ) ∈ A { i }    i = 1 , 2 , 4 G^{(i)}F^{(1)}\in A\{i\}\ \ i=1,2,4 G(i)F(1)A{i}  i=1,2,4
    证明: 当 i = 1 时 : F G G ( 1 ) F ( 1 ) F G = F G ( ∵ G G ( 1 ) = I r , F ( 1 ) F = I r ) 当i=1时:FGG^{(1)}F^{(1)}FG=FG(\because GG^{(1)}=I_r,F^{(1)}F=I_r) i=1FGG(1)F(1)FG=FG(GG(1)=Ir,F(1)F=Ir) i = 1 i=1 i=1时, G ( 1 ) F ( 1 ) ∈ A { 1 } G^{(1)}F^{(1)}\in A\{1\} G(1)F(1)A{1} 得证。 i = 2 时 : G ( 2 ) F ( 1 ) F G G ( 2 ) F ( 1 ) = G ( 2 ) G G ( 2 ) F ( 1 ) = G ( 2 ) F ( 1 ) i=2时:G^{(2)}F^{(1)}FGG^{(2)}F^{(1)}=G^{(2)}GG^{(2)}F^{(1)}=G^{(2)}F^{(1)} i=2G(2)F(1)FGG(2)F(1)=G(2)GG(2)F(1)=G(2)F(1) i = 2 i=2 i=2时, G ( 2 ) F ( 1 ) ∈ A { 2 } G^{(2)}F^{(1)}\in A\{2\} G(2)F(1)A{2} 得证。 i = 4 时 : ( G ( 4 ) F ( 1 ) F G ) H = ( G ( 4 ) G ) H = G ( 4 ) G ( 对 G 运 用 { 4 } 逆 的 定 义 ) = G ( 4 ) I n G = G ( 4 ) F ( 1 ) F G i=4时:(G^{(4)}F^{(1)}FG)^{H}=(G^{(4)}G)^{H}=G^{(4)}G(对G运用\{4\}逆的定义)=G^{(4)}I_nG=G^{(4)}F^{(1)}FG i=4(G(4)F(1)FG)H=(G(4)G)H=G(4)G(G{4})=G(4)InG=G(4)F(1)FG i = 4 i=4 i=4时, G ( 4 ) F ( 1 ) ∈ A { 4 } G^{(4)}F^{(1)}\in A\{4\} G(4)F(1)A{4} 得证。

  2. G ( 1 ) F ( i ) ∈ A { i }    i = 1 , 2 , 3 G^{(1)}F^{(i)}\in A\{i\}\ \ i=1,2,3 G(1)F(i)A{i}  i=1,2,3
    证明: 当 i = 1 时 : F G G ( 1 ) F ( 1 ) F G = F G 当i=1时:FGG^{(1)}F^{(1)}FG=FG i=1FGG(1)F(1)FG=FG i = 1 i=1 i=1时, G ( 1 ) F ( 1 ) ∈ A { 1 } G^{(1)}F^{(1)}\in A\{1\} G(1)F(1)A{1} 得证。 i = 2 时 : G ( 1 ) F ( 2 ) F G G ( 1 ) F ( 2 ) = G ( 1 ) F ( 2 ) F F ( 2 ) = G ( 1 ) F ( 2 ) i=2时:G^{(1)}F^{(2)}FGG^{(1)}F^{(2)}=G^{(1)}F^{(2)}FF^{(2)}=G^{(1)}F^{(2)} i=2G(1)F(2)FGG(1)F(2)=G(1)F(2)FF(2)=G(1)F(2) i = 2 i=2 i=2时, G ( 1 ) F ( 2 ) ∈ A { 2 } G^{(1)}F^{(2)}\in A\{2\} G(1)F(2)A{2} 得证。 i = 3 时 : ( F G G ( 1 ) F ( 3 ) ) H = ( F F ( 3 ) ) H = F F ( 3 ) ( 对 F 运 用 { 3 } 逆 的 定 义 ) = F I n F ( 3 ) = F G G ( 1 ) F ( 3 ) i=3时:(FGG^{(1)}F^{(3)})^{H}=(FF^{(3)})^{H}=FF^{(3)}(对F运用\{3\}逆的定义)=FI_nF^{(3)}=FGG^{(1)}F^{(3)} i=3(FGG(1)F(3))H=(FF(3))H=FF(3)(F{3})=FInF(3)=FGG(1)F(3) i = 3 i=3 i=3时, G ( 1 ) F ( 3 ) ∈ A { 3 } G^{(1)}F^{(3)}\in A\{3\} G(1)F(3)A{3} 得证。

  3. G ( 1 ) F + ∈ A { 1 , 2 , 3 } , G + F ( 1 ) ∈ A { 1 , 2 , 4 } G^{(1)}F^+\in A\{1,2,3\},G^+F^{(1)}\in A\{1,2,4\} G(1)F+A{1,2,3},G+F(1)A{1,2,4}
    由以上1,2的证明可以很自然地得到这个结论。

  4. A + = G + F ( 1 , 3 ) = G ( 1 , 4 ) F + A^+=G^+F^{(1,3)}=G^{(1,4)}F^+ A+=G+F(1,3)=G(1,4)F+
    证明: G + F ( 1 ) ∈ A { 1 , 2 , 4 } G^+F^{(1)}\in A\{1,2,4\} G+F(1)A{1,2,4} 由2可知 G + F ( 3 ) ∈ A { 3 } G^+F^{(3)}\in A\{3\} G+F(3)A{3}综上可以得到: A + = G + F ( 1 , 3 ) A^+=G^+F^{(1,3)} A+=G+F(1,3)同理可得第二个结论

  5. A + = G + F + = G H ( G G H ) − 1 ( F H F ) − 1 F H = G H ( F H A G H ) − 1 F H A^+=G^+F^+=G^H(GG^H)^{-1}(F^HF)^{-1}F^H=G^H(F^HAG^H)^{-1}F^H A+=G+F+=GH(GGH)1(FHF)1FH=GH(FHAGH)1FH
    证明: ∵ G H ( G G H ) − 1 ∈ G { 1 , 2 , 4 } , ( F H F ) − 1 F H ∈ F { 1 , 2 , 3 } \because G^H(GG^H)^{-1}\in G\{1,2,4\},(F^HF)^{-1}F^H\in F\{1,2,3\} GH(GGH)1G{1,2,4},(FHF)1FHF{1,2,3} ∴ 由 前 面 结 论 可 得 A + = G + F + = G H ( G G H ) − 1 ( F H F ) − 1 F H \therefore 由前面结论可得A^+=G^+F^+=G^H(GG^H)^{-1}(F^HF)^{-1}F^H A+=G+F+=GH(GGH)1(FHF)1FH 又 G H ( G G H ) − 1 ( F H F ) − 1 F H = G H ( F H F G G H ) − 1 F H = G H ( F H A G H ) − 1 F H 又 G^H(GG^H)^{-1}(F^HF)^{-1}F^H=G^H(F^HFGG^H)^{-1}F^H=G^H(F^HAG^H)^{-1}F^H GH(GGH)1(FHF)1FH=GH(FHFGGH)1FH=GH(FHAGH)1FH
    16广义逆的计算及应用_第2张图片

3.矩阵方程 A X B = D AXB=D AXB=D的相容性条件及通解

3.1 定理1

矩阵方程 A X B = D AXB=D AXB=D相容的充要条件: A A ( 1 ) D B ( 1 ) B = D , AA^{(1)}DB^{(1)}B=D, AA(1)DB(1)B=D,在相容条件下的矩阵方程通解为 { A ( 1 ) D B ( 1 ) + Y − A ( 1 ) A Y B B ( 1 ) ∣ Y 为 阶 数 合 适 的 任 意 矩 阵 } \{A^{(1)}DB^{(1)}+Y-A^{(1)}AYBB^{(1)}|Y为阶数合适的任意矩阵\} {A(1)DB(1)+YA(1)AYBB(1)Y}
证明:
充分性: 令 X = A ( 1 ) D B ( 1 ) , 显 然 A A ( 1 ) D B ( 1 ) B = A X B = D 令X=A^{(1)}DB^{(1)},显然AA^{(1)}DB^{(1)}B=AXB=D X=A(1)DB(1),AA(1)DB(1)B=AXB=D充分性得证
必要性:
已知 A X B = D AXB=D AXB=D有解,那么 D = A X B = A A ( 1 ) A X B B ( 1 ) B = A A ( 1 ) D B ( 1 ) B D=AXB=AA^{(1)}AXBB^{(1)}B=AA^{(1)}DB^{(1)}B D=AXB=AA(1)AXBB(1)B=AA(1)DB(1)B接下来证明通解。

通解有两个含义:1,解集中的任何元素为方程的解。2,方程的任意解均可以由集合中的元素表示出来
X = A ( 1 ) D B ( 1 ) + Y − A ( 1 ) A Y B B ( 1 ) X=A^{(1)}DB^{(1)}+Y-A^{(1)}AYBB^{(1)} X=A(1)DB(1)+YA(1)AYBB(1)带入可得 A X B = A ( A ( 1 ) D B ( 1 ) + Y − A ( 1 ) A Y B B ( 1 ) ) B = A A ( 1 ) D B ( 1 ) B + A Y B − A Y B = D AXB=A(A^{(1)}DB^{(1)}+Y-A^{(1)}AYBB^{(1)})B=AA^{(1)}DB^{(1)}B+AYB-AYB=D AXB=A(A(1)DB(1)+YA(1)AYBB(1))B=AA(1)DB(1)B+AYBAYB=D
再证明第二点: A X B = D AXB=D AXB=D X = A ( 1 ) D B ( 1 ) + X − A ( 1 ) D B ( 1 ) = A ( 1 ) D B ( 1 ) + X − A ( 1 ) A X B B ( 1 ) X=A^{(1)}DB^{(1)}+X-A^{(1)}DB^{(1)}=A^{(1)}DB^{(1)}+X-A^{(1)}AXBB^{(1)} X=A(1)DB(1)+XA(1)DB(1)=A(1)DB(1)+XA(1)AXBB(1)对应于 Y = X Y=X Y=X的情况。

3.2 推论

  • 线性方程组 A X = b AX=b AX=b有解的充要条件为: A A ( 1 ) b = b ( 取 B 为 I n ) , AA^{(1)}b=b(取B为I_n), AA(1)b=b(BIn),通解为 { A ( 1 ) b + y − A ( 1 ) A y ∣ y 为 列 向 量 } \{A^{(1)}b+y-A^{(1)}Ay|y为列向量\} {A(1)b+yA(1)Ayy}
  • A { 1 } ( A X A = A 的 解 ) A\{1\}(AXA=A的解) A{1}(AXA=A)为如下集合 { A ( 1 ) D A ( 1 ) + X − A ( 1 ) A X A A ( 1 ) } , \{A^{(1)}DA^{(1)}+X-A^{(1)}AXAA^{(1)}\}, {A(1)DA(1)+XA(1)AXAA(1)},四个 A { 1 } A\{1\} A{1}可互不相同
    16广义逆的计算及应用_第3张图片

极小范数解

16广义逆的计算及应用_第4张图片
16广义逆的计算及应用_第5张图片
16广义逆的计算及应用_第6张图片
对通解的理解还是有点欠缺
16广义逆的计算及应用_第7张图片

你可能感兴趣的:(矩阵论,通信,广义逆的计算)