请问哪位大侠知道求伪逆矩阵的方法

来源:百度知道 编辑:UC知道 时间:2024/05/31 10:45:16
不用matlab软件求,小弟想知道计算上有什么的技巧吗?请给个解题例子,感激不尽!

如果A列满秩,那么pinv(A)=(A'*A)^{-1}*A'。
如果A行满秩,那么pinv(A)=pinv(A')'。
如果秩亏损,那么只好先做奇异值分解A=UDV',U,V是正交阵,D是对角阵。
然后取对角阵S,如果D(i,i)=0,那么S(i,i)=0,如果D(i,i)<>0,那么S(i,i)=1/D(i,i)。于是pinv(A)=VSU'。