fm_mn.m

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

function fm_mn(flag)
% FM_MN defines Monomial Loads
%
% FM_MN(FLAG)
%       FLAG = 0 initialization
%       FLAG = 1 algebraic equations
%       FLAG = 2 algebraic Jacobians
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2006 Federico Milano

global Mn PQ DAE Bus Settings Varname

V = DAE.V(Mn.bus);

switch flag
 case 0
  i = find(Mn.con(:,8));
  for j = 1:length(i)
    k = i(j);
    idx = findbus(PQ,Mn.bus(k));
    if isempty(idx),
      fm_disp(['No PQ load found for initializing monomial ', ...
               'load at bus ',Varname.bus{Mn.bus(k)}])
    else
      P = PQ.P0(idx)*Mn.con(k,4)/100;
      Q = PQ.Q0(idx)*Mn.con(k,5)/100;
      PQ = pqsub(PQ,idx,P,Q);
      Mn.con(k,4)  = Mn.con(k,4)*PQ.P0(idx)/(V(k)^Mn.con(k,6))/100;
      Mn.con(k,5)  = Mn.con(k,5)*PQ.Q0(idx)/(V(k)^Mn.con(k,7))/100;
      PQ = remove(PQ,idx,'zero');
    end
  end
  fm_disp('Initialization of monomial loads completed.')
 case 1
  if Settings.init
    DAE.gp = DAE.gp + sparse(Mn.bus,1,Mn.con(:,4).*V.^Mn.con(:,6),Bus.n,1);
    DAE.gq = DAE.gq + sparse(Mn.bus,1,Mn.con(:,5).*V.^Mn.con(:,7),Bus.n,1);
  elseif Mn.init
    i = Mn.init;
    DAE.gp = DAE.gp + sparse(Mn.bus(i),1,Mn.con(i,4).*V(i).^Mn.con(i,6),Bus.n,1);
    DAE.gq = DAE.gq + sparse(Mn.bus(i),1,Mn.con(i,5).*V(i).^Mn.con(i,7),Bus.n,1);
  end
 case 2
  if Settings.init
    DAE.J12 = DAE.J12 + sparse(Mn.bus,Mn.bus,Mn.con(:,4).* ...
			       Mn.con(:,6).*V.^(Mn.con(:,6)-1),Bus.n,Bus.n);
    DAE.J22 = DAE.J22 + sparse(Mn.bus,Mn.bus,Mn.con(:,5).* ...
			       Mn.con(:,7).*V.^(Mn.con(:,7)-1),Bus.n,Bus.n);
  elseif Mn.init
    i = Mn.init;
    DAE.J12 = DAE.J12 + sparse(Mn.bus(i),Mn.bus(i),Mn.con(i,4).* ...
			       Mn.con(i,6).*V(i).^(Mn.con(i,6)-1),Bus.n,Bus.n);
    DAE.J22 = DAE.J22 + sparse(Mn.bus(i),Mn.bus(i),Mn.con(i,5).* ...
			       Mn.con(i,7).*V(i).^(Mn.con(i,7)-1),Bus.n,Bus.n);
  end
end