MATLAB解二次规划代码.pdfVIP

黑龙江八一农垦大学 数学建模选修课 MATLAB MATLAB 第八章 MMAATTLLAABB解二次非线性规划 一、数学模型 1 T T min x Hx+ f x x 2 A⋅ x≤ b Aeq⋅x= beq lb≤ x≤ ub 其中H为二次型矩阵，A、Aeq分别为不等式约束与等式约束系数矩阵，f,b,beq,lb,ub,x 为向量。 二、求解二次规划问题函数quadprog 调用格式： X= quadprog(H,f,A,b) X= quadprog(H,f,A,b,Aeq,beq) X= quadprog(H,f,A,b,Aeq,beq,lb,ub) X= quadprog(H,f,A,b,Aeq,beq,lb,ub,x0) X= quadprog(H,f,A,b,Aeq,beq,lb,ub,x0,options) [x,fval]= quadprog(…) [x,fval,exitflag]= quadprog(…) [x,fval,exitflag,output]= quadprog(…) [x,fval,exitflag,output,lambda]= quadprog(…) 说明：输入参数中，x0为初始点；若无等式约束或无不等式约束，就将相应的矩阵和向量 设置为空；options为指定优化参数。输出参数中，x是返回最优解；fval是返回解所对应的 目标函数值；exitflag是描述搜索是否收敛；output是返回包含优化信息的结构。Lambda是 返回解x入包含拉格朗日乘子的参数。 三、举例 例1：求解：二次规划问题 min f(x x)= -2x-6x+x -2xx+2x2 2 1, 2 1 2 1 1 2 2 s.t x+x≤2 1 2 -x+2x≤21 2 x≥0, x≥01 2 将其写成标准形式： T 2 − 2 x − 2 x 1 ⎡ ⎤⎡ ⎤ ⎡ ⎤ ⎡ ⎤ min f(x,x )= (x,x ) 1 + 1 1 2 1 2 ⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥ 2 − 2 4 x −6 x ⎣ ⎦ 2 ⎣ ⎦ 2 ⎣ ⎦ ⎣ ⎦ s.t x+x≤2 1 2 -x+2x≤21 2 x≥0, x≥01 2 程序： 编写M文件mainG.m如下： H=[1 -1;-1 2]; f=[-2;-6]; 1 黑龙江八一农垦大学 数学建