x1quadprog.m

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

% x1QuadProg.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  1/2(5x1^2 + 6x1x2 + 5x2^2)-95x1-105x2
% s.t. x1 + x2 <= 10
%      3x1 + x2 <= 15
%      2x1 + 3x2 <=30
%      -15x1 + 13x2 <=0
%      x1>=0, x2>=0
%
% 注意：原问题不是标准形式，要先把原问题转化为二次规划标准形式
%   min 0.5*x'Hx + f'x
%   s.t. Ax <= b
%        LB <= X <= UB
%   即可确定有关参数：H,f,A,b,lb,ub

clear all
clc

H = [5 3; 3 5];
f = [-95; -105];    
A = [1  2; 3  1; 2  3; -15  13];
b = [10; 15; 30; 0];
lb = [0; 0];
ub = [inf; inf];

[x,fval,exitflag,output,lambda] = quadprog(H,f,A,b,[],[],lb,ub)
% Results: x = [4.0000; 3.0000]；minimum value: -596.5000