diff --git a/run.m b/run.m index c510371..ec990da 100644 --- a/run.m +++ b/run.m @@ -59,7 +59,7 @@ while(k<=kmax && maxD> EPS) jaco=Jacobi(Balance,busNum,QGi,Vmf1,YdotSinVolt,YdotCosVolt,diag_Volt_YdotSin,diag_Volt_YdotCos);%雅克比矩阵 [dV, dVangle]=Solv(busNum,jaco,dP,dQ);%解出修正量 [Vmf1, Vaf1]=Modify(Vmf1,Vaf1,dV,dVangle,1); - fprintf('第 %d 次迭代, 最大不平衡量为 %f\n\n',k,full(maxD)); + fprintf('第 %d 次迭代, 最大不平衡量为 %f\n',k,full(maxD)); %转换为三相电压 VoltpABC=Tf2p*conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']); VoltpA=conj(VoltpABC(1,:)');