matlab 求解非线性规划（fmincon函数）

一、化标准型（注意min，还有<=）

其中X为n维变元向量，G(X)与Ceq(X)均为非线性函数组成的向量，其它变量的含义与线性规划、二次规划中相同.

二、matlab求解非线性规划的求解步骤

1）首先建立M文件fun.m,定义目标函数F(X):

function f = fun(X);
f = F(x)

2）若约束条件中有非线性约束：G(X)  0 或 Ceq(X) = 0，则建立M文件nonlcon.m定义函数G(X)与Ceq(X)：

function [G, Ceq] = nonlcon(X)
G = ...
Ceq = ...

3）建立主程序.非线性规划求解的函数是fmincon

 代码：

 x=fmincon(‘fun’,X0,A,b,Aeq,beq,VLB,VUB,’nonlcon’,options) 

三、例子

 1）先建立M文件 fun.m,定义目标函数

function f=fun(x);   
f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1);

2）再建立M文件mycon.m定义非线性约束

 function [g,ceq]=mycon(x)
 % 非线性不等式约束(记得标准化为<=,并且移项)
 g=[1.5+x(1)*x(2)-x(1)-x(2); % 一定要注意写法的规范，再次强调这里的x是一个向量！不能把x(1)写成x1
    -x(1)*x(2)-10];
 % 非线性等式约束
 ceq=[];

3）主程序main.m

x0=[-1;1]; % 根据约束条件随便取一个
A=[];b=[];
Aeq=[1 1];beq=[0];
vlb=[];vub=[];
[x,fval]=fmincon('fun',x0,A,b,Aeq,beq,vlb,vub,'mycon')