setx0.m

来自「用于电力系统的一个很好的分析软件」· M 代码 · 共 16 行

function a = setx0(a)

global DAE

if ~a.n, return, end

% variable initialization
DAE.x(a.q1) = 1;
DAE.y(a.q) = 1;

% pilot bus voltage reference
a.con(:,5) = DAE.y(a.vbus);

fm_disp('Initialization of Central Area Controllers completed.')