对于实际问题,把它公式化为非线性规划问题后,可以利用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)2-4x(1)x(2)+4x(2)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)2+x(2)<=0;
-x(1)-x(2)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