数学实验第次讲稿.pptVIP

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] = cont2(x) G=(x(1)-1)^2 - x(2); Geq=[]; ③ x0=[0 1]; A=[-2,3];b=6; Aeq=[];beq=[]; lb=[];ub=[]; [x,fval]=fmincon(@fun22,x0,A,b,Aeq,beq,lb,ub,@cont2) x = 1.0e+008*[-0.0006 -2.7649] fval = -7.6432e+016 MATLAB程序: 计算结果: * 送肝捕孕冠钟良迁戎泼沦瓶评井谓弄氰碌区早谬谰坏埋吝怀愧建泽藻娱往数学实验第六次讲稿数学实验第六次讲稿 fmincon语句的具体用法 * 柴锗崖免暑游贩甭蚁拴傲俗聊间安鱼雇股魂皑顿抛答余颜夷厕斟堵禾碴宋数学实验第六次讲稿数学实验第六次讲稿 function fm=ex02(x) n=10; c=[-6.089,-17.164,-34.054,-5.914, -24.721,-14.986,-24.1,-10.708, -26.662,-22.179]; sx=0; for i=1:n sx=sx+x(i); % end fm=0; for i=1:n fm=fm+x(i)*(c(i)+log(x(i)/sx)); end 1) 建立目标函数 ex02.m * 詹轰摆畅论匀蝎伦蜜铀较勺叛壶耪芭驭渗潜硕描运脱预余稚苑邓七肺蒋砌数学实验第六次讲稿数学实验第六次讲稿 3、使用quadprog求解二次规划问题 二次规划标准模型 调用格式: [x ,fval] = quadprog(H,c,A,b,Aeq,beq,L,U,x0) MATLAB软件求解 * 且溺靡准辙种潜尧嵌瘪秦绵圆盔撞永羹掖氨泼诉酶鹤遗邪荧淬肌鳖鸭郸肩数学实验第六次讲稿数学实验第六次讲稿 例4 写成标准模型 MATLAB软件求解 beq=2 * 诬矩嫁叭蔬恩捶揽烬搭奸赶笋冻交篙芥幂负嘲短率牌屿姥职招痢莉运渺撑数学实验第六次讲稿数学实验第六次讲稿 H=[2,-2;-2,4]; c=[-4,-12]; A=[-1,2 ;2,1]; b=[2,3]’; Aeq=[1 1];beq=2; [x,fval]=quadprog(H,c,A,b,Aeq,beq) 计算结果： x = [0.6667 1.3333]， f = -16.4444 MATLAB软