qiantui1.m

前推潮流算法

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