fxcall.m

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

  DAE.Fx = DAE.Fx - sparse(e2d52,a.delta(is52),b522.*M2(is52),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2d52,e2q52,b522.*a.c1(is52),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2d52,e2d52,b521+b522.*a.c2(is52)+(~a.u(is52)),DAE.n,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.p(is52),e2q52,Gp1(is52),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.p(is52),e2d52,Gp2(is52),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.q(is52),e2q52,Gq1(is52),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.q(is52),e2d52,Gq2(is52),DAE.m,DAE.n);
  DAE.Fy = DAE.Fy + sparse(e1q52,a.vf(is52),a526,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2q52,a.vf(is52),a523,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q52,bs52,a525.*M1(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q52,vs52,a525.*M3(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q52,bs52,a522.*M1(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q52,vs52,a522.*M3(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2d52,bs52,b522.*M2(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2d52,vs52,b522.*M4(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q52,a.p(is52),a526.*Kp(is52),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q52,a.p(is52),a523.*Kp(is52),DAE.n,DAE.m);
end

% Model V Type 3
if (~isempty(is53))
  Wn = rad*a.u(is53);
  d53 = a.delta(is53);
  o53 = a.omega(is53);
  e1q53 = a.e1q(is53);
  psid53 = a.psid(is53);
  psiq53 = a.psiq(is53);
  e1q(is53) = DAE.x(e1q53);
  psid(is53) = DAE.x(psid53);
  psiq(is53) = DAE.x(psiq53);
  a531 = (xd(is53)-xd1(is53))./(xd(is53)+xl(is53));
  a532 = a.u(is53)./(1-a531);
  a534 = a.u(is53)./Td10(is53);
  c534 = a.u(is53)./(xd(is53)+xl(is53));
  c535 = a.u(is53)./(xq(is53)+xl(is53));
  q1 = -2*ra(is53).*(e1q(is53)-psid(is53)).*c534.*c534.*iM(is53);
  q2 = -2*ra(is53).*psiq(is53).*c535.*c535.*iM(is53);
  q3 = -(a534+a531.*Wn.*ra(is53).*c534).*a532;
  q4 = a531.*Wn.*ra(is53).*c534.*a532;
  DAE.Fx = DAE.Fx + sparse(o53,e1q53,q1,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(o53,psiq53,q2,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(o53,psid53,q1,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q53,d53,a531.*Wn.*vg(is53).*cc(is53).*a532,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q53,o53,a531.*Wn.*psiq(is53).*a532,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1q53,o53,a532.*a534.*Kw(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1q53,e1q53,q3-(~a.u(is53)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q53,psiq53,a531.*Wn.*omega(is53).*a532,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1q53,psid53,q4,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq53,d53,Wn.*vg(is53).*ss(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq53,o53,Wn.*psid(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq53,psiq53,Wn.*ra(is53).*c535+(~a.u(is53)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq53,psid53,Wn.*omega(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid53,d53,Wn.*vg(is53).*cc(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid53,o53,Wn.*psiq(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid53,e1q53,Wn.*ra(is53).*c534,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid53,psiq53,Wn.*omega(is53),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psid53,psid53,Wn.*ra(is53).*c534+(~a.u(is53)),DAE.n,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.p(is53),e1q53,vg(is53).*c534.*ss(is53),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(a.p(is53),psiq53,vg(is53).*c535.*cc(is53),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(a.p(is53),psid53,vg(is53).*c534.*ss(is53),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.q(is53),e1q53,vg(is53).*c534.*cc(is53),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(a.q(is53),psiq53,vg(is53).*c535.*ss(is53),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(a.q(is53),psid53,vg(is53).*c534.*cc(is53),DAE.m,DAE.n);
  DAE.Fy = DAE.Fy + sparse(e1q53,a.vf(is53),a532.*a534,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q53,a.p(is53),a532.*a534.*Kp(is53),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e1q53,bs53,a531.*Wn.*vg(is53).*cc(is53).*a532,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q53,vs53,a531.*Wn.*ss(is53).*a532,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(psiq53,bs53,Wn.*vg(is53).*ss(is53),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(psiq53,vs53,Wn.*cc(is53),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(psid53,bs53,Wn.*vg(is53).*cc(is53),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(psid53,vs53,Wn.*ss(is53),DAE.n,DAE.m);
end

% Model VI
if(~isempty(is6))
  o6 = a.omega(is6);
  e1q6 = a.e1q(is6);
  e1d6 = a.e1d(is6);
  e2q6 = a.e2q(is6);
  e2d6 = a.e2d(is6);
  e1d(is6) = DAE.x(e1d6);
  e1q(is6) = DAE.x(e1q6);
  e2d(is6) = DAE.x(e2d6);
  e2q(is6) = DAE.x(e2q6);
  Taa = a.con(:,24);
  gd = xd2(is6)./xd1(is6).*Td20(is6)./Td10(is6).*(xd(is6)-xd1(is6));
  gq = xq2(is6)./xq1(is6).*Tq20(is6)./Tq10(is6).*(xq(is6)-xq1(is6));
  a1 = a.u(is6)./Td20(is6);
  a2 = a1.*(xd1(is6)-xd2(is6)+gd);
  a3 = a.u(is6).*Taa(is6)./Td10(is6)./Td20(is6);
  a4 = a.u(is6)./Td10(is6);
  a5 = a4.*(xd(is6)-xd1(is6)-gd);
  a6 = a4.*(1-Taa(is6)./Td10(is6));
  b1 = a.u(is6)./Tq20(is6);
  b2 = b1.*(xq1(is6)-xq2(is6)+gq);
  b3 = a.u(is6)./Tq10(is6);
  b4 = b3.*(xq(is6)-xq1(is6)-gq);
  p6 = a.p(is6);
  q6 = a.q(is6);
  DAE.Fx = DAE.Fx + sparse(o6,e2q6,N1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(o6,e2d6,N2(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1q6,a.delta(is6),a5.*M1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1q6,o6,a6.*Kw(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q6,e1q6,a4.*synsat(a,2,e1q(is6),is6)+(~a.u(is6)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q6,e2q6,a5.*a.c3(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q6,e2d6,a5.*a.c1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1d6,a.delta(is6),b4.*M2(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1d6,e1d6,b3+(~a.u(is6)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1d6,e2q6,b4.*a.c1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1d6,e2d6,b4.*a.c2(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2q6,a.delta(is6),a2.*M1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2q6,o6,a3.*Kw(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2q6,e1q6,a1,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2q6,e2q6,a1+a2.*a.c3(is6)+(~a.u(is6)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2q6,e2d6,a2.*a.c1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2d6,a.delta(is6),b2.*M2(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2d6,e1d6,b1,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2d6,e2q6,b2.*a.c1(is6),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2d6,e2d6,b1+b2.*a.c2(is6)+(~a.u(is6)),DAE.n,DAE.n);
  DAE.Gx = DAE.Gx + sparse(p6,e2q6,Gp1(is6),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(p6,e2d6,Gp2(is6),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(q6,e2q6,Gq1(is6),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(q6,e2d6,Gq2(is6),DAE.m,DAE.n);
  DAE.Fy = DAE.Fy + sparse(e1q6,a.vf(is6),a6,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2q6,a.vf(is6),a3,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q6,bs6,a5.*M1(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q6,vs6,a5.*M3(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e1d6,bs6,b4.*M2(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e1d6,vs6,b4.*M4(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q6,bs6,a2.*M1(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q6,vs6,a2.*M3(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2d6,bs6,b2.*M2(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2d6,vs6,b2.*M4(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q6,p6,a6.*Kp(is6),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q6,p6,a3.*Kp(is6),DAE.n,DAE.m);
end

% Model VIII
if ~isempty(is8)
  Wn = rad*a.u(is8);
  d8 = a.delta(is8);
  o8 = a.omega(is8);
  e1q8 = a.e1q(is8);
  e1d8 = a.e1d(is8);
  e2q8 = a.e2q(is8);
  e2d8 = a.e2d(is8);
  psid8 = a.psid(is8);
  psiq8 = a.psiq(is8);
  e1d(is8) = DAE.x(e1d8);
  e1q(is8) = DAE.x(e1q8);
  e2d(is8) = DAE.x(e2d8);
  e2q(is8) = DAE.x(e2q8);
  psid(is8) = DAE.x(psid8);
  psiq(is8) = DAE.x(psiq8);
  Taa = a.con(:,24);
  gd = xd2(is8)./xd1(is8).*Td20(is8)./Td10(is8).*(xd(is8)-xd1(is8));
  gq = xq2(is8)./xq1(is8).*Tq20(is8)./Tq10(is8).*(xq(is8)-xq1(is8));
  a18 = a.u(is8)./Td20(is8);
  a28 = a18.*(xd1(is8)-xd2(is8)+gd);
  a38 = a.u(is8).*Taa(is8)./Td10(is8)./Td20(is8);
  a48 = a.u(is8)./Td10(is8);
  a58 = a48.*(xd(is8)-xd1(is8)-gd);
  a68 = a48.*(1-Taa(is8)./Td10(is8));
  b18 = a.u(is8)./Tq20(is8);
  b28 = b18.*(xq1(is8)-xq2(is8)+gq);
  b38 = a.u(is8)./Tq10(is8);
  b48 = b38.*(xq(is8)-xq1(is8)-gq);
  xd2(is8) = xd2(is8)+xl(is8);
  xq2(is8) = xq2(is8)+xl(is8);
  p8 = a.p(is8);
  q8 = a.q(is8);
  q1 = 2*ra(is8).*(e2d(is8)+psiq(is8))./xq2(is8).*iM(is8)./xq2(is8);
  q2 = 2*ra(is8).*(e2q(is8)-psid(is8))./xd2(is8).*iM(is8)./xd2(is8);
  DAE.Fx = DAE.Fx - sparse(o8,e2d8,q1,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(o8,e2q8,q2,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(o8,psiq8,q1,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(o8,psid8,q2,DAE.n,DAE.n);
  
  DAE.Fx = DAE.Fx + sparse(e1q8,o8,a68.*Kw(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q8,e1q8,a48.*synsat(a,2,e1q(is8),is8)+(~a.u(is8)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1q8,e2q8,a58./xd2(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e1q8,psid8,a58./xd2(is8),DAE.n,DAE.n);
  
  DAE.Fx = DAE.Fx - sparse(e1d8,e1d8,b38+(~a.u(is8)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1d8,e2d8,b48./xq2(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e1d8,psiq8,b48./xq2(is8),DAE.n,DAE.n);
  
  DAE.Fx = DAE.Fx + sparse(e2q8,o8,a38.*Kw(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2q8,e1q8,a18,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2q8,e2q8,a18+a28./xd2(is8)+(~a.u(is8)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(e2q8,psid8,a28./xd2(is8),DAE.n,DAE.n);
  
  DAE.Fx = DAE.Fx + sparse(e2d8,e1d8,b18,DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2d8,e2d8,b18+b28./xq2(is8)+(~a.u(is8)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(e2d8,psiq8,b28./xq2(is8),DAE.n,DAE.n);
  
  DAE.Fx = DAE.Fx - sparse(psiq8,a.delta(is8),Wn.*vg(is8).*ss(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq8,o8,Wn.*psid(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq8,e2d8,Wn.*ra(is8)./xq2(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq8,psiq8,Wn.*ra(is8)./xq2(is8)+(~a.u(is8)),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psiq8,psid8,Wn.*omega(is8),DAE.n,DAE.n);
  
  DAE.Fx = DAE.Fx + sparse(psid8,a.delta(is8),Wn.*vg(is8).*cc(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid8,o8,Wn.*psiq(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid8,e2q8,Wn.*ra(is8)./xd2(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx + sparse(psid8,psiq8,Wn.*omega(is8),DAE.n,DAE.n);
  DAE.Fx = DAE.Fx - sparse(psid8,psid8,Wn.*ra(is8)./xd2(is8)+(~a.u(is8)),DAE.n,DAE.n);
  
  DAE.Gx = DAE.Gx + sparse(p8,e2q8,vg(is8)./xd2(is8).*ss(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(p8,e2d8,vg(is8)./xq2(is8).*cc(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(p8,psiq8,vg(is8)./xq2(is8).*cc(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(p8,psid8,vg(is8)./xd2(is8).*ss(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(q8,e2q8,vg(is8)./xd2(is8).*cc(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(q8,e2d8,vg(is8)./xq2(is8).*ss(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx + sparse(q8,psiq8,vg(is8)./xq2(is8).*ss(is8),DAE.m,DAE.n);
  DAE.Gx = DAE.Gx - sparse(q8,psid8,vg(is8)./xd2(is8).*cc(is8),DAE.m,DAE.n);
  
  DAE.Fy = DAE.Fy + sparse(e1q8,a.vf(is8),a68,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(e2q8,a.vf(is8),a38,DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(psiq8,bs8,Wn.*vg(is8).*ss(is8),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(psiq8,vs8,Wn.*cc(is8),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(psid8,bs8,Wn.*vg(is8).*cc(is8),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy + sparse(psid8,vs8,Wn.*ss(is8),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e1q8,p8,a68.*Kp(is8),DAE.n,DAE.m);
  DAE.Fy = DAE.Fy - sparse(e2q8,p8,a38.*Kp(is8),DAE.n,DAE.m);
end