qiantui3.m

前推潮流算法

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