低秩矩阵近似-奇异值分解
奇异值分解有广泛的用途,例如对于低秩矩阵近似(low-rank matrix approximation)问题,给定一个秩为。式(A.33)中的分解称为奇异值分解(Singular Value Decomposition, 简称 SVD),其中。的右奇异向量(right-singular vector),的左奇异向量(left-singular vector),的秩(rank)就等于非零奇异值的个
A.3 奇异值分解
任意实矩阵 A∈Rm×n\mathbf{A} \in \mathbb{R}^{m \times n}A∈Rm×n 都可分解为
A=UΣVT,(A.33) \mathbf{A} = \mathbf{U} \mathbf{\Sigma} \mathbf{V}^T, \tag{A.33} A=UΣVT,(A.33)
其中, U∈Rm×m\mathbf{U} \in \mathbb{R}^{m \times m}U∈Rm×m 是满足 UTU=I\mathbf{U}^T \mathbf{U} = \mathbf{I}UTU=I 的 mmm 阶酉矩阵(unitary matrix); V∈Rn×n\mathbf{V} \in \mathbb{R}^{n \times n}V∈Rn×n 是满足 VTV=I\mathbf{V}^T \mathbf{V} = \mathbf{I}VTV=I 的 nnn 阶酉矩阵; Σ∈Rm×n\mathbf{\Sigma} \in \mathbb{R}^{m \times n}Σ∈Rm×n 是 m×nm \times nm×n 的矩阵,其中 (Σ)ii=σi(\mathbf{\Sigma})_{ii} = \sigma_i(Σ)ii=σi 且其他位置的元素均为 0, σi\sigma_iσi 为非负实数且满足 σ1≥σ2≥…≥0\sigma_1 \geq \sigma_2 \geq \ldots \geq 0σ1≥σ2≥…≥0 。
式(A.33)中的分解称为奇异值分解(Singular Value Decomposition, 简称 SVD),其中 U\mathbf{U}U 的列向量 ui∈Rmu_i \in \mathbb{R}^mui∈Rm 称为 A\mathbf{A}A 的左奇异向量(left-singular vector), V\mathbf{V}V 的列向量 vi∈Rnv_i \in \mathbb{R}^nvi∈Rn 称为 A\mathbf{A}A 的右奇异向量(right-singular vector), σi\sigma_iσi 称为奇异值(singular value)。矩阵 A\mathbf{A}A 的秩(rank)就等于非零奇异值的个数。
奇异值分解有广泛的用途,例如对于低秩矩阵近似(low-rank matrix approximation)问题,给定一个秩为 rrr 的矩阵 A\mathbf{A}A ,欲求其最优 kkk 秩近似矩阵 A~\tilde{\mathbf{A}}A~ , k<rk < rk<r ,该问题可形式化为
minA~∈Rm×n∥A−A~∥Fs.t.rank(A~)=k.(A.34) \min_{\tilde{\mathbf{A}} \in \mathbb{R}^{m \times n}} \|\mathbf{A} - \tilde{\mathbf{A}}\|_F\\ s.t. \text{rank}(\tilde{\mathbf{A}}) = k . \tag{A.34} A~∈Rm×nmin∥A−A~∥Fs.t.rank(A~)=k.(A.34)
奇异值分解提供了上述问题的解析解:对矩阵 A\mathbf{A}A 进行奇异值分解后,将矩阵 Σ\mathbf{\Sigma}Σ 中的 r−kr-kr−k 个最小的奇异值置零获得矩阵 Σk\mathbf{\Sigma}_kΣk ,即仅保留最大的 kkk 个奇异值,则
Ak=UkΣkVkT(A.35) \mathbf{A}_k = \mathbf{U}_k \mathbf{\Sigma}_k \mathbf{V}_k^T \tag{A.35} Ak=UkΣkVkT(A.35)
就是式(A.34)的最优解,其中 Uk\mathbf{U}_kUk 和 Vk\mathbf{V}_kVk 分别是式(A.33)中的前 kkk 列组成的矩阵。这个结果称为 Eckart-Young-Mirsky 定理。
补充酉矩阵性质:
酉矩阵的性质
酉矩阵(Unitary Matrix)是复数域上满足以下条件的n×nn \times nn×n矩阵UUU:
U∗U=IU^* U = IU∗U=I
其中U∗U^*U∗表示UUU的共轭转置,III为单位矩阵。以下是其主要性质:
-
逆矩阵与共轭转置的关系
酉矩阵的逆矩阵等于其共轭转置:
U−1=U∗U^{-1} = U^*U−1=U∗ -
行列式的模长
酉矩阵的行列式绝对值为1:
∣det(U)∣=1|\det(U)| = 1∣det(U)∣=1 -
保持内积与向量长度
对任意向量x,y∈Cnx, y \in \mathbb{C}^nx,y∈Cn,酉矩阵保持内积和向量长度不变:
⟨Ux,Uy⟩=⟨x,y⟩,∥Ux∥=∥x∥\langle Ux, Uy \rangle = \langle x, y \rangle, \quad \|Ux\| = \|x\|⟨Ux,Uy⟩=⟨x,y⟩,∥Ux∥=∥x∥ -
特征值的性质
酉矩阵的特征值λ\lambdaλ满足∣λ∣=1|\lambda| = 1∣λ∣=1,即分布在复平面的单位圆上。 -
矩阵乘积的酉性
若UUU和VVV均为酉矩阵,则它们的乘积UVUVUV也是酉矩阵。 -
与正交矩阵的关系
在实数域中,酉矩阵退化为正交矩阵,满足Q⊤Q=IQ^\top Q = IQ⊤Q=I。 -
奇异值均为1
酉矩阵的奇异值分解(SVD)中所有奇异值为1,即:
U=UΣV∗ ⟹ Σ=IU = U \Sigma V^* \implies \Sigma = IU=UΣV∗⟹Σ=I
这一性质使其在信号处理和量子力学中尤为重要。
更多推荐

所有评论(0)