数建模培训最优化.ppt

x=fmincon(‘fun’,x0,A,b)x=fmincon(‘fun’,x0,A,b,Aeq,beq)x=fmincon(‘fun’,x0,A,b,Aeq,beq,lb,ub)x=fmincon(‘fun’,x0,A,b,Aeq,beq,lb,ub,nonlcon)x=FMINCON(FUN,X0,A,B,Aeq,Beq,LB,UB,NONLCON,OPTIONS)在Matlab优化工具箱中，fmincon函数是用SQP算法来解决一般的约束非线性规划的函数，它的命令格式为：上式中x为最优点；若将左端的x换为[x,f],则返回最优点x和最优值f。第59页,共86页，星期六，2024年，5月【例1】求解约束非线性规划：(初值为[1;1])解:首先建立一个m文件fun1.mfunctiony=fun1(x)y=-exp(x(1))*x(2)^2*(3-exp(x(1))-x(2)^2);存储为fun1.m首先将问题转化为matlab要求的格式;即求出fun,A,b,Aeq,Beq,X0,Lb,Ub第60页,共86页，星期六，2024年，5月function[g,cep]=fun2(x)g=[];%g为非线性不等式,且为g=0ceq=exp(x(1))+x(2)^2-3;%ceq为非线性等式然后存储为fun2.m建立主程序：A=[];b=[];Aeq=[];Beq=[];Lb=[];Ub=[];[x,f]=fmincon(‘fun1’,[1;1],[],[],[],[],[],[],’fun2’)-f建立非线性约束m-文件fun2.m运行结果为:x=0.88520.7592f=6.2043e-016最优点最优值第61页,共86页，星期六，2024年，5月【例2】求解约束非线性规划：解：首先建立一个m文件fun5.mfunctiony=fun5(x)y=(x(1)-1)^2+(x(2)-2)^2+(x(3)-3)^2+(x(4)-4)^2;存储为fun5.m文件.第62页,共86页，星期六，2024年，5月x0=[1;1;1;1];A=[1111;3321];B=[5;10];Aeq=[];Beq=[];Lb=[0;0;0;0];[x,g]=fmincon(‘fun5’,x0,A,B,Aeq,Beq,Lb)运行结果为:x=0.00000.66671.66652.6668g=6.3333建立主程序第63页,共86页，星期六，2024年，5月小结：用Matlab求解非线性规划问题，基本步骤：1.首先建立M文件fun.m,定义目标函数f（x）:functionf=fun(x);f=f（x）;2.若约束条件中有非线性约束:g(x)或Ceq(x)=0,则建立M文件nonlcon.m定义函数g(x)与Ceq(x):function[g,Ceq]=nonlcon(X)g=...Ceq=...3.建立主程序.并运行。第64页,共86页，星期六，2024年，5月s.t.例如：在对策论中：在最不利的条件下，寻求最有利的策略；在投资规划中要确定最大风险的最低限度；在城市规划中，要确定急救中心的位置，使其到所有地点最大距离为最小。最大最小化问题第65页,共86页，星期六，2024年，5月求解最大最小化问题的Matlab函数为fminimax.其调用格式如下：x=fminimax(F,x0,,A,b,Aeq,beq,lb,ub,nonlcon,options)或[x,fval,maxfval,exitflag,output]=fminimax(…)其中：x返回最优解；fval返回解x处的目标函数值；maxfval返回解x处的最大函数值；exitflag描述计算的退出条件；output返回包含优化信息的输出参数。第66页,共86页，星期六，2024年，5月例：求解下列最大最小化问题：首先编写一个M文件ff2.m，计算4个函数值。functionf=ff2(x)f(1)=3*x(1)^2+2*x(2)^2-12*x(1)+35;f(2)=5*x(1)*x(2)-4*x(2)+7;f(3)=x(1)^2+6*x(2);f(4)=4*x(1)^2+9*x(2)^2-12*x(1)*x(2)+20;第67页,共86页，星期六，2024年，5月然后，输入初值x0=(1,1)，并调用优化函