x2quadprog.m

来自「《实用化工计算机模拟：MATLAB在化学工程中的应用 》这本书光盘里的程序~」· M 代码 · 共 42 行

% x2QuadProg.m
% 二次规划（优化）问题的求解例子
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2002/05/31 $
%
% 优化问题：
% min  2.5x2^2 -35x1 - 56x2
% s.t. -3x1 + x2 = 7
%      8x1 + 3x2 <= 12
%      -3x1 + 4x2 >= -4
%
% 注意：原问题不是标准形式，要先把原问题转化为二次规划标准形式
% （将原第三个约束条件两边乘以-1）
%   min 0.5*x'Hx + f'x
%   s.t. Ax <= b
%        Aeq*x = beq
%        lb <= x <= ub
%   即：
%   min  1/2x'Hx+fx
%       X = [x1; x2]
%       H = [0 0; 0 5];
%       f = [-35; -56];
%      8x1 + 3x2 <= 12
%      3x1 - 4x2 <= 4

clear all
clc

H = [0 0; 0 5];
f = [-35; -56];    
A = [8 3; 3 -4];
b = [12; 4];
Aeq = [-3 1];
beq = [7];

[x,fval,exitflag,output,lambda] = quadprog(H,f,A,b,Aeq,beq)
% Results: x = [-0.5294 5.4118]；minimum value: -211.3114