作者celestialgod (天)
看板MATLAB
标题[问题] 解线性方程 \ vs pinv vs least-square
时间Wed Oct 12 21:18:39 2016
我遇到一个奇怪的问题，我看了一下 \ 的文件
遇到矩阵的时候，他会是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://webptt.com/cn.aspx?n=bbs/MATLAB/M.1476278323.A.B80.html
1F：推 sin55688: 你的例子是无穷多组解吧? 并不是无解 10/12 23:31
只是很好奇说 为什麽不会用有唯一的Moore–Penrose pseudoinverse做解
※ 编辑: celestialgod (36.235.43.110), 10/12/2016 23:54:00
2F：推 sin55688: 只能保证有解，但不一定有唯一解。 10/13 01:57
3F：→ sin55688: 剩下可能就要看他解least square是用什麽演算法... 10/13 01:58
4F：推 sunev: 反过来说，干嘛一定要用算得比较慢的pseudoinverse 10/13 08:12