fm_inilf.m

电力系统的psat

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-2005 Federico Milano
%
% This toolbox is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2.0 of the License, or
% (at your option) any later version.
%
% This toolbox is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANDABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
% General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this toolbox; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307,
% USA.

fm_var

% Deleting Self-generated functions
clear fm_call

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

% OPF, CPF, LIB & SNB variables
OPF.report = [];
OPF.uname = [];
OPF.fname = [];
OPF.varout = '';
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.fV = '';
Varname.uV = '';
Varname.fang = '';
Varname.uang = '';
Varname.fPm = '';
Varname.uPm = '';
Varname.fVf = '';
Varname.uVf = '';
Varname.fVref = '';
Varname.uVref = '';
Varname.fIf = '';
Varname.uIf = '';
Varname.fP = '';
Varname.uP = '';
Varname.fQ = '';
Varname.uQ = '';
Varname.fPflow = '';
Varname.uPflow = '';
Varname.fQflow = '';
Varname.uQflow = '';
Varname.fx = '';
Varname.ux = '';
Varname.fdet = '';
Varname.udet = '';
Varname.fautostate = '';
Varname.uautostate = '';
Varname.fautojlfr = '';
Varname.uautojlfr = '';
Varname.fautojlfdr = '';
Varname.uautojlfdr = '';
Varname.fpflf = '';
Varname.upflf = '';
Varname.fpflfd = '';
Varname.upflfd = '';

% 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.Fx  = [];
DAE.Fy = [];

% 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 = [];

% 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 = [];

% Slack variables
SW.con = [];
SW.n = 0;
SW.bus = [];

% PV variables
PV.con = [];
PV.n = 0;
PV.bus = [];
PV.pq.con = [];
PV.pq.n = 0;
PV.pq.bus = [];

% PQ  variables
PQ.con = [];
PQ.n = 0;
PQ.bus = [];
PQ.P0 = [];
PQ.Q0 = [];

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

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

% Frequency Dependent Load Variables
Fl.con = [];
Fl.n = 0;
Fl.bus = [];
Fl.a0 = [];
Fl.Dw = [];
Fl.x = [];

% 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.Gp = [];
Syn.Gq = [];
Syn.c1 = [];
Syn.c2 = [];
Syn.c3 = [];

% 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 = [];

% 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.V = [];
Fault.ang = [];
Fault.delta = 0;

% 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.con = [];
Tg.n = 0;
Tg.syn  = [];
Tg.bus  = [];
Tg.dat1 = [];
Tg.dat2 = [];
Tg.tg = [];
Tg.tg1 = [];
Tg.tg2  = [];
Tg.tg3 = [];

% 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 = [];

% Static Var Compensator variables
Svc.con = [];
Svc.bus = [];
Svc.n = 0;
Svc.bcv = [];
Svc.alpha = [];
Svc.vm = [];

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

% TCSC  variables
Tcsc.con = [];
Tcsc.n = 0;
Tcsc.bus1 = [];
Tcsc.bus2 = [];
Tcsc.Pe = [];
Tcsc.B = [];
Tcsc.SI = [];
Tcsc.x1 = [];
Tcsc.x2 = [];
Tcsc.x3 = [];

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

% StatCom variables
Statcom.con = [];
Statcom.n = 0;
Statcom.bus = [];
Statcom.Id = [];
Statcom.Iq = [];
Statcom.Vdc = [];
Statcom.alpha = [];
Statcom.m = [];
Statcom.Vmdc = [];
Statcom.Vmac = [];

% SSSC variables
Sssc.con = [];
Sssc.n = 0;
Sssc.bus1 = [];
Sssc.bus2 = [];
Sssc.dat = [];
Sssc.Vdc = [];
Sssc.beta = [];
Sssc.m = [];
Sssc.Pmac = [];
Sssc.Vmdc = [];

% UPFC variables
Upfc.con = [];
Upfc.n = 0;
Upfc.bus1 = [];
Upfc.bus2 = [];
Upfc.dat = [];
Upfc.Vmdc = [];
Upfc.Vmac = [];
Upfc.Vdc = [];
Upfc.alpha = [];
Upfc.m = [];
Upfc.x1 = [];
Upfc.x2 = [];
Upfc.Imd = [];
Upfc.Imq = [];

% RLC variables
RLC.con = [];
RLC.bus1 = [];
RLC.bus2 = [];
RLC.dat = [];
RLC.n = 0;
RLC.id = [];
RLC.iq = [];
RLC.vd = [];
RLC.vq = [];
RLC.check = 0;

% 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
Demand.con = [];
Demand.n = 0;
Demand.bus = [];

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

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;

% ----------------------------------------------------------- %
%                       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.syn = [];
Cluster.svc = [];
Cluster.exc = [];
Cluster.cac = [];
Cluster.Vs = [];

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 = [];

Jimma.con = [];
Jimma.n = 0;
Jimma.bus = [];
Jimma.dat = [];
Jimma.x = [];

Mixed.con = [];
Mixed.n = 0;
Mixed.bus = [];
Mixed.dat = [];
Mixed.x = [];
Mixed.y = [];