BUS+Hound 源代码

function Gycall(a) global Bus DAE Settings if ~a.n, return, end V = DAE.V(a.bus); P0 = a.con(:,2); Q0 = a.con(:,5); ap…

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 = zeros(Bus.n,1…

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

function a = add(a,data) global Bus newbus = Bus.int(round(data(1))); idx = findbus(a,newbus); if isempty(idx) a.n = …

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

function Fxcall(a) global DAE Bus Pod if ~a.n, return, end V1 = DAE.V(a.bus1); V2 = DAE.V(a.bus2); t1 = DAE.a(a.bus1); …

function gcall(a) global Bus DAE if ~a.n, return, end V1 = DAE.V(a.bus1); V2 = DAE.V(a.bus2); t1 = DAE.a(a.bus1); t2 =…

function p = base(p) global Line Bus Settings if ~p.n, return, end p.B = p.Cp./Line.con(p.line,9); p.y = (1-p.Cp)./Line.…

function a = gcall(a) global Bus DAE if ~a.n, return, end V1 = DAE.V(a.bus1); V2 = DAE.V(a.bus2); a1 = DAE.a(a.bus1); …

function p = base(p) global Line Bus if ~p.n, return, end fm_errv(p.con(:,4),'Upfc',p.bus1) Vb = Bus.con(p.bus1,2); p.c…