ex_infeas_denseformat.m

来自「HERE IS A GOOD PSO TOOL BOX」· M 代码 · 共 39 行

%----------------------------------------------------------------------
%
%       INFEASIBLE GP EXAMPLE
%
%----------------------------------------------------------------------
%
%   minimize:   x^-2 + y^-2
%   subject to: x^-1 <= 1
%               y^-1 <= 1
%               x + y <= 0.9
%
%   This problem has no feasible set.

clear all; close all;

%----------------------------------------------------------------------
%      PROBLEM DATA IN POSYNOMIAL FORM (DENSE FORMAT)
%----------------------------------------------------------------------

A   = [-2   0 ;...
        0  -2 ;...
       -1   0 ;...
        0  -1 ;...
        1   0 ;...
        0   1 ];
b   = [ 1 1 1 1 1/0.9 1/0.9]';
szs = [ 2 1 1 2]';

%----------------------------------------------------------------------
%      SOLVE THE PROBLEM IN POSYNOMIAL FORM
%----------------------------------------------------------------------
[x1,status1,lambda1,nu1] = gpposy(A,b,szs);

%----------------------------------------------------------------------
%      SOLVE THE PROBLEM IN CONVEX FORM
%----------------------------------------------------------------------
[x2,status2,lambda2,nu2,mu2] = gpcvx(A,log(b),szs);