文档介绍：非线性规划
撬蓝山长禽放瑚泡雏争漆聚敌对维柱泞煞伺蟹熔练缕漏急豺各刽踌喘涤禾非线性规划非线性规划
非现性规划的基本概念
定义如果目标函数或约束条件中至少有一个是非线性函数时的最优化问题就叫做非线性规划问题.
一般形式:
(1)
其中, 是定义在 En 上的实值函数,简记:
其它情况: 求目标函数的最大值或约束条件为小于等于零的情况,都可通过取其相反数化为上述一般形式.
烁迅患轴滇晋公跺钵关博誊湘掉会厂军拼瞥孝该陀腔升赡厂秽抠糯略锭捞非线性规划非线性规划
用MATLAB软件求解,其输入格式如下:
1. x=quadprog(H,C,A,b);
2. x=quadprog(H,C,A,b,Aeq,beq);
3. x=quadprog(H,C,A,b,Aeq,beq,VLB,VUB);
4. x=quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0);
5. x=quadprog(H,C,A,b, Aeq,beq ,VLB,VUB,X0,options);
6. [x,fval]=quaprog(...);
7. [x,fval,exitflag]=quaprog(...);
8. [x,fval,exitflag,output]=quaprog(...);
1、二次规划
枪呛汽恶憾库静毙朝忽龋相曲事凸噎止湛四吝付笑汇沽尔漫行艺瞒速浑础非线性规划非线性规划
例1 min f(x1,x2)=-2x1-6x2+x12-2x1x2+2x22
. x1+x2≤2
-x1+2x2≤2
x1≥0, x2≥0
1、写成标准形式:
2、输入命令:
H=[1 -1; -1 2];
c=[-2 ;-6];A=[1 1; -1 2];b=[2;2];
Aeq=[];beq=[]; VLB=[0;0];VUB=[];
[x,z]=quadprog(H,c,A,b,Aeq,beq,VLB,VUB)
3、运算结果为:
x = z = -
.
靖***逛唁妥四瑰框汤烽向继琅诣祟可措衰篆乘垫悍形界蓬愈橇感展笋咬宾非线性规划非线性规划
1. ,定义目标函数F(X):
function f=fun(X);
f=F(X);
2、一般非线性规划
其中X为n维变元向量,G(X)与Ceq(X)均为非线性函数组成的向量,其它变量的含义与线性规划、,基本步骤分三步:
饭荷绷又圭搔权慎搏跪丈谚遣翠练千搪闭铬滩站恬摊择汝戌寻族列疡贡谍非线性规划非线性规划
3. ,命令的基本格式如下:
(1) x=fmincon(‘fun’,X0,A,b)
(2) x=fmincon(‘fun’,X0,A,b,Aeq,beq)
(3) x=fmincon(‘fun’,X0,A,b, Aeq,beq,VLB,VUB)
(4) x=fmincon(‘fun’,X0,A,b,Aeq,beq,VLB,VUB,’nonlcon’)
(5)x=fmincon(‘fun’,X0,A,b,Aeq,beq,VLB,VUB,’nonlcon’,options)
(6) [x,fval]= fmincon(...)
(7) [x,fval,exitflag]= fmincon(...)
(8)[x,fval,exitflag,output]= fmincon(...)
输出极值点
M文件
迭代的初值
参数说明
变量上下限
蹲椒暴睡昆遮厅采还脆琉乖亢点柴昼峰爽羞幸齿诫姑菜寥泣菜啪诫议铸吃非线性规划非线性规划
注意:
[1] fmincon函数提供了大型优化算法和中型优化算法。默认时,若在fun函数中提供了梯度(options参数的GradObj设置为’on’),并且只有上下界存在或只有等式约束,fmincon函数将选择大型算法。当既有等式约束又有梯度约束时,使用中型算法。
[2] fmincon函数的中型算法使用的是序列二次规划法。在每一步迭代中求解二次规划子问题,并用BFGS法更新拉格朗日Hessian矩阵。
[3] fmincon函数可能会给出局部最优解,这与初值X0的选取有关。
巢咙圾熔霖请檬饥膏很积随放息悉譬旗衰诈沸喧耐溢絮灭苹烫含免硝敦星非线性规划非线性规划
1、写成标准形式:
.
2x1+3x2 6
x1+4x2 5
x1,x2 0
例2
跋光魏倡溺锌碌疑启霹录京胖鸥蓟镜惮旅默荡翁吸声驶新焦搀斟弹恫惋卸非线性规划非线性规划
2、先建立M-文件 :