From ecaed830b547d025369bd6137707f6a96856556a Mon Sep 17 00:00:00 2001 From: facat Date: Tue, 22 Apr 2014 09:29:28 +0800 Subject: [PATCH] =?UTF-8?q?=E4=BF=AE=E6=AD=A3=E5=8E=9F=E6=9D=A5=E5=A4=9A?= =?UTF-8?q?=E6=8D=A2=E4=BA=86=E4=B8=80=E8=A1=8C=E7=9A=84=E9=94=99=E8=AF=AF?= =?UTF-8?q?=E3=80=82?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: facat --- run.m | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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,:)');