fxcall.m

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

function Fxcall(p)

global DAE

if ~p.n, return, end

if p.ty1
  tin = p.dat1(:,16);
  u = tin < p.con(p.ty1,5) & tin > p.con(p.ty1,6); % windup limiter
  DAE.Fx = DAE.Fx - sparse(p.dat1(:,1),p.dat1(:,14),u.*p.dat1(:,7).*p.dat1(:,4),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat1(:,14),p.dat1(:,1),p.dat1(:,12).*p.dat1(:,15).*p.dat1(:,10),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat1(:,14),p.dat1(:,2),p.dat1(:,12).*p.dat1(:,15),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat1(:,14),p.dat1(:,3),p.dat1(:,15),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(p.dat1(:,1),p.dat1(:,1),p.dat1(:,7),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(p.dat1(:,2),p.dat1(:,2),p.dat1(:,8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(p.dat1(:,3),p.dat1(:,3),p.dat1(:,9),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat1(:,2),p.dat1(:,1),p.dat1(:,8).*p.dat1(:,11),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat1(:,3),p.dat1(:,2),p.dat1(:,9).*p.dat1(:,13),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat1(:,3),p.dat1(:,1),p.dat1(:,10).*p.dat1(:,13).*p.dat1(:,9),DAE.n,DAE.n);
end

if p.ty2
  tin = DAE.x(p.tg) - p.dat2(:,7).*(DAE.x(p.dat2(:,9)) - p.con(p.ty2,3)) + p.dat2(:,5);
  u = tin < p.dat2(:,3) & tin > p.dat2(:,4); % windup limiter
  DAE.Fx = DAE.Fx - sparse(p.dat2(:,9),p.dat2(:,9),u.*p.dat2(:,7).*p.dat2(:,10),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(p.dat2(:,9),p.dat2(:,1),u.*p.dat2(:,10),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(p.dat2(:,1),p.dat2(:,1),p.dat2(:,6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(p.dat2(:,1),p.dat2(:,9),p.dat2(:,6).*p.dat2(:,8),DAE.n,DAE.n);
end