setx0.m

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

function a = setx0(a)

global Bus DAE PQ Varname

if ~a.n, return, end

V = DAE.V(a.bus);

for i = 1:a.n
  idx = findbus(PQ,a.bus(i));
  if isempty(idx)
    fm_disp(['No PQ load found for initializing frequency ', ...
             'dependent load at bus ',Varname.bus{a.bus(i)}])
  else
    P = PQ.P0(idx)*a.con(i,2)/100;
    Q = PQ.Q0(idx)*a.con(i,5)/100;
    PQ = pqsub(PQ,idx,P,Q);
    a.con(i,2) = a.con(i,2)*PQ.P0(idx)/(V(i)^a.con(i,3))/100;
    a.con(i,5) = a.con(i,5)*PQ.Q0(idx)/(V(i)^a.con(i,6))/100;
    PQ = remove(PQ,idx,'zero');
  end
end
DAE.x(a.x) = 0;
a.a0 = DAE.a(a.bus);
a.Dw = zeros(a.n,1);

%check limits
fm_disp('Initialization of Frequency Dependent Loads completed.')