📄 psat.m
字号:
% Frequency Dependent LoadFl = struct('con',[], ... 'n',0, ... 'bus',[], ... 'a0',[], ... 'Dw',[], ... 'x',[]);% Transmission Line structureLines = struct('con',[],'n',0,'bus1',[],'bus2',[]);if ~command_line_psat set(hdlp, ... 'YData',410-(idx-25)-[20 40 40 20 20], ... 'FaceColor',[0.50 1.00-idx/600 0.50+idx/600]) idx = idx + 30; pause(0.1)end% Synchronous Machine structureSyn = struct('con',[], ... 'n',0, ... 'bus',[], ... 'Id',[], ... 'Iq',[], ... 'Pg',[], ... 'Qg',[], ... 'J11',[], ... 'J12',[], ... 'J21',[], ... 'J22',[], ... 'delta',[], ... 'omega',[], ... 'e1q',[], ... 'e1d',[], ... 'e2q',[], ... 'e2d',[], ... 'psiq',[], ... 'psid',[], ... 'pm',[], ... 'vf',[], ... 'Gp',[], ... 'Gq',[], ... 'c1',[], ... 'c2',[], ... 'c3',[]);% Excitation Control structureExc = struct('con',[], ... 'n',0, ... 'bus',[], ... 'syn',[], ... 'vrif',[], ... 'vrif0',[], ... 'vr1',[], ... 'vr2',[], ... 'vr3',[], ... 'vm',[], ... 'vf',[]);% Subtransmission Area Equivalents structureSAE1 = struct('con',[], ... 'n',0, ... 'm',[]);SAE2 = struct('con',[], ... 'n',0, ... 'm1',[], ... 'm2',[]);SAE3 = struct('con',[], ... 'n',0, ... 'm1',[], ... 'm2',[], ... 'm3',[]);% Fault structureFault = struct('con',[], ... 'n',0, ... 'bus',[], ... 'dat',[], ... 'V',[], ... 'ang',[], ... 'delta',0);Breaker = struct('con',[], ... 'n',0, ... 'bus',[]);% Tap changer & load structureTap = struct('con',[], ... 'n',0, ... 'bus',[], ... 'm',[]);% Asynchronous Motor structureMot = struct('con',[], ... 'n',0, ... 'bus',[], ... 'dat',[], ... 'slip',[], ... 'e1r',[], ... 'e1m',[], ... 'e2r',[], ... 'e2m',[]);% Thermostatically contrlled loadThload = struct('con',[], ... 'n',0, ... 'bus',[], ... 'T',[], ... 'G',[]);% Turbine Governor structureTg = struct('con',[], ... 'n',0, ... 'syn',[], ... 'bus',[], ... 'dat1',[], ... 'dat2',[], ... 'tg',[], ... 'tg1',[], ... 'tg2',[], ... 'tg3',[]);% On Load Tap Changer structureLtc = struct('con',[], ... 'n',0, ... 'bus1',[], ... 'bus2',[], ... 'dat',[], ... 'm',[]);if ~command_line_psat set(hdlp, ... 'YData',410-(idx-25)-[20 40 40 20 20], ... 'FaceColor',[0.50 1.00-idx/600 0.50+idx/600]) idx = idx + 30; pause(0.1)end% Over Excitation Limiter structureOxl = struct('con',[], ... 'n',0, ... 'bus',[], ... 'syn',[], ... 'exc',[], ... 'v',[], ... 'If',[]);% Power system stabilizer structurePss = struct('con',[], ... 'n',0, ... 'bus',[], ... 'syn',[], ... 'exc',[], ... 'v1',[], ... 'v2',[], ... 'v3',[], ... 'va',[], ... 'vss',[], ... 'Vs',[], ... 's1',[]);% Static Var Compensator structureSvc = struct('con',[], ... 'n',0, ... 'bus',[], ... 'bcv',[], ... 'alpha',[], ... 'vm',[], ... 'Be',[], ... 'Vref',[]);% HVDC structureHvdc = struct('con',[], ... 'n',0, ... 'bus1', [], ... 'bus2',[], ... 'dat',[], ... 'Id',[], ... 'xr',[], ... 'xi',[]);% TCSC structureTcsc = struct('con',[], ... 'n',0, ... 'bus1',[], ... 'bus2',[], ... 'Pe',[], ... 'B',[], ... 'SI',[], ... 'x1',[], ... 'x2',[], ... 'x3',[]);% StatCom structureStatcom = struct('con',[], ... 'n',0, ... 'bus',[], ... 'dat',[], ... 'Id',[], ... 'Iq',[], ... 'Vdc',[], ... 'alpha',[], ... 'm',[], ... 'Vmdc',[], ... 'Vmac',[]);% SSSC structureSssc = struct('con',[], ... 'n',0, ... 'bus1',[], ... 'bus2',[], ... 'dat',[], ... 'Vdc',[], ... 'beta',[], ... 'm',[], ... 'Pmac',[], ... 'Vmdc',[]);% UPFC structureUpfc = struct('con',[], ... 'n',0, ... 'bus1',[], ... 'bus2',[], ... 'dat',[], ... 'Vmdc',[], ... 'Vmac',[], ... 'Vdc',[], ... 'alpha',[], ... 'm',[], ... 'x1',[], ... 'x2',[], ... 'Imd',[], ... 'Imq',[]);if ~command_line_psat set(hdlp, ... 'YData',410-(idx-25)-[20 40 40 20 20], ... 'FaceColor',[0.50 1.00-idx/600 0.50+idx/600]) idx = idx + 30; pause(0.1)end% RLC structureRLC = struct('con',[], ... 'dat',[], ... 'bus1',[], ... 'bus2',[], ... 'n',0, ... 'vd',[], ... 'vq',[], ... 'id',[], ... 'iq',[]);% Mass structureMass = struct('con',[], ... 'syn',[], ... 'n',0, ... 'delta_HP',[], ... 'omega_HP',[], ... 'delta_IP',[], ... 'omega_IP',[], ... 'delta_LP',[], ... 'omega_LP',[], ... 'delta_EX',[], ... 'omega_EX',[]);% SSR structureSSR = struct('con',[], ... 'bus',[], ... 'n',0, ... 'Id',[], ... 'Iq',[], ... 'If',[], ... 'Edc',[], ... 'Eqc',[], ... 'Tm',[], ... 'Efd',[], ... 'delta_HP',[], ... 'omega_HP',[], ... 'delta_IP',[], ... 'omega_IP',[], ... 'delta_LP',[], ... 'omega_LP',[], ... 'delta', [], ... 'omega', [], ... 'delta_EX',[], ... 'omega_EX',[]);% Voltage stability & Market structureDemand = struct('con',[], ... 'n',0, ... 'bus',[]);Supply = struct('con',[], ... 'n',0, ... 'bus',[]);Rsrv = struct('con',[], ... 'n',0, ... 'bus',[]);Rmpg = struct('con',[], ... 'n',0, ... 'sup',[], ... 'bus',[]);Rmpl = struct('con',[], ... 'n',0, ... 'dem',[], ... 'bus',[]);Vltn = struct('con',[], ... 'n',0, ... 'bus',[]);Ypdp = struct('con',[], ... 'day',[], ... 'week',[], ... 'year',[], ... 'd',1, ... 'w',1, ... 'y',1, ... 'len', 0);if ~command_line_psat set(hdlp, ... 'YData',410-(idx-25)-[20 40 40 20 20], ... 'FaceColor',[0.50 1.00-idx/600 0.50+idx/600]) idx = idx + 30; pause(0.1) hold offendif failed, disp(' '), disp('PSAT is not properly initialized.')else, if ~command_line_psat, fm_main, endendif ~command_line_psat, close(hdl), endclear failed hdl hdlp idx a count fid nname sline tipi ncompclear psatver psatdate psatdir lipdir command_line_psat% ----------------------------------------------------------- %% 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 - Fuel Cell ModelSofc = struct('con',[],'n',0,'bus',[],'Ik',[],'Vk',[],'pH2',[],'pH20',[],'pO2',[],'qH2',[],'m',[]);% Cac - Central Area ControllerCac = struct('con',[],'n',0,'bus',[],'q1',[],'q',[]);% Cluster - Cluster ControllersCluster = struct('con',[],'n',0,'exc',[],'syn',[],'svc',[],'bus',[],'cac',[],'Vs',[]);% Exload - Exponential recovery loadExload = struct('con',[],'n',0,'dat',[],'bus',[],'xp',[],'xq',[]);% Phs - Phase Shifting TransformerPhs = struct('con',[],'n',0,'bus1',[],'bus2',[],'alpha',[],'Pm',[]);% Wind speedWind = struct('con',[],'n',0,'speed',[],'vwa',[],'vw',[]);% Constant speed wind turbineCswt = struct('con',[],'n',0,'bus',[],'wind',[],'dat',[],'omega_wr',[], ... 'omega_m',[],'gamma',[],'e1r',[],'e1m',[]);% Doubly fed induction generatorDfig = struct('con',[],'n',0,'bus',[],'wind',[],'dat',[], ... 'omega_m',[],'theta_p',[],'idr',[],'iqr',[]);% Direct drive synchronous generatorDdsg = struct('con',[],'n',0,'bus',[],'wind',[],'dat',[], ... 'omega_m',[],'theta_p',[],'ids',[],'iqs',[],'idc',[]);% Busfreq - Bus frequency measurementBusfreq = struct('con',[],'n',0,'dat',[],'bus',[],'x',[],'w',[]);% Pmu - Phasor measurement unitPmu = struct('con',[],'n',0,'bus',[],'dat',[],'vm',[],'thetam',[]);% Jimma - Jimma's loadJimma = struct('con',[],'n',0,'dat',[],'bus',[],'x',[]);% Mixed - Mixed loadMixed = struct('con',[],'n',0,'dat',[],'bus',[],'x',[],'y',[]);⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -