xfirst.m

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

function xfirst(a)

if ~a.n, return, end

global DAE Settings

I0 = a.con(:,26);
P0 = a.con(:,27);
V0 = a.con(:,28);
uI = a.u.*a.dat(:,9);
uP = a.u.*a.dat(:,10);
uV = a.u.*a.dat(:,11);

DAE.y(a.cosa) = 0.5*a.u.*(a.dat(:,3)+a.dat(:,4));
DAE.y(a.cosg) = a.u.*a.dat(:,6); % minimum gamma
DAE.y(a.Vrdc) = uI + uP + uV.*V0;
DAE.y(a.Vidc) = uI + uP + uV.*V0;
DAE.y(a.yr) = uI.*I0 + uP.*P0 + uV.*V0;
DAE.y(a.yi) = uI.*I0 + uP.*P0 + uV.*V0;
DAE.y(a.Sr) = 15*a.u.*a.con(:,3)/Settings.mva;
DAE.y(a.Si) = 15*a.u.*a.con(:,3)/Settings.mva;

DAE.x(a.Idc) = uI.*I0 + uP.*P0;
DAE.x(a.xr) = DAE.y(a.cosa);
DAE.x(a.xi) = DAE.y(a.cosg);