bus+communication 源代码

function a = base(a) if ~a.n, return, end global Bus Settings Vn = a.con(:,3); Vb = getkv(Bus,a.bus,1); fm_errv(Vn,'Fault',a.b…

function a = intervention(a,t) global DAE Bus Settings persistent angles voltages if ~a.n, return, end % do not repeat computat…

function a = setup(a) global Bus if isempty(a.con), return, end a.n = length(a.con(:,1)); [a.bus,a.vbus] = getbus(Bus,a.con(:,1…

function a = base(a) if ~a.n, return, end global Settings Bus fm_errv(a.con(:,3),'Induction Motor',a.bus); Vb2new = getkv(Bus,…

function p = base(p) global Bus Settings if ~p.n, return, end fm_errv(p.con(:,3),'PQ Bus',p.bus); Vb = getkv(Bus,p.bus,1); p.co…

function adjgen(a) global Bus DAE if ~a.n, return, end idx = find(a.gen); if isempty(idx), return, end p = a.P0(idx); q = a.Q…

function p = base(p) global Bus Settings if ~p.n, return, end V1 = getkv(Bus,p.bus1,1); V2 = getkv(Bus,p.bus2,1); VL1 = p.con(…

function [n,idx,data] = gams(a) global Bus n = int2str(a.n); idx = sparse(a.bus,[1:a.n],1,Bus.n,a.n); data = getzeros(Bus); if a…

function p = base(p) global Bus Settings if ~p.n, return, end fm_errv(p.con(:,3),'Slack Bus',p.bus); Vb = getkv(Bus,p.bus,1); p…

function [qmax,qmin] = swlim(a) global Bus if ~a.n qmax = []; qmin = []; return end qmax = find(Bus.Qg(a.bus) > a.con(:,6…