fcall.m

用于电力系统的一个很好的分析软件

function p = fcall(p)

global DAE

if ~p.n, return, end


if p.ty1
  tg1 = DAE.x(p.dat1(:,1));
  tg2 = DAE.x(p.dat1(:,2));
  tg3 = DAE.x(p.dat1(:,3));
  wref = DAE.y(p.wref(p.ty1));
  tin = p.dat1(:,6) + p.dat1(:,4).*(wref-DAE.x(p.dat1(:,14)));
  tin = max(tin,p.con(p.ty1,6));
  tin = min(tin,p.con(p.ty1,5));
  p.dat1(:,15) = tin;
  DAE.f(p.dat1(:,1)) = p.u(p.ty1).*p.dat1(:,7).*(-tg1 + tin);
  DAE.f(p.dat1(:,2)) = p.u(p.ty1).*p.dat1(:,8).*(-tg2 + p.dat1(:,11).*tg1);
  DAE.f(p.dat1(:,3)) = p.u(p.ty1).*p.dat1(:,9).*(-tg3 + p.dat1(:,13).*(tg2 + p.dat1(:,10).*tg1));
end

if p.ty2
  wref = DAE.y(p.wref(p.ty2));
  tg = DAE.x(p.dat2(:,1));
  DAE.f(p.dat2(:,1)) = -p.u(p.ty2).*p.dat2(:,6).*(tg + p.dat2(:,8).*(DAE.x(p.dat2(:,9)) - wref));
end