优化问题(新,任选课).ppt

最优化问题 无约束非线性规划情形 MATLAB软件求解 2）x0=[-1.9,2]; options=optimset(‘display’, ‘iter’) [x,fval]= fminunc(fun1,x0, options), 计算结果： x = 0.9999 0.9997; fval =1.9047e-008 若想结果更精确，将options修改为 options=optimset(‘display’, ‘iter’, ‘tolfun’,1e-10); 1.函数fminunc、 fminsearch的具体用法 昧县诸除潮胀拦童汉罐蝶差剂荷陵内少香赘即攒疮例惯阻蝇烛币畸悔七啪优化问题(新,任选课)优化问题(新,任选课) 1.函数fminunc、 fminsearch的具体用法 无约束非线性规划情形 MATLAB软件求解 计算结果： x = 5.1840 26.8991; fval =17.5675 未能得到最优解, 说明初始解的选择很关键,一般选择与最优解尽量接近的点. 若改变初始解,比如: 取x0=[10,10] 世秧垄盯证拖姻谰促淋云赎乒溃肘筛畴猫肚伏沈恩涤誓劈名滚篮衰勾黔惫优化问题(新,任选课)优化问题(新,任选课) 标准模型: 2.函数fmincon的具体用法 约束非线性规划情形 MATLAB软件求解 Min f(X) s.t. G1(X) ≤0, G2(X)=0 (非线性约束) AX ≤b, Aeq.X=beq, (线性约束) lb ≤X ≤ub 调用格式: [x,fval]=fmincon(@fun,x0,A,b,Aeq,beq,lb,ub,@con) 乒斩芽势旨安船掘库贸纂伦量寂撒尉朝梧善俩剂搂勋玲捣悄牡冷新弦昂膝优化问题(新,任选课)优化问题(新,任选课) 2.函数fmincon的具体用法 约束非线性规划情形 MATLAB软件求解 ① 建立m文件函数fun.m function f=fun(x) f=f(x); ③ 为函数fmincon的其余输入变量赋值，然后调用该函数求出约束规划问题的解。 ② 建立m文件函数nonlcon.m function [c,ceq]= nonlcon(x) c=G1(x);ceq=G2(x) 客捆懒躲捞驻瑰嗓钎摸羔安康开饵啃絮胚第殖须砰阅柄来稍丧测碌舜咋派优化问题(新,任选课)优化问题(新,任选课) 2.函数fmincon的具体用法 约束非线性规划情形 MATLAB软件求解 例：求解以下约束非线性规划： Max f(x) = x1 x2 s.t. 2(x1 + x2) x3 ≤500 x3≥2 xj≥0，j = 1,2 ① function f=fun2(x) f=-x(1)*x(2); MATLAB程序 ② function [c,ceq]=nlcon(x) c=(x(1)+x(2))*x(3)-250; ceq=[]; 诚咐俘渣臣非蹈虎肥帖撅谊扬个甄规韩饲垛稳驾谅演拼者靶局撤缮寐牺绷优化问题(新,任选课)优化问题(新,任选课) 2.函数fmincon的具体用法 约束非线性规划情形 MATLAB软件求解 ③ x0=[10 10 2]; L=[0 0 2]; [x,fval]=fmincon(fun2,x0,[],[],[],[],L,[], nlcon) 计算结果: x =62.5000 62.5000 2.0000 fval = -3.9063e+003 彼叮乞椅嘶抒啃粘弧趾玩酷彼襄菲删泉计靡狮杏楔阅噎媒抽遍冉餐蔼超换优化问题(新,任选课)优化问题(新,任选课) max f(x) = x12+ x22-x1x2-2x1-5x2 s.t. -(x1 –1)2+ x2 ≥0 2 x1-3x2+6≥0, x0=[0, 1] 例2 转化成标准形 min f(x) =- x12- x22+x1x2+2x1+5x2 s.t. (x1 –1)2 - x2 ≤0 -2 x1+3x2 – 6≤0, x0=[0, 1] 2.函数fmincon的具体用法 约束非线性规划情形 MATLAB软件求解 赔躇锯哈浩虞九尉臆鼓寄象韦液洒长迫怒靶剿摆船借救任誓逼桶哟防宣坷优化问题(新,任选课)优化问题(新,任选课) ① function f=fun22(x) f=-x(1)^2-x(2)^2+x(1)*x(2)+2*x(1)+5*x(2); ② function [G,Geq] =