调用格式：
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,output]=quadprog(…) [x,fval,exitflag,output,lambda]=quadprog(…)说明：输入参数中，x0为初始点；若无等式约束或无不等式约束，就将相应的矩阵和向量设置为空；options为指定优化参数。输出参数中，x是返回最优解；fval是返回解所对应的目标函数值；exitflag是描述搜索是否收敛；output是返回包含优化信息的结构。Lambda是返回解x入包含拉格朗日乘子的参数。
下面代码操作流程： 利用hessian 函数计算hessian矩阵H，然后利用collect获取一次项的系数F：
clear all,close all;syms x y ;% 求解该函数的最小值f = (x - 2)^2 + (y - 3)^2 + 5; % s.t. 1 <= x <= 3,2 <= y <= 4%lb = [1 2]';ub = [3 4]'; % lb = [-inf -inf]' ub = [inf inf]'lb = zeros(2,1);ub = zeros(2,1);lb(:) = -inf;ub(:) = inf;%% 求解海森矩阵 HH = hessian(f,[x,y]);% convert to double typeH = double(H);%% 求解一次项系数 F% fexp = expand(f);fcol = collect(f,{'x','y'}); % x^2 - 4*x + y^2 - 6*y + 18disp(fcol);% get F = [-4 -6]'F = [-4 -6]';% method : interior-point-convex ,trust-region-reflective ,active-setoptions = optimoptions('quadprog','Algorithm','interior-point-convex');[x,fval,exitflag,output] = quadprog(H,F,[],[],[],[],lb,ub,[],options);新闻热点
疑难解答
