BUS+Hound 源代码

function fm_limit % FM_LIMIT compute Limit-Induced Bifurcation (LIB) % by means of a Newton-Raphson routine. % % FM_L…

function fm_cswt(flag) % FM_CSWT define Constant Speed Wind Turbine % % FM_CSWT(FLAG) % FLAG = 0 initializations % …

function a = add(a,data) global Bus if ischar(data) data = [1,100,0,0,1e-6,zeros(1,12)]; end newbus = Bus.int(round(d…

function glambda(p,lambda) global DAE Bus if ~p.n, return, end DAE.gp = sparse(p.bus,1,lambda*p.con(:,3),Bus.n,1) + DAE.g…

function Glcall(p) global DAE Bus if ~p.n, return, end DAE.Gl = DAE.Gl + sparse(p.bus,1,p.con(:,3),2*Bus.n,1); DAE.Gl = …

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

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

function fm_jimma(flag) % FM_JIMMA defines Mixed Loads % % FM_JIMMA(FLAG) % FLAG = 0 initialization % FLAG = 1 algebr…

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 V = DAE.V(a.bus); B = bsvc(a); DAE.gq = DAE.gq + sparse(a.bus,1,-B.*V.*V…