fxcall.m

电力系统分析计算程序

function Fxcall(p)

global DAE 

if ~p.n, return, end

if p.ty1
  pm1 = p.pm(p.ty1);
  tg1 = p.dat1(:,1);
  tg2 = p.dat1(:,2);
  tg3 = p.dat1(:,3);
  u1 = p.u(p.ty1);
  tin = p.dat1(:,15);
  u = tin < p.con(p.ty1,5) & tin > p.con(p.ty1,6); % windup limiter
  DAE.Fx = DAE.Fx - sparse(tg1,p.dat1(:,14),u.*u1.*p.dat1(:,7).*p.dat1(:,4),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(tg1,tg1,p.dat1(:,7),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(tg2,tg2,p.dat1(:,8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(tg3,tg3,p.dat1(:,9),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(tg2,tg1,u1.*p.dat1(:,8).*p.dat1(:,11),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(tg3,tg2,u1.*p.dat1(:,9).*p.dat1(:,13),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(tg3,tg1,u1.*p.dat1(:,10).*p.dat1(:,13).*p.dat1(:,9),DAE.n,DAE.n);
  DAE.Fy = DAE.Fy + sparse(tg1,p.wref(p.ty1),u.*u1.*p.dat1(:,7).*p.dat1(:,4),DAE.n,DAE.m);
  DAE.Gx = DAE.Gx - sparse(pm1,tg3,u1,DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(pm1,tg2,u1.*p.dat1(:,12),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(pm1,tg1,u1.*p.dat1(:,10).*p.dat1(:,12),DAE.m,DAE.n);
end

if p.ty2
  tg = p.dat2(:,1);
  pm2 = p.pm(p.ty2);
  tin = DAE.y(p.pm(p.ty2));
  u = tin < p.dat2(:,3) & tin > p.dat2(:,4) & p.u(p.ty2); % windup limiter
  DAE.Fx = DAE.Fx - sparse(tg,tg,p.dat2(:,6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(tg,p.dat2(:,9),p.u(p.ty2).*p.dat2(:,6).*p.dat2(:,8),DAE.n,DAE.n);
  DAE.Fy = DAE.Fy + sparse(tg,p.wref(p.ty2),p.u(p.ty2).*p.dat2(:,6).*p.dat2(:,8),DAE.n,DAE.m);
  DAE.Gx = DAE.Gx - sparse(pm2,tg,u,DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(pm2,p.dat2(:,9),u.*p.dat2(:,7),DAE.m,DAE.n);
end