SVD详解

 矩阵的奇异值(singular value decompostion, SVD)在最优化问题、特征值问题、最小二乘方问题、广义逆矩阵问题及统计学方面都有重要的应用。本人在研究生物基因组规模的代谢网络和图像识别的相关算法时,都大量的涉及了这一方面的内容,因此开始对此进行初步的理解和研究。我们从基础的概念开始,逐渐深入。

1 奇异矩阵

对于一个方阵A,如果行列式|A|=0,则A是奇异矩阵;如果|A| !=0, 则是非奇异矩阵。

已知 |A| !=0  -> A可逆, 所以可逆矩阵是非奇异矩阵;

如果A是奇异矩阵,则 AX=0 有无穷多个解,AX=b有无穷多个解或者无解;(why?)

如果A是非奇异矩阵,则AX=0有零解,AX=b有唯一解。

A(n×n) is a singular matrix <=> Rank(A)<n; 

A(n×n) is a non-singular matrix <=> Rank(A)=n.

2 正交矩阵:如果实方阵Q满足 QTQ = I,或Q-1 = QT则称Q为正交矩阵。

 容易证明,正交矩阵的列向量是两两正交的单位向量。

  推论:1 正交向量是非奇异的。

              证明:det(QTQ)= detQT×detQ = (detQ)2 =  1 ≠ 0

                       所以|Q| ≠0

           2 正交矩阵的逆矩阵仍然是正交矩阵。

              证明

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在Python中,可以使用numpy包的linalg.svd()函数来计算奇异值分解(SVD)。下面是一个使用numpy的linalg.svd()函数求解SVD的示例代码: ```python import numpy as np A = np.array([[2, 4], [1, 3], [0, 0], [0, 0]]) u, s, v = np.linalg.svd(A) print(u, s, v) ``` 这段代码将矩阵A进行了奇异值分解,分别得到了矩阵U、奇异值数组S和矩阵V的转置。其中,U是一个形状为(m, m)的矩阵,S是一个形状为(k)的一维数组,V是一个形状为(n, n)的矩阵。这里的m、n和k分别表示矩阵A的行数、列数和奇异值的个数。 另外,如果你想自定义一个svd()函数来进行奇异值分解,可以使用numpy自带的svd函数。以下是一个自定义的svd()函数的示例代码: ```python import numpy as np def svd(M): u, s, v = np.linalg.svd(M) return u, s, v # 调用自定义的svd()函数 A = np.array([[2, 4], [1, 3], [0, 0], [0, 0]]) u, s, v = svd(A) print(u, s, v) ``` 这段代码定义了一个名为svd()的函数,接受一个形状为(m, n)的矩阵作为输入,并返回奇异值分解得到的矩阵U、奇异值数组S和矩阵V的转置。使用时,你可以将矩阵A作为参数传入svd()函数,然后得到相应的结果。 希望这个回答对你有所帮助!<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [奇异值分解 SVD 的数学解释](https://blog.csdn.net/fan_fan_feng/article/details/79668052)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT0_1"}}] [.reference_item style="max-width: 50%"] - *2* [【SVD(奇异值分解)】详解及python-Numpy实现](https://blog.csdn.net/weixin_43821215/article/details/126003959)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT0_1"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值