MATLAB中求非方陣的廣義逆矩陣

如果不是方陣的話，會有個廣義逆矩陣（偽逆），matlab裡面用pinv()實現。
基本文法為X=pinv(A),X=pinv(A,tol),其中tol為誤差,pinv為pseudo-inverse的縮寫：max(size(A))*norm(A)*eps。函數返回一個與A的轉置矩陣A' 同型的矩陣X，並且滿足：AXA=A,XAX=X.此時，稱矩陣X為矩陣A的偽逆，也稱為廣義逆矩陣。

1. >> A=rand(2,3)
2. A =
3.     0.9572    0.8003    0.4218
4.     0.4854    0.1419    0.9157
5. >> B=pinv(A)
6. B =
7.     0.5985   -0.0838
8.     0.7633   -0.5430
9.    -0.4355    1.2206
10. >> A*B*A
11. ans =
12.     0.9572    0.8003    0.4218
13.     0.4854    0.1419    0.9157

複製代碼