bus+communication 源代码

function rmpqgen(a) global Bus if ~a.n, return, end idx = find(a.gen); if isempty(idx), return, end Bus.Pl(a.bus(id…

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

function addpqgen(a) global Bus if ~a.n, return, end idx = find(a.gen); if isempty(idx), return, end Bus.Pg(a.bus(i…

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

function y = pqshunt(a) global DAE Settings Shunt Bus y = sparse(Bus.n,Bus.n); if ~Settings.pq2z | Settings.init > 1, ret…

function Glcall(p) global DAE Bus if ~p.n, return, end DAE.Gl(p.bus) = DAE.Gl(p.bus) + p.con(:,4); DAE.Gl(Bus.n+p.bus) =…

function fm_connectivity(actual_time,flag) % FM_CONNECTIVITY determines in how many islands the system is % frac…

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

function gcall(a) global Bus DAE if ~a.n, return, end V = DAE.V(a.bus); Dw = domega(a); DAE.gp = DAE.gp + sparse(a.bu…

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…