BUS+Hound 源代码

function fm_ncomp % FM_NCOMP search components used in the current data % file and initializes fields used for power %…

function Fxcall(a) global DAE Bus if ~a.n, return, end iTf = 1./a.con(:,5); Kv = a.con(:,12); DAE.Fx = DAE.Fx + spars…

function gcall(a) global Bus DAE if ~a.n, return, end x = DAE.x(a.x); V1 = DAE.V(a.bus); Tf = a.con(:,5); Plz = a.con(…

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 V1 = DAE.V(a.bus); Tf = a.con(:,5); Plz = a.con(:,6); Pli = a.c…

function fm_pl(flag) % FM_PL define Polinomial Loads % % FM_PL(FLAG) % FLAG -> 0 initialization % FLAG -> 1 algeb…

function fex_nrloop(tol,itermax) % FEX_NRLOOP -> called by FEX_ABCD % Perform Newton Raphson for non linear algebraic equation …

function gcall(a) global Bus DAE if ~a.n, return, end DAE.gq = DAE.gq - sparse(a.bus,1,DAE.x(a.ist).* DAE.V(a.bus),Bus.n,…

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

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…