setx0.m

这是一个很适合研究和学习用的电力系统仿真软件

function a = setx0(a)

global DAE Line

if ~a.n, return, end

ty3 = a.con(:,2) == 3;
typwr = a.con(:,10) == 1;
tyang = a.con(:,10) == 2;

kp = a.Cp./(1-a.Cp);

DAE.x(a.vcs) = kp.*ssscden(a);
DAE.x(a.vpi) = DAE.x(a.vcs(find(ty3)));

% reset transmission line reactance and admittance matrix
Line.con(a.line,9) = Line.con(a.line,9) + a.xcs;
fm_y;

[Ps0,Qs0,Pr0,Qr0,Is0,Ir0] = fex_lineflows;
a.Pref = Ps0(a.line) + ty3.*tyang.*Pr0(a.line);

vcs_max = a.con(:,8);
vcs_min = a.con(:,9);

idx = find(DAE.x(a.vcs) > vcs_max);
if idx, ssscwarn(a,idx,' Vs is over its max limit.'), end
idx = find(DAE.x(a.vcs) < vcs_min);
if idx, ssscwarn(a,idx,' Vs is under its min limit.'), end
DAE.x(a.vcs) = max(DAE.x(a.vcs),vcs_min);
DAE.x(a.vcs) = min(DAE.x(a.vcs),vcs_max);

% reference voltage signal
a.V0 = DAE.x(a.vcs);

fm_disp('Initialization of SSSC completed.')