ch20-最大最小化.doc

第20章 最大最小化 算法 通常我们遇到的都是目标函数的最大和最小化问题，但是在某些情况下，则要求使最 大值最小化才有意义。例如城市规划中需要确定急救中心，消防中心的位置，可取的目标函数应该是到最远点的距离为最小，而不是到所有目的地的距离和为最小，这是两种完全不同的准则，在控制理论、逼近论、决策论中也使用最大最小化原则。 最大最小化问题的数学模型为 式中x,b,beq,lb和ub为向量，A和Aeq为矩阵，c(x),ceq(x)和F(x)为函数，返回向量。 F(x),c(x)和ceq(x)可以是非线性函数。 MATLAB优化工具箱中采用序列二次规划法求解最大最小化问题。 有关函数介绍 fminmax使多目标函数中的最坏情况达到最小化。给定初值估计，改值必须服从一定的约束条件。其调用格式为： x=fminimax(fun,x0)初值为x0，找到fun函数的最大最小化解 。 x=fminimax（fun,x0,A,b）给定线性不等式A*x=b,求解最大最小化问题。 x=fminimax(fun,x,A,b,Aeq,beq)给定线性等式，Aeq*x=beq,求解最大最小化问题。如果没有不等式存在，则设置A=[ ],b=[ ]. x=fminimax(fun,x,A,b,Aeq,beq,lb,ub),为设计变量定义一系列下限lb和上限ub,使得总有lb=x=ub. x=fminimax(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)在nonlcon参数中给定非线性不等式约束c(x)或等式约束ceq(x)，fminimax函数要求c(x)=0且ceq(x)=0. 若没有边界存在，则设置lb=[ ]和(或)ub=[ ]. x=fminimax(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options)用options给定的参数进行优化。 x=fminimax(fun,x0,A,b,Aeq,beq,lb,ub,nonicon,options,P1,P2,…)将问题参数P1,P2等直接传递给函数fun和nonlcon。如果不需要变量A,b，Aeq，beq，lb，ub，nonlcon和options则将它们设置为空矩阵。 [x,fval]=fminimax(…)返回解x处的目标函数值。 [x,fval,maxfval]=fminimax(…)返回解x处的最大函数值。 [x,fval,maxfval,exitflag]=fminimax(…)返回exitflag参数，描述函数计算的退出条件。 [x,fval,maxfval,exitflag,output]=fminimax(…)返回描述优化信息的结构输出output参数。 [x,fval,maxfval,exitflag,output,lambda]=fminimax(…)返回包含解x处拉格朗日乘子的lambda参数。 调用格式中的maxfval变量为解x处函数值的最大值，即，maxfval=max{fun(x)}. fminimax函数使用序列二次规划法（SQp)进行计算。对一维搜索法和Hess矩阵的计算进行了修改。在一维搜索中，将精确目标函数和另外目标函数一起使用。当有一个目标函数不在发生改善时，一维搜索终止。修改的Hess矩阵借助于本问题的结构，也被采用。 使用fminmax函数时需要注意下面几个问题： （1）在options,MinAbsMax中设置F最大绝对值最小化了的目标数。该目标应该放到F的第一个元素中去。例如，考虑上面(??)的问题，需要找到x值，使下式的最大绝对值最小化; 通过调用fminimax函数来求解 x0=[0.1;0.1]; % 设置初值 options=optimset(‘MinAbsMax’,5); % 使最大绝对值最小化 [x,fval]=fminmax(fun,x0,[],[] ,[],[],[],[],[],options) 经过7次迭代以后，得到问题的解 ?? （2）当提供了等式约束并且在二次子问题中发现并剔除了因变等式时，则在过程标题中打印 字样（当输出选项设置为 ）。因变等式只是在等式连续的情况下才被剔除。若系统不连续，则子问题不可行并且在过程标题中打印 字样。 另外，要求目标函数必须连续，否限fminimax函数可能给出局部最优解。 应用实例 【例20-1】 定位问题 设某城市有某种物品的10个需求点，第i个需求点的坐标为（），道路网与坐标轴平行，彼此正交。现打算建一个该物品的供应中心，且由于受到城市某些条件的限制，该供应中心只能设在x界于[5，8]， y界于[5，8]的范围内。问该中心应建在何处为好？ 点的坐标为：