minopf.m

一种基于MATLAB的多目标最优化仿真工具箱

% MINOPF: Optimal Power Flow using MINOS; don't call directly, use mopf instead
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u);
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u,
%             mpopt);
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u,
%             mpopt, jacoption);
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u,
%             mpopt, N, FPARM, H, Cw);
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u,
%             mpopt, N, FPARM, H, Cw, jacoption);
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u,
%             mpopt, N, FPARM, H, Cw, z0, zl, zu);
%
% [bus, gen, branch, f, info, g, jac, xr, pimul] = ...
%             minopf(baseMVA, bus, gen, gencost, branch, col, A, l, u,
%             mpopt, N, FPARM, H, Cw, z0, zl, zu, jacoption);
%
% MINOPF is a MEX file that attempts to solve the optimal power flow problem
% specified in matrices baseMVA, bus, gen, branch, and gencost.  Matrices
% A, l and u impose additional linear constraints on the optimization variables,
% so that l <= A*x <= u .  Vector col is a vector of indexes into the columns
% of the data tables, organized as follows:
%
% position
%  in
%  col   description                                data type
% -------------------------------------------------------------------
%  1     PQ: value that identifies a PQ bus          bus classification
%  2     PV: value that identifies a PV bus          bus classification
%  3     REF: value that identifies reference bus    bus classification
%  4     NONE: value for unknown type of bus         bus classification
%  5     BUS_I: index of column holding bus numbers  index into bus table
%  6     BUS_TYPE: idx of col. holding bus types     index into bus table
%  7     PD: index of column holding bus load        index into bus table
%  8     QD: index of column holding bus Q load      index into bus table
%  9     GS: index of column holding bus shunt GS    index into bus table
% 10     BS: index of column holding bus shunt BS    index into bus table
% 11     BUS_AREA: idx of col. holding area          index into bus table
% 12     VM: index of column holding voltage magnit  index into bus table
% 13     BASE_KV: index of column holding base volta index into bus table
% 14     ZONE: index of column holding zone number   index into bus table
% 15     VMAX
% 16     VMIN
% 17     LAM_P
% 18     LAM_Q
% 19     MU_VMAX
% 20     MU_VMIN
% 21     GEN_BUS
% 22     PG
% 23     QG
% 24     QMAX
% 25     QMIN
% 26     VG
% 27     MBASE
% 28     GEN_STATUS
% 29     PMAX
% 30     PMIN
% 31     MU_PMAX
% 32     MU_PMIN
% 33     MU_QMAX
% 34     MU_QMIN
% 35     PW_LINEAR
% 36     POLYNOMIAL
% 37     MODEL
% 38     STARTUP
% 39     SHUTDOWN
% 40     NCOST
% 41     COST
% 42     F_BUS
% 43     T_BUS
% 44     BR_R
% 45     BR_X
% 46     BR_B
% 47     RATE_A
% 48     RATE_B
% 49     RATE_C
% 50     TAP
% 51     SHIFT
% 52     BR_STATUS
% 53     PF
% 54     QF
% 55     PT
% 56     QT
% 57     MU_SF
% 58     MU_ST

%   MINOPF for MATPOWER
%   $Id: minopf.m,v 1.5 2007/07/05 19:55:31 ray Exp $
%   by Carlos E. Murillo-Sanchez, PSERC Cornell & Universidad Autonoma de Manizales
%   Copyright (c) 2000-2004 by Power System Engineering Research Center (PSERC)
%   See http://www.pserc.cornell.edu/minopf/ for more info.

c1 = [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
      VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN];
c2 = [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
        GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN];
c3 = [PW_LINEAR, POLYNOMIAL, MODEL, STARTUP, SHUTDOWN, NCOST, COST];
c4 = [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, ...
        RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST];
col = [ c1 c2 c3 c4];


%[buso, geno, brancho, f, info, g, jac] = minopf(baseMVA, bus, gen, ...
%                                     gencost, branch, col, A, l, u, option);