diff --git a/run.m b/run.m
index 4396098..52e6db3 100644
--- a/run.m
+++ b/run.m
@@ -2,7 +2,7 @@ clc
 clear
 lineZ=readLineZ('.\feeder13\lineParameter.txt');
 [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
-    phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY]=dataRead(lineZ,'.\feeder13\data1.txt');
+    phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY]=dataRead(lineZ,'.\feeder13\data2.txt');
 a=exp(1j*2*pi/3);
 Tp2f=1/3*[1 1 1;
     1 a a^2;
@@ -12,16 +12,19 @@ Tf2p=inv(Tp2f);
 fsY1amp=abs(fsY1);
 [r,c,fsY1ang]=find(fsY1);
 fsY1ang=angle(fsY1ang);
+%%参考电压
+refV0=1;
+refV1=0;
+refV2=0;
+%refV1=1;
 Pabc=phaseASpotLoadP+phaseBSpotLoadP+phaseCSpotLoadP;
 Qabc=phaseASpotLoadQ+phaseBSpotLoadQ+phaseCSpotLoadQ;
 busNum=length(phaseASpotLoadP);
 %给序电压赋初值
 Vmf1=sparse(ones(busNum,1));
-Vmf2=sparse(ones(busNum,1));
-Vmf0=sparse(busNum,1);
-Vaf1=sparse(busNum,1);
-Vaf2=sparse(busNum,1);
-Vaf0=sparse(busNum,1);
+Vaf1=sparse(zeros(busNum,1));
+Vf1=Vmf1.*exp(1j*Vaf1);
+%Vaf1=sparse(zeros(busNum,1));
 %先求解正序的
 PQi=nodeNum;
 PG=sparse(busNum,1);
@@ -38,6 +41,7 @@ fsY00=fsY0;
 fsY22=fsY2;
 Vf2=sparse(busNum,1);
 If2=sparse(busNum,1);
+% Vf0=sparse(busNum,1);
 Vf0=sparse(busNum,1);
 If0=sparse(busNum,1);
 %准备序矩阵
@@ -49,14 +53,6 @@ fsY0(:,Balance)=0;
 fsY0=fsY0+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
 while(k<=kmax && maxD> EPS)
     k=k+1;
-    [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
-        PQi,PG,QG,QGi,PD,QD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
-    maxD=max(abs([dP;dQ;]));
-    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));
-    %%
     %转换为三相电压
     VoltpABC=Tf2p*conj([ Vf0'; (Vmf1.*exp(1j*Vaf1))'; Vf2']);
     VoltpA=conj(VoltpABC(1,:)');
@@ -73,16 +69,21 @@ while(k<=kmax && maxD> EPS)
     If2(Balance)=0;
     Vf0=fsY0\If0;
     Vf2=fsY2\If2;
+%     Vf2=Vf2+refV2;
+%     Vf0=Vf0+refV0;
+    [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ...
+        PQi,PG,QG,QGi,PD,QD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量
+    maxD=max(abs([dP;dQ;]));
+    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));
 end
 Vf1=Vmf1.*exp(1j*Vaf1);
 %%
-% V=[1.01*exp(1j*0);1.01*exp(-1j*pi/360*120);1.01*exp(-1j*pi/360*240)];
-% abs(Tp2f*V)
-% Tf2p*Tp2f*V
-% Vf0=Vf0+refV0;
-% Vf2=Vf2+refV2;
-% Vf0=Vf0*exp(-1j*(-1.0471975511966));
-% Vf1=Vf1*exp(-1j*(1.0471975511966));
+V=[1.0*exp(1j*0);1.0*exp(-1j*pi/180*120);1.0*exp(1j*pi/180*120)];
+abs(Tp2f*V);
+Tf2p*Tp2f*V;
 %%
 (Vf0.*conj(fsY00*Vf0)+Vf1.*conj(fsY11*Vf1)+Vf2.*conj(fsY22*Vf2))*3;
 conj(Tf2p*[If0(2);If1(2);If2(2)]).*(Tf2p*[Vf0(2);Vf1(2);Vf2(2)]);