“minFunc” MATLAB 软件包（由 Mark Schmidt 编写）是最好的软件包之一： http ://www.di.ens.fr/~mschmidt/Software/minFunc.html

这个包相当成熟，已经在 ML 社区中广泛使用。它包含一个稳定的 L-BFGS 实现以及相关的牛顿和准牛顿优化方法。

有限内存 BFGS 的一个很好的 Matlab 实现是 Tim Kelley 的《优化迭代方法》一书（SIAM，1999；可从出版商网站免费下载的 PDF）。你可以在这里找到他的 Matlab 代码。

关于您的第二个问题：我建议您自己实现该算法，因为这样更容易适应您的特定问题的特殊结构（参数调整，利用衍生信息等）。