Re: [問題] Conjugate Gradient Method in optim
※ 引述《wheado (principal component QQ)》之銘言:
: [程式範例]:
: http://ideone.com/xCSQ8s
: [環境敘述]:
: R version 3.3.2 (2016-10-31)
: [關鍵字]:
: nonlinear conjugate gradient for minimization
: 感謝大家的指導跟建議
: 這真的好困難QQ
我目前測出來的結果是...應該是演算法本身的問題
那篇教學裡有說非線性的conjugate gradient只有在起始點接近local minima時才會收斂
我在[-4,4]^2裡用你的code跑了10,000個起始點,
黑色的點代表以該點為起始點有收斂到正確的極值 f(x)=-1
http://i.imgur.com/74MDxkO.png
看得出來只有在 (pi/2, -pi/2) modulo pi 的附近才會收斂到正確的點
其他地方只能聽天由命
也可能是程式裡有bug 但我看了半天還是看不出來
不是numDeriv的問題, 換成analytical solution也是一樣
P.S. 你忘記update i=i+1
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 68.52.13.65
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1479182086.A.A4D.html
※ 編輯: abc2090614 (68.52.13.65), 11/15/2016 12:08:52
推
11/15 13:08, , 1F
11/15 13:08, 1F
→
11/15 13:08, , 2F
11/15 13:08, 2F
→
11/15 13:09, , 3F
11/15 13:09, 3F
→
11/15 13:09, , 4F
11/15 13:09, 4F
→
11/15 14:03, , 5F
11/15 14:03, 5F
→
11/15 14:03, , 6F
11/15 14:03, 6F
推
11/15 14:08, , 7F
11/15 14:08, 7F
→
11/15 14:08, , 8F
11/15 14:08, 8F
→
11/15 14:08, , 9F
11/15 14:08, 9F
→
11/15 14:08, , 10F
11/15 14:08, 10F
→
11/15 14:11, , 11F
11/15 14:11, 11F
推
11/15 14:19, , 12F
11/15 14:19, 12F
→
11/15 15:18, , 13F
11/15 15:18, 13F
推
11/15 16:13, , 14F
11/15 16:13, 14F
→
11/16 02:22, , 15F
11/16 02:22, 15F
→
11/16 02:22, , 16F
11/16 02:22, 16F
→
11/16 02:22, , 17F
11/16 02:22, 17F
推
11/16 08:04, , 18F
11/16 08:04, 18F
→
11/16 08:04, , 19F
11/16 08:04, 19F
→
11/16 08:07, , 20F
11/16 08:07, 20F
→
11/16 08:07, , 21F
11/16 08:07, 21F
→
11/16 08:07, , 22F
11/16 08:07, 22F
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 2 篇):