Question

请教一个关于matlab函数fmincon的问题

在一篇论文中，提出了10个约束条件。其中Rb,e,Rr是设计变量，h，β1, β2是常量。 y, y' 是θ的函数。在C1中，max的意思是θ取[0,β1], C2中，θ取[0,β2]，也就是约束条件含有参数θ。
作者提到，用的是fmincon这个函数，但是这个函数怎么处理含参约束呢？即C1,C2,C4,C5这五个约束怎么处理？

Answer 1

根据Matlab帮助的说明，fmincon的约束里可以给c(x)<=0这种形式的约束。这里的c(x)是任意的函数。你可以写一个function文件，返回值是max|phi|-30，这样就相当于是约束C1。fmincon会在求解过程中告诉你的function具体的Rb,e,Rr值，然后h、beta是你事先知道的，你在function里面只需要改变theta的值，求出|phi|的最大值即可。

追问

你可以写一个function文件，返回值是max|phi|-30。
我的问题是：|phi|这个函数，是包含Rb,e,Rr,θ四个变量。而Rb,e,Rr三个变量是设计变量，而θ不是。如何编写一个函数文件，使返回值是max|phi|-30。

追答

这里要用到的是nonlcon这个参数，Matlab帮助里面的说明是：

举个例子，假设phi=Rb-(theta-pi/2)^2

那么你的function文件内容就是：

function [c,ceq]=mycon(x)
%假设x(1)代表Rb
     beta=pi;    %这里设定参数β1
     dt=0.001;   %这里设定theta的搜索步长
     theta=0:dt:beta;
     c=max(abs(x(1)-(theta-pi/2).^2))-pi/6;
     ceq=[];
end function

Answer 2

这个是个带约束的非线性规划函数，我用过。我可以解决。