qiantui2.m

前推潮流算法

dPkt(t)=(r(t)*(P(t)^2+Q(t)^2)/(abs(U(t))^2))/1000
dQkt(t)=(x(t)*(P(t)^2+Q(t)^2)/(abs(U(t))^2))/1000
Pkt(t)=P(t)+dPkt(t)
Qkt(t)=Q(t)+dQkt(t)