SVD奇异值分解

本文详细介绍了矩阵的特征值分解(EVD)和奇异值分解(SVD),包括它们的定义、应用场景以及相互之间的联系。EVD适用于实对称矩阵,SVD则对任意矩阵有效。SVD可以视为EVD的一种推广,当A为方阵时,EVD的特征值即为SVD中奇异值的平方。通过EVD和SVD,可以进行矩阵的谱分析、低秩近似和数据分析等任务。

一、特征值分解EVD

其中:

1、Q为标准正交阵,Σ为对角矩阵,λ称为特征值,q是Q(特征矩阵)的列向量,称为特征向量。

2、A为方阵,是m*m的实对称矩阵

二、奇异值分解

其中:

1、A是m*n的矩阵

2、U是m*m的标准正交阵,V是n*n的标准正交阵。

3、\Sigma是m*n的矩阵

三、如何做SVD,如何联系SVD和EVD

EVD是SVD的特殊形式,EVD要求更严格,需要A为方阵。

其中:

1、AA^{T}是m*m的,A^{T}A是n*n的

2、而对于A,其特征值λ(奇异值)是AA^{T}A^{T}A的特征值开根号得到。

 

 

 四、举例

1、计算AA^{T}A^{T}A

 

 2、计算A^{T}A的特征值、特征向量,从而得到V、\Sigma ^{T}\Sigma

 

 3、 计算AA^{T}的特征值、特征向量U、\Sigma \Sigma ^{T}

4、 计算A的奇异值

\delta = \sqrt{3}  

\delta = 1

5、最终奇异值分解

参考:

1、 SVD(奇异值分解)小结 - EndlessCoding - 博客园 (cnblogs.com)

2、奇异值分解(SVD) - 知乎 (zhihu.com)

SVD奇异值分解是一种常用的矩阵分解方法,可用于信号和图像的降噪。其原理和方法如下: ### 原理 奇异值分解可以将一个矩阵分解为三个矩阵的乘积,基本思想是把矩阵分解为U、Σ和V三个部分,其中U和V是正交矩阵,Σ是对角矩阵。通常噪声对应较小的奇异值,通过对Σ中的奇异值进行适当处理,如丢弃小奇异值并重构矩阵,就可以达到一定程度的降噪效果[^2][^4]。 ### 方法 在数字信号处理中,使用SVD进行降噪的具体步骤如下: 1. **构建矩阵**:将待降噪的信号构成一个矩阵[^2]。 2. **SVD分解**:对该矩阵进行SVD分解,得到U、Σ和V [^2]。 3. **选择奇异值**:根据降噪的需求,选择保留的奇异值个数,将其余的奇异值置零。较小的奇异值通常对应噪声,去除这些值可以减少噪声影响 [^2][^3][^4]。 4. **重构矩阵**:根据保留的奇异值,重新构建矩阵 [^2]。 5. **逆SVD运算**:对重构后的矩阵进行逆SVD运算,得到降噪后的信号 [^2]。 在图像降噪中,通过调整保留的奇异值数量k,可以在噪声去除和细节保留之间取得平衡。较小的k值(如5或10)可以有效去除噪声,但会丢失图像细节;较大的k值(如50或100)能保留更多细节,但降噪效果会降低 [^3]。 ### 代码示例(Python) 以下是一个简单的Python示例,使用NumPy库进行SVD降噪: ```python import numpy as np # 生成含噪声的信号 signal = np.linspace(0, 10, 100) noise = np.random.normal(0, 1, 100) noisy_signal = signal + noise # 构建矩阵 matrix = noisy_signal.reshape(-1, 1) # SVD分解 U, S, V = np.linalg.svd(matrix, full_matrices=False) # 选择保留的奇异值个数 k = 1 S[k:] = 0 # 重构矩阵 reconstructed_matrix = np.dot(U * S, V) # 提取降噪后的信号 denoised_signal = reconstructed_matrix.flatten() print("降噪前信号:", noisy_signal) print("降噪后信号:", denoised_signal) ```
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值