fm_inilf.m

这是一个很适合研究和学习用的电力系统仿真软件

function fm_inilf
% FM_INILF initialize all system and component variables
%          for power flow computations
%
% FM_INILF
%
%see also FM_SLF
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Update:    22-Feb-2003
%Update:    09-Jul-2003
%Version:   1.0.2
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2006 Federico Milano

fm_var

% Deleting Self-generated functions
clear fm_call

% General Variables
DAE.V = [];
DAE.a = [];
DAE.kg = 0;
Settings.iter = 0;

% OPF, CPF, LIB & SNB variables
OPF.report = [];
OPF.init = 0;
OPF.line = 0;
CPF.init = 0;
LIB.lambda = 0;
LIB.dldp = [];
LIB.bus = [];
LIB.init = 0;
SNB.lambda = 0;
SNB.dldp = [];
SNB.bus = [];
SNB.init = 0;

% Variable Name Structure
Varname.bus = '';
Varname.uvars = '';
Varname.fvars = '';
Varname.nvars = 0;
Varname.idx = [];
Varout.t = [];
Varout.vars = [];
Varout.idx = [];

% State variables
DAE.t = -1;
DAE.x = [];
DAE.f = [];
DAE.n= 0;
DAE.npf = 0;

% power flow vectors & matrices
DAE.g = [];
DAE.gp = [];
DAE.glfp = [];
DAE.gq = [];
DAE.glfq = [];
DAE.J11 = [];
DAE.J12 = [];
DAE.J21 = [];
DAE.J22 = [];
DAE.Jlf = [];
DAE.Jlfv = [];
DAE.Jlfd = [];
DAE.Gx = [];
DAE.Gl = [];
DAE.Gk = [];
DAE.Fx = [];
DAE.Fy = [];
DAE.Fl = [];
DAE.Fk = [];

% small signal stability analysis parameters
SSSA.pf = [];
SSSA.eigs = [];
SSSA.report = [];

% Phasor Measurment Units
PMU.number = 0;
PMU.voltage = '';
PMU.angle = '';
PMU.location = '';
PMU.report = '';
PMU.measv = 0;
PMU.measc = 0;
PMU.pseudo = 0;
PMU.noobs = 0;

% do not clear component variables if the file is not read again
if ~clpsat.readfile, return, end

% Bus variables
Bus.con = [];
Bus.n = 0;
Bus.int = [];
Bus.Pl = [];
Bus.Ql = [];
Bus.Pg = [];
Bus.Qg = [];
Bus.island = [];

% Area variables
Area.con = [];
Area.n = 0;
Area.slack = [];

% Line variables
Line.con= [];
Line.n = 0;
Line.Y = [];
Line.Bp = [];
Line.Bpp = [];
Line.from= [];
Line.to= [];

% Shunt variables
Shunt.con = [];
Shunt.bus = [];
Shunt.g = [];
Shunt.b = [];

% Three-Winding Transformers
Twt.con = [];

% initialize components
SW = init(SW);
PV = init(PV);
PQ = init(PQ);
Fl = init(Fl);
Demand = init(Demand);
Supply = init(Supply);

% Polinomial Load  variables
Pl.con = [];
Pl.n = 0;
Pl.bus = [];
Pl.init = [];

% Monomial Load  variables
Mn.con = [];
Mn.n = 0;
Mn.bus = [];
Mn.init = [];

% Synchronous Machine  variables
Syn.con = [];
Syn.n = 0;
Syn.bus = [];
Syn.Id = [];
Syn.Iq = [];
Syn.Pg = [];
Syn.Qg = [];
Syn.J11 = [];
Syn.J12 = [];
Syn.J21 = [];
Syn.J22 = [];
Syn.delta = [];
Syn.omega = [];
Syn.e1q = [];
Syn.e1d = [];
Syn.e2q = [];
Syn.e2d = [];
Syn.psiq = [];
Syn.psid = [];
Syn.pm = [];
Syn.vf = [];
Syn.Pg0 = [];
Syn.Gp = [];
Syn.Gq = [];
Syn.c1 = [];
Syn.c2 = [];
Syn.c3 = [];
Syn.cluster = [];

% Mass variables
Mass.con = [];
Mass.n = 0;
Mass.syn = [];
Mass.delta_HP = [];
Mass.omega_HP = [];
Mass.delta_IP = [];
Mass.omega_IP = [];
Mass.delta_LP = [];
Mass.omega_LP = [];
Mass.delta_EX = [];
Mass.omega_EX = [];

% Excitation Control  variables
Exc.con = [];
Exc.n = 0;
Exc.syn = [];
Exc.bus = [];
Exc.vrif = [];
Exc.vrif0 = [];
Exc.vr1  = [];
Exc.vr2 = [];
Exc.vr3 = [];
Exc.vm = [];
Exc.vf = [];
Exc.cluster = [];

% Subtransmission Area Equivalents  variables
SAE1.con = [];
SAE1.m = [];
SAE2.con = [];
SAE3.con = [];
SAE1.n = 0;
SAE2.n = 0;
SAE3.n = 0;
SAE2.m1 = [];
SAE2.m2 = [];
SAE3.m1 = [];
SAE3.m2 = [];
SAE3.m3 = [];

% Fault  variables
Fault.con = [];
Fault.dat = [];
Fault.n = 0;
Fault.bus = [];
Fault.ang = [];

% Breaker  variables
Breaker.con = [];
Breaker.n = 0;
Breaker.line = [];
Breaker.bus = [];

% Tap changer & load  variables
Tap.con = [];
Tap.n = 0;
Tap.bus = [];
Tap.m = [];

% Asynchronous Motor  variables
Mot.con = [];
Mot.n = 0;
Mot.bus = [];
Mot.slip = [];
Mot.e1r = [];
Mot.e1m = [];
Mot.e2r = [];
Mot.e2m = [];

% Thermostatically controlled load variables
Thload.con = [];
Thload.n = 0;
Thload.bus = [];
Thload.T = [];
Thload.G = [];

% Turbine Governor  variables
Tg = init(Tg);

% On Load Tap Changer  variables
Ltc.con = [];
Ltc.bus1 = [];
Ltc.bus2 = [];
Ltc.n = 0;
Ltc.dat = [];
Ltc.m = [];

% Over Excitation Limiter  variables
Oxl.con = [];
Oxl.n = 0;
Oxl.bus = [];
Oxl.syn = [];
Oxl.exc = [];
Oxl.v = [];
Oxl.If = [];

% Power system stabilizer  variables
Pss.con = [];
Pss.n = 0;
Pss.v1 = [];
Pss.v2 = [];
Pss.v3 = [];
Pss.va = [];
Pss.vss = [];
Pss.Vs = [];
Pss.s1 = [];
Pss.bus = [];
Pss.syn = [];
Pss.exc = [];

% FACTS
Svc  = init(Svc);
Statcom = init(Statcom);
Tcsc = init(Tcsc);
Sssc = init(Sssc);
Upfc = init(Upfc);

% HVDC  variables
Hvdc.con = [];
Hvdc.n = 0;
Hvdc.bus1 = [];
Hvdc.bus2 = [];
Hvdc.dat = [];
Hvdc.Id = [];
Hvdc.xr = [];
Hvdc.xi = [];

% Lines variables
Lines.con = [];
Lines.n = 0;
Lines.bus1 = [];
Lines.bus2 = [];

% SSR variables
SSR.con = [];
SSR.bus = [];
SSR.n = 0;
SSR.Id = [];
SSR.Iq = [];
SSR.If = [];
SSR.Edc = [];
SSR.Eqc = [];
SSR.Tm = [];
SSR.Efd = [];
SSR.delta_HP = [];
SSR.omega_HP = [];
SSR.delta_IP = [];
SSR.omega_IP = [];
SSR.delta_LP = [];
SSR.omega_LP = [];
SSR.delta = [];
SSR.omega = [];
SSR.delta_EX = [];
SSR.omega_EX = [];

% Voltage stability & Market variables
Rsrv.con = [];
Rsrv.n = 0;
Rsrv.bus = [];

Rmpg.con = [];
Rmpg.n = 0;
Rmpg.sup = [];
Rmpg.bus = [];

Rmpl.con = [];
Rmpl.n = 0;
Rmpl.dem = [];
Rmpl.bus = [];

Vltn.con = [];
Vltn.n = 0;
Vltn.bus = [];

Ypdp.con = [];
Ypdp.day = [];
Ypdp.week = [];
Ypdp.year = [];
Ypdp.d = 1;
Ypdp.w = 1;
Ypdp.y = 1;
Ypdp.len = 1;

Pod.con = [];
Pod.n = [];
Pod.svc = [];
Pod.upfc = [];
Pod.tcsc = [];
Pod.sssc = [];
Pod.statcom = [];
Pod.idx = [];
Pod.type = [];
Pod.v1 = [];
Pod.v2 = [];
Pod.v3 = [];
Pod.Vs = [];
Pod.u = [];

% ----------------------------------------------------------- %
%                       W A R N I N G                         %
% ----------------------------------------------------------- %
% Following lines have been written by the UDM build utility. %
% This utility requires you do NOT change anything beyond     %
% this point in order to be able to correctly install and     %
% uninstall UDMs.                                             %
% ----------------------------------------------------------- %

Sofc.con = [];
Sofc.n = 0;
Sofc.bus = [];
Sofc.Ik = [];
Sofc.Vk = [];
Sofc.pH2 = [];
Sofc.pH20 = [];
Sofc.pO2 = [];
Sofc.qH2 = [];
Sofc.m = [];

Cac.con = [];
Cac.n = 0;
Cac.bus = [];
Cac.q1 = [];
Cac.q = [];

Cluster.con = [];
Cluster.n = 0;
Cluster.bus = [];
Cluster.svc = [];
Cluster.exc = [];
Cluster.cac = [];
Cluster.Vs = [];
Cluster.dVsdQ = [];
Cluster.u = [];

Exload.con = [];
Exload.n = 0;
Exload.bus = [];
Exload.dat = [];
Exload.xp = [];
Exload.xq = [];

Phs.con = [];
Phs.n = 0;
Phs.bus1 = [];
Phs.bus2 = [];
Phs.alpha = [];
Phs.Pm = [];

Wind.con = [];
Wind.n = 0;
Wind.speed = [];
Wind.vwa = [];
Wind.vw = [];

Cswt.con = [];
Cswt.n = 0;
Cswt.bus = [];
Cswt.dat = [];
Cswt.wind = [];
Cswt.omega_wr = [];
Cswt.omega_m = [];
Cswt.gamma = [];
Cswt.e1r = [];
Cswt.e1m = [];

Dfig.con = [];
Dfig.n = 0;
Dfig.bus = [];
Dfig.dat = [];
Dfig.wind = [];
Dfig.omega_m = [];
Dfig.theta_p = [];
Dfig.idr = [];
Dfig.iqr = [];

Ddsg.con = [];
Ddsg.n = 0;
Ddsg.bus = [];
Ddsg.dat = [];
Ddsg.wind = [];
Ddsg.omega_m = [];
Ddsg.theta_p = [];
Ddsg.ids = [];
Ddsg.iqs = [];
Ddsg.idc = [];

Busfreq.con = [];
Busfreq.n = 0;
Busfreq.bus = [];
Busfreq.dat = [];
Busfreq.x = [];
Busfreq.w = [];

Pmu.con = [];
Pmu.n = 0;
Pmu.bus = [];
Pmu.dat = [];
Pmu.vm = [];
Pmu.thetam = [];

Mixload = init(Mixload);
Jimma = init(Jimma);