Matlab求解非线性规划问题（用fmincon函数求解NLP）

对于实际问题，把它公式化为非线性规划问题后，可以利用matlab中的fmincon函数求解。首先需要将非线性规划问题的数学模型写成标准形式：
min f(x)
s.t.
Ax<=B;%线性不等式约束
Aeq*x==Beq；%线性等式约束
LB<=x<=UB；%变量上下限约束
con(x)<=0；%非线性不等式约束
coneq(x)==0；%非线性等式约束

matlab中求解非线性规划的函数：
[X,Y] = FMINCON(FUN,X0,A,B,Aeq,Beq,LB,UB,NONLCON,OPTIONS)
其中，
返回值是最优解X和对应的目标值Y；
FUN是目标函数，用f(x)的.m文件名表示；
X0是给定的变量初始值；
NONLCON是非线性约束，要用[con(x),coneq(x)]的.m文件名表示；
某一约束若没有，则对应的位置赋值为[]，OPTIONS设定求解方法，可以使用默认值.
示例如下：
min f(x)=2x(1)²-4x(1)x(2)+4x(2)²-6x(1)-3x(2)
s.t.
x(1)+x(2)<=3;
4*x(1)+x(2)<=9;
x(1),x(2)>=0;
-x(1)²+x(2)<=0;
-x(1)-x(2)²+2==0;

编写目标函数的.m文件：

function objvalue=eg4_f1(x)
objvalue=2*x(1)^2-4*x(1)*x(2)+4*x(2)^2-6*x(1)-3*x(2);
end

编写非线性约束的.m文件：

function [cong,conh]=eg4_f2(x)
cong=-x(1)^2+x(2);
conh=-x(1)-x(2)^2+2;
end

命令行输入：
[x,y]=fmincon(‘eg4_f1’,rand(2,1),[],[],[],[],zeros(2,1),[],‘eg4_f2’)
求解结果如下：

[x,y]=fmincon('eg4_f1',rand(2,1),[],[],[],[],zeros(2,1),[],'eg4_f2')

Local minimum found that satisfies the constraints.

Optimization completed because the objective function is non-decreasing in 
feasible directions, to within the default value of the optimality tolerance,
and constraints are satisfied to within the default value of the constraint tolerance.

<stopping criteria details>


x =

    1.5283
    0.6868


y =

   -8.8706