最优化拟牛顿法问题！！请帮忙！

apolloxht

目标函数是光滑函数还是多峰的，约束条件为x>=0，采用罚函数法(外点法)，用拟牛顿法作为优化方法；罚因子M=5好像在非可行域是截至的，就是当优化到边界时即x=0时它的方向就没什么变化；当M取很小时，如0.0001到是能够优化几个循环，但结果偏差很大。
主要原因是：到边界时不能找到更好的方向；那位能告诉我如何解决这个问题??
感激不尽！