BUS+Hound 源代码

function y = pqshunt(a,idx) global DAE Shunt Bus g = a.con(idx,4)./(DAE.V(a.bus(idx)).^2); b = -a.con(idx,5)./(DAE.V(a.bus(idx))…

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

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

function fm_fl(flag) % FM_FL define General Exponential Voltage Frequency Dependent Loads % %FM_FL(FLAG) % FLAG = 0 initiali…

function [n,idx,data] = gams(a) global Bus n = int2str(a.n); idx = sparse(a.bus,[1:a.n],1,Bus.n,a.n); data = zeros(Bus.n,1); if …

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

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

function [qmax,qmin] = swlim(a,q) global Bus if ~a.n qmax = []; qmin = []; return end qmax = find(Bus.Qg(a.bus) > a.con(:,…

function fm_mixed(flag) % FM_MIXED defines Mixed Loads % % FM_MIXED(FLAG) % FLAG = 0 initialization % FLAG = 1 algebr…

function Fxcall(a) global DAE Bus Pod if ~a.n, return, end V1 = DAE.V(a.bus1); V2 = DAE.V(a.bus2); t1 = DAE.a(a.bus1); t2 = DAE…