BUS+Hound 源代码

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

function p = base(p) global Line Bus Settings if ~p.n, return, end p.B = p.Cp./Line.con(p.line,9); p.y = (1-p.Cp)./Line.con(p.l…

clear all, clc, close all % PF data in MATLAB format % IEEE 14 Bus Test System Sb = 100; econv = 0.001; mxitr = 10; % …

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

function fm_tcsc(flag) % FM_TCSC defines Synchronous Machines % % FM_TCSC(FLAG) % FLAG = 0 initializations % FLAG = 1…

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 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 fm_eigen(type) % FM_EIGEN compute eigenvalues of static and dynamic Jacobian % matrices % % FM_EIGEN(TYPE,REPORT…