NLopt中的无导数算法


(本文所涉及各量皆为C风格,即包含的头文件是 nlopt.h而非 nlopt.hpp
首先说一下NLopt中算法的命名规则,C风格的命名规则为 NLOPT_{G,L}{N,D}_xxxx,其中 G/L代表全局最优与局部最优, N/D代表无导数与有导数。例如 NLOPT_LN_COBYLA则表示无导数局部最优的COBYLA算法。
这里的导数是我们算完之后填入

if (grad)
{
   
	/* 各变量的导数 */
}

中的,但是在很多情况下,当问题很复杂时,很难求解出其导数,这个时候就需要使用无导数的算法。

无导数全局最优算法

这种算法的命名为NLOPT_GN_xxxx。最简单的方法就是在编译器中输入NLOPT_GN_,然后它会自己弹出来有哪些算法,一个一个试,看哪个满足自己的要求即可:)。

DIRECT和DIRECT-L

DIRECT指的是Dividing Rectangles算法。
首先是基本的NLOPT_GN_DIRECTNLOPT_GN_DIRECT_L。其中NLOPT_GN_DIRECT将更多的精力放在了全局,而NLOPT_GN_DIRECT_L算法在局部的表现会更好,其适用于没有过多局部最小值的函数。
其次是NLOPT_GN_DIRECT_L_RAND,该变体加入了一点随机性以帮助确定将哪个维度二等分。
由于DIRECT算法将各约束缩放为超立方体(与正方体类似,只是正方体为三维,而超立方体为指定维度),因此所有维度的搜索权重是相同的。如果你要搜索的空间与超立方体相差甚远,则更适合用DIRECT的变体NLOPT_GNL_DIRECT_NOSCAL, NLOPT_GN_DIRECT_L_NOSCALNLOPT_GN_DIRECT_L_RAND_NOSCAL
同时,还可以使用基于原始的Fortran语言写的NLopt算法NLOPT_GN_ORIG_DIRECT, NLOPT_GN_ORIG_DIRECT_L,在不同的情况下,这两个的表现可能会更好,不过也可能会变差,需要实际尝试。

  • 4
    点赞
  • 26
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
NLOPT是一个用于非线性优化的开源库。它支持多种优化算法和约束条件,并可以在C语言进行使用。下面是一个简单的示例代码,演示了如何在C使用NLOPT进行优化: ```c #include <stdio.h> #include <math.h> #include <nlopt.h> // 目标函数 double myfunc(unsigned n, const double *x, double *grad, void *my_func_data) { if (grad) { // 计算梯度,如果需要的话 grad[0] = 0.0; grad[1] = 0.5 / sqrt(x[1]); } return sqrt(x[1]); } int main() { double lb[2] = {0.0, 0.01}; // 变量下界 double ub[2] = {1.0, 1.0}; // 变量上界 double x[2] = {0.5, 0.5}; // 初始变量值 nlopt_opt opt; opt = nlopt_create(NLOPT_LD_MMA, 2); // 选择优化算法和变量数量 nlopt_set_lower_bounds(opt, lb); // 设置变量下界 nlopt_set_upper_bounds(opt, ub); // 设置变量上界 nlopt_set_min_objective(opt, myfunc, NULL); // 设置目标函数 double minf; // 最小目标函数值 if (nlopt_optimize(opt, x, &minf) < 0) { printf("Optimization failed!\n"); } else { printf("Found minimum at f(%g,%g) = %0.10g\n", x[0], x[1], minf); } nlopt_destroy(opt); // 销毁优化对象 return 0; } ``` 在上面的示例,使用了NLOPT的LD_MMA算法进行优化。你可以根据你的具体需求选择其他的优化算法。首先,需要定义一个目标函数`myfunc`,然后创建一个优化对象`opt`,并设置变量的上下界、目标函数等。最后,通过调用`nlopt_optimize`函数进行优化,得到最小目标函数值。 请注意,这只是NLOPT的简单示例,你可以根据你的具体问题进行相应的修改和扩展。你可以在NLOPT的官方文档找到更多关于其用法和参数设置的信息。希望这对你有所帮助!如果你有任何其他问题,请随时提问。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值