L33-左右逆和伪逆-习题集

linear-algebra

L33-左右逆和伪逆-习题集

参考

问题 32.1

求矩阵 A=[101010]A=\begin{bmatrix} 1 & 0 & 1 \\ 0 & 1 & 0 \end{bmatrix} 的右逆 right inverse

解答

矩阵 AA 行满秩,存在右逆,根据公式可以求出它的右逆 Aright1=AT(AAT)1A_{right}^{-1}=A^{T}(AA^{T})^{-1}

Aright1=AT(AAT)1=[100110]([101010][100110])1=[100110]([2001])1\begin{aligned} A_{right}^{-1}&=A^{T}(AA^{T})^{-1} \\ &= \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 1 & 0 \end{bmatrix} (\begin{bmatrix} 1 & 0 & 1 \\ 0 & 1 & 0 \end{bmatrix} \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 1 & 0 \end{bmatrix})^{-1} \\ &= \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 1 & 0 \end{bmatrix} (\begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix})^{-1} \end{aligned}

可以通过构建增广矩阵的算法来求解矩阵 [2001]\begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix} 的逆矩阵,也可以通过逆矩阵的代数表达式 A1=1detACTA^{-1}=\cfrac{1}{detA}C^{T} 来求解逆矩阵

可得 ([2001])1=[1/2001](\begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix})^{-1}=\begin{bmatrix} 1/2 & 0 \\ 0 & 1 \end{bmatrix}

所以矩阵 AA 的右逆为

Aright1=AT(AAT)1=[100110][1/2001]\begin{aligned} A_{right}^{-1}&=A^{T}(AA^{T})^{-1} \\ &= \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} 1/2 & 0 \\ 0 & 1 \end{bmatrix} \end{aligned}

可以通过 AAright1=IAA_{right}^{-1}=I 验证运算结果

问题 32.2

矩阵 A=[4386]A=\begin{bmatrix} 4 & 3 \\ 8 & 6 \end{bmatrix} 是否存在左逆 left inverse?是否存在右逆 right inverse?是否存在伪逆 pseudoinverse?如果存在请求出具体的矩阵。

解答

由于矩阵 AA 的第一行向量与第二行向量线性相关 2row1=row22row1=row2 所以矩阵 AA 的秩为 r=1r=1,即矩阵 AA 行不满秩,且列不满秩,则矩阵没有右逆和左逆。

提示

也可以通过行列式 detAdetA 来判断矩阵是否为可逆矩阵

detA=4×63×8=2424=0detA=4 \times 6 - 3 \times 8=24-24=0

假如存在右逆 Aright1A_{right}^{-1},则根据行列式的特性 9 detCdetD=det(CD)detC \cdot detD=det(CD) 可得 detAdet(Aright1)=det(AAright1)=detIdetA \cdot det(A_{right}^{-1})=det(AA_{right}^{-1})=detIdetA=0detA=0,则可得一个矛盾的等式 0=detI0=detI,所以矩阵 AA 不存在右逆

证明左逆不存在的步骤类似

矩阵 AA 具有伪逆 A+=VΣ+UTA^{+}=V \Sigma^{+} U^{T}

L30-奇异值分解课堂中对矩阵 AA 进行了 SVD 奇异值分解,所求得的各部分的矩阵如下

A=UΣVT=[1/52/52/51/5][125000][0.80.60.60.8]\begin{aligned} A&=U \Sigma V^{T} \\ &= \begin{bmatrix} 1/\sqrt{5} & 2/\sqrt{5} \\ 2/\sqrt{5} & -1/\sqrt{5} \end{bmatrix} \begin{bmatrix} \sqrt{125} & 0 \\ 0 & 0 \end{bmatrix} \begin{bmatrix} 0.8 & 0.6 \\ 0.6 & -0.8 \end{bmatrix} \end{aligned}

在矩阵 AA 的伪逆表达式中矩阵 Σ+\Sigma^{+} 是对角矩阵 Σ\Sigma 的伪逆,可以根据矩阵 Σ\Sigma 得出

Σ+=[1/125000]\Sigma^{+}= \begin{bmatrix} 1/\sqrt{125} & 0 \\ 0 & 0 \end{bmatrix}

所以矩阵 AA 的伪逆为

A+=VΣ+UT=[0.80.60.60.8][1/125000][1/52/52/51/5]=[0.80.60.60.8][1/55000]1/5[1221]=[0.80.60.60.8][1/25000][1221]=1125[4836]\begin{aligned} A^{+}&=V \Sigma^{+} U^{T} \\ &=\begin{bmatrix} 0.8 & 0.6 \\ 0.6 & -0.8 \end{bmatrix} \begin{bmatrix} 1/\sqrt{125} & 0 \\ 0 & 0 \end{bmatrix} \begin{bmatrix} 1/\sqrt{5} & 2/\sqrt{5} \\ 2/\sqrt{5} & -1/\sqrt{5} \end{bmatrix} \\ &=\begin{bmatrix} 0.8 & 0.6 \\ 0.6 & -0.8 \end{bmatrix} \begin{bmatrix} 1/5\sqrt{5} & 0 \\ 0 & 0 \end{bmatrix} 1/\sqrt{5}\begin{bmatrix} 1 & 2 \\ 2 & -1 \end{bmatrix} \\ &=\begin{bmatrix} 0.8 & 0.6 \\ 0.6 & -0.8 \end{bmatrix} \begin{bmatrix} 1/25 & 0 \\ 0 & 0 \end{bmatrix} \begin{bmatrix} 1 & 2 \\ 2 & -1 \end{bmatrix} \\ &=\cfrac{1}{125} \begin{bmatrix} 4 & 8 \\ 3 & 6 \end{bmatrix} \end{aligned}

Copyright © 2024 Ben

Theme BlogiNote

Icons from Icônes