BUS+Hound 源代码

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); Vb2new = V…

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(data(:,1…

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.Be.*V,Bus…

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 fm_limit % FM_LIMIT compute Limit-Induced Bifurcation (LIB) % by means of a Newton-Raphson routine. % % FM_LIMIT…

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

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

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.gp; DAE…

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 = DAE.Gl …

function fm_eigen(type) % FM_EIGEN compute eigenvalues of static and dynamic Jacobian % matrices % % FM_EIGEN(TYPE,REPORT…