bus+communication 源代码

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

function Fxcall(a) global DAE Bus Settings if ~a.n, return, end Kpf = a.con(:,5); Kpv = a.con(:,6); Tpv = a.con(:,8); …

function gcall(a) global Bus DAE Settings if ~a.n, return, end x = DAE.x(a.x); y = DAE.x(a.y); V1 = DAE.V(a.bus); t1 =…

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 Settings if ~a.n, return, end V1 = DAE.V(a.bus); Kpf = a.con(:,5); Kpv = a.con(:,6); …

function fm_spf % FM_SPF solve standard power flow by means of the NR method % and fast decoupled power flow (XB and BX va…

function fm_thload(flag) % FM_THLOAD defines Thermostatically Controlled Loads % % FM_THLOAD(FLAG) % FLAG = 0 initializ…

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

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);

function fm_pareto % FM_PARETO determine a Pareto set of the multi-objective % OPF Problem % % FM_PARETO % %see al…