[问题] Conjugate Gradient Method in optim

楼主: wheado (principal component QQ)   2016-11-14 15:23:44
[问题类型]:
程式咨询(我想用R 做某件事情,但是我不知道要怎么用R 写出来)
[软件熟悉度]:
入门(写过其他程式,只是对语法不熟悉)
[问题叙述]:
我想自己练习写一个Nonlinear Conjugate Gradient Method,可是我参考文献
http://www.cs.cmu.edu/%7Equake-papers/painless-conjugate-gradient.pdf
中的 "52页 B4-table" 所写的程式码却跟R内建optim找到的不一样,
很明显不是最小值,试了很久还是不懂哪里出问题,我把我的程式贴在下面。
我想写出来的程式是想跟optim一样只需要
Input 函数 起始点 迭代次数 误差
就可以 output 一个解使得函数是local min
我使用numDeriv套件找到该点的 梯度 跟 hessian
[程式范例]:
http://ideone.com/xCSQ8s
[环境叙述]:
R version 3.3.2 (2016-10-31)
[关键字]:
nonlinear conjugate gradient for minimization
感谢大家的指导跟建议
这真的好困难QQ

Links booklink

Contact Us: admin [ a t ] ucptt.com