diff --git a/Modification.m b/Modification.m
index cbe7d32..977785d 100644
--- a/Modification.m
+++ b/Modification.m
@@ -17,6 +17,8 @@ t=deltX(size(PVi,1)+size(PGi,1)+1:ContrlCount)';
 t(Busnum+Balance)=0;
 %Volt=Volt+AlphaP*t(2:2:2*Busnum);暂时改一下20111227
 %UAngel=UAngel+AlphaP*t(1:2:2*Busnum);暂时改一下20111227
+balVolt=Volt(Balance);
 Volt=Volt+AlphaP*t(1:Busnum);
+Volt(Balance)=balVolt;
 UAngel=UAngel+AlphaP*t(Busnum+1:2*Busnum);
 end
\ No newline at end of file
diff --git a/OPF.m b/OPF.m
index 340ce87..902bf09 100644
--- a/OPF.m
+++ b/OPF.m
@@ -1,12 +1,12 @@
 tic
 clear
 %[kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,LineLimti,LineLimtj,LinePLimt,PG,QG,PD,QD,CenterA,LineCount,PGi,PVQU,PVQL]=pf('5sj.txt');
-[kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL]=pf('ieee3001.dat');
+[kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL]=pf('ieee10471.dat');
 GB=full(GB);
 %PVi电压节点序号
 %PVu电压节点电压标幺值
-Volt
-UAngel*180/3.1415926
+Volt;
+UAngel*180/3.1415926;
 %sprintf('%f\n',Volt);
 %sprintf('%f\n',Angel);
 %% 初值
diff --git a/OPF_Init.m b/OPF_Init.m
index 6a80fb1..f533a3c 100644
--- a/OPF_Init.m
+++ b/OPF_Init.m
@@ -2,8 +2,8 @@ function [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount]=O
 RestraintCount=size(PVi,1)+size(PGi,1)+Busnum; %约束条件数
 t_Bal_volt=Volt(Balance);
 Volt=ones(1,Busnum);
-%Volt(Balance)=t_Bal_volt;
-Volt(Balance)=1;
+Volt(Balance)=t_Bal_volt;
+%Volt(Balance)=1;
 UAngel=zeros(1,Busnum);
 Init_Z=ones(1,RestraintCount);
 Init_W=-1*ones(1,RestraintCount);
diff --git a/SolveIt.m b/SolveIt.m
index e7e8fe2..98405a4 100644
--- a/SolveIt.m
+++ b/SolveIt.m
@@ -8,15 +8,19 @@ aa=[
     deltH',zeros(size(Init_Y,2));
     ];
 yy=[LxComa;-Ly];
-% t=size(PVi,1)+size(PGi,1);
-% aa(t+2*Balance-1,:)=0;
-% aa(:,t+2*Balance-1)=0;
-% aa(t+2*Balance-1,t+2*Balance-1)=1;
-%暂时改一下
+%% 平衡节点电压不变 
+t=size(PVi,1)+size(PGi,1);
+aa(t+Balance,:)=0;
+aa(:,t+Balance)=0;
+aa(t+Balance,t+Balance)=1;
+deltG(t+Balance,:)=0;
+%%
 t=size(PVi,1)+size(PGi,1)+Busnum;
 aa(t+Balance,:)=0;
 aa(:,t+Balance)=0;
 aa(t+Balance,t+Balance)=1;
+deltG(t+Balance,:)=0;
+%%
 dxdy=aa\yy;
 dX=dxdy(1:ContrlCount);
 dY=dxdy(ContrlCount+1:ContrlCount+2*Busnum);