fxcall.m

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

function Fxcall(p)

global DAE

if ~p.n, return, end

V1 = DAE.y(p.v1);
V2 = DAE.y(p.v2);
a1 = DAE.y(p.bus1);
a2 = DAE.y(p.bus2);
ss = sin(a1-a2);
cc = cos(a1-a2);
den = ssscden(p);
kp = p.Cp./(1-p.Cp);
vcs = DAE.x(p.vcs);
Tr = p.con(:,7);
vcs_max = p.con(:,8);
vcs_min = p.con(:,9);

Ms = p.y./den;

P1vs = Ms.*V1.*V2.*ss;  
P2vs = -P1vs;
Q1vs = Ms.*V1.*(V1-V2.*cc);
Q2vs = Ms.*V2.*(V2-V1.*cc);

u = (vcs <= vcs_max & vcs >= vcs_min & p.u);
ty3 = p.con(:,2) == 3;

DAE.Fx = DAE.Fx - sparse(p.vcs,p.vcs,1./Tr,DAE.n,DAE.n);
if ~isempty(find(ty3))
  DAE.Fx = DAE.Fx - sparse(p.vpi,p.vpi,~p.u,DAE.n,DAE.n);
end

DAE.Fy = DAE.Fy + sparse(p.vcs,p.v0,u./Tr,DAE.n,DAE.m);

DAE.Gx = DAE.Gx + sparse(p.bus1,p.vcs,u.*P1vs,DAE.m,DAE.n);  
DAE.Gx = DAE.Gx + sparse(p.bus2,p.vcs,u.*P2vs,DAE.m,DAE.n);
DAE.Gx = DAE.Gx + sparse(p.v1,p.vcs,u.*Q1vs,DAE.m,DAE.n);
DAE.Gx = DAE.Gx + sparse(p.v2,p.vcs,u.*Q2vs,DAE.m,DAE.n);    

a = find(p.u.*u.*ty3);
if a
  
  c = find(p.con(a,10) ~= 2);
  b = a(c);
  Kpr = p.con(a,11);
  Kin = p.con(a,12);
  c1 = p.y(a).*(1+vcs(a)./den(a));
  c2 = p.y(a).*vcs(a)./(den(a).^3);    
  Ms = p.y(a)./den(a);
  M2 = c2.*V1(a).*V2(a).*ss(a);
  M1 = -M2;
  M3 = c2.*(V2(a).*cc(a)-V1(a));
  M4 = c2.*(V1(a).*cc(a)-V2(a));              
  Jps1 = V1(a).*V2(a).*(M1.*ss(a)+c1.*cc(a));
  Jps2 = V1(a).*V2(a).*(M2.*ss(a)-c1.*cc(a));
  Jps3 = V2(a).*ss(a).*(M3.*V1(a)+c1);
  Jps4 = V1(a).*ss(a).*(M4.*V2(a)+c1);
  Jps5 = Ms.*V1(a).*V2(a).*ss(a);

  DAE.Gx = DAE.Gx - sparse(p.v0(b),p.vcs(b),Jps5(c).*Kpr(b),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(p.v0(a),p.vpi(a),1,DAE.m,DAE.n);

  DAE.Fx = DAE.Fx - sparse(p.vpi(b),p.vcs(b),Jps5(c).*Kin(b),DAE.n,DAE.n);
  
  DAE.Fy = DAE.Fy + sparse(p.vpi(a),p.pref(a),Kin(a),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(p.vpi(b),p.bus1(b),Jps1(c).*Kin(b),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(p.vpi(b),p.bus2(b),Jps2(c).*Kin(b),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(p.vpi(b),p.v1(b),Jps3(c).*Kin(b),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(p.vpi(b),p.v2(b),Jps4(c).*Kin(b),DAE.n,DAE.m);
  
end