bdmatrix.m

电力系统分析计算程序

function bdmatrix(a)

global LA DAE 

LA.b_tcsc = [];  
LA.d_tcsc = [];  

if ~a.n
  fm_disp('* * * No TCSC device found')
  return
end

ix = find(a.con(:,3) == 1); % index - operation mode x constant
ip = find(a.con(:,3) == 2); % index - operation mode P constant

Fu = sparse(DAE.n,a.n);
Gu = sparse(a.x0(ix),ix,a.u(ix),DAE.m,a.n) + ...
    sparse(a.pref(ip),ip,a.u(ip),DAE.m,a.n);

% B & D matrix for TCSC references
LA.d_tcsc = -full(DAE.Gy\Gu);
LA.b_tcsc = full(Fu + DAE.Fy*LA.d_tcsc);