setpm.m

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

function setpm(p)

global DAE Syn

if ~p.n, return, end

if p.ty1
  Syn.pm(p.syn(p.ty1)) = DAE.x(p.tg3) + ...
      p.dat1(:,12).*(DAE.x(p.tg2) + p.dat1(:,10).*DAE.x(p.tg1));
end

if p.ty2
  Syn.pm(p.syn(p.ty2)) = DAE.x(p.tg) - ...
      p.dat2(:,7).*(DAE.x(p.dat2(:,9)) - p.con(p.ty2,3)) + p.dat2(:,5);
  Syn.pm(p.syn(p.ty2)) = max(Syn.pm(p.syn(p.ty2)), p.dat2(:,4));
  Syn.pm(p.syn(p.ty2)) = min(Syn.pm(p.syn(p.ty2)), p.dat2(:,3));
end