※ 引述《Warlock (覡)》之銘言：
: ※ 引述《chanbruce (非線性)》之銘言：
: : 我嘗試用pinv(A)
: : A=[1 2 3
: : 4 5 6
: : 7 8 9]
: : s=pinv(A)
: : s*A
: : ans =
: : 0.8333 0.3333 -0.1667
: : 0.3333 0.3333 0.3333
: : -0.1667 0.3333 0.8333
: : ?什麼答案不是 I
: 那個A是 Singular 的
: 意即根本沒有矩陣可以與A相乘 結果乘出來是 I的
: 看determinant就知道了....
: 若 存在 B 使得 A*B=I
: (A,B 為 方陣)
: 則因為
: det(A*B) = det(A) * det (B)
: 且 A*B=I, det(I)=1, det(A)=0
: det(A*B)=det(I)=1=det(A)*det(B)=0*det(B)
: 故 det(B) * 0 = 1
: 根本沒有B這種東西~~~
: pseudo-inverse 的定義要好好看一看~~~
用pinv(A)來表示A的pseudo-inverse
pinv(A)*A*pinv(A)=pinv(A)
A*pinv(A)*A=A
(pinv(A)*A)^-1=pinv(A)*A
(A*pinv(A))^-1=A*pinv(A)
以上是pseudo-inverse的定義
所以你用pinv(A)*A不會得到I
但是如果在乘上一個pinv(A)會得到pinv(A)
pinv(A)*A得到的只是一個很像I的東西....假的I
因為你用的是假的inverse.....
--
因為我....
交友不慎...行為差池....
所以才會導致這樣的結果....
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.25.179