交直流程序

function [Bbus Bf Pbusinj Pfinj] = makeBdc（baseMVA bus branch）
%MAKEBDC   Builds the B matrices and phase shift injections for DC power flow.
%   [BBUS BF PBUSINJ PFINJ] = MAKEBDC（MPC）
%   [BBUS BF PBUSINJ PFINJ] = MAKEBDC（baseMVA BUS BRANCH）
%
%   Returns the B matrices and phase shift injection vectors needed for
%   a DC power flow. The bus real power injections are related to bus
%   voltage angles by 
%       P = BBUS * Va + PBUSINJ
%   The real power flows at the from end the lines are related to the bus
%   voltage angles by
%       Pf = BF * Va + PFINJ
%   Does appropriate conversions to p.u.
%   Bus numbers must be consecutive beginning at 1 （i.e. internal ordering）.
%


%   Example:
%       [Bbus Bf Pbusinj Pfinj] = makeBdc（baseMVA bus branch）;
%
%   See also DCPF.

%   MATPOWER
%   Copyright （c） 1996-2016 Power Systems Engineering Research Center （PSERC）
%   by Carlos E. Murillo-Sanchez PSERC Cornell & Universidad Nacional de Colombia
%   and Ray Zimmerman PSERC Cornell
%
%   This file is part of MATPOWER.
%   Covered by the 3-clause BSD License （see LICENSE file for details）.
%   See http://www.pserc.cornell.edu/matpower/ for more info.

%% extract from MPC if necessary
if nargin < 3
    mpc     = baseMVA;
    baseMVA = mpc.baseMVA;
    bus     = mpc.bus;
    branch  = mpc.branch;
end

%% constants
nb = size（bus 1）;          %% number of buses
nl = size（branch 1）;       %% number of lines

%% define named indices into bus branch matrices
[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] = idx_bus;
[F_BUS T_BUS BR_R BR_X BR_B RATE_A RATE_B RATE_C ...
    
    TAP SHIFT BR_STATUS PF QF ~ QT MU_SF MU_ST ...
    ANGMIN ANGMAX MU_ANGMIN MU_ANGMAX]