实践Yalmip使用Knitro的一些总结

Posted 2020-11-20 landiljy

Yalmip使用Knitro的一些总结

1.软件

　　Knitro 11.0.1  

　　Win64（包含安装包和确定机器ID的软件）：链接：https://pan.baidu.com/s/14IfxlAdo3myYiAeCGzfSQQ 密码：2k9q

　　

图1  Knitro 11.0.1 win64 下载地址

　　官网获取地址：https://www.artelys.com/en/optimization-tools/knitro/downloads

　　通过学术邮箱可以申请一个月的完整版使用权。

图2 Knitro 的申请

　　Yalmip的获取教程见https://www.cnblogs.com/landiljy/p/6776434.html

 

2.Yalmip使用Knitro

• Yalmip和Knitro都安装好以后，命令行窗口输入yalmiptest，看下是否找到knitro求解器
• Yalmip中调用knitro。见下面代码

1 options = sdpsettings(‘verbose‘,2,‘solver‘,‘knitro‘,‘savesolveroutput‘,1);%verbose设置显示信息的程度；savesolveroutput默认为0，表示不保存求解结果，设为1表示保存求解结果。
2 sol = optimize(Constraints,F,options);%Constraints为约束；F为目标函数；options为属性。

• 按照上面步骤设置savesolveroutput=1后，报错。
• 然后查看了Yalmip调用Knitro求解器的函数文件，发现了问题。方法如下：
• 按照路径F:MATLAB2014AddToolYALMIP-mastersolverscallknitro.m（红色部分为软件里的路径部分）,打开callknitro.m文件。
• 可以发现图4中保存计算结果的语句中fmin语句有错（没有fmin这个变量）。图3中返回的是fval。
• 其中这里可以参考：http://blog.sina.com.cn/s/blog_50c154510102wusb.html

图3  Yalmip调用Knitro的Matlab版本接口

图4 Yalmip调用Knitro求解器的函数

• 设置Knitro的参数方法如下：
• Knitro的Matlab版本的Nlp例子路径：C:Program FilesArtelysKnitro 11.0.1knitromatlab

[x,fval,exitflag,output,lambda,grad,hess] = knitromatlab(@objfun,x0,A,b,Aeq,beq,lb,ub,@constfun,extendedFeatures,[],‘nlp2options.opt‘);%红色字体是参数设置文档

• 图5中exampleNLP2.m文件是Knitro给的NLP的例子，其中nlp2options.opt是Knitro参数设置文档，具体格式见图6。
• 图6中加了‘#’的语句，则‘#’后面的语句不执行（即蓝色框）。红色框为一个参数的设置语句(参数详见：https://www.gams.com/latest/docs/S_KNITRO.html#KNITROopttol)，后半句为说明这个参数对应Yalmip中的参数名称（即图7中蓝色框部分）.
• 那么我们将nlp2options.opt文件拷贝到我们自己项目的路径下，然后修改callknitro.m文件，具体修改见最后的代码红色字体部分。
• 那么我们以后要修改Knitro的参数，只需要在callknitro.m文件中修改即可。

 

图5  Knitro的Matlab版本的例子

图6 Knitro的参数设置文档

图7 Yalmip中Knitro的参数列表

 

1 solvertime = tic;
2 %自己加的对options的设置---start----2018.06.27
3 model.options.knitro.optionsfile = ‘nlp2options.opt‘;  %参数设置文档
4 %自己加的对options的设置---end----
5 [xout,fval,exitflag,output,lambda] = knitromatlab_mip(funcs.objective,model.x0,model.A,full(model.b),model.Aeq,full(model.beq),model.lb,model.ub,funcs.constraints,model.xType,model.objFnType,model.cineqFnType,model.extendedFeatures,model.options.knitro,model.options.knitro.optionsfile);
6 solvertime = toc(solvertime);