[問題] 解線性方程 \ vs pinv vs least-square

看板MATLAB作者 (天)時間8年前 (2016/10/12 21:18), 8年前編輯推噓3(301)
留言4則, 2人參與, 最新討論串1/1
我遇到一個奇怪的問題,我看了一下 \ 的文件 遇到矩陣的時候,他會是least-square的解 但是我剛好碰到一個例子是下面這樣,解就會不一樣 X = [13.5, 6.75, 6.75;6.75,6.75,0;6.75,0,6.75]; y = [0.62;0.31;0.31]; 1. X\y % [-0.0791, 0.125, 0.125] 2. pinv(X) * y % [0.0306, 0.0153, 0.0153] 3. pinv(X * X') * X' * y % [0.0306, 0.0153, 0.0153] 我的Matlab是R2016a 想問問看有沒有其他人有一樣問題,或是其他版本沒錯的... 還是我對\有所誤會 -- R資料整理套件系列文: magrittr #1LhSWhpH (R_Language) http://tinyurl.com/j3ql84c data.table #1LhW7Tvj (R_Language) http://tinyurl.com/hr77hrn dplyr(上) #1LhpJCfB (R_Language) http://tinyurl.com/jtg4hau dplyr(下) #1Lhw8b-s (R_Language) tidyr #1Liqls1R (R_Language) http://tinyurl.com/jq3o2g3 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.235.43.110 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1476278323.A.B80.html

10/12 23:31, , 1F
你的例子是無窮多組解吧? 並不是無解
10/12 23:31, 1F
只是很好奇說 為什麼不會用有唯一的Moore–Penrose pseudoinverse做解 ※ 編輯: celestialgod (36.235.43.110), 10/12/2016 23:54:00

10/13 01:57, , 2F
只能保證有解,但不一定有唯一解。
10/13 01:57, 2F

10/13 01:58, , 3F
剩下可能就要看他解least square是用什麼演算法...
10/13 01:58, 3F

10/13 08:12, , 4F
反過來說,幹嘛一定要用算得比較慢的pseudoinverse
10/13 08:12, 4F
文章代碼(AID): #1N_ZWpk0 (MATLAB)
文章代碼(AID): #1N_ZWpk0 (MATLAB)