case36.m.bak

比较经典的matlab的电力系统潮流计算程序

function [baseMVA, bus, gen, branch, area, gencost] = case
%CASE    Defines the power flow data in a format similar to PTI.
%   [baseMVA, bus, gen, branch, area, gencost] = case
%   The format for the data is similar to PTI format except where noted.
%   An item marked with (+) indicates that it is included in this data
%   but is not part of the PTI format. An item marked with (-) is one that
%   is in the PTI format but is not included here.
%
%   Bus Data Format
%       1   bus number (1 to 29997)
%       2   bus type
%               PQ bus          = 1
%               PV bus          = 2
%               reference bus   = 3
%               isolated bus    = 4
%       3   Pd, real power demand (MW)
%       4   Qd, reactive power demand (MVAR)
%       5   Gs, shunt conductance (MW (demanded?) at V = 1.0 p.u.)
%       6   Bs, shunt susceptance (MVAR (injected?) at V = 1.0 p.u.)
%       7   area number, 1-100
%       8   Vm, voltage magnitude (p.u.)
%       9   Va, voltage angle (degrees)
%   (-)     (bus name)
%       10  baseKV, base voltage (kV)
%       11  zone, loss zone (1-999)
%   (+) 12  maxVm, maximum voltage magnitude (p.u.)
%   (+) 13  minVm, minimum voltage magnitude (p.u.)
%
%   Generator Data Format
%       1   bus number
%   (-)     (machine identifier, 0-9, A-Z)
%       2   Pg, real power output (MW)
%       3   Qg, reactive power output (MVAR)
%       4   Qmax, maximum reactive power output (MVAR)
%       5   Qmin, minimum reactive power output (MVAR)
%       6   Vg, voltage magnitude setpoint (p.u.)
%   (-)     (remote controlled bus index)
%       7   mBase, total MVA base of this machine, defaults to baseMVA
%   (-)     (machine impedance, p.u. on mBase)
%   (-)     (step up transformer impedance, p.u. on mBase)
%   (-)     (step up transformer off nominal turns ratio)
%       8   status, 1 - machine in service, 0 - machine out of service
%   (-)     (% of total VARS to come from this gen in order to hold V at
%               remote bus controlled by several generators)
%       9   Pmax, maximum real power output (MW)
%       10  Pmin, minimum real power output (MW)
%
%   Branch Data Format
%       1   f, from bus number
%       2   t, to bus number
%   (-)     (circuit identifier)
%       3   r, resistance (p.u.)
%       4   x, reactance (p.u.)
%       5   b, total line charging susceptance (p.u.)
%       6   rateA, MVA rating A (long term rating)
%       7   rateB, MVA rating B (short term rating)
%       8   rateC, MVA rating C (emergency rating)
%       9   ratio, transformer off nominal turns ratio ( = 0 for lines )
%           (taps at 'from' bus, impedance at 'to' bus, i.e. ratio = Vf / Vt)
%       10  angle, transformer phase shift angle (degrees)
%   (-)     (Gf, shunt conductance at from bus p.u.)
%   (-)     (Bf, shunt susceptance at from bus p.u.)
%   (-)     (Gt, shunt conductance at to bus p.u.)
%   (-)     (Bt, shunt susceptance at to bus p.u.)
%       11  initial branch status, 1 - in service, 0 - out of service
%
% (+) Area Data Format
%       1   i, area number
%       2   price_ref_bus, reference bus for that area
% 
% (+) Generator Cost Data Format
%       NOTE: If gen has n rows, then the first n rows of gencost contain
%       the cost for active power produced by the corresponding generators.
%       If gencost has 2*n rows then rows n+1 to 2*n contain the reactive
%       power costs in the same format.
%       1   model, 1 - piecewise linear, 2 - polynomial
%       2   startup, startup cost in US dollars
%       3   shutdown, shutdown cost in US dollars
%       4   n, number of cost coefficients to follow for polynomial
%           (or data points for piecewise linear) total cost function
%       5 and following, cost data, piecewise linear data as:
%                   x0, y0, x1, y1, x2, y2, ...
%           and polynomial data as, e.g.:
%                   c2, c1, c0
%           where the polynomial is c0 + c1*P + c2*P^2
%
% << this file created [97-Aug-26 12:29:15] by PB::System version 1.3 >>

%%-----  Power Flow Data  -----%%
%% system MVA base
baseMVA = 100.0000;

%% bus data
bus = [
	1	3	0.0	0.0	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	2	2	0.0	0.0	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	3	2	0.0	0.0	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	4	1	0.0	0.0	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	5	1	90.0000	30.0000	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	6	1	0.0	0.0	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	7	1	100.0000	35.0000	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	8	1	0.0	0.0	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
	9	1	125.0000	50.0000	0.0	0.0	1	1.0000	0.0000	345.0000	1	1.1000	0.9000;
];

%% generator data
gen = [
1	0.0	0.0	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;
2	6.0	3.6	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;	
3	3.1	0.0	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;
4	1.6	0.7	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;
5	4.3	3.34	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;	
6	-0.01	0.0	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;
7	2.25	0.0	300.0000	-300.0000	1.0000	100.0000	1	250.0000	10.0000;
];

%% branch data
branch = [
	1	24	0.0	0.015	0.0	220.0000	220.0000	220.0000	1.075	0	1;
	2	9	0.0	0.0217	0.0	220.0000	220.0000	220.0000	1.075	0	1;
	3	22	0.0	0.0124	0.0	220.0000	220.0000	220.0000	1.075	0	1;
	4	19	0.0	0.064	0.0	220.0000	220.0000	220.0000	1.025   0	1;
	5	18	0.0	0.0375	0.0	220.0000	220.0000	220.0000	1.075   0	1;
	6	17	0.0	0.0337	0.0	220.0000	220.0000	220.0000	1.0     0	1;
	7	30	0.0	0.0438	0.0	220.0000	220.0000	220.0000	1.025   0	1;
	8	31	0.0	0.0328	0.0	220.0000	220.0000	220.0000	1.025   0	1;
	10	9	0.0	-0.002	0.0	220.0000	220.0000	220.0000	1.0	0	1;
	10	11	0.0	0.018	0.0	220.0000	220.0000	220.0000	1.025   0	1;
	15	12	0.0	0.018	0.0	220.0000	220.0000	220.0000	1.025   0	1;
	17	13	0.0	0.01	0.0	220.0000	220.0000	220.0000	1.025   0	1;
	15	14	0.0	-0.002	0.0	220.0000	220.0000	220.0000	1.0	0	1;
	17	16	0.0	0.001	0.0	220.0000	220.0000	220.0000	1.027	0	1;
	9	22	0.0559	0.218	0.1954	220.0000	220.0000	220.0000	0	0	1;
	9	23	0.0034  0.0131	0.0	220.0000	220.0000	220.0000	0	0	1;
	9	24	0.0147	0.104	0.0	220.0000	220.0000	220.0000	0	0	1;
	11	25	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	12	26	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	12	27	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	13	28	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	14	19	0.0034	0.02	0.0	220.0000	220.0000	220.0000	0	0	1;
	16	18	0.0033	0.0333	0.0	220.0000	220.0000	220.0000	0	0	1;
	16	19	0.0578	0.218	0.1887	220.0000	220.0000	220.0000	0	0	1;
	16	20	0.0165	0.0662	0.2353	220.0000	220.0000	220.0000	0	0	1;
	16	21	0.0374	0.178	0.164	220.0000	220.0000	220.0000	0	0	1;
	16	29	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	19	21	0.0114	0.037	0.0	220.0000	220.0000	220.0000	0	0	1;
	19	30	0.0196	0.854	0.08	220.0000	220.0000	220.0000	0	0	1;	
	20	22	0.0214	0.0859	0.3008	220.0000	220.0000	220.0000	0	0	1;
	21	22	0.015	0.0607	0.2198	220.0000	220.0000	220.0000	0	0	1;
	22	23	0.0537	0.19	0.1653	220.0000	220.0000	220.0000	0	0	1;
	23	24	0.0106	0.074	0.0	220.0000	220.0000	220.0000	0	0	1;
	25	26	0.0033	0.0343	1.88	220.0000	220.0000	220.0000	0	0	1;
	27	28	0.00245	0.0255	1.395	220.0000	220.0000	220.0000	0	0	1;
	16	0	0.0	0.0	1.0	220.0000	220.0000	220.0000	0	0	1;
	25	0	0.0	0.0	-1.3665	220.0000	220.0000	220.0000	0	0	1;
	26	0	0.0	0.0	-1.3665	220.0000	220.0000	220.0000	0	0	1;
	27	0	0.0	0.0	-1.3665	220.0000	220.0000	220.0000	0	0	1;
	30	31	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	31	33	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	29	32	0.0	0.0001	0.0	220.0000	220.0000	220.0000	0	0	1;
	32	33	0.01	0.04	0.20	220.0000	220.0000	220.0000	0	0	1;
];

%%-----  OPF Data  -----%%
%% area data
area = [
	1	5;
];

%% generator cost data
gencost = [
	2	1500.00	0.00	3	0.11	5	150;
	2	2000.00	0.00	3	0.085	1.2	600;
	2	3000.00	0.00	3	0.1225	1	335;
];

return;