BUS+Hound 源代码

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

function Gycall(a) global Bus DAE if ~a.n, return, end DAE.J22 = DAE.J22 - sparse(a.bus,a.bus,DAE.x(a.ist),Bus.n,Bus.n);

function gcall(a) global Bus DAE if ~a.n, return, end V = DAE.V(a.bus); B = bsvc(a); DAE.gq = DAE.gq + sparse(a.bus,1…

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

function a = add(a,data) global Bus a.n = a.n + length(data(1,:)); a.con = [a.con; data]; a.bus = [a.bus; Bus.int(round(da…

function Gycall(a) global Bus DAE if ~a.n, return, end V = DAE.V(a.bus); DAE.J22 = DAE.J22 + sparse(a.bus,a.bus,-2*a.B…

function q = getq(a,idx) global Bus if ~a.n, return, end B = bsvc(a); V = DAE.V(a.bus); q = B(idx).*V(idx).*V(idx);

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

function 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); a2 =…

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